--- DarthVader/OrbitalInfo/src/igrf_sub.for	2012/05/15 15:00:51	1.7
+++ DarthVader/OrbitalInfo/src/igrf_sub.for	2014/08/29 07:06:41	1.10
@@ -1,5 +1,5 @@
-        subroutine igrf_sub(xlat,xlong,year,height,
-     &          xl,icode,dip,dec)
+c        subroutine igrf_sub(xlat,xlong,year,height,
+c     &          xl,icode,dip,dec)
 c----------------------------------------------------------------
 c   INPUT:
 c     xlat     geodatic latitude in degrees
@@ -14,27 +14,27 @@
 c     dip     geomagnetic inclination in degrees
 c     dec     geomagnetic declination in degress
 c----------------------------------------------------------------
-
-      REAL              LATI,LONGI
-      COMMON/GENER/     UMR,ERA,AQUAD,BQUAD
-      SAVE /GENER/
-C
-      CALL INITIZE
-      ibbb=0
-      ALOG2=ALOG(2.)
-      ISTART=1
-        lati=xlat
-        longi=xlong
+c
+c      REAL              LATI,LONGI
+c      COMMON/GENER/     UMR,ERA,AQUAD,BQUAD
+c      SAVE /GENER/
+C
+c      CALL INITIZE
+c      ibbb=0
+c      ALOG2=ALOG(2.)
+c      ISTART=1
+c        lati=xlat
+c        longi=xlong
 c
 C----------------CALCULATE PROFILES-----------------------------------
 c
-        CALL FELDCOF(YEAR,DIMO)
-        CALL FELDG(LATI,LONGI,HEIGHT,BNORTH,BEAST,BDOWN,BABS)
-        CALL SHELLG(LATI,LONGI,HEIGHT,DIMO,XL,ICODE,BAB1)
-        DIP=ASIN(BDOWN/BABS)/UMR
-      DEC=ASIN(BEAST/SQRT(BEAST*BEAST+BNORTH*BNORTH))/UMR
-      RETURN
-      END
+c        CALL FELDCOF(YEAR,DIMO)
+c        CALL FELDG(LATI,LONGI,HEIGHT,BNORTH,BEAST,BDOWN,BABS)
+c        CALL SHELLG(LATI,LONGI,HEIGHT,DIMO,XL,ICODE,BAB1)
+c        DIP=ASIN(BDOWN/BABS)/UMR
+c      DEC=ASIN(BEAST/SQRT(BEAST*BEAST+BNORTH*BNORTH))/UMR
+c      RETURN
+c      END
 c
 c
 C SHELLIG.FOR, Version 2.0, January 1992
@@ -677,7 +677,7 @@
          F = GHA(J) * 1.D-5
          F0 = F0 + F * F
  1234 CONTINUE
-      DIMO = DSQRT(F0)
+      DIMO = REAL(DSQRT(F0))
       
       GH1(1) =  0.0
       I=2          
@@ -693,13 +693,13 @@
          IF(IS.EQ.0) F0 = F0 * (2.D0 * X - 1.D0) / X
          F = F0 * 0.5D0                                    
          IF(IS.EQ.0) F = F * SQRT2
-         GH1(I) = GHA(I-1) * F0
+         GH1(I) = GHA(I-1) * REAL(F0)
          I = I+1                                         
          DO 9 M=1,N                                    
             F = F * (X + M) / (X - M + 1.D0)                 
             IF(IS.EQ.0) F = F * DSQRT((X - M + 1.D0) / (X + M))
-            GH1(I) = GHA(I-1) * F
-            GH1(I+1) = GHA(I) * F
+            GH1(I) = GHA(I-1) * REAL(F)
+            GH1(I+1) = GHA(I) * REAL(F)
             I=I+2
  9       CONTINUE                                          
       RETURN
@@ -947,7 +947,7 @@
 C ASTRONOMICAL UNION .
 C-----------------------------------------------------------------
         INTEGER TL1,TL2,TL3
-        CHARACTER (len=258) TP1,TP2,TP3
+        CHARACTER (len=*) :: TP1,TP2,TP3
         INTEGER L1,L2,L3
         CHARACTER *258 P1,P2,P3
         COMMON/PPATH/ L1,L2,L3,P1, P2, P3