求CAPE的算法或gs

BDRUSH

我也想知道咧！同求好人分享

o_O_小岛

下载fnl资料···里面有cape计算值···可以少计算

sukeyangel

去下了一个。
!----------计算CAPE
total=0.0
do pp=lfcp,elp,-0.5
        do i=1,z-1
        if((pp-p(i))*(pp-p(i+1))<=0)then
        tt=(t(i)-t(i+1))*(pp-p(i+1))/(p(i)-p(i+1))+t(i+1)
        endif
        enddo
        call ttp(pp,thse,tt,tpp)
  tppv=(tpp+273.16)*(1+0.378*(6.11*10**(a*tpp/(b+tpp)))/pp)
  call tenvironment(pp,z,t,td,p,tce,tde)
  tenv=(tce+273.16)*(1+0.378*(6.11*10**(a*tde/(b+tde)))/pp)
  cape=(rd*(tppv-tenv)/pp)*0.5
  total=total+cape
end do
100 write(*,*)"total_cape=",total,"J/KG"
stop
end

这一段就够了吧

司空摘月

风之痕 发表于 2011-9-10 14:17
**** 作者被禁止或删除 内容自动屏蔽 ****

嘿嘿 我最近正好要算那个cape值，我也看了下plotskew.gs里面的量有很多。不知道你把有关cape的gs整理了没，能不能偷懒用下你的啊。谢谢

司空摘月

司空摘月 发表于 2014-9-9 15:21
嘿嘿 我最近正好要算那个cape值，我也看了下plotskew.gs里面的量有很多。不知道你把有关cape的gs整理了没 ...

”风之痕“，你发的内容被屏蔽了，你要是解决了用gs画cape值的问题，可以发我邮箱么？谢谢

百年孤独

{:eb502:}{:eb502:}

Anonymous888

我也想要知道。

kyoslsl

{:5_213:}

anica

同求                        

赫小猫

anica 发表于 2016-8-22 18:02
同求

请问你现在有计算cape的相关程序吗？