如何判断零度层临近上下两层的气压层

绯宜物语

数据是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