- 积分
- 9006
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-3-17
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2014-9-18 20:30:27
|
显示全部楼层
本帖最后由 绯宜物语 于 2014-9-18 20:31 编辑
数据是1000hpa到500hpa共16层的温度和geopotential数据。geopotential这个应该是气压层高度吧?就按照你说的方法判断并计算出了零度层高度,输出grd文件,ctl文件的缺测值也有设置。判断、计算是这样的:
if( tem(ix,iy,iz,it)*tem(ix,iy,iz+1,it).LE. 0 ) then
h(ix,iy,iz,it)=(geo(ix,iy,iz+1,it)-geo(ix,iy,iz,it))/(tem(ix,iy,iz+1,it)-tem(ix,iy,iz,it))*(-tem(ix,iy,iz,it))+geo(ix,iy,iz,it)
write(3) h(ix,iy,iz,it)
endif
ctl:
dset f:\tg.grdtitle tg
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 16 levels 1000 975 950 925 900 875
850 825 800 775 750 700 650 600 550 500
tdef 1 linear 00z01JUL1979 1mo
vars 1
h 16 t,z,y,x heightendvars
|
|