C	-*-fortran-*- CRANE'S GLOBAL MODEL
C                REF: TP 1770, PP.27-30
C    ************************************************
      REAL R,H,L,G,THETA,D,ALPHA,A,B,D1,X,Y,U,TERM1,TERM2
      REAL TERM3,LEFF,ATT,KE

C                      INPUTS
      PARAMETER(PI=3.1415927)
      PARAMETER(E=2.7182818)
      WRITE(6,*) 'ENTER GROUND HEIGHT (KM), ELEVATION (DEG)'
      READ(5,*)  G,THETA
      WRITE(6,*) 'A*R^B COEFFICIENTS A,B'
      READ(5,*)  A,B
      WRITE(6,*) ' 0 RAIN RATE ENDS PROGRAM'
      WRITE(6,*) '  '
110      WRITE(6,*) 'RAIN-RATE R=? RAIN HEIGHT H=?'
      READ(5,*)  R,H
      IF (R .EQ. 0.0) THEN 
        GOTO 22
      END IF
C   ************************************************
      L=(H-G)/SIN(THETA*PI/180)
      D=L*COS(THETA*PI/180)
      WRITE(6,*)  'HORIZONTAL DISTANCE D=',D
      WRITE(6,*)  'SLANT PATH DISTANCE L=',L
      ALPHA=A*R**B
      D1=3.8-.6*LOG(R)
      IF (D1 .GT. 22.5) THEN 
        GOTO 350
      END IF
      X=2.3*R**(-.17)
      Y=.026-.03*LOG(R)
      U=(1/D1)*(LOG(X)+Y*D1)
      TERM1=(E**(U*B*D1)-1)/(U*B)
      TERM2=(X**B*E**(Y*B*D1))/(Y*B)
      TERM3=(X**B*E**(Y*B*D))/(Y*B)
      IF (D .LT. D1) THEN 
         TERM1=(E**(U*B*D)-1)/(U*B)
         TERM2=0
         TERM3=0
      END IF 
      LEFF=(TERM1-TERM2+TERM3)/COS(THETA*PI/180)
      ATT=ALPHA*LEFF
      KE=LEFF/L
      WRITE(6,*)  'R=',R,'   ALPHA=',ALPHA,'   KE=',KE,'   ATT=',ATT
      GOTO 110
350      WRITE(6,*)  'RAIN RATE OUT OF RANGE'
      GOTO 110
22      END