荒城旧梦
* PHI < 90 => 向前散射辐射
* PHI > 90 => 向后散射辐射
例如，如果设定在西边，传输到东南方向的辐射有45度的相对方位角。
说明：SBDART目前被设置成模拟最多40个计算的天顶角和40个方位角的辐射。虽
然这些限制可能会扩大，要知道，运行很大数量的SBDART模式将大大增加运行时间和内存要求。在对DEC Alpha进行的测试中，运行时间大约与NSTR^2一致，对于NSTR小于40。代码的内存使用情况也大约与NSTR^2一致。
说明：太阳方位角的默认值，SAZA=180。该SAZA的值将导致前向散射太阳辐射在
PHI=0附近出现。
辐射边界条件 ============
IBCND：
=0：一般情况下：边界层是下列情况的任何组合： * 来自顶部的光束 (参见FBEAM)
* 来自顶部的各向同性光束 (参见FISOT)
* 来自顶部的热辐射 (参见TEMIS,TTEMP) * 内部热辐射源 (参见TEMPER)
* 底部的反射 (参见LAMBER,ALBEDO,HL) * 来自底部的热辐射 (参见BTEMP)
=1：从顶部到底部的各向同性光束，为了得到反照率和整个介质与入射光束角的透
射率；
在该情况下，唯一的输入变量是NLYR, DTAUC, SSALB, PMOM, NSTR, USRANG,
NUMU, UMU, ALBEDO ,DELTAM, PRNT, HEADER, 和数组维数。
NOPLNK,LAMBER都被假设为真 ，边界层底部可以有任何反照率值。唯一的输
出是ALBMED,TRNMED。在这种情况下，UMU被解释为光束的角度阵列。
如果USRANG = TRUE，它们必须为正和单调递增，并且通过以下方式返回：然
而，在内部，对于UMU的否定在增加，因此MAXUMU必须至少为2*NUMU。
如果USRANG = FALSE，UMU作为NSTR/2的角度余弦值求积分的方式返回，并且
单调递增。
FISOT： 上边界各向同性光束强度。（如果热辐射为正，单位是w/sq m，否则是任意单位）。
相应的入社通量是pi（3.14159...）时间‘FISOT’。
TEMIS： 顶层的发射率，如果NOTHRM=1不使用。
BTEMP：以开尔文为单位的表面温度。如果未设置，那么表面温度降被设置为模式大气的底
层的温度。表面反射率是通过反照率计算得到（参见ISALB和ALBCON）。当NOTHRM=1，BTEMP的输入值忽略。当SPOWDER=.true.，BTEMP指定了仅低于颗粒层水平的表面
31 / 32
温度。
离散纵标法特定的输出选项： ==========================
PRNT(k)：k
1 输入变量（除了PMOM） 2 通量
3 用户水平和角度的强度
4 作为应变量太阳天顶角的平面透射率和平面反照率（IBCND = 1） 5 用于每一层的相函数矩PMOM
（仅当PRNT(1) = TRUE，并且只用于有散射作业的云层）
===============================================================================
32 / 32
