
      SUBROUTINE SCATT
C     SUBROUTINE TO COMPUTE SCATTER PARAMETERS
C
      COMMON /M/F,D,NS,A,DH,DHS,S,E,POL,KM
      COMMON /MP/ H1E,H2E,H1G,H2G,DLS1,DLS2,DL1,DL2,DL,DLS,TE1,TE2,TE,KL
      COMMON /MAR14/D3,D4,T5
      COMMON /MLDS/ AG,AD,AS,ACR,AED,MD,AH50,AH5,D5,MS,AES,DX,H5
      REAL NS,MD,MDO,MS,MSS,MDS,K1,K2,K3,K4
      KK=0
10    KK=KK+1
      D5=DL+ 200.
      D6=DL+ 400.
11    T5=TE+D5/A
      T6=TE+D6/A
      H5=AMIN1(((1./H1E +1./H2E)/(T5*F*ABS(.007-.058*T5))), (15.))
      H6=AMIN1(((1./H1E+1./H2E)/(T6*F*ABS(.007-.058*T6))),(15.))
      S5=H5+10.*ALOG10(F*T5**4)-.1*(NS-301.)*EXP(-T5*D5/40.)
      S6=H6+10.*ALOG10(F*T6**4)-.1*(NS-301.)*EXP(-T6*D6/40.)
      IF(T5*D5 .LE.10.) AH5=S5+103.4+.332*T5*D5-10.*ALOG10(T5*D5)
      IF(T6*D6 .LE.10.)AH6=S6+103.4+.332*T6*D6-10.*ALOG10(T6*D6)
      IF(T5*D5 .GT. 10. .AND.T5*D5.LE. 70.) AH5=S5+97.1+.212*T5*D5-2.5*
     CALOG10(T5*D5)
      IF(T6*D6 .GT. 10. .AND.T6*D6 .LE. 70.) AH6=S6+97.1+.212*T6*D6-2.5
     C*ALOG10(T6*D6)
      IF(T5*D5 .GT. 70.) AH5=S5+86.8+.157*T5*D5+5.*ALOG10(T5*D5)
      IF(T6*D6 .GT. 70.) AH6=S6+86.8+.157*T6*D6+5.*ALOG10(T6*D6)
      MS=(AH6-AH5)/(D6-D5)
      IF (KK .EQ. 2) GO TO 25
      IF (H5 .EQ. 10) GO TO 30
      IF (KK .EQ. 1) GO TO 20
25    MS=MSS
      AH50=AH5
      AH5=AH5S
      D5=D5S
      GO TO 30
20    DH=0.
C!!!  DL1=DLS1*EXP(-.07*SQRT(DH/H1E))
C!!!  DL2=DLS2*EXP(-.07*SQRT(DH/H2E))
      DL1=DLS1*EXP(-.07*SQRT(DH/AMAX1(5.,H1E)))
      DL2=DLS2*EXP(-.07*SQRT(DH/AMAX1(5.,H2E)))
      DL=DL1+DL2
C!!!  TE1=(.00065/DLS1)*((DLS1/DL1-1.)*DH-3.077*H1G)
C!!!  TE2=(.00065/DLS2)*((DLS2/DL2-1.)*DH-3.077*H2G)
      TE1=(.00065/DLS1)*((DLS1/DL1-1.)*DH-3.077*H1E)
      TE2=(.00065/DLS2)*((DLS2/DL2-1.)*DH-3.077*H2E)
      TE=AMAX1((TE1+TE2),(-DL/A))
      T=TE+D/A
      AH5S=AH5
      MSS=MS
      D5S=D5
      GO TO 10
30    CONTINUE
      RETURN
      END

