This document will go through the steps to run a wrf model for South Idaho for a 3 day forecast. Other wrf runs will be similar, but not the same.

Transcription

1 Running WRF on R1 Introduction Web Docs This document will go through the steps to run a wrf model for South Idaho for a 3 day forecast. Other wrf runs will be similar, but not the same. Runtime Notes Note: the daily modeling process has been scripted and can found in cmsharedscriptsbuildswrfgcc3.6run_scripts Setting up your environment module load gcc4.8.4 openmpigcc pbspro export WRF_BASE_DIR=PATHTOWRFINSTALL export WRF_VERSION=3.7 export WRF_RUN_DIR=$WRF_BASE_DIR$WRF_VERSION source $WRF_RUN_DIRWPSWRF_Set_Env.bash export Run_Date=`date +"%Y%m%d"` export Log_File=$WRF_RUN_DIRDAILY_LOGSWRF_$Run_Date.log export WORKING_DATE=`date +"%Y-%m-%d"` Getting GFS data cd $WRF_RUN_DIRWPS_GEOG mkdir $WRF_RUN_DIRWPS_GEOGS-Idaho cd S-Idaho cmsharedscriptsbuildswrfgcc3.6run_scriptswrf_g et_gfs_parallel.sh ${Run_Date} this will call a script to get the GFS data in parallel Setting up your WPS and namelist.wps if you havent already done so, back up your WPS directory as a template cd $WRF_RUN_DIR cp -r WPS WPS-Template cd $WRF_RUN_DIRWPS cp cmsharedscriptsbuildswrfgcc3.6namelistsnameli st.wps-2-domain-72hr-template namelist.wps

2 this will copy a template namelist.wps that has the necessary settings to do a daily forecast. Below is a finished namelist.wps, the areas you will need to edit are highlighted. &share wrf_core = 'ARW', max_dom = 2, start_date = ' _06:00:00',' _06:00:00', end_date = ' _06:00:00',' _06:00:00', interval_seconds = 10800, io_form_geogrid = 2, debug_level = 0 &geogrid parent_id = 0,1, parent_grid_ratio = 1,3, i_parent_start = 1,30, j_parent_start = 1,18, e_we = 120,211, e_sn = 80,142, geog_data_res = '30s','30s' dx = 9000 dy = 9000 map_proj ='lambert', ref_lat = , ref_lon = , truelat1 = , truelat2 = , stand_lon = , ref_x = 60.0, ref_y = 40.0, geog_data_path = 'cmhomeggrentzwrf3.7wps_geog', opt_geogrid_tbl_path = 'cmhomeggrentzwrf3.7wpsgeogrid' &ungrib out_format prefix = 'WPS', = 'GFS1' &metgrid fg_name = 'GFS1' io_form_metgrid = 2, opt_metgrid_tbl_path = 'cmhomeggrentzwrf3.7wpsmetgrid', WPS: now you are ready to run geogrid.geogrid.exe link the ungrib data ln -sf ungribvariable_tablesvtable.gfs Vtable.link_grib.csh $WRF_RUN_DIRWPS_GEOGS-Idaho run ungrib.ungrib.exe run metgrid.metgrid.exe WRFV3: if you havent already, back up your WRFV3run directory

3 cp -r $WRF_RUN_DIRWRFV3run $WRF_RUN_DIRWRFV3run- Template cd $WRF_RUN_DIRWRFV3run now you will create a namelist.input file. Below is an example of an edited cmsharedscriptsbuildswrfgcc3.6namelistsnamelist.input-2-domain- 72hr-Template file and the areas that need editing are highlighted. &time_control run_days = 3, run_hours = 0, run_minutes = 0, run_seconds = 0, start_year = 2015, 2015, start_month = 07, 07, start_day = 23, 23, start_hour = 06, 06, start_minute = 00, 00, start_second = 00, 00, end_year = 2015, 2015, end_month = 07, 07, end_day = 26, 26, end_hour = 06, 06, end_minute = 00, 00, end_second = 00, 00, interval_seconds = 10800, input_from_file =.true.,.true., history_interval = 60, 60, frames_per_outfile = 1000, 1000, restart =.false., restart_interval = 1440, io_form_history = 2, io_form_restart = 2, io_form_input = 2, io_form_boundary = 2, debug_level = 0, &domains time_step = 54, time_step_fract_num = 0, time_step_fract_den = 1, max_dom = 2, e_we = 120,211, e_sn = 80,142, e_vert = 40,40, p_top_requested = 5000, num_metgrid_levels = 27, num_metgrid_soil_levels = 4, dx = 9000, 3000, dy = 9000, 3000, grid_id = 1, 2, parent_id = 1, 1, i_parent_start = 1, 30, j_parent_start = 1, 18, parent_grid_ratio = 1, 3, parent_time_step_ratio = 1, 3, feedback = 0, smooth_option = 0, &physics mp_physics = 8, 8, ra_lw_physics = 3, 3, ra_sw_physics = 3, 3, radt = 9, 9,

4 sf_sfclay_physics = 2, 2, sf_surface_physics = 2, 2, bl_pbl_physics = 2, 2, bldt = 0, 0, cu_physics = 1, 0, cudt = 5, 5, surface_input_source = 1, num_soil_layers = 4, sf_urban_physics = 0, 0, &dynamics w_damping = 0, diff_opt = 1, 1, km_opt = 4, 4, diff_6th_opt = 0, 0, diff_6th_factor = 0.12, 0.12, base_temp = 290., damp_opt = 0, zdamp = 5000., 5000., dampcoef = 0.2, 0.2, khdif = 0, 0, kvdif = 0, 0, non_hydrostatic =.true.,.true., moist_adv_opt = 1, 1, scalar_adv_opt = 1, 1, &bdy_control spec_bdy_width = 5, spec_zone = 1, relax_zone = 4, specified =.true.,.false., nested =.false.,.true., &grib2 &namelist_quilt nio_tasks_per_group = 0, nio_groups = 1, Link Metgrid data and run real.exe ln -sf...wpsmet_em.d0*..real.exe now you are ready to run wrf.exe. Below is a pbs submit script for running wrf.exe on the R1 cluster #!binbash # Submit command is: qsub sge-hello-world-gcc.bash # Set the job name #PBS -N WRFDaily # Specify a queue #PBS -q cpu_amd_smp #PBS -j oe # Specify the resources you would like #PBS -l select=4:ncpus=16:mpiprocs=16:mem=24gb

5 ##PBS -l select=64 # Use the current working directory cd $PBS_O_WORKDIR #PBS -S binbash # Load Source System Default Modules echo "Sourcing (loading) default modules \n". etcprofile.dmodules.sh source ~WRF3.7WPSWRF_Set_Env.bash # Add Modules You Need echo "Adding required job modules \n" module add shared module add pbspro module add gcc module add openmpigcc #PBS -v PATH,LD_LIB_PATH # Output How Many Slots The Job Used echo "=========================================================== =====" echo "PE_HOSTFILE file:" cat $PE_HOSTFILE echo "$TMPDIRmachines file:" #cat $TMPDIRmachines echo "=========================================================== =====" echo JOB_ID=$PBS_JOBID echo QUEUE=$PBS_QUEUE echo PBS_CWD_PATH=$PBS_O PATH echo PBS_O_HOME=$PBS_O_HOME echo PBS_O_LOGNAME=$PBS_O_LOGNAME echo PBS_O_PATH=$PBS_O_PATH echo PBS_O_SHELL=$PBS_O_SHELL echo PBS_O_HOST=$PBS_O_HOST echo PBS_WORKDIR=$PBS_O_WORKDIR echo PBS_ENVIRONMENT=$PBS_ENVIRONMENT echo MPI_HOME=$MPI_HOME echo LOADEDMODULES=$LOADEDMODULES # Use MPIRUN To Process The Job echo "***** Starting the MPI process \n" mpirun -v wrf.exe Post-Processing submit the script to the cluster: qsub PBS_WRF-cpu.bash Now you should have two wrfout files in your $WRF_RUN_DIRWRFV3run directory. With these you can generate images of the model.

6 unset NCARG_ROOT cp cmdata1afloreswrf- ScriptsRun_Scriptsncl_plotting_scriptswrf_Temperat ure.ncl.template.wrf_temperature.ncl sed -i 's Working-Date '$WORKING_DATE' g' wrf_temperature.ncl sed -i 's _06: _'$Offset': g' wrf_temperature.ncl mkdir output_$working_date mkdir output_$working_datetemp module load ncarg ffmpeg ncl wrf_temperature.ncl here you can stitch together the images created to make a video cd output_$working_datetemp ffmpeg -r 8 -i frame%05d.png -c:v libx264 -r 8 -pix_fmt yuv420p -s 1920x1080 output.mp4 ffmpeg -r 8 -i frame%05d.png -c:v libtheora -r 8 -pix_fmt yuv420p -s 1920x1080 -q:v 5 output.ogg mv -f output.mp4 cmdata1afloress-idaho3daywrf_$ {WORKING_DATE}_Temperature.mp4 mv -f output.ogg cmdata1afloress-idaho3daywrf_$ {WORKING_DATE}_Temperature.ogg now you should have an mp4 and ogg movie file for the temperature of south Idaho for the next 3 days.