--- DarthVader/OrbitalInfo/src/igrf_sub.for	2009/08/04 14:01:27	1.3
+++ DarthVader/OrbitalInfo/src/igrf_sub.for	2009/08/11 12:58:22	1.4
@@ -2,17 +2,17 @@
      &          xl,icode,dip,dec)
 c----------------------------------------------------------------
 c   INPUT:
-c	xlat	geodatic latitude in degrees
-c	xlong	geodatic longitude in degrees
-c	year	decimal year (year+month/12.0-0.5 or year+day-of-year/365 
-c		or 366 if leap year) 
-c	height	height in km
+c     xlat     geodatic latitude in degrees
+c     xlong    geodatic longitude in degrees
+c     year     decimal year (year+month/12.0-0.5 or year+day-of-year/365 
+c          or 366 if leap year) 
+c     height     height in km
 c   OUTPUT:
-c	xl	L value
-c	icode	=1  L is correct; =2  L is not correct;
-c		=3  an approximation is used
-c	dip	geomagnetic inclination in degrees
-c	dec	geomagnetic declination in degress
+c     xl     L value
+c     icode     =1  L is correct; =2  L is not correct;
+c          =3  an approximation is used
+c     dip     geomagnetic inclination in degrees
+c     dec     geomagnetic declination in degress
 c----------------------------------------------------------------
 
       REAL              LATI,LONGI
@@ -20,7 +20,7 @@
       SAVE /GENER/
 C
       CALL INITIZE
-	ibbb=0
+      ibbb=0
       ALOG2=ALOG(2.)
       ISTART=1
         lati=xlat
@@ -32,7 +32,7 @@
         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
+      DEC=ASIN(BEAST/SQRT(BEAST*BEAST+BNORTH*BNORTH))/UMR
       RETURN
       END
 c
@@ -342,17 +342,17 @@
 C-- because 1E-38 is the minimal allowable arg. for ALOG in our envir.
 C-- D. Bilitza, Nov 87.
 C
-11    	FI=0.5*ABS(FI)/SQRT(B0)+1E-12                       
+11         FI=0.5*ABS(FI)/SQRT(B0)+1E-12                       
 C
 C*****COMPUTE L FROM B AND I.  SAME AS CARMEL IN INVAR.  
 C
 C-- Correct dipole moment is used here. D. Bilitza, Nov 87.
 C
       DIMOB0=DIMO/B0
-	arg1=alog(FI)
-	arg2=alog(DIMOB0)
-c	arg = FI*FI*FI/DIMOB0
-c	if(abs(arg).gt.88.0) arg=88.0
+      arg1=alog(FI)
+      arg2=alog(DIMOB0)
+c    arg = FI*FI*FI/DIMOB0
+c    if(abs(arg).gt.88.0) arg=88.0
       XX=3*arg1-arg2
       IF(XX.GT.23.0) GOTO 776   
       IF(XX.GT.11.7) GOTO 775