《WRF模式运行及后处理简介-于恩涛.ppt》由会员分享，可在线阅读，更多相关《WRF模式运行及后处理简介-于恩涛.ppt（47页珍藏版）》请在三一文库上搜索。
1、WRF模式运行及后处理 简介 竺可桢-南森国际研究中心 于恩涛 2010-9-9 WRFWRF模式模式简简介介 模式模式编译编译 1 模式模式设设置及运行置及运行 2 后后处处理理 3 4 报告提纲 WRF: Weather Research and Forecasting Model 可用于业务和科研 由包括NCEP、NCAR、NOAA等多家政府机构、 科研机构及大学合作开发与更新 采用F90语言，模块化代码编写，便于新模块耦 合与更新 WRFWRF 是是？ 哪些人使用WRF? 大气科学家(动力过程，物理过程 ，天气气候研究) 天气气候业务预报人员 气象应用(空气质量，水文) WRF应用领域
2、： ARW min = 60 , 10 , sf_sfclay_physics Surface layer = 1 , 1 , sf_surface_physics Surface = 1 , 1 , bl_pbl_physics Boundary layer = 1 , 1 , bldt Boundary layer time step; min= 0, 0, cu_physics cumulus scheme = 1 , 0 , cudt cumulus time step; min = 5 , isfflx = 1, ifsnow = 0, icloud = 1, surface_inp
3、ut_source = 1, num_soil_layers = 5 , mp_zero_out = 0 , Microphysics (mp) scheme每个区域必须一致，但cumulus (cu) schemes各个区域可以不同 不同的组合结果可能不同，但也可能一些组合模式无法运行，需要 多多尝试 real mpirun -np 2 real.exe tail f rsl.out.0000 - extrapolating TEMPERATURE near sfc: i,j,psfc, p target d01 2007-09-03_12:00:00 forcing artificial
4、silty clay loam LAND CHANGE = 0 WATER CHANGE = 0 d01 2007-09-03_12:00:00 Timing for processing 0 s. LBC valid between these times 2007-09-03_09:00:00.0000 2007-09- 03_12:00:00 d01 2007-09-03_12:00:00 Timing for output 0 s. d01 2007-09-03_12:00:00 Timing for loop # 13 = 0 s. d01 2007-09-03_12:00:00 r
5、eal_em: SUCCESS COMPLETE REAL_EM INIT wrf.exe real.exe 生成文件： wrfbdy_d01 and wrfinput_d01 (NetCDF files) 多个区域 wrfinput会有多个 SST update则会有wrflowinp mpirun -np 4 wrf.exe & 结果文件： wrfout_d01* rsl.out/rsl.error* WRF后处理 http:/www.mm m.ucar.edu/wrf /OnLineTutoria l/Graphics/NC L/index.html 推荐 NCL load “$NCAR
6、G_ROOT/lib/ncarg/nclscripts/ csm/gsn_code.ncl load “$NCARG_ROOT/lib/ncarg/nclscripts/ wrf/WRFUserARW.ncl begin a = addfile(./geo_em.d01.nc,r) wks = gsn_open_wks(pdf,plt_ter5) res = True resMainTitle = GEOGRID FIELDS” pltres = True mpres = True ter = wrf_user_getvar(a,HGT_M,0) rescnFillOn = True resC
7、ontourParameters = (/0.,1000.,50./) contour = wrf_contour(a,wks,ter,res) plot = wrf_map_overlays(a,wks,(/contour/), pltres,mpres) end Special WRF NCL Functions wrf_user_getvar Get native and diagnostic variables wrf_contour / wrf_vector Create line/shaded & vector plots wrf_map_overlays / wrf_overla
8、ys Overlay plots created with wrf_contour and wrf_vector wrf_user_intrp3d / wrf_user_intrp2d Interpolate horizontally to a given pressure/height (3d data only) Interpolate vertically along a given line wrf_user_ll_to_ij / wrf_user_ij_to_ll Convert: lat/lon ij wrf_user_list_times Get list if times av
9、ailable in input file wrf_user_unstagger Unstaggers an array Special WRF NCL Functions wrf_user_getvar Get fields from input file ter = wrf_user_getvar(a,”HGT”,0) ter=a-HGT(0,:,:) t2 = wrf_user_getvar(a,”T2”,-1) t2=a-T2 avo/pvo: Absolute/Potential Vorticity, cape_2d: 2D mcape/mcin/lcl/lfc, !cape_3d:
10、 3D cape/cin, dbz/mdbz: Reflectivity (3D and max), geopt/geopotential: Geopotential, p/pres/pressure: Pressure, !rh/rh2: Relative Humidity (3D and 2m), slp: Sea Level Pressure, !td/td2: Dew Point Temperature (3D and 2m), tc/tk: Temperature (C and F), th/theta: Potential Temperature, z/height: Height
11、, ua/va/wa: wind on mass points, uvmet/uvmet10: wind rotated to earth coordinates (3D and 10m) WRF NCL Functions example slp = wrf_user_getvar(a,slp,5) t2 = wrf_user_getvar(a,T2,5) u10 = wrf_user_getvar(a,U10,5) v10 = wrf_user_getvar(a,V10,5) oscnLineColor = NavyBlue c_slp = wrf_contour(a,wks,slp,os
12、) otcnFillOn = True c_tc = wrf_contour(a,wks,t2,ot) ovNumVectors = 47 vec = wrf_vector(a,wks,u10,v10,ov) plot = wrf_map_overlays(a, wks, (/c_tc,c_slp,vec/), pltres, mpres) ARWpost 作用：将WRF结果转换为grads或vis5d文件 如果需要转换为vis5d文件，需要vis5d库文件 下载与编译： 编译时需要WRF库文件 ./configure ./compile namelist.ARWpost namelist.A
13、RWpost namelist.ARWpost namelist.ARWpost 实例 2000-2009中国地区夏季降水模拟 实例 2000-2009中国地区夏季降水模拟 web资源 WRF model users site http:/www.mmm.ucar.edu/wrf/users/user_main.html ARW users guide http:/www.mmm.ucar.edu/wrf/users/docs/user_guide/contents.html WRF-ARW/WPS online tutorial http:/www.mmm.ucar.edu/wrf/OnLineTutorial/index.htm WRF namelist description http:/www.mmm.ucar.edu/wrf/users/docs/user_guide/users_guide_chap 5.html#Nml Tutorial presentations http:/www.mmm.ucar.edu/wrf/users/tutorial/tutorial_presentation.htm 谢谢！