#
# A simple test of rt4 using a 2 mm/hr rain layer (spherical drops) at 85 GHz.
#
rm -f testr.lay testr.sca testr.rts testr.out testr.out.check
#   Output the layer file
cat >testr.lay <<EOF
 1.00  280.0   0.03  'testr.rts'
 0.00  285.0   0.00  ' '
EOF
#   Output the scattering file
cat >testr.sca <<EOF
C  WAVELENGTH=   .3506000E+00
C  COATED_SPHERE=N
C  REFRACT_INDEX=   .3170000E+01  -.1750000E+01
C  RADIUS_RANGE=    .3000000E-02    .3000000E+00
C  NUMBER_INTEGRATION_STEPS= 100
C  GAMMA_DIST_A=   .1600000E+05
C  GAMMA_DIST_B=   .7100000E+02
C  GAMMA_DIST_ALPHA=   .0000000E+00
C  GAMMA_DIST_GAMMA=   .1000000E+01
   .5069279E+00    EXINCTION
   .2169107E+00    SCATTERING
   .4278926    SINGLE SCATTERING ALBEDO
   13    DEGREE OF LEGENDRE SERIES
   0    1.00000000  -.32726057   .16964847  -.33956743  1.00000000   .16964847
   1     .54995612  -.11363125  1.47401795  -.15818725   .54995612  1.47401795
   2     .57262352   .27346170   .55901253   .27545005   .57262352   .55901253
   3     .19675512   .09886793   .17743348   .13869054   .19675512   .17743348
   4     .08139809   .04933270   .04076688   .05860883   .08139809   .04076688
   5     .02272976   .01365347   .00933083   .01807309   .02272976   .00933083
   6     .00664139   .00418458   .00173696   .00515890   .00664139   .00173696
   7     .00166838   .00104457   .00031744   .00134169   .00166838   .00031744
   8     .00041854   .00026656   .00004724   .00033132   .00041854   .00004724
   9     .00009750   .00006206   .00000545   .00007808   .00009750   .00000545
  10     .00002228   .00001438  -.00000012   .00001759   .00002228  -.00000012
  11     .00000477   .00000311  -.00000032   .00000372   .00000477  -.00000032
  12     .00000096   .00000063  -.00000012   .00000072   .00000096  -.00000012
  13     .00000017   .00000011  -.00000003   .00000013   .00000017  -.00000003
EOF
#
#   Use scatcnv to make a radiative transfer scattering file for rt4
scatcnv <<EOF >/dev/null
testr.sca
testr.rts
8
G
0
EOF
#
#   Run rt4 for a water surface
rt4 <<EOF >/dev/null
2
8
G
testr.lay
285
F
(3.17,-1.75)
3.0
3510
T
VH
2
1 2
testr.out
EOF
#
#
#   Compare output with results below
cat >testr.out.check <<EOF
C  NUMMU=  8  NUMAZI=  1  AZIORDER=  0  NSTOKES=2
C  LAYER_FILE=testr.lay                       
C  SRC_CODE=2   QUAD_TYPE=GAUSSIAN        
C  GROUND_TEMP=  285.00   GROUND_TYPE=FRESNEL         
C  GROUND_INDEX=   3.1700  -1.7500   SKY_TEMP=    3.00
C  WAVELENGTH= .351000E+04
C  UNITS=KELVINS - EBB               OUTPUT_POLARIZATION=VH
C    Z       MU    FLUX/RADIANCE (V,H)
   1.000 -2.00000  242.02  225.56
   1.000  2.00000    3.00    3.00
   1.000  -.98940  225.90  224.94
   1.000  -.94458  229.15  224.20
   1.000  -.86563  234.69  223.20
   1.000  -.75540  241.77  222.61
   1.000  -.61788  248.71  223.53
   1.000  -.45802  252.11  227.11
   1.000  -.28160  246.92  231.53
   1.000  -.09501  231.01  224.90
   1.000   .09501    3.00    3.00
   1.000   .28160    3.00    3.00
   1.000   .45802    3.00    3.00
   1.000   .61788    3.00    3.00
   1.000   .75540    3.00    3.00
   1.000   .86563    3.00    3.00
   1.000   .94458    3.00    3.00
   1.000   .98940    3.00    3.00
    .000 -2.00000  247.63  210.72
    .000  2.00000  140.52  137.63
    .000  -.98940  214.61  213.12
    .000  -.94458  219.16  211.35
    .000  -.86563  227.43  208.47
    .000  -.75540  239.41  205.22
    .000  -.61788  254.55  203.38
    .000  -.45802  270.48  206.83
    .000  -.28160  280.64  222.67
    .000  -.09501  276.36  246.81
    .000   .09501  254.67  243.13
    .000   .28160  210.67  203.14
    .000   .45802  167.88  163.23
    .000   .61788  139.39  136.60
    .000   .75540  120.87  119.30
    .000   .86563  108.93  108.14
    .000   .94458  101.60  101.30
    .000   .98940   97.83   97.77
EOF
#
echo "Differences in output listed below:"
diff testr.out.check testr.out | more
