WRF模式運行及后處理簡介-于恩濤

WRF模式運行及后處理簡介竺可楨-南森國際研究中心于恩濤2010-9-91WRF模式簡介2模式編譯3模式設置及運行4后處理報告提綱WRF:WeatherResearchandForecastingModel可用于業(yè)務和科研由包括NCEP、NCAR、NOAA等多家政府機構、科研機構及大學合作開發(fā)與更新采用F90語言，模塊化代碼編寫，便于新模塊耦合與更新WRF

是？哪些人使用WRF??大氣科學家(動力過程，物理過程，天氣氣候研究)?天氣氣候業(yè)務預報人員?氣象應用(空氣質(zhì)量，水文)WRF應用領域：

?ARW&NMM–物理過程/參數(shù)化研究–個例模擬–天氣預報系統(tǒng)研究–數(shù)據(jù)同化?ARW–區(qū)域氣候模擬–耦合大氣化學–全球模擬–Idealized模擬實例2000-2009中國地區(qū)夏季降水模擬Download&compilemodelDownloadsourcecodeCompileWRFCompileWPSWPS&WRF下載與解壓1.下載和2.設置netcdf路徑〔netcdf庫必須由本機編譯器編譯〕3解壓>>2.ExportNETCDF=/usr/local…3.>tarzxfv

>tarzxfvWRF

configure生成configure.wrf文件WRF

compile(1)WRF

compile(2)如編譯成功，那么會生成可執(zhí)行文件Tips:1.編譯首先從最簡單開始，即最好先從串行開始編譯，編譯成功再嘗試編譯并行。2.重新編譯時，./clean–a去除所有，再開始重新編譯。WPStip:WPS使用WRF庫函數(shù)，因此WPS必須放在WRF源代碼同一個目錄下。WPS

compileSETUP&RUNMODELWPS流程(1)設定domain，準備靜態(tài)數(shù)據(jù)geogrid.exe(2)解壓強迫場資料(GFS,NAM...)ungrib.exe(3)將靜態(tài)數(shù)據(jù)和強迫場資料轉為WRF需要格式metgrid.exe所有設定通過namelist.wpsnamelist.wps&sharewrf_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:00:00','2007-09-03_12:00:00',interval_seconds=21600,io_form_geogrid=2,/

start_date,end_date

一個區(qū)域一列

interval_seconds是強迫場資料的時間間隔namelist.wps

(續(xù))

&geogridparent_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,199geog_data_res='2m','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'/thereismore…geogrid–設置區(qū)域，生成靜態(tài)數(shù)據(jù)>./geogrid.exe *生成

geo_em.d01.nc

*“Successfulcompletionofgeogrid.”運行成功

Ungrib準備強迫場資料強迫場可以使用在分析或全球模式結果.需要指定適宜的Vtable來解壓grib1&grib2數(shù)據(jù)“Successfulcompletionofungrib.”說明運行成功>metgrid.exe“Successfulcompletionofmetgrid.”...Output...met_em.d01.2007-09-02_00:00:00.ncmet_em.d01.2007-09-03_00:00:00.nc …….metgrid:///DomainWizard.htmlWRF運行流程運行

real.exe(生成WRF運行所需的輸入資料)運行

wrf.exe都使用namelist.input來控制namelist.input&time_controlrun_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

=09,09,end_day=03,03,end_hour=12,12,end_minute

=00,00,end_second

=00,00,Forstart_*,end_*,onecolumnperdomainnamelist.input

(續(xù))

interval_seconds=10800input_from_file=.true.,.true.,history_interval=60,60,frames_per_outfile

=6,6,restart=.false.,restart_interval

=5000,interval_seconds

同

namelist.wpshistory_interval–結果文件輸出頻率frames_per_outfile

–每個結果文件保存多少時次

restart

=.true.是否重啟運行restart_interval輸出重啟文件的間隔

namelist.input

(續(xù))

&domainstime_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_sn=1,1,1,e_sn=76,199,91,s_vert=1,1,1,e_vert=31,31,31,num_metgrid_levels

=27dx=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,3,namelist.input

(續(xù))

&physicsmp_physics[Microphysics]

=1,1,ra_lw_physics[Longwaverad]=1,1,ra_sw_physics[Shortwaverad]=1,1,radt[Radiationtimestep;min]=60,10,sf_sfclay_physics[Surfacelayer]

=1,1,sf_surface_physics[Surface]

=1,1,bl_pbl_physics[Boundarylayer]

=1,1,bldt[Boundarylayertimestep;min]=0,0,cu_physics[cumulusscheme]=1,0,cudt[cumulustimestep;min]=5,isfflx=1,ifsnow

=0,icloud

=1,surface_input_source

=1,num_soil_layers

=5,mp_zero_out

=0,Microphysics(mp)scheme每個區(qū)域必須一致，但cumulus(cu)schemes各個區(qū)域可以不同不同的組合結果可能不同，但也可能一些組合模式無法運行，需要多多嘗試real>mpirun-np2real.exe>tail–f-->extrapolatingTEMPERATUREnearsfc:i,j,psfc,ptargetd012007-09-03_12:00:00forcingartificialsiltyclayloamLANDCHANGE=0WATERCHANGE=0d012007-09-03_12:00:00Timingforprocessing0s.LBCvalidbetweenthesetimes2007-09-03_09:00:00.00002007-09-03_12:00:00d012007-09-03_12:00:00Timingforoutput0s.d012007-09-03_12:00:00Timingforloop#13=0s.d012007-09-03_12:00:00real_em:SUCCESSCOMPLETEREAL_EMINITwrf.exereal.exe生成文件：wrfbdy_d01andwrfinput_d01(NetCDFfiles)多個區(qū)域wrfinput會有多個SSTupdate那么會有wrflowinp>mpirun-np4wrf.exe& 結果文件：wrfout_d01*rsl.out/rsl.error*WRF后處理推薦NCLload“$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"load“$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"begina=addfile("./geo_em.d01.nc","r")wks=gsn_open_wks("pdf","plt_ter5")res=Trueres@MainTitle="GEOGRIDFIELDS”pltres=Truempres=Trueter=wrf_user_getvar(a,"HGT_M",0)res@cnFillOn=Trueres@ContourParameters=(/0.,1000.,50./)contour=wrf_contour(a,wks,ter,res)plot=wrf_map_overlays(a,wks,(/contour/),\pltres,mpres)endSpecialWRFNCLFunctionswrf_user_getvar Getnativeanddiagnosticvariableswrf_contour/wrf_vector Createline/shaded&vectorplotswrf_map_overlays/wrf_overlays Overlayplotscreatedwithwrf_contourandwrf_vectorwrf_user_intrp3d/wrf_user_intrp2d Interpolatehorizontallytoagivenpressure/height(3ddataonly) Interpolateverticallyalongagivenline

wrf_user_ll_to_ij/wrf_user_ij_to_ll Convert:lat/lonijwrf_user_list_times Getlistiftimesavailableininputfilewrf_user_unstagger UnstaggersanarraySpecialWRFNCLFunctionswrf_user_getvar

Getfieldsfrominputfileter=wrf_user_getvar(a,”HGT”,0){ter=a->HGT(0,:,:)}t2=wrf_user_getvar(a,”T2”,-1){t2=a->T2}avo/pvo:Absolute/PotentialVorticity,cape_2d:2Dmcape/mcin/lcl/lfc,!cape_3d:3Dcape/cin,dbz/mdbz:Reflectivity(3Dandmax),geopt/geopotential:Geopotential,p/pres/pressure:Pressure,!rh/rh2:RelativeHumidity(3Dand2m),slp:SeaLevelPressure,!td/td2:DewPointTemperature(3Dand2m),tc/tk:Temperature(CandF),th/theta:PotentialTemperature,z/height:Height,ua/va/wa:windonmasspoints,uvmet/uvmet10:windrotatedtoearthcoor