times = wrf_user_list_times(a) minLat = min(lat) maxLat = max(lat) minLon = min(lon) maxLon = max(lon) ;calculate the total precipitation rain1 = wrf_user_getvar(a,"RAINNC",-1) rain2 = wrf_user_getvar(a,"RAINC",-1) rain3 = wrf_user_getvar(a,"RAINSH",-1) dim_size...
filename = "wrfout_d01_2023-07-01_00_00_00" a = addfile(filename,"r") ;---Read several WRF variables at first time step it = 0 tc = wrf_user_getvar(a,"tc",it) ; 3D temperature u = wrf_user_getvar(a,"ua",it) ; 3D U at mass points v = wrf_user_getvar(a,"va",...
以下是一个示例NCL脚本,用于处理WRFOUT文件并绘制地表温度图。 1. NCL脚本的基本结构 一个NCL脚本通常包含以下几个部分: 加载必要的库:这些库提供了处理和可视化数据所需的函数。 读取WRFOUT文件:使用addfile或addfiles函数读取WRFOUT文件。 提取变量:使用wrf_user_getvar等函数从WRFOUT文件中提取所需的变量。
pm_ctrl_plane = wrf_user_intrp3d(pm_ctrl,z,"v",plane,angle,False) pblh_ctrl = wrf_user_getvar(a, "PBLH",0) ;二维的 pblh_plane_ctrl = wrf_user_intrp2d(pblh_ctrl,plane,angle,False) ;画图函数 contour_pm= wrf_contour(a,wks,pm_plane(0:zmax_pos,lat_min:lat_max),opts_cn)...
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为选中区域,其他区域会设为缺测 借助此方法,可以仅绘制或计算部分区域。
u = wrf_user_getvar(a,"U10",-1) output_u=new((/station_num,filenum+5/), float) output_u(:,0)=station output_u(:,1)=lon output_u(:,2)=lat do istation=0,station_num-1 loc = wrf_user_ll_to_ij(a,lon(istation), lat(istation), True)-1 ...
Time = wrf_user_getvar(a,"Times",-1) n = dimsizes(Time) strs_raw = new(85,string) time_num = ispan(0,85-1,1)+8 ;时间间隔8小时 str1 = charactertostring(Time(0,0:9)) str2 = charactertostring(Time(0,11:18)) time_num@units = "hours since "+str1+" "+str2 utc_date ...
USTCACELab 便捷处理方法 NCL的wrf_user_getvar()函数可以将WRFOutput直接插值到高度层或气压层，实现垂直坐标转换，便于NCL继续处理。Example:USTCACELab Compare H=wrf_user_getvar(fi,”height”,time);直接读取计算层（N层）H=(“PH”+”PHB”)/9.81 H=fi->height ;通过位势计算 网格层(N+1层)
NCL的wrf_user_getvar()函数可以将WRFOutput直接插值到高度层或气 压层,实现垂直坐标转换,便于NCL继续处理。 便捷处理方法 Example: USTCACELab Compare H=wrf_user_getvar(fi,”height”,time);直接读取计算层(N层) H=(“PH”+”PHB”)/9.81;通过位势计算网格层(N+1层) ...
USTCACELab 便捷处理方法 NCL的wrf_user_getvar()函数可以将WRFOutput直接插值到高度层或气压层，实现垂直坐标转换，便于NCL继续处理。Example:USTCACELab Compare H=wrf_user_getvar(fi,”height”,time);直接读取计算层（N层）H=(“PH”+”PHB”)/9.81 H=fi->height ;通过位势计算 网格层(N+1层)