wrfrun 0.1.7__py3-none-any.whl

wrfrun/res/extension/plotgrids.ncl

@@ -0,0 +1,216 @@
+
+;   Script display location of model domains
+;   Only works for ARW domains
+;   Only works for NCL versions 6.2 or later
+;   Reads namelist file directly
+
+load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
+load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
+
+begin
+;
+
+; Check the version of NCL
+  version = systemfunc("ncl -V")
+  if(version.lt.6.2) then
+    print("You need NCL V6.2 or later to run this script. Try running plotgrids.ncl. Stopping now...")
+    return
+  end if
+
+; We generate plots, but what kind do we prefer?
+; type = "x11"
+; type = "pdf"
+; type = "ps"
+; type = "ncgm"
+  type = "png"
+  wks = gsn_open_wks(type,"wps_show_dom")
+
+; read the following namelist file
+  filename = "namelist.wps"
+
+; Set the colors to be used
+  colors = (/"white","black","White","ForestGreen","DeepSkyBlue","Red","Blue"/)
+  gsn_define_colormap(wks, colors)  
+
+
+; Set some map information ; line and text information
+  mpres = True
+  mpres@mpFillOn = True
+  mpres@mpFillColors  = (/"background","DeepSkyBlue","ForestGreen","DeepSkyBlue", "transparent"/)
+  mpres@mpDataBaseVersion           = "Ncarg4_1"
+  mpres@mpGeophysicalLineColor      = "Black"
+  mpres@mpGridLineColor             = "Black"
+  mpres@mpLimbLineColor             = "Black"
+  mpres@mpNationalLineColor         = "Black"
+  mpres@mpPerimLineColor            = "Black"
+  mpres@mpUSStateLineColor          = "Black"
+;  mpres@mpOutlineBoundarySets       = "AllBoundaries"
+  ;mpres@mpGridSpacingF              = 45
+  mpres@tiMainString                = " WPS Domain Configuration  "
+
+  lnres = True 
+  lnres@gsLineThicknessF = 2.5
+  lnres@domLineColors    = (/ "white", "Red" , "Red" , "Blue" /)
+
+  txres = True
+  txres@txFont = "helvetica-bold"
+  ;txres@txJust = "BottomLeft"
+  txres@txJust = "TopLeft"
+  txres@txPerimOn = False
+  txres@txFontHeightF = 0.015
+
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+; Do not change anything between the ";;;;;" lines
+
+  maxdom = 21
+  nvar = 19
+  parent_idn = new (maxdom,integer)
+  parent_grid_ration = new (maxdom,integer)
+  i_parent_startn = new (maxdom,integer)
+  j_parent_startn = new (maxdom,integer)
+  e_wen = new (maxdom,integer)
+  e_snn = new (maxdom,integer)
+  plotvar = new((/maxdom,nvar/),float)
+  plotvar@_FillValue = -999.0
+
+  plotvar = wrf_wps_read_nml(filename)
+
+  mpres@max_dom = floattointeger(plotvar(0,0))
+  mpres@dx = plotvar(0,1)
+  mpres@dy = plotvar(0,2)
+  if (.not.ismissing(plotvar(0,3))) then
+    mpres@ref_lat = plotvar(0,3)
+  else
+    mpres@ref_lat = 0.0
+  end if
+  if (.not.ismissing(plotvar(0,4))) then
+    mpres@ref_lon = plotvar(0,4)
+  else
+    mpres@ref_lon = 0.0
+  end if
+  if (.not.ismissing(plotvar(0,5))) then
+    mpres@ref_x = plotvar(0,5)
+  end if
+  if (.not.ismissing(plotvar(0,6))) then
+    mpres@ref_y = plotvar(0,6)
+  end if
+  mpres@truelat1 = plotvar(0,7)
+  mpres@truelat2 = plotvar(0,8)
+  mpres@stand_lon = plotvar(0,9)
+  mproj_int = plotvar(0,10)
+  mpres@pole_lat = plotvar(0,11)
+  mpres@pole_lon = plotvar(0,12)
+
+  do i = 0,maxdom-1
+    parent_idn(i) = floattointeger(plotvar(i,13))
+    parent_grid_ration(i) = floattointeger(plotvar(i,14))
+    i_parent_startn(i) = floattointeger(plotvar(i,15))
+    j_parent_startn(i) = floattointeger(plotvar(i,16))
+    e_wen(i) = floattointeger(plotvar(i,17))
+    e_snn(i) = floattointeger(plotvar(i,18))
+  end do
+
+
+  if(mpres@max_dom .gt. 1) then
+    do i = 1,mpres@max_dom-1
+
+      ;Making sure edge is nested grid is at least 5 grid points from mother domain.
+      if(i_parent_startn(i) .lt. 5) then
+        print("Warning: Western edge of grid must be at least 5 grid points from mother domain!")
+      end if
+      if(j_parent_startn(i) .lt. 5) then
+        print("Warning: Southern edge of grid must be at least 5 grid points from mother domain!")
+      end if
+      pointwe = (e_wen(i)-1.)/parent_grid_ration(i)
+      pointsn = (e_snn(i)-1.)/parent_grid_ration(i)
+      gridwe = e_wen(parent_idn(i)-1)-(pointwe+i_parent_startn(i))
+      gridsn = e_snn(parent_idn(i)-1)-(pointsn+j_parent_startn(i))
+      if(gridwe .lt. 5) then
+        print("Warning: Eastern edge of grid must be at least 5 grid points from mother domain!")
+      end if
+      if(gridsn .lt. 5) then
+        print("Warning: Northern edge of grid must be at least 5 grid points from mother domain!")
+      end if
+
+      ;Making sure nested grid is fully contained in mother domain.
+      gridsizewe = (((e_wen(parent_idn(i)-1)-4)-i_parent_startn(i))*parent_grid_ration(i))-(parent_grid_ration(i)-1)
+      gridsizesn = (((e_snn(parent_idn(i)-1)-4)-j_parent_startn(i))*parent_grid_ration(i))-(parent_grid_ration(i)-1)
+      if(gridwe .lt. 5) then
+        print("Warning: Inner nest (domain = " + (i+1) + ") is not fully contained in mother nest (domain = " + parent_idn(i) + ")!")
+        print("For the current setup of mother domain = " + parent_idn(i) + ", you can only have a nest of size " + gridsizewe + "X" + gridsizesn + ". Stopping Program!")
+        exit
+      end if
+      if(gridsn .lt. 5) then
+        print("Warning: Inner nest (domain = " + (i+1) + ") is not fully contained in mother nest (domain = " + parent_idn(i) + ")!")
+        print("For the current setup of mother domain = " + parent_idn(i) + ", you can only have a nest of size " + gridsizewe + "X" + gridsizesn + ". Stopping Program!")
+        exit
+      end if
+
+      ;Making sure the nest ends at a mother grid domain point.
+      pointwetrunc = decimalPlaces(pointwe,0,False)
+      pointsntrunc = decimalPlaces(pointsn,0,False)
+      if((pointwe-pointwetrunc) .ne. 0.) then
+        nest_we_up = (ceil(pointwe)*parent_grid_ration(i))+1
+        nest_we_dn = (floor(pointwe)*parent_grid_ration(i))+1
+        print("Nest does not end on mother grid domain point. Try " + nest_we_dn + " or " + nest_we_up + ".")
+      end if
+      if((pointsn-pointsntrunc) .ne. 0.) then
+        nest_sn_up = (ceil(pointsn)*parent_grid_ration(i))+1
+        nest_sn_dn = (floor(pointsn)*parent_grid_ration(i))+1
+        print("Nest does not end on mother grid domain point. Try " + nest_sn_dn + " or " + nest_sn_up + ".")
+      end if
+
+    end do
+  end if
+
+  mpres@parent_id = parent_idn(0:mpres@max_dom-1)
+  mpres@parent_grid_ratio = parent_grid_ration(0:mpres@max_dom-1)
+  mpres@i_parent_start = i_parent_startn(0:mpres@max_dom-1)
+  mpres@j_parent_start = j_parent_startn(0:mpres@max_dom-1)
+  mpres@e_we = e_wen(0:mpres@max_dom-1)
+  mpres@e_sn = e_snn(0:mpres@max_dom-1)
+
+  if(mproj_int .eq. 1) then
+    mpres@map_proj = "lambert"
+    mpres@pole_lat = 0.0
+    mpres@pole_lon = 0.0
+  else if(mproj_int .eq. 2) then
+    mpres@map_proj = "mercator"
+    mpres@pole_lat = 0.0
+    mpres@pole_lon = 0.0
+  else if(mproj_int .eq. 3) then
+    mpres@map_proj = "polar"
+    mpres@pole_lat = 0.0
+    mpres@pole_lon = 0.0
+  else if(mproj_int .eq. 4) then
+    mpres@map_proj = "lat-lon"
+  end if
+  end if
+  end if
+  end if
+
+; Deal with global wrf domains that don't have dx or dy
+
+  if (mpres@dx.lt.1e-10 .and. mpres@dx.lt.1e-10) then
+    mpres@dx = 360./(mpres@e_we(0) - 1)
+    mpres@dy = 180./(mpres@e_sn(0) - 1)
+    mpres@ref_lat = 0.0
+    mpres@ref_lon = 180.0
+  end if
+
+  mp = wrf_wps_dom (wks,mpres,lnres,txres)
+
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+
+; Now you can add some information to the plot. 
+; Below is an example of adding a white dot over the DC location.
+  ;pmres = True
+  ;pmres@gsMarkerColor = "White"
+  ;pmres@gsMarkerIndex = 16
+  ;pmres@gsMarkerSizeF = 0.01
+  ;gsn_polymarker(wks,mp,-77.26,38.56,pmres)
+
+
+  frame(wks)           ; lets frame the plot - do not delete
+
+end

wrfrun/res/job_scheduler/pbs.template

@@ -0,0 +1,6 @@
+#PBS -q batch
+#PBS -N wrfrun
+#PBS -o {STDOUT_LOG_PATH}
+#PBS -e {STDERR_LOG_PATH}
+#PBS -l nodes={NODE_NUM}:ppn={CORE_NUM}
+#PBS -l walltime=9999:00:00

wrfrun/res/job_scheduler/slurm.template

@@ -0,0 +1,6 @@
+#SBATCH -J wrfrun                   %Job's name
+#SBATCH -N{NODE_NUM}                %node num
+#SBATCH -n{CORE_NUM}                %core num
+#SBATCH --time=9999:00:00           %maximum run time
+#SBATCH --output={STDOUT_LOG_PATH}  %stdout output log path
+#SBATCH --error={STDERR_LOG_PATH}   %error output log path

wrfrun/res/namelist/namelist.input.da_wrfvar.template

@@ -0,0 +1,261 @@
+&wrfvar1
+    print_detail_grad = .false.
+/
+
+&wrfvar2
+    calc_w_increment = .false.
+    dt_cloud_model = .false.
+    write_mod_filtered_obs = .false.
+    wind_sd = .true.
+    wind_sd_buoy = .true.
+    wind_sd_synop = .true.
+    wind_sd_ships = .true.
+    wind_sd_metar = .true.
+    wind_sd_sound = .true.
+    wind_sd_pilot = .true.
+    wind_sd_airep = .true.
+    wind_sd_qscat = .true.
+    wind_sd_tamdar = .true.
+    wind_sd_geoamv = .true.
+    wind_sd_mtgirs = .true.
+    wind_sd_polaramv = .true.
+    wind_sd_profiler = .true.
+    wind_stats_sd = .false.
+    qc_rej_both = .false.
+/
+
+&wrfvar3
+    ob_format = 2
+    ob_format_gpsro = 1
+/
+
+&wrfvar4
+    use_synopobs = .true.
+    use_shipsobs = .true.
+    use_soundobs = .true.
+    use_qscatobs = .true.
+    use_metarobs = .true.
+    use_pilotobs = .true.
+    use_airepobs = .true.
+    use_tamdarobs = .true.
+    use_geoamvobs = .true.
+    use_polaramvobs = .true.
+    use_bogusobs = .true.
+    use_buoyobs = .true.
+    use_profilerobs = .true.
+    use_satemobs = .true.
+    use_gpspwobs = .true.
+    use_gpsrefobs = .true.
+    use_ssmiretrievalobs = .false.
+    use_airsretobs = .true.
+    use_radarobs = .false.
+    use_radar_rv = .false.
+    use_radar_rf = .false.
+    use_amsuaobs = .false.
+    use_amsubobs = .false.
+    use_eos_amsuaobs = .false.
+    use_iasiobs = .false.
+    use_mhsobs = .false.
+    use_seviriobs = .false.
+    use_airsobs = .false.
+    use_atmsobs = .false.
+    use_ssmisobs = .false.
+    use_hirs3obs = .false.
+    use_hirs4obs = .false.
+    use_msuobs = .false.
+    use_mwtsobs = .false.
+    use_mwhsobs = .false.
+    use_amsr2obs = .false.
+    use_obs_errfac = .false.
+/
+
+&wrfvar5
+    check_max_iv = .true.
+    max_error_t = 5.0
+    max_error_uv = 5.0
+    max_error_pw = 5.0
+    max_error_ref = 5.0
+    max_error_rh = 5.0
+    max_error_q = 5.0
+    max_error_p = 5.0
+    max_error_tb = 5.0
+    max_error_thickness = 5.0
+    max_error_rv = 5.0
+    max_error_rf = 5.0
+/
+
+&wrfvar6
+    max_ext_its = 1
+    ntmax = 150
+    eps = 0.0001, 0.0001, 0.0001, 0.0001, 0.0001
+    orthonorm_gradient = .true.
+/
+
+&wrfvar7
+    cv_options = 3
+    len_scaling1 = 1.0
+    len_scaling2 = 1.0
+    len_scaling3 = 1.0
+    len_scaling4 = 1.0
+    len_scaling5 = 1.0
+    var_scaling1 = 1.0
+    var_scaling2 = 1.0
+    var_scaling3 = 1.0
+    var_scaling4 = 1.0
+    var_scaling5 = 1.0
+/
+
+&wrfvar8
+/
+
+&wrfvar9
+    trace_use = .false.
+/
+
+&wrfvar10
+/
+
+&wrfvar11
+    check_rh = 0
+    sfc_assi_options = 1
+    calculate_cg_cost_fn = .true.
+/
+
+&wrfvar12
+/
+
+&wrfvar13
+/
+
+&wrfvar14
+    rtminit_nsensor = 22
+    rtminit_platform = 1, 1, 1, 1, 9, 10, 1, 1, 1, 10, 1, 1, 10, 12, 12,
+                       9, 17, 2, 1, 1, 1, 1
+    rtminit_satid = 15, 16, 18, 19, 2, 2, 15, 16, 17, 2, 18, 19, 2, 2, 3,
+                    2, 0, 16, 15, 16, 17, 18
+    rtminit_sensor = 3, 3, 3, 3, 3, 3, 4, 4, 4, 16, 15, 15, 15, 21, 21,
+                     11, 19, 10, 0, 0, 0, 0
+    thinning = .true.
+    thinning_mesh = 45.0, 45.0, 45.0, 45.0, 45.0, 45.0, 45.0, 45.0, 45.0,
+                    45.0, 45.0, 45.0, 45.0, 45.0, 45.0, 45.0, 45.0, 45.0
+    qc_rad = .true.
+    write_iv_rad_ascii = .false.
+    write_oa_rad_ascii = .true.
+    use_antcorr = .true., .true., .true., .true., .true., .true., .true.,
+                  .true., .true., .true., .true., .true., .true., .true.,
+                  .true., .true., .true., .true.
+    rtm_option = 1
+    only_sea_rad = .true.
+    use_varbc = .true.
+    varbc_factor = 1.0
+    freeze_varbc = .true.
+    varbc_nobsmin = 10
+    use_crtm_kmatrix = .true.
+    use_rttov_kmatrix = .true.
+    rttov_emis_atlas_ir = 1
+    rttov_emis_atlas_mw = 1
+/
+
+&wrfvar15
+    num_pseudo = 0
+    pseudo_x = 129.0
+    pseudo_y = 26.0
+    pseudo_z = 1.0
+    pseudo_val = -500.0
+    pseudo_err = 1.0
+/
+
+&wrfvar16
+/
+
+&wrfvar17
+/
+
+&wrfvar18
+    analysis_date = '2022-05-19_18:00:00.0000'
+/
+
+&wrfvar19
+    pseudo_var = 'p'
+/
+
+&wrfvar20
+/
+
+&wrfvar21
+    time_window_min = '2022-05-19_17:00:00.0000'
+/
+
+&wrfvar22
+    time_window_max = '2022-05-19_19:00:00.0000'
+/
+
+&wrfvar23
+/
+
+&time_control
+    start_year = 2022
+    start_month = 5
+    start_day = 19
+    start_hour = 18
+    start_minute = 0
+    end_year = 2022
+    end_month = 5
+    end_day = 19
+    end_hour = 18
+    end_minute = 0
+/
+
+&dfi_control
+/
+
+&domains
+    e_we = 110
+    e_sn = 120
+    e_vert = 44
+    dx = 30000
+    dy = 30000
+    nproc_x = -1
+    nproc_y = -1
+    numtiles = 6
+    interp_theta = .true.
+    hypsometric_opt = 2
+/
+
+&physics
+    mp_physics = 2
+    progn = 0
+    no_mp_heating = 0
+    ra_lw_physics = 4
+    ra_sw_physics = 4
+    radt = 15
+    sf_sfclay_physics = 1
+    sf_surface_physics = 2
+    bl_pbl_physics = 1
+    num_land_cat = 21
+    num_soil_layers = 4
+    cu_physics = 1
+    cudt = 0
+    isfflx = 1
+    mp_zero_out = 2
+/
+
+&fdda
+/
+
+&dynamics
+    km_opt = 4
+    damp_opt = 3
+    dampcoef = 0.2
+/
+
+&bdy_control
+/
+
+&grib2
+/
+
+&namelist_quilt
+    nio_tasks_per_group = 0
+    nio_groups = 1
+/