大气净辐射通量计算

一颗梅子在南京

有谁算过大气的净辐射通量值吗？我找到的计算方法如下图：
file:///C:/Users/cgm/AppData/Roaming/Tencent/Users/811090821/QQ/WinTemp/RichOle/GSUF1%601_MU_WWZ$64%7D$1CQ4.jpg
然后我编写的fortran程序是：
do i=1,12
      j=1
   BD1(i,j)=BD(i,j)*ltrans(i,j)
   do j=2,ny
      BD1(i,j)=(BD1(i,j-1)+BD(i,j))*ltrans(i,j)
   enddo
   enddo
   do j=1,ny
      do i=1,12
      Flup(j)=Flup(j)+BDG(i)*ltrans1(i,j)+BD1(i,j)
   enddo
   enddo
   do i=1,12
      ffdown(i)=0
   enddo
   do i=1,12
      j=ny
      BD2(i,j)=BD(i,j)*ltrans(i,j)
   do j=ny-1,1,-1
      BD2(i,j)=(BD2(i,j+1)+BD(i,j))*ltrans(i,j)
   enddo
   enddo
   do j=ny,1,-1
      do i=1,12
      Fldown(j)=Fldown(j)+ffdown(i)*ltrans(i,j)/ltrans1(i,j)+BD2(i,j)
   enddo
   enddo
   do i=1,ny
      Fn(i)=Fldown(i)-Flup(i)
   enddo
   do i=1,ny
      write(*,*)Fn(i)
   enddo
其中：i为波段数，j为垂直分层,ltrans(i,j)为每一层i波段的透过率，ltrans1(i,j)为波段i从底层到j层的透过率，其余量都与公式中的量相对应！
但是净辐射通量的结果Fn(i)=Fldown(i)-Flup(i)有点问题
file:///C:/Users/cgm/AppData/Roaming/Tencent/Users/811090821/QQ/WinTemp/RichOle/7FYAZOSAWO%60]P%25YY9P_PR5O.jpg
理论上底层应该是辐射降温的，净辐射通量值不应该是正值的！
不知道问题出在哪了，求牛人们帮忙解决！谢谢

一颗梅子在南京

图片没显示，重新贴上：

没有放弃过

浮躁  还是浮躁

一颗梅子在南京

没有放弃过 发表于 2013-3-7 22:24

                               
登录/注册后可看大图

浮躁  还是浮躁

这不是一个让你吐槽的帖子，想吐槽的请绕道，谢谢