anchor.diff:          GO TO 50
      50 CONTINUE
         CCOR=EXP(CCOR)
  -      END FUNCTION  CCOR
  +       RETURN
  +      END
   C-----------------------------------------------------------------------
         FUNCTION  CCOR2(ALT, R,H1,ZH,H2)
   C       O&O2 CHEMISTRY/DISSOCIATION CORRECTION FOR MSIS MODELS
  
positive.code:       ALAST=ALT
        RETURN
    100 FORMAT(1X,'MASS', I5, '  NOT VALID')
        END Subroutine GTS7
  C-----------------------------------------------------------------------
        SUBROUTINE METERS(METER)
        Implicit None
  C      Convert outputs to Kg & Meters if METER true
        logical,Intent(In) :: METER
        Integer IMR
        COMMON/METSEL/IMR
        SAVE
        IMR=0
        IF(METER) IMR=1
        END SUBROUTINE METERS
  C-----------------------------------------------------------------------
        Real FUNCTION SCALH(ALT,XM,TEMP)
        Implicit None
        real, intent(in) :: alt, xm, temp
        real gsurf,re, g
        real,parameter :: RGAS=831.4
  C      Calculate scale height (km)
        COMMON/PARMB/GSURF,RE
        SAVE
        G=GSURF/(1.+ALT/RE)**2
        SCALH=RGAS*TEMP/(G*XM)
        END FUNCTION SCALH
  C-----------------------------------------------------------------------
        Real FUNCTION GLOBE7(YRD,SEC,LAT,LONG,TLOC,F107A,F107,AP,P)
  C       CALCULATE G(L) FUNCTION
  C       Upper Thermosphere Parameters
        Real,Intent(IN) :: YRD,SEC,LAT, LONG,TLOC,F107A,F107,AP(*)
  !     inout since P modified in function
        Real,Intent(inout) :: P(*)
  
  
        Real SV(25)
        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,
  
negative.code:       DATA MN2/4/,ZN2/72.5,55.,45.,32.5/
        DATA ZMIX/62.5/,ALAST/99999./,MSSL/-999/
        DATA SV/25*1./
        IF(ISW.NE.64999) CALL TSELEC(SV)
  C      Put identification data into common/datime/
        DO 1 I=1,3
  