diff --git a/geopack/src/Modif.for b/geopack/src/Modif.for
index 85b7243..d7a8198 100644
--- a/geopack/src/Modif.for
+++ b/geopack/src/Modif.for
@@ -7,4 +7,15 @@
       CALL TRACE_08 (XI,YI,ZI,DIR,DSMAX,ERR,RLIM,R0,IOPT,PARMOD,
      * T96_01,IGRF_GSW_08,XF,YF,ZF,XX,YY,ZZ,L,LMAX)
       
-      END
\ No newline at end of file
+      END
+
+      SUBROUTINE T96_01_MODIFIED (IOPT,PARMOD,X,Y,Z,BX,BY,BZ)
+
+      DIMENSION PARMOD(10)
+      EXTERNAL T96_01
+
+      COMMON /GEOPACK1/ PSI
+
+      CALL T96_01 (IOPT,PARMOD,PSI,X,Y,Z,BX,BY,BZ)
+
+      END
diff --git a/geopack/src/include/geopack.hh b/geopack/src/include/geopack.hh
index 8f32228..5ddecc0 100644
--- a/geopack/src/include/geopack.hh
+++ b/geopack/src/include/geopack.hh
@@ -22,8 +22,13 @@ extern "C" {
         float* FP_GSW_X_i, float* FP_GSW_Y_i, float* FP_GSW_Z_i,
         float XX[], float YY[], float ZZ[], int* L, int* LMAX);
 
+    void t96_01_modified_ (int* IOPT, float PARMOD[], float* sat_pos_GSW_X_i, float* sat_pos_GSW_Y_i,
+        float* sat_pos_GSW_Z_i, float* Bx_GSW, float* By_GSW, float* Bz_GSW);
+
+    void igrf_gsw_08_ (float* pos_GSW_X_i, float* pos_GSW_Y_i, float* pos_GSW_Z_i, float* HX_GSW, float* HYGSW, float* HZGSW);
+
 #ifdef __cplusplus
 }
 #endif
 
-#endif
\ No newline at end of file
+#endif
--
libgit2 0.21.2