matlab 温图,NCL画假相当位温图

; 假相当位温

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)</