anchor.diff:     10 CONTINUE
         GOTO 90
      50 CONTINUE
  -      DD=DENSM(ALT,1.,0,TZ,MN3,ZN3,TN3,TGN3,MN2,ZN2,TN2,TGN2)                
  +      DD=DENSM(ALT,1.,0.,TZ,MN3,ZN3,TN3,TGN3,MN2,ZN2,TN2,TGN2)                
         T(2)=TZ
      90 CONTINUE
         ALAST=ALT
  
positive.code: C        D(9) - Anomalous oxygen NUMBER DENSITY(CM-3)
  C        T(1) - EXOSPHERIC TEMPERATURE
  C        T(2) - TEMPERATURE AT ALT
  C
        DIMENSION D(9),T(2),AP(7),DS(9),TS(2)
        COMMON/METSEL/IMR
        CALL GTD7(IYD,SEC,ALT,GLAT,GLONG,STL,F107A,F107,AP,MASS,D,T)
  C       TOTAL MASS DENSITY
  
negative.code: C
  C        To get current values of SW: CALL TRETRV(SW)
  C
        Real,Intent(OUT)   :: D(9), T(2)
        CHARACTER(len=4) :: ISDATE(3),ISTIME(2),NAME(2),
       & ISD(3),IST(2),NAM(2)
        
        DIMENSION AP(7),DS(9),TS(2)
        DIMENSION ZN3(5),ZN2(4),SV(25)
        COMMON/GTS3C/TLB,S,DB04,DB16,DB28,DB32,DB40,DB48,DB01,ZA,T0,Z0
       & ,G0,RL,DD,DB14,TR12
  