anchor.diff:        IF(SW(7).EQ.0) GOTO 200
         T71 = P(12)*PLG(3,2)*CD14*SWC(5)
         T72 = P(13)*PLG(3,2)*CD14*SWC(5)
  -      T(7) = 
  +      T(7) =
        1 ((P(4)*PLG(2,2) + P(5)*PLG(4,2)
        2 + T71)*CTLOC
        4 + (P(7)*PLG(2,2) + P(8)*PLG(4,2)
  
positive.code: C
          IF(MASS.EQ.48) THEN
           D(6) = 1.66E-24*(4.*D(1)+16.*D(2)+28.*D(3)+32.*D(4)+40.*D(5)+
       &       D(7)+14.*D(8))  
           IF(IMR.EQ.1) D(6)=D(6)/1000.
           ENDIF
           T(2)=TZ
     10 CONTINUE
        GOTO 90
     50 CONTINUE
        DD=DENSM(ALT,1.,0,TZ,MN3,ZN3,TN3,TGN3,MN2,ZN2,TN2,TGN2)                
        T(2)=TZ
     90 CONTINUE
        ALAST=ALT
  
negative.code: C        D(2) - O NUMBER DENSITY(CM-3)
  C        D(3) - N2 NUMBER DENSITY(CM-3)
  C        D(4) - O2 NUMBER DENSITY(CM-3)
  C        D(5) - AR NUMBER DENSITY(CM-3)
  C        D(6) - TOTAL MASS DENSITY(GM/CM3) [includes anomalous oxygen]
  C        D(7) - H NUMBER DENSITY(CM-3)
  C        D(8) - N NUMBER DENSITY(CM-3)
  