#!/bin/csh
# Runs a general test of RT3 with multiple layers and solar and thermal sources
#
rm -f testa.lay rayleigh.sca mietest.sca testa.out.check
#   Output the layer file
cat >testa.lay <<EOF
15.00  200.0   0.02  'rayleigh.sca'
 5.00  270.0   0.05  'mietest.sca'
 0.00  300.0   0.00  ' '
EOF
#
#   Output the scattering files
cat >rayleigh.sca <<EOF
   .0500000E+00    EXINCTION
   .0500000E+00    SCATTERING
  1.0000000    SINGLE SCATTERING ALBEDO
    2    DEGREE OF LEGENDRE SERIES
   0    1.00000000  -.50000000   .00000000   .00000000  1.00000000   .00000000
   1     .00000000   .00000000  1.50000000   .00000000   .00000000  1.50000000
   2     .50000000   .50000000   .00000000   .00000000   .50000000   .00000000
EOF
#
cat >mietest.sca <<EOF
   .1000000E+01    EXINCTION
   .9900000E+00    SCATTERING
   .9900000    SINGLE SCATTERING ALBEDO
   11    DEGREE OF LEGENDRE SERIES
   0    1.00000000  -.32071711   .71206342  -.01882245  1.00000000   .71206342
   1    1.45529318  -.20350675  1.76014119  -.04725108  1.45529318  1.76014119
   2    1.05402631   .24638948  1.06682431   .00894436  1.05402631  1.06682431
   3     .39758994   .18605748   .39651104   .04505815   .39758994   .39651104
   4     .11659302   .07124848   .09576412   .00958275   .11659302   .09576412
   5     .02387477   .01700757   .01765088   .00215761   .02387477   .01765088
   6     .00395010   .00302534   .00261549   .00029195   .00395010   .00261549
   7     .00053888   .00043592   .00032713   .00003502   .00053888   .00032713
   8     .00006372   .00005326   .00003583   .00000337   .00006372   .00003583
   9     .00000667   .00000572   .00000351   .00000029   .00000667   .00000351
  10     .00000063   .00000055   .00000031   .00000002   .00000063   .00000031
  11     .00000006   .00000005   .00000003   .00000000   .00000006   .00000003
EOF
#
#   Run rt3 
rt3 <<EOF >&! rt3.log
4           Nstokes
4           Nquad
G
4           Nazimuth
testa.lay
N           no delta-M scaling
3           solar and thermal source
5.0         solar flux
60.0        solar angle
300.0       ground temp
L
0.25        ground albedo
0.0         sky temp
3.0         wavelength (microns)
W
IQ
3           number of output levels
1 2 3       output level numbers
3           number of output azimuths
testa.out
EOF
#
#
#   Compare output with results below
cat >testa.out.check <<EOF
C  NUMMU=  4  NUMAZI=  9  AZIORDER=  4  NSTOKES=4
C  LAYER_FILE=testa.lay                          DELTA-M=N
C  SRC_CODE=3   QUAD_TYPE=GAUSSIAN        
C  DIRECT_FLUX= .50000E+01   DIRECT_MU= .500000
C  GROUND_TEMP=  300.00   GROUND_TYPE=LAMBERTIAN      
C  GROUND_ALBEDO=  .25000   SKY_TEMP=     .00
C  WAVELENGTH= .300000E+01
C  UNITS=WATTS/(M^2 MICRON STER)     OUTPUT_POLARIZATION=IQ
C    Z      PHI     MU    FLUX/RADIANCE (I,Q,U,V)
  15.000    .0 -2.00000   .173477E+01   .103900E+00   .000000E+00   .000000E+00
  15.000    .0  2.00000   .500000E+01   .000000E+00   .000000E+00   .000000E+00
  15.000    .0  -.96029   .415400E+00  -.188827E+00   .000000E+00   .000000E+00
  15.000    .0  -.79667   .451745E+00  -.219329E+00   .000000E+00   .000000E+00
  15.000    .0  -.52553   .579206E+00  -.203552E+00   .000000E+00   .000000E+00
  15.000    .0  -.18343   .799029E+00  -.137430E+00   .000000E+00   .000000E+00
  15.000    .0   .18343   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000    .0   .52553   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000    .0   .79667   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000    .0   .96029   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000  90.0  -.96029   .443776E+00   .145556E+00   .508178E-01   .100812E-03
  15.000  90.0  -.79667   .465895E+00   .155671E+00   .123671E+00   .172961E-03
  15.000  90.0  -.52553   .506812E+00   .180072E+00   .215800E+00   .119817E-03
  15.000  90.0  -.18343   .572843E+00   .230011E+00   .347868E+00   .250606E-04
  15.000  90.0   .18343   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000  90.0   .52553   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000  90.0   .79667   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000  90.0   .96029   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000 180.0  -.96029   .496458E+00  -.882229E-01   .203017E-07  -.120093E-10
  15.000 180.0  -.79667   .608726E+00  -.173126E-01   .124362E-07  -.290620E-10
  15.000 180.0  -.52553   .750658E+00   .267501E-01   .488537E-09  -.297141E-10
  15.000 180.0  -.18343   .917189E+00  -.951650E-02  -.209440E-07  -.424919E-11
  15.000 180.0   .18343   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000 180.0   .52553   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000 180.0   .79667   .000000E+00   .000000E+00   .000000E+00   .000000E+00
  15.000 180.0   .96029   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   5.000    .0 -2.00000   .124860E+01   .450068E-01   .000000E+00   .000000E+00
   5.000    .0  2.00000   .241529E+01   .820379E-01   .000000E+00   .000000E+00
   5.000    .0  -.96029   .316952E+00  -.571097E-01   .000000E+00   .000000E+00
   5.000    .0  -.79667   .378777E+00  -.746844E-01   .000000E+00   .000000E+00
   5.000    .0  -.52553   .528770E+00  -.771516E-01   .000000E+00   .000000E+00
   5.000    .0  -.18343   .836870E+00  -.524060E-01   .000000E+00   .000000E+00
   5.000    .0   .18343   .640948E+00   .851421E-02   .000000E+00   .000000E+00
   5.000    .0   .52553   .556200E+00   .221156E-01   .000000E+00   .000000E+00
   5.000    .0   .79667   .420970E+00  -.145090E-01   .000000E+00   .000000E+00
   5.000    .0   .96029   .319709E+00  -.692199E-01   .000000E+00   .000000E+00
   5.000  90.0  -.96029   .310467E+00   .406900E-01   .158511E-01   .199327E-03
   5.000  90.0  -.79667   .336355E+00   .480606E-01   .404127E-01   .389393E-03
   5.000  90.0  -.52553   .386372E+00   .658319E-01   .752772E-01   .381569E-03
   5.000  90.0  -.18343   .460334E+00   .100630E+00   .127413E+00   .373166E-04
   5.000  90.0   .18343   .438575E+00   .148283E+00   .202021E+00   .428581E-04
   5.000  90.0   .52553   .362782E+00   .134683E+00   .157325E+00   .241136E-04
   5.000  90.0   .79667   .301335E+00   .118329E+00   .922815E-01   .131270E-04
   5.000  90.0   .96029   .274216E+00   .111656E+00   .382001E-01   .529811E-05
   5.000 180.0  -.96029   .313108E+00  -.204385E-01   .540970E-08  -.240508E-10
   5.000 180.0  -.79667   .346697E+00   .166653E-02   .343395E-08  -.676079E-10
   5.000 180.0  -.52553   .398128E+00   .151647E-01   .481160E-09  -.106247E-09
   5.000 180.0  -.18343   .453772E+00   .758518E-02  -.501378E-08  -.967635E-10
   5.000 180.0   .18343   .566832E+00  -.656012E-01  -.231479E-07  -.374676E-11
   5.000 180.0   .52553   .390841E+00  -.143244E+00  -.278118E-07  -.210807E-11
   5.000 180.0   .79667   .273935E+00  -.161544E+00  -.256513E-07  -.114760E-11
   5.000 180.0   .96029   .246343E+00  -.142586E+00  -.223435E-07  -.463174E-12
    .000    .0 -2.00000   .262221E+00   .000000E+00   .000000E+00   .000000E+00
    .000    .0  2.00000   .510073E+00  -.474961E-02   .000000E+00   .000000E+00
    .000    .0  -.96029   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.79667   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.52553   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.18343   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000    .0   .18343   .124286E+00  -.409200E-02   .000000E+00   .000000E+00
    .000    .0   .52553   .151508E+00  -.428797E-02   .000000E+00   .000000E+00
    .000    .0   .79667   .177364E+00  -.560435E-02   .000000E+00   .000000E+00
    .000    .0   .96029   .191069E+00  -.970279E-02   .000000E+00   .000000E+00
    .000  90.0  -.96029   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.79667   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.52553   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.18343   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000  90.0   .18343   .121056E+00  -.102552E-02   .196468E-02  -.126150E-04
    .000  90.0   .52553   .146406E+00   .221397E-02   .336247E-02  -.516552E-04
    .000  90.0   .79667   .170769E+00   .742230E-02   .428931E-02  -.834948E-04
    .000  90.0   .96029   .187139E+00   .119689E-01   .276823E-02  -.572090E-04
    .000 180.0  -.96029   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.79667   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.52553   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.18343   .825162E-01   .000000E+00   .000000E+00   .000000E+00
    .000 180.0   .18343   .122545E+00  -.628812E-02  -.457769E-09  -.108266E-10
    .000 180.0   .52553   .146764E+00  -.923285E-02  -.103255E-08  -.957088E-11
    .000 180.0   .79667   .168759E+00  -.129317E-01  -.181443E-08  -.255211E-11
    .000 180.0   .96029   .184590E+00  -.147828E-01  -.235030E-08   .252312E-11
EOF
#
echo "Differences in output listed below:"
diff testa.out.check testa.out | more


