--- DarthVader/OrbitalInfo/src/igrf_sub.for	2011/11/23 21:19:34	1.5
+++ DarthVader/OrbitalInfo/src/igrf_sub.for	2012/05/15 15:00:51	1.7
@@ -191,7 +191,7 @@
 C          B0           MAGNETIC FIELD STRENGTH IN GAUSS
 C-----------------------------------------------------------------------
       DIMENSION         V(3),U(3,3),P(8,100),SP(3)
-      COMMON            X(3),H(144)
+      COMMON            X(3),H(196)
       COMMON/FIDB0/     SP
       SAVE /FIDB0/
       COMMON/GENER/     UMR,ERA,AQUAD,BQUAD
@@ -394,7 +394,7 @@
 C* CALLS ENTRY POINT FELDI IN GEOMAGNETIC FIELD SUBROUTINE FELDG   *
 C*******************************************************************
       DIMENSION         P(7),U(3,3)
-      COMMON            XI(3),H(144)
+      COMMON            XI(3),H(196)
 C*****XM,YM,ZM  ARE GEOMAGNETIC CARTESIAN INVERSE CO-ORDINATES          
       ZM=P(3)                                                           
       FLI=P(1)*P(1)+P(2)*P(2)+1E-15
@@ -476,11 +476,11 @@
 C                 POINTING IN THE TANGENTIAL PLANE TO THE NORTH, EAST
 C                 AND DOWNWARD.   
 C-----------------------------------------------------------------------
-      DIMENSION         V(3),B(3),G(144)
+      DIMENSION         V(3),B(3),G(196)
       CHARACTER*258      NAME
       INTEGER NMAX
       REAL TIME
-      COMMON            XI(3),H(144)
+      COMMON            XI(3),H(196)
       COMMON/MODEL/     G,NMAX,TIME,NAME
       SAVE/MODEL/
       COMMON/GENER/     UMR,ERA,AQUAD,BQUAD
@@ -581,8 +581,8 @@
       CHARACTER*258    FIL1, FIL2           
       CHARACTER*258    FILMOD
 C     ### FILMOD, DTEMOD arrays +1
-c     DIMENSION       GH1(144),GH2(120),GHA(144),FILMOD(14),DTEMOD(14)
-      DIMENSION       GH1(144),GH2(120),GHA(144),FILMOD(3),DTEMOD(3)
+c     DIMENSION       GH1(144),GH2(120),GHA(144),FILMOD(14),DTEMOD(14) 
+      DIMENSION       GH1(196),GH2(196),GHA(196),FILMOD(3),DTEMOD(3)
       DOUBLE PRECISION X,F0,F 
       INTEGER L1,L2,L3
       INTEGER NMAX
@@ -615,9 +615,9 @@
 c     FILMOD(1) = 'OrbitalInfo/src/dgrf00.dat'
 c     FILMOD(2) = 'OrbitalInfo/src/igrf05.dat'
 c     FILMOD(3) = 'OrbitalInfo/src/igrf05s.dat'
-c     WRITE(*,*) FILMOD(1)
-c     WRITE(*,*) FILMOD(2)
-c     WRITE(*,*) FILMOD(3)
+c      WRITE(*,*) FILMOD(1)
+c      WRITE(*,*) FILMOD(2)
+c      WRITE(*,*) FILMOD(3)
 c      DATA   FILMOD / 'dgrf00.dat', 'igrf05.dat', 'igrf05s.dat'/
       DATA   DTEMOD / 2005., 2010., 2015./
 c     
@@ -652,10 +652,13 @@
       FIL1 = FILMOD(L)   
       DTE2 = DTEMOD(L+1) 
       FIL2 = FILMOD(L+1) 
+c      WRITE(*,*) FIL1
+c      WRITE(*,*) FIL2
 c      print *, "que"
 C--   GET IGRF COEFFICIENTS FOR THE BOUNDARY YEARS
       CALL GETSHC (IU, FIL1, NMAX1, ERAD, GH1, IER)  
       IF (IER .NE. 0) STOP                           
+c      print *, "quessss"
       CALL GETSHC (IU, FIL2, NMAX2, ERAD, GH2, IER)  
       IF (IER .NE. 0) STOP                    
 c      print *, "quj"
@@ -944,7 +947,7 @@
 C ASTRONOMICAL UNION .
 C-----------------------------------------------------------------
         INTEGER TL1,TL2,TL3
-        CHARACTER *258 TP1,TP2,TP3
+        CHARACTER (len=258) TP1,TP2,TP3
         INTEGER L1,L2,L3
         CHARACTER *258 P1,P2,P3
         COMMON/PPATH/ L1,L2,L3,P1, P2, P3