#!/bin/csh
# Runs the Mie test described in "A New Polarized Atmospheric Radiative 
#   Transfer Model", K. F. Evans and G. L. Stephens, J. Quant. Spectrosc.
#    Radiative Transfer, Vol 46, No. 5, pp 413-423, 1991.
#
rm -f mietest.lay mietest.sca mietest.out mietest.out.check
#   Output the layer file
cat >mietest.lay <<EOF
  1.00   0.00    0.000000  'mietest.sca'
  0.00   0.00    0.000000  '           '
EOF
#
#   Output the scattering file
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 and output at 0, 90, and 180 degrees azimuth
rt3 <<EOF  >&! rt3.log
4              Nstokes
8              Nquad
G 
8              Nazimuth
mietest.lay
N              no delta-M scaling
1              solar source
0.628318531    solar flux
78.46304097    solar angle
0.0            ground temp
L
0.10           surface albedo
0.0            sky temp
0.951
W
IQ
2           number of output levels
1 2         output level numbers
3           number of output azimuths
mietest.out
EOF
#
#   Compare output with results below
cat >mietest.out.check <<EOF
C  NUMMU=  8  NUMAZI= 17  AZIORDER=  8  NSTOKES=4
C  LAYER_FILE=mietest.lay                        DELTA-M=N
C  SRC_CODE=1   QUAD_TYPE=GAUSSIAN        
C  DIRECT_FLUX= .62832E+00   DIRECT_MU= .200000
C  GROUND_TEMP=     .00   GROUND_TYPE=LAMBERTIAN      
C  GROUND_ALBEDO=  .10000   SKY_TEMP=     .00
C  WAVELENGTH= .951000E+00
C  UNITS=WATTS/(M^2 MICRON STER)     OUTPUT_POLARIZATION=IQ
C    Z      PHI     MU    FLUX/RADIANCE (I,Q,U,V)
   1.000    .0 -2.00000   .365133E+00   .489285E-01   .000000E+00   .000000E+00
   1.000    .0  2.00000   .628319E+00   .000000E+00   .000000E+00   .000000E+00
   1.000    .0  -.98940   .611250E-01  -.254753E-01   .000000E+00   .000000E+00
   1.000    .0  -.94458   .759512E-01  -.301313E-01   .000000E+00   .000000E+00
   1.000    .0  -.86563   .101371E+00  -.337486E-01   .000000E+00   .000000E+00
   1.000    .0  -.75540   .142650E+00  -.352959E-01   .000000E+00   .000000E+00
   1.000    .0  -.61788   .208195E+00  -.334187E-01   .000000E+00   .000000E+00
   1.000    .0  -.45802   .312367E+00  -.266875E-01   .000000E+00   .000000E+00
   1.000    .0  -.28160   .484328E+00  -.143500E-01   .000000E+00   .000000E+00
   1.000    .0  -.09501   .807490E+00   .153739E-03   .000000E+00   .000000E+00
   1.000    .0   .09501   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000    .0   .28160   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000    .0   .45802   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000    .0   .61788   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000    .0   .75540   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000    .0   .86563   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000    .0   .94458   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000    .0   .98940   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000  90.0  -.98940   .553862E-01   .219558E-01   .181689E-02   .204991E-04
   1.000  90.0  -.94458   .585238E-01   .237238E-01   .429232E-02   .441628E-04
   1.000  90.0  -.86563   .645430E-01   .271399E-01   .708604E-02   .611213E-04
   1.000  90.0  -.75540   .741562E-01   .326561E-01   .104047E-01   .674601E-04
   1.000  90.0  -.61788   .885759E-01   .410670E-01   .145388E-01   .595037E-04
   1.000  90.0  -.45802   .109780E+00   .537914E-01   .199176E-01   .334282E-04
   1.000  90.0  -.28160   .140755E+00   .735436E-01   .272017E-01  -.161735E-04
   1.000  90.0  -.09501   .184828E+00   .106090E+00   .375727E-01  -.100117E-03
   1.000  90.0   .09501   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000  90.0   .28160   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000  90.0   .45802   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000  90.0   .61788   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000  90.0   .75540   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000  90.0   .86563   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000  90.0   .94458   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000  90.0   .98940   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000 180.0  -.98940   .514661E-01  -.175194E-01   .353333E-08  -.241748E-11
   1.000 180.0  -.94458   .511342E-01  -.122200E-01   .334104E-08  -.702084E-11
   1.000 180.0  -.86563   .544662E-01  -.681647E-02   .323669E-08  -.125962E-10
   1.000 180.0  -.75540   .612331E-01  -.129271E-02   .316845E-08  -.180448E-10
   1.000 180.0  -.61788   .716076E-01   .436260E-02   .306566E-08  -.218610E-10
   1.000 180.0  -.45802   .861062E-01   .100360E-01   .283014E-08  -.217545E-10
   1.000 180.0  -.28160   .105191E+00   .151938E-01   .232724E-08  -.133365E-10
   1.000 180.0  -.09501   .127362E+00   .171386E-01   .139210E-08   .139267E-10
   1.000 180.0   .09501   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000 180.0   .28160   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000 180.0   .45802   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000 180.0   .61788   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000 180.0   .75540   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000 180.0   .86563   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000 180.0   .94458   .000000E+00   .000000E+00   .000000E+00   .000000E+00
   1.000 180.0   .98940   .000000E+00   .000000E+00   .000000E+00   .000000E+00
    .000    .0 -2.00000   .274800E-01   .000000E+00   .000000E+00   .000000E+00
    .000    .0  2.00000   .273970E+00   .317170E-01   .000000E+00   .000000E+00
    .000    .0  -.98940   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.94458   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.86563   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.75540   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.61788   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.45802   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.28160   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000    .0  -.09501   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000    .0   .09501   .934986E-01   .567385E-02   .000000E+00   .000000E+00
    .000    .0   .28160   .161225E+00   .969435E-02   .000000E+00   .000000E+00
    .000    .0   .45802   .184881E+00   .668130E-02   .000000E+00   .000000E+00
    .000    .0   .61788   .171583E+00  -.814586E-03   .000000E+00   .000000E+00
    .000    .0   .75540   .145162E+00  -.920166E-02   .000000E+00   .000000E+00
    .000    .0   .86563   .117483E+00  -.162311E-01   .000000E+00   .000000E+00
    .000    .0   .94458   .936345E-01  -.209165E-01   .000000E+00   .000000E+00
    .000    .0   .98940   .753632E-01  -.231434E-01   .000000E+00   .000000E+00
    .000  90.0  -.98940   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.94458   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.86563   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.75540   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.61788   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.45802   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.28160   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000  90.0  -.09501   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000  90.0   .09501   .561758E-01   .155300E-01   .949841E-02   .140749E-04
    .000  90.0   .28160   .764336E-01   .267527E-01   .127407E-01  -.974867E-05
    .000  90.0   .45802   .815112E-01   .307685E-01   .128460E-01  -.168168E-04
    .000  90.0   .61788   .788949E-01   .301894E-01   .111981E-01  -.116161E-04
    .000  90.0   .75540   .744413E-01   .282078E-01   .896557E-02  -.330049E-05
    .000  90.0   .86563   .702974E-01   .261542E-01   .657858E-02   .281815E-05
    .000  90.0   .94458   .672162E-01   .245569E-01   .418008E-02   .474342E-05
    .000  90.0   .98940   .654567E-01   .236249E-01   .181520E-02   .287773E-05
    .000 180.0  -.98940   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.94458   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.86563   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.75540   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.61788   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.45802   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.28160   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000 180.0  -.09501   .872073E-02   .000000E+00   .000000E+00   .000000E+00
    .000 180.0   .09501   .488879E-01   .841337E-03  -.132059E-08  -.213737E-10
    .000 180.0   .28160   .600969E-01   .584603E-03  -.215189E-08  -.129592E-10
    .000 180.0   .45802   .597792E-01  -.317858E-02  -.295357E-08  -.429956E-11
    .000 180.0   .61788   .555423E-01  -.804907E-02  -.352726E-08   .214458E-12
    .000 180.0   .75540   .518829E-01  -.128188E-01  -.387878E-08   .153983E-11
    .000 180.0   .86563   .504972E-01  -.170131E-01  -.405942E-08   .118486E-11
    .000 180.0   .94458   .522332E-01  -.203676E-01  -.412493E-08   .381463E-12
    .000 180.0   .98940   .577439E-01  -.226181E-01  -.412231E-08  -.802475E-13
EOF

echo "Differences in output listed below:"
diff mietest.out.check mietest.out | more

