; 假相当位温
begin
fi = addfile("./new.grib2", "r")
time = fi->initial_time0_hours
;
YYYYMMDDHH = tostring(cd_calendar(time, -3))
dims = dimsizes(YYYYMMDDHH)
cst=new((dims),string)
cst1=new((dims),string)
do i = 0,dims-1
utc = str_get_cols(YYYYMMDDHH(i), 0, 7)
hr = str_get_cols(YYYYMMDDHH(i), 8, 9)
cst(i) = systemfunc("date +'%Y%m%d%H' -d " + "'" + utc + " " + hr + " " + 8 + " hour'")
cst1(i) = str_get_cols(cst(i), 2, 9)
end do
;读取数据
rh=fi->RH_P0_L100_GLL0(:,10:30,:,:)
;print(rh(0,0,0)) ;rh=17.2%
lev=fi->lv_ISBL0(10:30)
levv=lev/100
;print(levv)
dims_lev = dimsizes(levv)
;print(levv)</