function wrf_user_getvar (nc_file, fld, it) Usage: ter = wrf_user_getvar (a, “HGT”, 0) nc_file:wrfout的nc文件名称 fld:变量名 it:时间设置。假如it=-1,则所有时间均输出。 例如:tc=wrf_user_getvar(f,"tc",-1) === ;*** load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code...
1.z = wrf_user_getvar(a_1, "z",-1) ; grid point height 2.pt = wrf_user_getvar(a_1, "th",-1) ; potential temperature (theta) 位温 th/theta是位温,tc是温度(摄氏度),tk是温度(开尔文) 3.floatSWDOWN(Time, south_north, west_east) ;地表高度的向下的短波辐射通量 SWDOWN:description...
wrf_user_getvar wrf_user_getvar function wrf_user_getvar (nc_file, fld, it)Usage: ter = wrf_user_getvar (a, “HGT”, 0)Get fields from a netCDF file for 1. Any given time by setting it to the time required.2. For all times in the input file(s), by setting it = -1 3...
z = wrf_user_getvar(a,"height",-1) qvapor = wrf_user_getvar(a,"QVAPOR",-1) ; Qv qv = qvapor*1000. qv@units ="g/kg"g = 9.80665 ; m/s2 g@units ="m/s2"u_500 = wrf_user_intrp3d(u,p,"h",500.,0.,False) v_500 = wrf_user_intrp3d(v,p,"h",500.,0.,False) ...
slp = wrf_user_getvar(a,"slp",-1) 1. 2. 3. 4. 5. 6. 7. 8. wrf-python os.chdir('F:/wrfout/wrfoutnew/') file_list = os.listdir(filepath) files = [os.path.join(filepath,x) for x in file_list] wrflist=[Dataset(d) for d in files] ...
(0,:)time=wrf_user_getvar(a,"XTIME",-1)u10=wrf_user_getvar(a,"U10",0)v10=wrf_user_getvar(a,"V10",0)slp=wrf_user_getvar(a,"slp",0)t2=wrf_user_getvar(a,"T2",0)td=wrf_user_getvar(a,"td",0)rainc=wrf_user_getvar(a,"RAINC",0)rainnc=wrf_user_getvar(a,"RAINNC"...
wrf_user_getvar (ucar.edu) RH2.png 风:10米处的风 U10和V10.png 气压:海平面气压 slp.png NCL: Writing CSV (comma-separated values) files (ucar.edu) write_csv_5.ncl 这个脚本非常实用,我贴在这里方便学习。 技术点包括: ①读多个文件 ②根据多个站点经纬度找格点位置 ③将结果写成csv文件输出 ;...
【NCL】 · 4篇 a=addfile("我是路径/wrfout.nc","r")t=wrf_user_getvar(a,"tc",0), land=wrf_user_getvar(a,"LANDMASK",0);land 中海洋区域数据为0, 陆地区域为1mt=mask(t,land,0);此时0为选中区域,其他区域会设为缺测 借助此方法,可以仅绘制或计算部分区域。
(ifil)+.nc,r) slp=wrf_user_getvar(a,PMSL,0) ;Getseallevelpressurefrom file slp = slp * 0.01 ; We want slp in Pa wrf_smooth_2d( slp, 3 ) slp@units = hPa ; Make sure the units match contour = wrf_contour(a,wks,slp,opts) plot = wrf_map_overlays(a,wks,(/contour/),plt...
16、es = True resMainTitle = GEOGRID FIELDS” pltres = True mpres = True ter = wrf_user_getvar(a,HGT_M,0) rescnFillOn = True resContourParameters = (/0.,1000.,50./) contour = wrf_contour(a,wks,ter,res) plot = wrf_map_overlays(a,wks,(/contour/), pltres,mpres) end,32,...