荒城旧梦5. 对每个波长的记录。输出格式：
write(*,*) '\用于后处理程序)
do m=1,nw
write(*,*)
& wl,ffv,topdn,topup,topdir,botdn,botup,botdir
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
其中，
WL =波长 （微米） FFV =过滤函数值
TOPDN =在ZOUT(2)的总向下流量km （w/m2/micron） TOPUP =在ZOUT(2)的总向上流量km （w/m2/micron） TOPDIR=在ZOUT(2)的直接向下流量km （w/m2/micron） BOTDN =在ZOUT(1)的总向下流量km （w/m2/micron） BOTUP =在ZOUT(1)的总向上流量km （w/m2/micron） BOTDIR=在ZOUT(1)的直接向下流量km （w/m2/micron） NPHI =用户的方位角数 NZEN =用户的天顶角数
PHI =用户指定的方位角数 （degrees） UZEN =用户指定的天顶角数 （degrees） VAEN =用户指定的天底角 （degrees） UURS =海拔ZOUT(2)（顶部) (w/m2/um/str)
用户角度的辐射
说明：由SBDART输出的辐射代表了散射辐射。其中不包括太阳直射。同样，记住
UURS代表了用户指定样本方向的辐射。因此，通过UURS的角度整合计算出的辐照度将由于直接光束的缺少而不会产生BOTDN，并且它将可能由于取样不会产生（BOTDN-BOTDIR）。
说明：如果设置了IDAY，那么PHI将是辐射传输的指示方向
6. 除了辐射是在ZOUT(1)海拔高度（底部）之外，和IOUT=5一样
7. 每个波长在每一层的辐射通量。该输出项可以产出大量的输出结果如果使用了很
多的样本点
write(*,*) '\用于后处理程序)
write(*,*) nz ; z水平的数量
26 / 32
write(*,*) nw ; 波长的数值
do j=1,nw
write(*,*) wl write(*,*)
& (Z(i),i=nz,1,-1), ; altitude (km)
& (fdird(i),i=1,nz), ; downward direct flux (w/m2/um) & (fdifd(i),i=1,nz), ; downward diffuse flux (w/m2/um) & (flxdn(i),i=1,nz), ; total downward flux (w/m2/um) & (flxup(i),i=1,nz) ; total upward flux (w/m2/um) Enddo
10. 每运行一次输出一个记录，通过波长积分得到。 输出量是，（采用梯形积分）
WLINF,WLSUP,FFEW,TOPDN,TOPUP,TOPDIR,BOTDN,BOTUP,BOTDIR
WLINF = 较低波长限制 (microns) WLSUP = 较高波长限制 (microns) FFEW = 滤波函数等效宽度 (microns) TOPDN =在ZOUT(2)的总向下流量km (w/m2) TOPUP =在ZOUT(2)的总向上流量km (w/m2) TOPDIR=在ZOUT(2)的直接向下流量km (w/m2) BOTDN =在ZOUT(1)的总向下流量km (w/m2) BOTUP =在ZOUT(1)的总向上流量km (w/m2) BOTDIR=在ZOUT(1)的直接向下流量km (w/m2)
说明：为了得到光谱通量密度（w/m2/micron），将FFEW中的任何这些量分开
处理
11. 通过波长积分得到的每个大气层的辐射通量。输出格式：
write(*,*) nz,phidw
do i=1,nz
write(*,*) zz,pp,fxdn(i),fxup(i),fxdir(i),dfdz,heat endd
其中，nz = 大气层数
phidw = 滤波函数等效宽度 (um) zz = 水平海拔高度 (km) pp = 水平压力 (mb) fxdn = 向下通量(直接+传播) (W/m2) fxup = 向上通量 (W/m2) fxdir = 向下通量, 只有直接 (W/m2) dfdz = 辐射能通量散度 (mW/m3) heat = 加热速率 (K/day)
27 / 32
说明：dfdz(i)和heat(i)定义在每层中心，即，在i-1层和i层之间。
20. 在ZOUT(2)公里的辐射输出
输出格式：
write(*,*) wlinf,wlsup,ffew,topdn,topup,topdir,
& botdn,botup,botdir write(*,*) nphi,nzen
write(*,*) (phi(i),i=1,nphi) write(*,*) (uzen(j),j=1,nzen)
write(*,*) ((r(i,j),i=1,nphi),j=1,nzen)
输出的第一项记录与格式IOU=10相同
(WLINF,WLSUP,FFEW,TOPDN,TOPUP,TOPDIR,BOTDN,BOTUP,BOTDIR) 额外的记录包括：
NPHI = 用户的方位角数
NZEN = 用户的天顶角数
PHI = 用户相对方位角 (nphi 数值) UZEN = 用户天顶角 (nzen 数值)
R = 辐射阵列nphi,nzen) (W/m2/sr)
说明：如果设置了IDAY，那么PHI将是辐射传输的指示方向
21. 除了辐射是在ZOUT(1)km海拔高度之外，和IOUT=20一样
22. 通过对波长积分得到的每个大气层的辐射和通量。 输出格式：
write(*,*) nphi,nzen,nz,ffew
write(*,*) (phi(i),i=1,nphi) write(*,*) (uzen(j),j=1,nzen) write(*,*) (z(k),k=nz,1,-1) write(*,*) (fxdn(k),k=1,nz) write(*,*) (fxup(k),k=1,nz) write(*,*) (fxdir(k),k=1,nz)
write(*,*) (((uurl(i,j,k),i=1,nphi),j=1,nzen),k=1,nz)
其中，nphi = 用户指定的方位角数
nzen = 用户指定的天顶角数 nz = 大气层数
ffew = 滤波函数等效宽度 (um)
phi = 用户指定方位角 (degrees) uzen = 用户指定天顶角 (degrees) z = 大气层海拔高度 (km)
28 / 32
fxdn = 向下通量 (直接+传输) (W/m2) fxup = 向上通量 (W/m2) fxdir =向下通量, 只有直接 (W/m2) uurl = 每层辐射 (W/m2/str)
说明：如果设置了IDAY，那么PHI将是辐射传输的指示方向
23. 和IOUT=20一致除了
对应于ZOUT(1)的下半球辐射输出 对应于ZOUT(2)的上半球辐射输出
使用此输出格式来确定在散射层上方和下方的辐射。例如，如果ZCLOUD=1
并且TCLOUD=10，你可以通过设置IOUT=23, ZOUT=1,2来得到在云的上方和下方区域的散射辐射
说明：如果设置了IDAY，那么PHI将是辐射传输的指示方向
===============================================================================
离散纵坐标选项 ==============
DELTAM：如果设置为真，那么使用delta-m方法（参见Wiscombe, 1977）。 该方法本质上是用于多个辐射流的delta-Eddingto近似值。
一般来说，对于一个给定的数流，强度和通量将会更准确，对于有很大峰值的相函
数如果“DELTAM”设置为真。在前向散射方向上10度左右范围内的强度经常不是非常准确，然而，因此，当主要的兴趣在这个所谓的“光环区”时，应设置DELTAM为假。（默认值=真）
NSTR： 用于计算天顶角的数量。NSTR必须能被2整除。使用NSTR=4将减少计算通量需要的
时间，和NSTR=16相比以因子5下降，其对于准确度几乎没有影响（大约0.5%的改变当设置DELTAM为真）。辐射的计算需要更多的流束。NSTR的默认值取决于IOUT的值。计算通量使用的默认值（IOUT=1，7，10，11）的NSTR=4。计算辐射使用的默认值（IOUT=5，6，20，21，22，23）的NSTR=20。
CORINT：当设置为真，对于delta-M的尺度效应，强度均正确（参见Nakajima and Tanaka, 1988）。
当为假，强度均不正确。
一般来说，如果束流源存在，DELTAM为真，并且待解决问题中包含散射，那么
CORINT应设置为真。然而，CORINT=FALSE程序运行得更快，并且会在光环区外产生相当正确的强度结果。如果CORINT=TRUE并且相函数的矩被PMAER(IAER=5)指定，或者在aerosol.dat (IAER=-1)中设置，那么指定足够数量的相函数矩来适当地模拟前向单次散射辐射就变得重要。否则，如果提供的矩太少，强度可能更准确当CORINT=FALSE。
默认值CORINT=.false.
29 / 32
CORINT的输入值可以忽略当：
1）在辐射模式中，即ioutb不等于（5,6,20,21,22） 2）没有束流源（FBEAM = 0），或者 3）没有散射（对于所有层SSALB=0）
辐射输出 ========
NZEN; 用户的天顶角数。如果该参数被指定，SBDART将输出在NZEN天顶角处的辐射量，
均匀地间隔输入数组UZEN的前两个值。例如，
nzen=9,
uzen=0,80 将产生在天顶角0,10,20,30,40,50,60,70,80处的输出。
UZEN： 用户天顶角。如果NZEN被指定，那么UZEN被解释为天顶角范围的限制，并且只要
求前两个元素满足。如果NZEN没有被设定，那么最多到NSTR的UZEN的值可以被指定。如果NZEN和UZEN均未被指定，并且要求计算辐射量（IOUT=5,6,20,21,22,23），那么使用天顶角的默认设置，其取决于如下的IOUT值：
* IOUT=5 or 20: NZEN=18, UZEN=0,85
* IOUT=6 or 21: NZEN=18, UZEN=95,180 * IOUT=22 or 23: NZEN=18, UZEN=0,180
说明：UZEN指定辐射传输处的天顶角：
UZEN = 0 => 直接辐射传输
UZEN < 90 =>上半球辐射 UZEN > 90 => 下半球辐射
UZEN = 180 => 直接向下辐射传输
VZEN： 用户天底角。这只是指定用户辐射角度方向的一种替代方法，借用uzen=180-vzen
NPHI： 用户方位角数。如果该参数被指定，SBDART将输出在NPHI方位角处的辐射量，均
匀地间隔输入数组NPHI的前两个值。例如，nphi=7，phi=0,180将会产生在天顶角0,30,60,90,120,150,180处的输出
PHI： 用户相对方位角。如果NPHI被指定，那么PHI被解释为方位角范围的限制，并且只
要求前两个元素满足。如果NPHI没有被设定，那么最多到NSTR的PHI的值可以被指定。如果NPHI和PHI均未被指定，并且要求计算辐射量（IOUT=5,6,20,21,22,23），那么使用方位角的默认设置，等价于NPHI=19, PHI=0,180的情况。
说明：如果设置了SAZA，那么PHI=0就代表往正北方向传输的辐射。PHI按照俯视地
球表面的顺时针方向增加。
如果为设置SAZA，那么PHI是前向散射方向的相对方位角。
30 / 32
