User Guide SBDART Input
User Guide SBDART Input
User Guide SBDART Input
This file documents input parameters for SBDART, (Santa Barbara DISORT
Atmospheric Radiative Transfer). SBDART is a software tool that computes
plane-parallel radiative transfer in clear and cloudy conditions within the
Earth's atmosphere and at the surface. For a general description and review
of the program please refer to Ricchiazzi et al 1998. (Bulletin of the
American Meteorological Society, October 1998).
SBDART's main input file is called INPUT. This file contains a single
NAMELIST input block also named INPUT. A significant advantage of NAMELIST
input is that not all elements of an input block need be specified by the
user. Since most of the code inputs have been initialized with reasonable
default values, new users can start by specifying just a few interesting
input parameters. The default state of input parameters may be determined
by removing INPUT from the current working directory. When SBDART detects
the absence of file INPUT, it will print the default settings of all input
parameters. This output may be redirected to a file for editing.
The default configuration of INPUT is as follows:
==========================================================================
&INPUT
idatm
wlinf
sza
nf
alat
pbar
uo3
xrsc
xco2
xco
xnh3
xo4
sc
lwp
krhclr
taerst
rhaer
abaer
pmaer
nothrm
zgrid1
zout
lamber
prnt
temis
uzen
phi
ttemp
idb
/
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
4
0.550
0.0
2
-64.7670
-1.0
-1.0
1.0
-1.0
-1.0
-1.0
1.0
1.0,3*0.0
5*0.0
0
5*0.0
-1.0
-1.0
940*0.0
-1
0.0
0.0,100.0
t
7*f
0.0
20*-1.0
20*-1.0
-1.0
20*0
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
amix
wlsup
csza
iday
alon
sclh2o
o3trp
xn2
xch4
xno2
xno
isalb
zcloud
nre
jaer
iaer
wlbaer
wbaer
zbaer
nosct
zgrid2
iout
ibcnd
ipth
nstr
vzen
imomc
btemp
= 0.0
= 0.550
= -1.0
= 0
= -64.0670
= -1.0
= -1.0
= -1.0
= -1.0
= -1.0
= -1.0
= 0
= 5*0.0
= 5*8.0
= 5*0
= 0
= 47*0.0
= 47*0.950
= 50*-1.0
= 0
= 30.0
= 10
= 0
= 1
= 4
= 20*90
= 3
= -1.0
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
isat
wlinc
solfac
time
zpres
uw
ztrp
xo2
xn2o
xso2
xhno3
albcon
tcloud
rhcld
zaer
vis
tbaer
gbaer
dbaer
kdist
ngrid
deltam
saza
fisot
nzen
nphi
imoma
spowder
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
0
0.0
1.0
16.0
-1.0
-1.0
0.0
-1.0
-1.0
-1.0
-1.0
0.0
5*0.0
-1.0
5*0.0
23.0
47*0.0
47*0.70
50*-1.0
3
50
t
180.0
0.0
0
0
3
f
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
=======================================================================
-------
atmospheric profile
aerosol information
spectral surface albedo
sensor filter function
solar spectrum
cloud vertical profile
(see
(see
(see
(see
(see
(see
input
input
input
input
input
input
quantity
quantity
quantity
quantity
quantity
quantity
IDATM)
IAER)
ISALB)
ISAT)
NF)
TCLOUD)
NF:
100
read(13,*,end=100) (wlsun(i),sun(i),i=1,5000)
continue
where, wlsun
sun
.25 to 4 micron
ISAT:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
(default)
the values of
effect if they
default values
Otherwise the
as specified in
METEO
GOES(EAST)
GOES(WEST)
AVHRR1(NOAA8)
AVHRR2(NOAA8)
AVHRR1(NOAA9)
AVHRR2(NOAA9)
AVHRR1(NOAA10)
AVHRR2(NOAA10)
AVHRR1(NOAA11)
AVHRR2(NOAA11)
GTR-100 ch1
GTR-100 ch2
GTR-100 410nm channel
GTR-100 936nm channel
MFRSR 415nm channel
MFRSR 500nm channel
MFRSR 610nm channel
MFRSR 665nm channel
MFRSR 862nm channel
MFRSR 940nm channel
AVHRR3 (nominal)
AVHRR4 (nominal)
AVHRR5 (nominal)
Biological action spectra for DNA damage by UVB radiation
AIRS1 380-460nm
AIRS2 520-700nm
AIRS3 670-975nm
AIRS4 415-1110nm
filt
WLINF:
WLSUP:
microns)
NOTE:
If ISAT eq -2, a rectangular filter (constant with
wavelength) is used, with central wavelength at WLINF
and an equivalent width of WLSUP (full width = WLSUP)
If ISAT eq -3, a triangular filter function is used with
the central wavelength at WLINF and an equivalent width of
WLSUP (full width = 2*WLSUP) (filter function is zero at end
points, and one at WLINF).
If ISAT eq -4, a gaussian filter function is used with
the central wavelength at WLINF and an equivalent width of
WLSUP (full width = 4*WLSUP)
If output is desired at a single wavelength, set
WLINF=WLSUP and ISAT=0. In this case, SBDART will
set WLINC=1 (the user specified value of WLINC is ignored)
and the output will be in units of (W/m2/um) for
irradiance and (W/m2/um/sr) for radiance.
WLINC:
ever is smaller.
wlinc < 0
1 >= wlinc > 0
wlinc > 1
SOLAR GEOMETRY
==============
SZA:
CSZA:
SAZA:
IDAY:
SOLFAC:
(jan 2)
4 -sea water
5 -sand
6 -vegetation
ALBCON:
ISALB=0 --
SC:
SC(2): Wind speed over ocean in m/s. Wind speed affects the
sunglint and foam contributions to the ocean brdf
(Default=5 m/s)
SC(3): Oceanic salinity in parts per thousand. Salinity affects
the Fresnel reflectivity of the ocean surface (Default=34.3)
NOTE: The oceanic brdf model is nearly the same as the
one used in the 6s radiative transfer model. The only
difference is that the brdf in SBDART does not depend
on the wind direction.
The brdf dependence on wind
direction was eliminated by averaging the 6s reflectance
model over 360 degrees of wind azimuth.
Note: The 6s model does not include any provisions for
wave facet shading, so results for shallow
NOTE:
---------- isalb=9 ---------SC specifies the Ross-thick Li-sparse brdf model parameters
SC(1)=isotropic coefficient
SC(2)=volumetric coefficient
SC(3)=geometric shadowing coefficient
SC(4)=hot spot magnitude
SC(5)=hot spot width
(Ross, 1981; Li and Strahler, 1992).
---------- isalb=10 ---------SC specifies mixing ratios snow, ocean, sand and vegetation
SC(1):
SC(2):
SC(3):
SC(4):
NOTE:
10
MODEL ATMOSPHERES
=================
IDATM:
ATMOSPHERIC PROFILE:
-----------------0 User Specified
1 TROPICAL
2 MID-LATITUDE SUMMER
3 MID-LATITUDE WINTER
4 SUB-ARCTIC SUMMER
5 SUB-ARCTIC WINTER
6 US62
-n List to standard out
default
ozone(atm-cm)
water vapor (g/cm2) total below_10km
------------------- ------------4.117
2.924
0.854
2.085
0.418
1.418
0.253
0.324
0.403
0.350
0.486
0.349
.0216
.0325
.0336
.0346
.0340
.0252
10
read(13,*) nn
do 10 i=1,nn
read(13,*) z(i),p(i),t(i),wh(i),wo(i)
continue
where nn is the number atmospheric layers
nn should be less or equal to than MXLY,
a parameter used in SBDART (see params.f)
to set the maximum number of levels in the
vertical grids.
z
11
UW:
UO3:
O3TRP:
ZTRP:
XN2:
XO2:
XCO2:
XCH4:
XN2O:
XCO:
XNH3:
XSO2:
XNO:
XHNO3:
XNO2:
volume
volume
volume
volume
volume
volume
volume
volume
volume
volume
volume
mixing
mixing
mixing
mixing
mixing
mixing
mixing
mixing
mixing
mixing
mixing
ratio
ratio
ratio
ratio
ratio
ratio
ratio
ratio
ratio
ratio
ratio
of
of
of
of
of
of
of
of
of
of
of
N2
O2
CO2
CH4
N2O
CO
NH3
SO2
NO
HNO3
NO2
(PPM,
(PPM,
(PPM,
(PPM,
(PPM,
(PPM,
(PPM,
(PPM,
(PPM,
(PPM,
(PPM,
default
default
default
default
default
default
default
default
default
default
default
= 781000.00 )
= 209000.00 )
=
360.00 )
=
1.74 )
=
0.32 )
=
0.15 )
=
5.0e-4)
=
3.0e-4)
=
3.0e-4)
=
5.0e-5)
=
2.3e-5)
12
XO4
PBAR:
ZPRES:
SCLH2O:
13
CLOUD PARAMETERS
================
ZCLOUD:
z(k)
...
2.5
14
5
4
3
2
1
2.0
1.5
1.0
0.5
0.0
6-9 km
4-6 km
1-4 km
total
cloud
optical
depth
6
24
30
60
opacity
(km-1)
2
12
10
scattering
parameters
mie(20)
[mie(10)*10+mie(20)*2]/12
mie(10)
15
For example,
ZCLOUD= 1 ,-5,0,0,0
TCLOUD= 10, 1,0,0,0
NRE=
10,20
#
#
#
#
#
#
#
#
#
#
#
#
#
#
16
NRE:
(up to 5 values)
re
(g/m2)
(um)
17
fwp
rei
cldfrac
#
#
#
#
#
IMOMC:
LWP:
18
tau =
3 Q(wl) * LWP
------------4 RHO * NRE
Let
tau = 30
RHCLD:
KRHCLR:
19
may drive the water vapor in clear layers to zero and still
be unable produce a given WVP. This parameter has no
effect when KDIST<0
ZAER:
TAERST:
moma
wl(k)
20
dtau(i,k)
waer(i,k)
21
tauaero(0.55um) = 3.912*(1.05*w+1.51*(1-w))/vis
where w is a weighting factor between the two
extremes and is given by
(1/vis-1/23)
w = ---------(1/5-1/23)
w = 1
vis < 5
w = 0
vis > 23
DBAER:
22
QBAER:
23
WLBAER:
WBAER:
GBAER:
PMAER:
/
pmaer(i) = | f(mu) P(i,mu) d mu
24
/
/
/
| f(mu) d mu
IMOMA:
SPOWDER:
isotropic scattering
rayleigh scattering
henyey_greenstein(g(re))
<-- default
haze L as specified by garcia/siewert
cloud c.1 as specified by garcia/siewert
spowder=t
iout=1
25
nre=
-100,
10
/
Similarly, the scattering properties of the surface and
atmosphere may be read from aerosol.dat with IAER=-1.
To model a semi-infinite granular surface layer the optical
depth of the bottom layer should be made very large, e.g.,
10000, as indicated in the example. However, a smaller
optical depth may also be specified in conjuntion with a
given value of sub-surface albedo selected with ISALB.
Thus, in the previous example the effect of a thin snow
layer covering a grass field may be modeled by setting,
tcloud=100,10 and isalb=6
NOTE: When SPOWDER is set, BTEMP no longer sets surface
skin temperature. Rather it sets the temperature of the
sub-surface level below the granular material. At present
there is no way to set the surface temperature (above the
graunual layer) to something other than the atmospheric
temperature of the bottom level.
=============================================================
NOTHRM:
nothrm=1
KDIST:
26
27
28
29
ZOUT:
IOUT:
value
----0.
1.
wavelength
filter function value
total downward flux at ZOUT(2) km
total upward flux at ZOUT(2) km
direct downward flux at ZOUT(2) km
total downward flux at ZOUT(1) km
total upward flux at ZOUT(1) km
direct downward flux at ZOUT(1) km
(microns)
(w/m2/micron)
(w/m2/micron)
(w/m2/micron)
(w/m2/micron)
(w/m2/micron)
(w/m2/micron)
5.
30
write(*,*) nphi,nzen
write(*,*) (phi(j),j=1,nphi)
write(*,*) (uzen(j),j=1,nzen)
do i=nzen,1,-1
write(*,*) (uurs(i,k),k=1,nphi)
enddo
enddo
where,
WL
=
FFV
=
TOPDN =
TOPUP =
TOPDIR=
BOTDN =
BOTUP =
BOTDIR=
NPHI =
NZEN =
PHI
=
UZEN =
VZEN =
UURS =
wavelength
filter function value
total downward flux at ZOUT(2) km
total upward flux at ZOUT(2) km
direct downward flux at ZOUT(2) km
total downward flux at ZOUT(1) km
total upward flux at ZOUT(1) km
direct downward flux at ZOUT(1) km
number of user azimuth angles
number of user zenith angles
user specified azimuth angles
user specified zenith angles
user specified nadir angles
radiance at user angles at
altitude ZOUT(2) (top)
(microns)
(w/m2/micron)
(w/m2/micron)
(w/m2/micron)
(w/m2/micron)
(w/m2/micron)
(w/m2/micron)
(degrees)
(degrees)
(degrees)
(w/m2/um/str)
7.
do j=1,nw
write(*,*) wl
write(*,*)
&
(Z(i),i=nz,1,-1),
&
(fdird(i),i=1,nz),
31
; altitude
; downward direct flux
(km)
(w/m2/um)
&
(fdifd(i),i=1,nz),
&
(flxdn(i),i=1,nz),
&
(flxup(i),i=1,nz)
enddo
10.
(microns)
(microns)
(microns)
(w/m2)
(w/m2)
(w/m2)
(w/m2)
(w/m2)
(w/m2)
nz
phidw
zz
pp
fxdn
fxup
fxdir
dfdz
heat
=
=
=
=
=
=
=
=
=
wlinf,wlsup,ffew,topdn,topup,topdir,
botdn,botup,botdir
nphi,nzen
(phi(i),i=1,nphi)
(uzen(j),j=1,nzen)
((r(i,j),i=1,nphi),j=1,nzen)
32
(WLINF,WLSUP,FFEW,TOPDN,TOPUP,TOPDIR,BOTDN,BOTUP,BOTDIR)
addition records contain:
NPHI
NZEN
PHI
UZEN
R
=
=
=
=
=
22.
nphi,nzen,nz,ffew
(phi(i),i=1,nphi)
(uzen(j),j=1,nzen)
(z(k),k=nz,1,-1)
(fxdn(k),k=1,nz)
(fxup(k),k=1,nz)
(fxdir(k),k=1,nz)
(((uurl(i,j,k),i=1,nphi),j=1,nzen),k=1,nz)
nphi
nzen
nz
ffew
phi
uzen
z
fxdn
fxup
fxdir
uurl
=
=
=
=
=
=
=
=
=
=
=
33
==========================================================================
DISORT options
==============
DELTAM:
NSTR:
CORINT:
34
NZEN:
nzen=9,
uzen=0,80
will cause output at zenith angles 0,10,20,30,40,50,60,70,80.
UZEN:
=
<
>
=
0
90
90
180
=>
=>
=>
=>
radiation
radiation
radiation
radiation
propagates directly up
in upper hemisphere
in lower hemisphere
propagates directly down
VZEN:
NPHI:
PHI:
35
* PHI < 90
* PHI > 90
36
TEMIS:
BTEMP:
PRNT(k):
=========================================================================
37