anchor.diff:        DENSM=D0
         IF(ALT.GT.ZN2(1)) GOTO 50
   C      STRATOSPHERE/MESOSPHERE TEMPERATURE
  -      Z=MAX(ALT,ZN2(MN2))
  +      Z=AMAX1(ALT,ZN2(MN2))
         MN=MN2
         Z1=ZN2(1)
         Z2=ZN2(MN)
  
positive.code:       PLG(7,3)=(11.*C*PLG(6,3)-7.*PLG(5,3))/4.
        PLG(8,3)=(13.*C*PLG(7,3)-8.*PLG(6,3))/5.
        PLG(4,4) = 15.*S2*S
        PLG(5,4) = 105.*S2*S*C
        PLG(6,4)=(9.*C*PLG(5,4)-7.*PLG(4,4))/2.
        PLG(7,4)=(11.*C*PLG(6,4)-8.*PLG(5,4))/3.
        XL=LAT
  
negative.code:      *  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
       *  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
       *  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
  C          S PARAM  
        DATA PJ1/
       *  9.56827E-01, 6.20637E-02, 3.18433E-02, 0.00000E+00, 0.00000E+00,
       *  3.94900E-02, 0.00000E+00, 0.00000E+00,-9.24882E-03,-7.94023E-03,
  