anchor.diff:        DATA RGAS/831.4/
         G=GSURF/(1.+ALT/RE)**2
         SCALH=RGAS*TEMP/(G*XM)
  -      RETURN
  -      END
  +      END FUNCTION SCALH
   C-----------------------------------------------------------------------
         FUNCTION GLOBE7(YRD,SEC,LAT,LONG,TLOC,F107A,F107,AP,P)
   C       CALCULATE G(L) FUNCTION 
   C       Upper Thermosphere Parameters
         REAL LAT, LONG
  -      DIMENSION P(1),SV(25),AP(1)
  +      Real P(*),SV(25),AP(*)
         COMMON/TTEST/TINF,GB,ROUT,T(15)
         COMMON/CSW/SW(25),ISW,SWC(25)
         COMMON/LPOLY/PLG(9,4),CTLOC,STLOC,C2TLOC,S2TLOC,C3TLOC,S3TLOC,
  
positive.code:       DO 50 I = 1,NSW
     50 TINF = TINF + ABS(SW(I))*T(I)
        GLOBE7 = TINF
        RETURN
        END
  C-----------------------------------------------------------------------
        SUBROUTINE TSELEC(SV)
  C        SET SWITCHES
  
negative.code:       FUNCTION VTST7(IYD,SEC,GLAT,GLONG,STL,F107A,F107,AP,IC)
  C       Test if geophysical variables or switches changed and save
  C       Return 0 if unchanged and 1 if changed
        Real AP(*)
        DIMENSION IYDL(2),SECL(2),GLATL(2),GLL(2),STLL(2)
        DIMENSION FAL(2),FL(2),APL(7,2),SWL(25,2),SWCL(25,2)
        COMMON/CSW/SW(25),ISW,SWC(25)
        SAVE
  