SBDART模式aerosol.dat文件

爱说笑LQT

想求助一下大家有谁用SBDART模式计算过气溶胶的辐射强迫呀，想问一下aerosol.dat文件要怎么去写呀，我怎么写结果都不对

trz

楼主你输入条件中，像函数那一项用到的legendre moments，怎么给的？我之前也想用来着，但是不会给这一项就没用

  IAER:         Boundary layer aerosol typeselector

               -1-read aerosol opticaldepth and scattering parameters
                  from aerosol.dat. Seesubroutine AEREAD.
               
                  the file format isreadable by the following Fortran code:

                  read(11,*) nn,moma
                  do k=1,huge(0)
                    read(11,*,end=100)wl(k)
                    do i=nz-nn+1,nz
                      read(11,*)dtau(i,k),waer(i,k),pmom(1:moma,i,k)
                    enddo
                  enddo
              100 continue

爱说笑LQT

trz 发表于 2021-11-13 15:03
楼主你输入条件中，像函数那一项用到的legendre moments，怎么给的？我之前也想用来着，但是不会给这一项就 ...

我设置为0，但是我输出的结果全是NAN，不知道这是为什么

trz

爱说笑LQT 发表于 2021-11-14 19:33
我设置为0，但是我输出的结果全是NAN，不知道这是为什么

我觉得是这一块的问题，legendre moments 表示的是像函数形状，可以用不对称因子g结合Henyey-Greenstein近似求出来（我不会），设为0可能不太对，说明好像有个例子                    
nstr=6
                     wlbaer=.400,.500,.600,.700
                     wbaer =0.88,0.89,0.91,0.92
                     pmaer= 0.80,0.70,0.60,0.50,
                            0.64,0.49,0.36,0.25,
                            0.51,0.34,0.22,0.12,
                            0.41,0.24,0.13,0.06,
                            0.33,0.17,0.08,0.03,
                            0.26,0.12,0.05,0.02
楼主学会了教教我

kkkppd

请问您知道怎么定义aerosol输入了吗，能否请教一下