WRF模式運(yùn)行及后處理簡介-于恩濤PPT演示課件

1、WRF模式運(yùn)行及后處理簡介,竺可楨-南森國際研究中心 于恩濤 2010-9-9 ,報(bào)告提綱,2,WRF: Weather Research and Forecasting Model 可用于業(yè)務(wù)和科研 由包括NCEP、NCAR、NOAA等多家政府機(jī)構(gòu)、科研機(jī)構(gòu)及大學(xué)合作開發(fā)與更新 采用F90語言，模塊化代碼編寫，便于新模塊耦合與更新,WRF 是,3,哪些人使用WRF? 大氣科學(xué)家(動(dòng)力過程，物理過程，天氣氣候研究) 天氣氣候業(yè)務(wù)預(yù)報(bào)人員 氣象應(yīng)用(空氣質(zhì)量，水文,4,WRF應(yīng)用領(lǐng)域： ARW&NMM 物理過程/參數(shù)化研究 個(gè)例模擬 天氣預(yù)報(bào)系統(tǒng)研究 數(shù)據(jù)同化 ARW 區(qū)域氣候模擬 耦合大氣化學(xué)

2、 全球模擬 Idealized模擬,5,實(shí)例,2000-2009中國地區(qū)夏季降水模擬,6,Download&compile model,Download source code Compile WRF Compile WPS,7,8,WPS&WRF下載與解壓,1. 下載WRFV3.2.1.TAR.gz 和 WPSV3.2.1.TAR.gz 2. 設(shè)置netcdf路徑（netcdf庫必須由本機(jī)編譯器編譯） 3 解壓,wget /wrf/src/WRFV3.2.1.TAR.gz wget /wrf/src/WP

3、SV3.2.1.TAR.gz 2. Export NETCDF=/usr/local 3. tar zxfv WRFV3.2.1.TAR.gz tar zxfv WPSV3.2.1.TAR.gz,9,WRF configure,生成configure.wrf文件,10,WRF compile(1,11,WRF compile(2,如編譯成功，則會(huì)生成可執(zhí)行文件,Tips: 1. 編譯首先從最簡單開始，即最好先從串行開始編譯，編譯成功再嘗試編譯并行。 2. 重新編譯時(shí)，./clean a 清除所有，再開始重新編譯。 WPS tip: WPS使用WRF庫函數(shù)，因此WPS必須放在WRF源代碼同一個(gè)目

4、錄下,12,WPS compile,13,SETUP & RUN MODEL,14,WPS,流程 (1) 設(shè)定domain，準(zhǔn)備靜態(tài)數(shù)據(jù) geogrid.exe (2) 解壓強(qiáng)迫場(chǎng)資料(GFS, NAM.) ungrib.exe (3) 將靜態(tài)數(shù)據(jù)和強(qiáng)迫場(chǎng)資料轉(zhuǎn)為WRF需要格式 metgrid.exe 所有設(shè)定通過namelist.wps,15,namelist.wps,share wrf_core = ARW, max_dom = 1, start_date = 2007-09-02_00:00:00,2007-09-02_00:00:00, end_date = 2007-09-03_12

5、:00:00,2007-09-03_12:00:00, interval_seconds = 21600, io_form_geogrid = 2,start_date, end_date 一個(gè)區(qū)域一列 interval_seconds是強(qiáng)迫場(chǎng)資料的時(shí)間間隔,16,namelist.wps (續(xù),geogrid parent_id = 1, 1, parent_grid_ratio = 1, 3, i_parent_start = 1, 53, j_parent_start = 1, 65, e_we = 100, 259, e_sn = 76, 199 geog_data_res = 2m,

6、2m, dx = 60000, dy = 60000, map_proj = lambert, ref_lat = 35, ref_lon = 105.5, truelat1 = 30.0, truelat2 = 60.0, stand_lon = 105.0, geog_data_path = geog,there is more,17,geogrid 設(shè)置區(qū)域，生成靜態(tài)數(shù)據(jù),./geogrid.exe * 生成 geo_em.d01.nc * “Successful completion of geogrid.” 運(yùn)行成功,18,Ungrib 準(zhǔn)備強(qiáng)迫場(chǎng)資料,強(qiáng)迫場(chǎng)可以使用在分析或全球模式

7、結(jié)果. 需要指定合適的Vtable來解壓grib1&grib2數(shù)據(jù),Successful completion of ungrib.” 說明運(yùn)行成功,19,metgrid.exe “Successful completion of metgrid.” .Output . met_em.d01.2007-09-02_00:00:00.nc met_em.d01.2007-09-03_00:00:00.nc,metgrid,20,/DomainWizard.html,21,WRF運(yùn)行,流程 運(yùn)行 real.exe (生成WRF運(yùn)行所需的輸入資料) 運(yùn)行 wrf

8、.exe 都使用namelist.input來控制,22,namelist.input,time_control run_days = 0, run_hours = 0, run_minutes = 0, run_seconds = 0, start_year = 2007 , 2007 , start_month = 09 , 09 , start_day = 02 , 02 , start_hour = 00 , 00 , start_minute = 00, 00, start_second = 00, 00, end_year = 2007 , 2007 , end_month = 0

9、9 , 09 , end_day = 03 , 03 , end_hour = 12 , 12 , end_minute = 00, 00, end_second = 00, 00,For start_*, end_*, one column per domain,23,namelist.input (續(xù),interval_seconds = 10800 input_from_file = .true., .true., history_interval = 60, 60, frames_per_outfile = 6, 6, restart = .false., restart_interv

10、al = 5000,interval_seconds 同 namelist.wps history_interval 結(jié)果文件輸出頻率 frames_per_outfile 每個(gè)結(jié)果文件保存多少時(shí)次 restart = .true. 是否重啟運(yùn)行 restart_interval 輸出重啟文件的間隔,24,namelist.input (續(xù),domains time_step = 360, time_step_fract_num = 0, time_step_fract_den = 1, max_dom = 1, s_we = 1, 1, 1, e_we = 100, 259, 94, s_s

11、n = 1, 1, 1, e_sn = 76, 199, 91, s_vert = 1, 1, 1, e_vert = 31, 31, 31, num_metgrid_levels = 27 dx = 60000, 12000, 333, dy = 60000, 12000, 333, grid_id = 1, 2, 3, parent_id = 0, 1, 2, i_parent_start = 0, 53, 30, j_parent_start = 0, 65, 30, parent_grid_ratio = 1, 3, 3, parent_time_step_ratio = 1, 3,

12、3,25,namelist.input (續(xù),min = 5 , isfflx = 1, ifsnow = 0, icloud = 1, surface_input_source = 1, num_soil_layers = 5 , mp_zero_out = 0,Microphysics (mp) scheme每個(gè)區(qū)域必須一致，但cumulus (cu) schemes各個(gè)區(qū)域可以不同 不同的組合結(jié)果可能不同，但也可能一些組合模式無法運(yùn)行，需要多多嘗試,26,real,mpirun -np 2 real.exe tail f rsl.out.0000 - extrapolating TEMP

13、ERATURE near sfc: i,j,psfc, p target d01 2007-09-03_12:00:00 forcing artificial 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

14、. d01 2007-09-03_12:00:00 Timing for loop # 13 = 0 s. d01 2007-09-03_12:00:00 real_em: SUCCESS COMPLETE REAL_EM INIT,27,wrf.exe,real.exe 生成文件： wrfbdy_d01 and wrfinput_d01 (NetCDF files) 多個(gè)區(qū)域 wrfinput會(huì)有多個(gè) SST update則會(huì)有wrflowinp mpirun -np 4 wrf.exe & 結(jié)果文件： wrfout_d01* rsl.out/rsl.error,28,WRF后處理,29,h

15、ttp://wrf/OnLineTutorial/Graphics/NCL/index.html,推薦,30,NCL,load “$NCARG_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 FIELD

16、S” 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,31,Special WRF NCL Functions,wrf_user_getvar Get native and diagnostic variables

17、 wrf_contour / wrf_vector Create line/shaded & vector plots wrf_map_overlays / wrf_overlays 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_

18、ll_to_ij / wrf_user_ij_to_ll Convert: lat/lon ij wrf_user_list_times Get list if times available in input file wrf_user_unstagger Unstaggers an array,32,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”,

19、-1) t2=a-T2 avo/pvo: Absolute/Potential Vorticity, cape_2d: 2D mcape/mcin/lcl/lfc, !cape_3d: 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

20、(3D and 2m), tc/tk: Temperature (C and F), th/theta: Potential Temperature, z/height: Height, ua/va/wa: wind on mass points, uvmet/uvmet10: wind rotated to earth coordinates (3D and 10m,33,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) 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/), pltr