计算大气热源Q1遇到的问题

LiuJiawei

现在网上流传的（http://bbs.06climate.com/forum.php?mod=viewthread&tid=18894）计算大气热源的程序基本是同一个。
师兄师姐们使用之后觉得量级差不多，但也没有细看程序。
我用NCL改编的时候发现了两个问题：
1.
平流计算中
real   ,parameter :: RAD=PI/180.0,DD=RR*RAD*2.5,
dx=cos((-87.5+2.5*(j-1))*RAD)
dy=2*dd

dx明显连单位都没有
我改为：dx=cos((-87.5+2.5*(j-1))*RAD)*RR*RAD*2.5*2（非边界处）
2.
REAL    :: PP(NL)=(/1000., 925., 850., 700., 600., 500., 400., 300., 250., 200., 150., 100.,70.,50.,30.,20.,10./)
pp的单位是hPa
而计算位温垂直输送时
dp=pp(nl)-pp(nl-1)
没有化成Pa，与上下程序不一致。
我改为：dp=pp(nl)-pp(nl-1)*100
但是，计算完成作图以后发现与国外文献比较，计算值大约偏大3倍。

不知道大家还有没有发现这个程序的问题，或者我哪里改的不对？

zongqing123

楼主可以把修改后的程序上传吗？谢谢，菜鸟一枚，希望从正确的程序开始学起

LiuJiawei

xqxsummer 发表于 2016-1-25 16:37
看了很多中文 单位都不一样呀 请问视热源英文是 apparent heating rate吗

常用的一般是整层积分的<Q1>，单位是w/m^2
和加热率Q1/cp，单位是K/s或者换算为K/d

建议你看以下两篇文献：
Yanai M, Tomita T. Seasonal and Interannual Variability of Atmospheric Heat Sources and Moisture Sinks as Determined from NCEP-NCAR Reanalysis[J]. Journal of Climate, 1998, 11(3):463-482.
赵煊, 徐海明, 徐蜜蜜,等. 春季中国东海黑潮区大气热源异常对中国东部降水的影响[J]. 气象学报, 2015, 73(2):263-275.

freekiller

贴图贴脚本      

LiuJiawei

freekiller 发表于 2015-12-10 08:42
贴图贴脚本

文献中该地区1980-1994年冬季平均的高值也就200W+，低值-100w左右。

freekiller

我看了一下，对于dx的计算我是赞同你的修改，cos函数需要转换为弧度后，还需要计算周长

freekiller

你再仔细看一下，期待你的结果

freekiller

另外ncl 中有函数计算任意两点间的曲面距离，你可以方便地计算不同纬度上的dx

LiuJiawei

freekiller 发表于 2015-12-10 10:27
另外ncl 中有函数计算任意两点间的曲面距离，你可以方便地计算不同纬度上的dx

好的，谢谢。这个程序是从FORTRAN的程序改写的，所以显得有些生硬。NCL中差分好像也有函数。我再改改。

LiuJiawei

我的程序中，还存在一处错误，是我改写时犯的错误，原程序中没有，求垂直输送项时时（p/p0）^ck与位温计算时的（p0/p）^ck不同。计算出来之后单独年的高值还是很高，但是经过多年平均的气候态时高值会回到参考文献相当的水平。

有需要的朋友，上面的程序请务必修改正确以后使用！

luzirs

问渠那得清如许，为有源头活水来

LiuJiawei

突然发现贴出版本的程序中spsp地表气压的读入时次也需要修改。
另外，多年平均之后的高值虽然快速降低，但是导师还是认为偏高。
如果发现其它问题，希望能够告知