Skip to content

Adding the WRF-Solar EPS model #1547

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
davegill merged 22 commits into from
Jan 24, 2022
Merged

Adding the WRF-Solar EPS model #1547

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
22 commits
Select commit Hold shift + click to select a range
f99755d
Adding WRF-Solar EPS
pedro-jm Aug 27, 2021
050edde
Perturbs MYNN variables in the PBL driver
pedro-jm Sep 8, 2021
327229b
Perturbs microphysics variables in the driver
pedro-jm Sep 9, 2021
6123639
Perturbs cldfra variables in the radiation driver
pedro-jm Sep 9, 2021
5e285b0
solar trigger reggie sans NMM
davegill Sep 16, 2021
4c0e463
Frequency perturbation updates controllod from the namelist
pedro-jm Sep 21, 2021
303d8cf
Merge remote-tracking branch 'myfork/wrfsolar_eps_devel' into wrfsola…
pedro-jm Sep 21, 2021
aa93147
Perturbs shcu variables in the driver
pedro-jm Sep 22, 2021
d535eb6
Improvements in the microphysics perturbations
pedro-jm Sep 22, 2021
d1ff02d
Impro pbl perturbations
pedro-jm Sep 22, 2021
2961741
Perturbations of radiation in the driver
pedro-jm Sep 23, 2021
b0e9daf
Perturbatons of the lsm in the driver
pedro-jm Sep 23, 2021
321de60
Impro lsm
pedro-jm Sep 23, 2021
5817938
Adds to README.namelist spdt
pedro-jm Sep 24, 2021
7c63ccc
trigger jenkins with no NMM HWRF requirement
davegill Sep 27, 2021
e6bfd42
Clean up, trigger jenkins to find failures
davegill Nov 17, 2021
ec932ae
EPS 17:15
davegill Dec 9, 2021
9b0812d
Fix MOD(A,B), where B is defined as zero
davegill Dec 19, 2021
258db09
STOCHPERT.TBL needs to be linked and removed
davegill Dec 19, 2021
e081bc0
num_pert_3d not dash 17:47
davegill Dec 20, 2021
ea74bfe
PEDRO: CHECK THIS! Need to switch some indexing for the 3d stochastic…
davegill Dec 20, 2021
4b71248
Impro to fix the MOD(A,B) statement
pedro-jm Jan 19, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions Makefile
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -584,6 +584,7 @@ em_real : wrf
ln -sf ../../run/ishmael-qi-qc.bin . ; \
ln -sf ../../run/ishmael-qi-qr.bin . ; \
ln -sf ../../run/BROADBAND_CLOUD_GODDARD.bin . ; \
ln -sf ../../run/STOCHPERT.TBL . ; \
if [ $(RWORDSIZE) -eq 8 ] ; then \
ln -sf ../../run/ETAMPNEW_DATA_DBL ETAMPNEW_DATA ; \
ln -sf ../../run/ETAMPNEW_DATA.expanded_rain_DBL ETAMPNEW_DATA.expanded_rain ; \
Expand Down Expand Up @@ -659,6 +660,7 @@ em_real : wrf
ln -sf ../../run/ishmael-qi-qc.bin . ; \
ln -sf ../../run/ishmael-qi-qr.bin . ; \
ln -sf ../../run/BROADBAND_CLOUD_GODDARD.bin . ; \
ln -sf ../../run/STOCHPERT.TBL . ; \
if [ $(RWORDSIZE) -eq 8 ] ; then \
ln -sf ../../run/ETAMPNEW_DATA_DBL ETAMPNEW_DATA ; \
ln -sf ../../run/ETAMPNEW_DATA.expanded_rain_DBL ETAMPNEW_DATA.expanded_rain ; \
Expand Down
1 change: 1 addition & 0 deletions Registry/Registry.EM_COMMON
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2547,6 +2547,7 @@ rconfig integer alevsiz_cu namelist,physics 1 1
rconfig integer aercu_opt namelist,physics 1 0 - "aercu_opt" "aerosol input option for multiscale KF" ""
rconfig real aercu_fct namelist,physics 1 1.0 - "aercu_fct" "aerosol multiplication factor" ""
rconfig integer aercu_used derived 1 0 - "aercu_used" "derived nml for packaging" ""
rconfig logical couple_farms namelist,physics 1 .false. - "couple_farms" "coupling FARMS radiation to land surface model" ""

#BSINGH - added shallowcu_forced_ra, numBins, thBinSize, rBinSize, minDeepFreq, minShallowFreq, shcu_aerosols_opt for CuP scheme

Expand Down
1 change: 1 addition & 0 deletions Registry/registry.dimspec
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10,6 +10,7 @@ ifdef EM_CORE=1

ifdef DA_CORE=0
dimspec stoclev 2 namelist=num_stoch_levels z num_stoch_levels
dimspec pertn3d 4 namelist=num_pert_3d z num_pert_3d
dimspec j 3 standard_domain y south_north
dimspec k 2 standard_domain z bottom_top
endif
Expand Down
79 changes: 79 additions & 0 deletions Registry/registry.stoch
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -83,6 +83,82 @@ state real ALPH_RAND3 - misc 1 - - "A
state real ALPH_RAND4 - misc 1 - - "ALPH_RAND4" "autoregressive coeff. for generic rand. pert." ""
state logical did_stoch - misc 1 - r "DID_STOCH" "Logical to tell us that we already did the initialization for dom 1" ""


# Variables for multi-pert table input
state real gridpt_stddev_mult3d {pertn3d} misc 1 - r "pert3d_std" "gridpoint standard deviation of random perturbations in mult3d" ""
state real stddev_cutoff_mult3d {pertn3d} misc 1 - r "pert3d_cut" "cutoff tails of pdf above this threshold standard deviation" ""
state real lengthscale_mult3d {pertn3d} misc 1 - r "pert3d_len" "Correlation length scale in meters for mult3d" "m"
state real timescale_mult3d {pertn3d} misc 1 - r "pert3d_tim" "Decorrelation time scale in s for mult3d" "s"
state integer mult3d_vertstruc {pertn3d} misc 1 - r "pert3d_ivs" "vertical structure for mult3d: 0=constant, 1=random phase" ""
state integer iseed_mult3d {pertn3d} misc 1 - r "pert3d_isd" "RANDOM SEED FOR mult3d " ""
state integer ISEEDARR_mult3d k{pertn3d} misc 1 - r "ISEEDARR_MULT3D" "Array to hold seed for restart, MULT3D_PERT" "" ""
state integer stepsp - misc 1 - r "stepsp" "Update perturbation every STEPSP time steps"

# Namelist variables for multi-pert option
rconfig integer multi_perturb namelist,stoch max_domains 0 - "stochastic forcing option: 0=none, 1=multi-perturbation option"
rconfig logical pert_farms namelist,stoch max_domains .false. - "adding perturbations to FARMS scheme"
rconfig real pert_farms_albedo namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_farms_aod namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_farms_angexp namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_farms_aerasy namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_farms_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_farms_qc namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_farms_qs namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig logical pert_deng namelist,stoch max_domains .false. - "adding perturbations to Deng scheme"
rconfig real pert_deng_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_deng_qc namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_deng_t namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_deng_w namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig logical pert_mynn namelist,stoch max_domains .false. - "adding perturbations to MYNN scheme"
rconfig real pert_mynn_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_mynn_qc namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_mynn_t namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_mynn_qke namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig logical pert_noah namelist,stoch max_domains .false. - "adding perturbations to Noah scheme"
rconfig real pert_noah_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_noah_t namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_noah_smois namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_noah_tslb namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig logical pert_thom namelist,stoch max_domains .false. - "adding perturbations to Thompson scheme"
rconfig real pert_thom_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_thom_qc namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_thom_qi namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_thom_qs namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_thom_ni namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig logical pert_cld3 namelist,stoch max_domains .false. - "adding perturbations to CLD3 scheme"
rconfig real pert_cld3_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real pert_cld3_t namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
rconfig real spdt namelist,stoch max_domains -1.0 - "Mins to update perts. If spdt<=0 => every time step."

# Perturbation arrays for multi-pert option
state real SPFORCC3d i{pertn3d}j misc 1 XY r "SPFORCC3d" "real spect. coeff. of random perturbation field" "" ""
state real SPFORCS3d i{pertn3d}j misc 1 XY r "SPFORCS3d" "imag. spect. coeff. of random perturbation field" "" ""
state real SP_AMP3d i{pertn3d}j misc 1 - r "SP_AMP3d" "amplitude of random perturbation field" "" ""
state real ALPH_RAND3d {pertn3d} misc 1 - - "ALPH_RAND3d" "autoregressive coeff. for generic rand. pert." ""
state real VERTSTRUCC3d ikj{pertn3d} dyn_em 1 - r "VERTSTRUCC3d" "vertical structure for stoch. pert. forcing " ""
state real VERTSTRUCS3d ikj{pertn3d} dyn_em 1 - r "VERTSTRUCS3d" "vertical structure for stoch. pert. forcing " ""
state real VERTAMPT3d k{pertn3d} misc 1 - r "VERTAMPT3d" "vert. amplitude of stoch. temperature perturb." "" ""

state real - i{stoclev}jf pert3d 1 - -
state real palbedo i{stoclev}jf pert3d 1 - h "PALBEDO" "" ""
state real paod i{stoclev}jf pert3d 1 - h "PAOD" "" ""
state real pangstrom i{stoclev}jf pert3d 1 - h "PANGSTROM" "" ""
state real passymfac i{stoclev}jf pert3d 1 - h "PASSYMFAC" "" ""
state real pqvapor i{stoclev}jf pert3d 1 - h "PQVAPOR" "" ""
state real pqcloud i{stoclev}jf pert3d 1 - h "PQCLOUD" "" ""
state real pqice i{stoclev}jf pert3d 1 - h "PQICE" "" ""
state real pqsnow i{stoclev}jf pert3d 1 - h "PQSNOW" "" ""
state real pni i{stoclev}jf pert3d 1 - h "PNI" "" ""
state real pth i{stoclev}jf pert3d 1 - h "PTH" "" ""
state real ptke i{stoclev}jf pert3d 1 - h "PTKE" "" ""
state real psmois i{stoclev}jf pert3d 1 - h "PSMOIS" "" ""
state real ptsoil i{stoclev}jf pert3d 1 - h "PTSOIL" "" ""
state real pw i{stoclev}jf pert3d 1 - h "PW" "" ""
state real p15 i{stoclev}jf pert3d 1 - h "P15" "" ""

# Namelist parameters for multi-pert option
rconfig integer num_pert_3d namelist,stoch 1 15 - "num_pert_3d" "number of 3d random fields in STOCHPERT.TBL plue one" ""

# Namelist parameters for random number streams
rconfig integer nens namelist,stoch 1 1 - "random number seed for ensemble members " "" ""

Expand Down Expand Up @@ -188,3 +264,6 @@ package skebs_perturb skebs_on==1 - state:ru_tendf_stoch,rv_tendf_
package random_perturb rand_perturb_on==1 - state:rand_pert,SPFORCS,SPFORCC,SP_AMP,VERTSTRUCC,VERTSTRUCS,VERTAMPT,RAND_REAL,RAND_IMAG,RAND_REAL_xxx,RAND_REAL_yyy,RAND_IMAG_xxx,RAND_IMAG_yyy

package stoch_param_perturb spp_on==1 - state:pattern_spp_conv,field_conv,SPFORCS2,SPFORCC2,SP_AMP2,pattern_spp_pbl,field_pbl,SPFORCS3,SPFORCC3,SP_AMP3,pattern_spp_lsm,field_sf,SPFORCS4,SPFORCC4,SP_AMP4,VERTSTRUCC,VERTSTRUCS,VERTAMPT,RAND_REAL,RAND_IMAG,RAND_REAL_xxx,RAND_REAL_yyy,RAND_IMAG_xxx,RAND_IMAG_yyy

package multi_perturb multi_perturb==1 - state:SPFORCS3d,SPFORCC3d,SP_AMP3d,VERTSTRUCC3d,VERTSTRUCS3d,VERTAMPT3d,RAND_REAL,RAND_IMAG,RAND_REAL_xxx,RAND_REAL_yyy,RAND_IMAG_xxx,RAND_IMAG_yyy,stepsp;pert3d:palbedo,paod,pangstrom,passymfac,pqvapor,pqcloud,pqice,pqsnow,pni,pth,ptke,psmois,ptsoil,pw

2 changes: 1 addition & 1 deletion clean
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -54,7 +54,7 @@ if ( "$arg" == '-a' || "$arg" == '-aa' ) then
( cd test ; rm -f */*.exe */ETAMPNEW_DATA* */GENPARM.TBL */LANDUSE.TBL */README.namelist */README.physics_files \
*/RRTM_DATA */SOILPARM.TBL */VEGPARM.TBL */MPTABLE.TBL */URBPARM.TBL */URBPARM_LCZ.TBL */grib2map.tbl \
*/CAM_ABS_DATA */CAM_AEROPT_DATA \
*/CCN_ACTIVATE.BIN \
*/CCN_ACTIVATE.BIN */STOCHPERT.TBL \
*/CAMtr_volume_mixing_ratio.RCP4.5 */CAMtr_volume_mixing_ratio.RCP6 */CAMtr_volume_mixing_ratio.RCP8.5 \
*/CAMtr_volume_mixing_ratio.A1B */CAMtr_volume_mixing_ratio.A2 */CAMtr_volume_mixing_ratio \
*/CLM_*DATA */RRTMG_LW_DATA */RRTMG_SW_DATA \
Expand Down
56 changes: 56 additions & 0 deletions dyn_em/module_first_rk_step_part1.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -320,7 +320,30 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
& ,swdown2=grid%swdown2, swddni2=grid%swddni2 & ! FARMS coupling
& ,swddif2=grid%swddif2, swddir2=grid%swddir2 & ! FARMS coupling
& ,swdownc2=grid%swdownc2, swddnic2=grid%swddnic2 & ! FARMS coupling
& ,couple_farms=config_flags%couple_farms &
#endif
! WRF-Solar EPS: FARMS
& ,multi_perturb=config_flags% multi_perturb &
& ,pert_farms=config_flags%pert_farms &
& ,perts_albedo=grid%pert3d(:,:,:,P_PALBEDO) &
& ,perts_aod=grid%pert3d(:,:,:,P_PAOD) &
& ,perts_angstrom=grid%pert3d(:,:,:,P_PANGSTROM) &
& ,perts_assymfac=grid%pert3d(:,:,:,P_PASSYMFAC) &
& ,perts_qvapor=grid%pert3d(:,:,:,P_PQVAPOR) &
& ,perts_qcloud=grid%pert3d(:,:,:,P_PQCLOUD) &
& ,perts_qsnow=grid%pert3d(:,:,:,P_PQSNOW) &
& ,pert_farms_albedo=config_flags%pert_farms_albedo &
& ,pert_farms_aod=config_flags%pert_farms_aod &
& ,pert_farms_angexp=config_flags%pert_farms_angexp &
& ,pert_farms_aerasy=config_flags%pert_farms_aerasy &
& ,pert_farms_qv=config_flags%pert_farms_qv &
& ,pert_farms_qc=config_flags%pert_farms_qc &
& ,pert_farms_qs=config_flags%pert_farms_qs &
! WRF-Solar EPS: cld3
& ,pert_cld3=config_flags%pert_cld3 &
& ,perts_th=grid%pert3d(:,:,:,P_PTH) &
& ,pert_cld3_qv=config_flags%pert_cld3_qv &
& ,pert_cld3_t=config_flags%pert_cld3_t &
!BSINGH - For WRFCuP scheme
& ,CU_PHYSICS=config_flags%cu_physics & !CuP, wig 5-Oct-2006
& ,SHALLOWCU_FORCED_RA=config_flags%shallowcu_forced_ra & !CuP, wig
Expand Down Expand Up @@ -968,6 +991,17 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
& ,spp_lsm=config_flags%spp_lsm,pattern_spp_lsm=grid%pattern_spp_lsm & !SPP
& ,field_sf=grid%field_sf & !SPP
& ,spp_pbl=config_flags%spp_pbl,pattern_spp_pbl=grid%pattern_spp_pbl & !SPP
! WRF-Solar EPS: Noah LSM
& ,multi_perturb=config_flags% multi_perturb &
& ,pert_noah=config_flags%pert_noah &
& ,perts_qvapor=grid%pert3d(:,:,:,P_PQVAPOR) &
& ,perts_th=grid%pert3d(:,:,:,P_PTH) &
& ,perts_smois=grid%pert3d(:,:,:,P_PSMOIS) &
& ,perts_tsoil=grid%pert3d(:,:,:,P_PTSOIL) &
& ,pert_noah_qv=config_flags%pert_noah_qv &
& ,pert_noah_t=config_flags%pert_noah_t &
& ,pert_noah_smois=config_flags%pert_noah_smois &
& ,pert_noah_tslb=config_flags%pert_noah_tslb &
& )

#ifdef WRF_HYDRO
Expand Down Expand Up @@ -1146,6 +1180,17 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
! For Wind Turbine Drag Parameterizations
& ,phb=grid%phb &
& ,XLAT_U=grid%xlat_u,XLONG_U=grid%xlong_u &
! WRF-Solar EPS: MYNN
& ,multi_perturb=config_flags%multi_perturb &
& ,pert_mynn=config_flags%pert_mynn &
& ,perts_qvapor=grid%pert3d(:,:,:,P_PQVAPOR) &
& ,perts_qcloud=grid%pert3d(:,:,:,P_PQCLOUD) &
& ,perts_th=grid%pert3d(:,:,:,P_PTH) &
& ,perts_tke=grid%pert3d(:,:,:,P_PTKE) &
& ,pert_mynn_qv=config_flags%pert_mynn_qv &
& ,pert_mynn_qc=config_flags%pert_mynn_qc &
& ,pert_mynn_t=config_flags%pert_mynn_t &
& ,pert_mynn_qke=config_flags%pert_mynn_qke &
!Variables required for camuwpbl scheme
& ,Z_AT_W=grid%z_at_w,CLDFRA_OLD_MP=grid%cldfra_old_mp &
& ,CLDFRA=grid%cldfra &
Expand Down Expand Up @@ -1507,6 +1552,17 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
& ,qke=grid%qke &
& ,PBLHAVG=grid%PBLHAVG, TKEAVG=grid%TKEAVG &
& ,BL_PBL_PHYSICS=config_flags%bl_pbl_physics &
! WRF-Solar EPS
& ,multi_perturb=config_flags%multi_perturb &
& ,pert_deng=config_flags%pert_deng &
& ,perts_qvapor=grid%pert3d(:,:,:,P_PQVAPOR) &
& ,perts_qcloud=grid%pert3d(:,:,:,P_PQCLOUD) &
& ,perts_th=grid%pert3d(:,:,:,P_PTH) &
& ,perts_w=grid%pert3d(:,:,:,P_PW) &
& ,pert_deng_qv=config_flags%pert_deng_qv &
& ,pert_deng_qc=config_flags%pert_deng_qc &
& ,pert_deng_t=config_flags%pert_deng_t &
& ,pert_deng_w=config_flags%pert_deng_w &
& )

#if (WRF_CHEM == 1)
Expand Down
Loading