Skip to content

Commit fed10f4

Browse files
pedro-jmpedro-jm
authored
Adding the WRF-Solar EPS model (wrf-model#1547)
TYPE: New feature KEYWORDS: WRF-Solar, ensemble, stochastic SOURCE: Pedro A. Jimenez, Ju Hye Kim, Jimy Dudhia (NCAR) and Jaemo Yang (NREL) DESCRIPTION OF CHANGES: Problem: The WRF-Solar model is a deterministic model tailored for solar energy applications. Having a probabilistic component tailored for solar energy application is missing. Solution: WRF-Solar Ensemble Prediction System (WRF-Solar EPS, Kim et al. 2021) was designed to provide probabilistic forecasts tailored for solar energy applications. We have selected the most relevant variables from 6 WRF packages using tangent linear methodologies (Yang et al. 2021). Then we have used the available stochastic infrastructure in WRF to perturb these variables at runtine. The scheme is activated setting `multi_perturb = 1` in the &stoch block of namelist.input. The characteristics of the perturbations are specified in a new file in the run directory called STOCHPERT.TBL. Then the perturbations to each module and variable are turned on/off via namelist settings (see run/README.namelist). More information can be found in this very preliminar website: https://ral.ucar.edu/solutions/products/wrf-solar-eps Link to Yang et al. 2021 article: https://www-sciencedirect-com.cuucar.idm.oclc.org/science/article/pii/S0038092X21002322?via%3Dihub Link to Kim et al. 2021 article: https://ieeexplore-ieee-org.cuucar.idm.oclc.org/document/9580552 LIST OF MODIFIED FILES: M Registry/Registry.EM_COMMON M Registry/registry.dimspec M Registry/registry.stoch M dyn_em/module_first_rk_step_part1.F M dyn_em/module_first_rk_step_part2.F M dyn_em/module_stoch.F M dyn_em/solve_em.F M phys/module_microphysics_driver.F M phys/module_pbl_driver.F M phys/module_ra_farms.F M phys/module_radiation_driver.F M phys/module_sf_noahdrv.F M phys/module_shallowcu_driver.F M phys/module_shcu_deng.F M phys/module_surface_driver.F M run/README.namelist A run/STOCHPERT.TBL M share/module_check_a_mundo.F TESTS CONDUCTED: 1. We have run multiple years WRF-Solar EPS and analyzed the characteristics of the probabilistic forecasts as well as calibrated forecasts. Our analysis do not show any potential problems at this point. 2. We ran a simulation with 10 ensemble members and calculated the spread of the ensemble (standard deviation) as a function of the lead time. The spread (this PR compared to our other tests) is very similar. Below is plot with the ensemble spread at given lead time original (left) and current mods (right). The area covered is CONUS. ![Screen Shot 2022-01-19 at 5 19 52 PM](https://user-images.githubusercontent.com/14111759/150239976-503980be-fbba-44b6-88f4-f473fcfed198.png) 3. Jenkins testing is all PASS. RELEASE NOTE: WRF-Solar was expanded to have a stochastic ensemble prediction system (WRF-Solar EPS) tailored for solar energy applications (Yang et al. 2021, Kim et al. 2022). The stochastic perturbations can be introduced into variables of six parameterizations controlling cloud and radiation processes. A more detailed description of the model is provided in the WRF-Solar EPS website: https://ral.ucar.edu/solutions/products/wrf-solar-eps Kim, J.-H., P.A. Jimenez, M. Sengupta, J. Yang, J. Dudhia, S. Alessandrini and Y. Xie, 2022: The WRF-Solar Ensemble Prediction System to provide solar irradiance probabilistic forecasts. IEEE J. of Photovoltaics (In press.) Yang, J., J.-H. Kim, P. A. Jimenez, M. Sengupta, J. Dudhia, Y. Xie, A. Golnas, R. Giering: An Efficient method to identify uncertainties of WRF-Solar variables in forecasting solar irradiance using a tangent linear sensitivity analysis. Solar Energy, 220, 509-522.
1 parent 0bda5e0 commit fed10f4

19 files changed

+1234
-36
lines changed

‎Makefile‎

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -582,6 +582,7 @@ em_real : wrf
582582
ln -sf ../../run/ishmael-qi-qc.bin . ; \
583583
ln -sf ../../run/ishmael-qi-qr.bin . ; \
584584
ln -sf ../../run/BROADBAND_CLOUD_GODDARD.bin . ; \
585+
ln -sf ../../run/STOCHPERT.TBL . ; \
585586
if [ $(RWORDSIZE) -eq 8 ] ; then \
586587
ln -sf ../../run/ETAMPNEW_DATA_DBL ETAMPNEW_DATA ; \
587588
ln -sf ../../run/ETAMPNEW_DATA.expanded_rain_DBL ETAMPNEW_DATA.expanded_rain ; \
@@ -663,6 +664,7 @@ em_real : wrf
663664
ln -sf ../../run/ishmael-qi-qc.bin . ; \
664665
ln -sf ../../run/ishmael-qi-qr.bin . ; \
665666
ln -sf ../../run/BROADBAND_CLOUD_GODDARD.bin . ; \
667+
ln -sf ../../run/STOCHPERT.TBL . ; \
666668
if [ $(RWORDSIZE) -eq 8 ] ; then \
667669
ln -sf ../../run/ETAMPNEW_DATA_DBL ETAMPNEW_DATA ; \
668670
ln -sf ../../run/ETAMPNEW_DATA.expanded_rain_DBL ETAMPNEW_DATA.expanded_rain ; \

‎Registry/Registry.EM_COMMON‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -2548,6 +2548,7 @@ rconfig integer alevsiz_cu namelist,physics 1 1
25482548
rconfig integer aercu_opt namelist,physics 1 0 - "aercu_opt" "aerosol input option for multiscale KF" ""
25492549
rconfig real aercu_fct namelist,physics 1 1.0 - "aercu_fct" "aerosol multiplication factor" ""
25502550
rconfig integer aercu_used derived 1 0 - "aercu_used" "derived nml for packaging" ""
2551+
rconfig logical couple_farms namelist,physics 1 .false. - "couple_farms" "coupling FARMS radiation to land surface model" ""
25512552

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

‎Registry/registry.dimspec‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -10,6 +10,7 @@ ifdef EM_CORE=1
1010

1111
ifdef DA_CORE=0
1212
dimspec stoclev 2 namelist=num_stoch_levels z num_stoch_levels
13+
dimspec pertn3d 4 namelist=num_pert_3d z num_pert_3d
1314
dimspec j 3 standard_domain y south_north
1415
dimspec k 2 standard_domain z bottom_top
1516
endif

‎Registry/registry.stoch‎

Lines changed: 79 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -83,6 +83,82 @@ state real ALPH_RAND3 - misc 1 - - "A
8383
state real ALPH_RAND4 - misc 1 - - "ALPH_RAND4" "autoregressive coeff. for generic rand. pert." ""
8484
state logical did_stoch - misc 1 - r "DID_STOCH" "Logical to tell us that we already did the initialization for dom 1" ""
8585

86+
87+
# Variables for multi-pert table input
88+
state real gridpt_stddev_mult3d {pertn3d} misc 1 - r "pert3d_std" "gridpoint standard deviation of random perturbations in mult3d" ""
89+
state real stddev_cutoff_mult3d {pertn3d} misc 1 - r "pert3d_cut" "cutoff tails of pdf above this threshold standard deviation" ""
90+
state real lengthscale_mult3d {pertn3d} misc 1 - r "pert3d_len" "Correlation length scale in meters for mult3d" "m"
91+
state real timescale_mult3d {pertn3d} misc 1 - r "pert3d_tim" "Decorrelation time scale in s for mult3d" "s"
92+
state integer mult3d_vertstruc {pertn3d} misc 1 - r "pert3d_ivs" "vertical structure for mult3d: 0=constant, 1=random phase" ""
93+
state integer iseed_mult3d {pertn3d} misc 1 - r "pert3d_isd" "RANDOM SEED FOR mult3d " ""
94+
state integer ISEEDARR_mult3d k{pertn3d} misc 1 - r "ISEEDARR_MULT3D" "Array to hold seed for restart, MULT3D_PERT" "" ""
95+
state integer stepsp - misc 1 - r "stepsp" "Update perturbation every STEPSP time steps"
96+
97+
# Namelist variables for multi-pert option
98+
rconfig integer multi_perturb namelist,stoch max_domains 0 - "stochastic forcing option: 0=none, 1=multi-perturbation option"
99+
rconfig logical pert_farms namelist,stoch max_domains .false. - "adding perturbations to FARMS scheme"
100+
rconfig real pert_farms_albedo namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
101+
rconfig real pert_farms_aod namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
102+
rconfig real pert_farms_angexp namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
103+
rconfig real pert_farms_aerasy namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
104+
rconfig real pert_farms_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
105+
rconfig real pert_farms_qc namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
106+
rconfig real pert_farms_qs namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
107+
rconfig logical pert_deng namelist,stoch max_domains .false. - "adding perturbations to Deng scheme"
108+
rconfig real pert_deng_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
109+
rconfig real pert_deng_qc namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
110+
rconfig real pert_deng_t namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
111+
rconfig real pert_deng_w namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
112+
rconfig logical pert_mynn namelist,stoch max_domains .false. - "adding perturbations to MYNN scheme"
113+
rconfig real pert_mynn_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
114+
rconfig real pert_mynn_qc namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
115+
rconfig real pert_mynn_t namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
116+
rconfig real pert_mynn_qke namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
117+
rconfig logical pert_noah namelist,stoch max_domains .false. - "adding perturbations to Noah scheme"
118+
rconfig real pert_noah_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
119+
rconfig real pert_noah_t namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
120+
rconfig real pert_noah_smois namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
121+
rconfig real pert_noah_tslb namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
122+
rconfig logical pert_thom namelist,stoch max_domains .false. - "adding perturbations to Thompson scheme"
123+
rconfig real pert_thom_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
124+
rconfig real pert_thom_qc namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
125+
rconfig real pert_thom_qi namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
126+
rconfig real pert_thom_qs namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
127+
rconfig real pert_thom_ni namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
128+
rconfig logical pert_cld3 namelist,stoch max_domains .false. - "adding perturbations to CLD3 scheme"
129+
rconfig real pert_cld3_qv namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
130+
rconfig real pert_cld3_t namelist,stoch max_domains 0.0 - "turn on/off the perturbation"
131+
rconfig real spdt namelist,stoch max_domains -1.0 - "Mins to update perts. If spdt<=0 => every time step."
132+
133+
# Perturbation arrays for multi-pert option
134+
state real SPFORCC3d i{pertn3d}j misc 1 XY r "SPFORCC3d" "real spect. coeff. of random perturbation field" "" ""
135+
state real SPFORCS3d i{pertn3d}j misc 1 XY r "SPFORCS3d" "imag. spect. coeff. of random perturbation field" "" ""
136+
state real SP_AMP3d i{pertn3d}j misc 1 - r "SP_AMP3d" "amplitude of random perturbation field" "" ""
137+
state real ALPH_RAND3d {pertn3d} misc 1 - - "ALPH_RAND3d" "autoregressive coeff. for generic rand. pert." ""
138+
state real VERTSTRUCC3d ikj{pertn3d} dyn_em 1 - r "VERTSTRUCC3d" "vertical structure for stoch. pert. forcing " ""
139+
state real VERTSTRUCS3d ikj{pertn3d} dyn_em 1 - r "VERTSTRUCS3d" "vertical structure for stoch. pert. forcing " ""
140+
state real VERTAMPT3d k{pertn3d} misc 1 - r "VERTAMPT3d" "vert. amplitude of stoch. temperature perturb." "" ""
141+
142+
state real - i{stoclev}jf pert3d 1 - -
143+
state real palbedo i{stoclev}jf pert3d 1 - h "PALBEDO" "" ""
144+
state real paod i{stoclev}jf pert3d 1 - h "PAOD" "" ""
145+
state real pangstrom i{stoclev}jf pert3d 1 - h "PANGSTROM" "" ""
146+
state real passymfac i{stoclev}jf pert3d 1 - h "PASSYMFAC" "" ""
147+
state real pqvapor i{stoclev}jf pert3d 1 - h "PQVAPOR" "" ""
148+
state real pqcloud i{stoclev}jf pert3d 1 - h "PQCLOUD" "" ""
149+
state real pqice i{stoclev}jf pert3d 1 - h "PQICE" "" ""
150+
state real pqsnow i{stoclev}jf pert3d 1 - h "PQSNOW" "" ""
151+
state real pni i{stoclev}jf pert3d 1 - h "PNI" "" ""
152+
state real pth i{stoclev}jf pert3d 1 - h "PTH" "" ""
153+
state real ptke i{stoclev}jf pert3d 1 - h "PTKE" "" ""
154+
state real psmois i{stoclev}jf pert3d 1 - h "PSMOIS" "" ""
155+
state real ptsoil i{stoclev}jf pert3d 1 - h "PTSOIL" "" ""
156+
state real pw i{stoclev}jf pert3d 1 - h "PW" "" ""
157+
state real p15 i{stoclev}jf pert3d 1 - h "P15" "" ""
158+
159+
# Namelist parameters for multi-pert option
160+
rconfig integer num_pert_3d namelist,stoch 1 15 - "num_pert_3d" "number of 3d random fields in STOCHPERT.TBL plue one" ""
161+
86162
# Namelist parameters for random number streams
87163
rconfig integer nens namelist,stoch 1 1 - "random number seed for ensemble members " "" ""
88164

@@ -188,3 +264,6 @@ package skebs_perturb skebs_on==1 - state:ru_tendf_stoch,rv_tendf_
188264
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
189265

190266
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
267+
268+
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
269+

‎clean‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -54,7 +54,7 @@ if ( "$arg" == '-a' || "$arg" == '-aa' ) then
5454
( cd test ; rm -f */*.exe */ETAMPNEW_DATA* */GENPARM.TBL */LANDUSE.TBL */README.namelist */README.physics_files \
5555
*/RRTM_DATA */SOILPARM.TBL */VEGPARM.TBL */MPTABLE.TBL */URBPARM.TBL */URBPARM_LCZ.TBL */grib2map.tbl \
5656
*/CAM_ABS_DATA */CAM_AEROPT_DATA \
57-
*/CCN_ACTIVATE.BIN \
57+
*/CCN_ACTIVATE.BIN */STOCHPERT.TBL \
5858
*/CAMtr_volume_mixing_ratio.RCP4.5 */CAMtr_volume_mixing_ratio.RCP6 */CAMtr_volume_mixing_ratio.RCP8.5 \
5959
*/CAMtr_volume_mixing_ratio.A1B */CAMtr_volume_mixing_ratio.A2 */CAMtr_volume_mixing_ratio \
6060
*/CAMtr_volume_mixing_ratio.SSP119 */CAMtr_volume_mixing_ratio.SSP126 \

‎dyn_em/module_first_rk_step_part1.F‎

Lines changed: 56 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -320,7 +320,30 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
320320
& ,swdown2=grid%swdown2, swddni2=grid%swddni2 & ! FARMS coupling
321321
& ,swddif2=grid%swddif2, swddir2=grid%swddir2 & ! FARMS coupling
322322
& ,swdownc2=grid%swdownc2, swddnic2=grid%swddnic2 & ! FARMS coupling
323+
& ,couple_farms=config_flags%couple_farms &
323324
#endif
325+
! WRF-Solar EPS: FARMS
326+
& ,multi_perturb=config_flags% multi_perturb &
327+
& ,pert_farms=config_flags%pert_farms &
328+
& ,perts_albedo=grid%pert3d(:,:,:,P_PALBEDO) &
329+
& ,perts_aod=grid%pert3d(:,:,:,P_PAOD) &
330+
& ,perts_angstrom=grid%pert3d(:,:,:,P_PANGSTROM) &
331+
& ,perts_assymfac=grid%pert3d(:,:,:,P_PASSYMFAC) &
332+
& ,perts_qvapor=grid%pert3d(:,:,:,P_PQVAPOR) &
333+
& ,perts_qcloud=grid%pert3d(:,:,:,P_PQCLOUD) &
334+
& ,perts_qsnow=grid%pert3d(:,:,:,P_PQSNOW) &
335+
& ,pert_farms_albedo=config_flags%pert_farms_albedo &
336+
& ,pert_farms_aod=config_flags%pert_farms_aod &
337+
& ,pert_farms_angexp=config_flags%pert_farms_angexp &
338+
& ,pert_farms_aerasy=config_flags%pert_farms_aerasy &
339+
& ,pert_farms_qv=config_flags%pert_farms_qv &
340+
& ,pert_farms_qc=config_flags%pert_farms_qc &
341+
& ,pert_farms_qs=config_flags%pert_farms_qs &
342+
! WRF-Solar EPS: cld3
343+
& ,pert_cld3=config_flags%pert_cld3 &
344+
& ,perts_th=grid%pert3d(:,:,:,P_PTH) &
345+
& ,pert_cld3_qv=config_flags%pert_cld3_qv &
346+
& ,pert_cld3_t=config_flags%pert_cld3_t &
324347
!BSINGH - For WRFCuP scheme
325348
& ,CU_PHYSICS=config_flags%cu_physics & !CuP, wig 5-Oct-2006
326349
& ,SHALLOWCU_FORCED_RA=config_flags%shallowcu_forced_ra & !CuP, wig
@@ -996,6 +1019,17 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
9961019
& ,spp_lsm=config_flags%spp_lsm,pattern_spp_lsm=grid%pattern_spp_lsm & !SPP
9971020
& ,field_sf=grid%field_sf & !SPP
9981021
& ,spp_pbl=config_flags%spp_pbl,pattern_spp_pbl=grid%pattern_spp_pbl & !SPP
1022+
! WRF-Solar EPS: Noah LSM
1023+
& ,multi_perturb=config_flags% multi_perturb &
1024+
& ,pert_noah=config_flags%pert_noah &
1025+
& ,perts_qvapor=grid%pert3d(:,:,:,P_PQVAPOR) &
1026+
& ,perts_th=grid%pert3d(:,:,:,P_PTH) &
1027+
& ,perts_smois=grid%pert3d(:,:,:,P_PSMOIS) &
1028+
& ,perts_tsoil=grid%pert3d(:,:,:,P_PTSOIL) &
1029+
& ,pert_noah_qv=config_flags%pert_noah_qv &
1030+
& ,pert_noah_t=config_flags%pert_noah_t &
1031+
& ,pert_noah_smois=config_flags%pert_noah_smois &
1032+
& ,pert_noah_tslb=config_flags%pert_noah_tslb &
9991033
& )
10001034
10011035
#ifdef WRF_HYDRO
@@ -1175,6 +1209,17 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
11751209
! For Wind Turbine Drag Parameterizations
11761210
& ,phb=grid%phb &
11771211
& ,XLAT_U=grid%xlat_u,XLONG_U=grid%xlong_u &
1212+
! WRF-Solar EPS: MYNN
1213+
& ,multi_perturb=config_flags%multi_perturb &
1214+
& ,pert_mynn=config_flags%pert_mynn &
1215+
& ,perts_qvapor=grid%pert3d(:,:,:,P_PQVAPOR) &
1216+
& ,perts_qcloud=grid%pert3d(:,:,:,P_PQCLOUD) &
1217+
& ,perts_th=grid%pert3d(:,:,:,P_PTH) &
1218+
& ,perts_tke=grid%pert3d(:,:,:,P_PTKE) &
1219+
& ,pert_mynn_qv=config_flags%pert_mynn_qv &
1220+
& ,pert_mynn_qc=config_flags%pert_mynn_qc &
1221+
& ,pert_mynn_t=config_flags%pert_mynn_t &
1222+
& ,pert_mynn_qke=config_flags%pert_mynn_qke &
11781223
!Variables required for camuwpbl scheme
11791224
& ,Z_AT_W=grid%z_at_w,CLDFRA_OLD_MP=grid%cldfra_old_mp &
11801225
& ,CLDFRA=grid%cldfra &
@@ -1536,6 +1581,17 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
15361581
& ,qke=grid%qke &
15371582
& ,PBLHAVG=grid%PBLHAVG, TKEAVG=grid%TKEAVG &
15381583
& ,BL_PBL_PHYSICS=config_flags%bl_pbl_physics &
1584+
! WRF-Solar EPS
1585+
& ,multi_perturb=config_flags%multi_perturb &
1586+
& ,pert_deng=config_flags%pert_deng &
1587+
& ,perts_qvapor=grid%pert3d(:,:,:,P_PQVAPOR) &
1588+
& ,perts_qcloud=grid%pert3d(:,:,:,P_PQCLOUD) &
1589+
& ,perts_th=grid%pert3d(:,:,:,P_PTH) &
1590+
& ,perts_w=grid%pert3d(:,:,:,P_PW) &
1591+
& ,pert_deng_qv=config_flags%pert_deng_qv &
1592+
& ,pert_deng_qc=config_flags%pert_deng_qc &
1593+
& ,pert_deng_t=config_flags%pert_deng_t &
1594+
& ,pert_deng_w=config_flags%pert_deng_w &
15391595
& )
15401596

15411597
#if (WRF_CHEM == 1)

0 commit comments

Comments
 (0)