Skip to content

Add RRTMG cloud overlap option 5#1591

Merged
weiwangncar merged 5 commits intowrf-model:developfrom
weiwangncar:rrtmg_cloud
Dec 17, 2021
Merged

Add RRTMG cloud overlap option 5#1591
weiwangncar merged 5 commits intowrf-model:developfrom
weiwangncar:rrtmg_cloud

Conversation

@weiwangncar
Copy link
Collaborator

@weiwangncar weiwangncar commented Dec 7, 2021

TYPE: enhancement

KEYWORDS: RRTMG, cloud overlap option

SOURCE: Mike Iacono, John Henderson, AER

DESCRIPTION OF CHANGES:

This PR adds a new cloud overlap option for RRTMG, exponential-random (cldovrlp=5). Also added is an optional decorrelation length option - either fixed at 2500 m (idcor=0), or a function of latitude (idcor=1). The decorrelation option works for cloud overlapping option 4 (exponential) too. Tests from one 15 km, east Asian domain in summer, suggests that compared to the default option, maximum-random, this exponential-random produces a slightly cooler atmosphere, as shown below.

Fig.1 SWDOWN difference at 2 pm local time.
diff_swdown_2pm

Fig.2 LWDNB difference at 11 pm local time.
diff_lwdnb_11pm

Fig. 3. Domain-averaged temp and height differences at 2 pm and 11 pm. The unit for T is K, and m for H.
pr1598

For more information, please see this report.

LIST OF MODIFIED FILES:
M Registry/Registry.EM_COMMON
M dyn_em/module_first_rk_step_part1.F
M phys/module_ra_rrtmg_lw.F
M phys/module_ra_rrtmg_sw.F
M phys/module_radiation_driver.F
M run/README.namelist
M share/output_wrf.F

TESTS CONDUCTED:

  1. Tested in a 15 km east Asian domain.
  2. The Jenkins tests are all passing.

RELEASE NOTE: This PR adds a new cloud overlap option for RRTMG, exponential-random set by namelist cldovrlp=5. A decorrelation length option, idcor, is also added. It is set either fixed at 2500 m (idcor=0), or a function of latitude (idcor=1). The decorrelation option works for cloud overlapping option 4 (exponential) too.

@weiwangncar
Copy link
Collaborator Author

Test Type              | Expected  | Received |  Failed
= = = = = = = = = = = = = = = = = = = = = = = =  = = = =
Number of Tests        : 19           17
Number of Builds       : 48           45
Number of Simulations  : 160           155        0
Number of Comparisons  : 101           100        0

Failed Simulations are: 
None
Which comparisons are not bit-for-bit: 
None

@davegill
Copy link
Contributor

davegill commented Dec 7, 2021

@weiwangncar
Wei,
Would you add a quick comment about the associated decorrelation length scale options associated with the new cloud overlap, both in the description and in the release notes.

@davegill
Copy link
Contributor

davegill commented Dec 7, 2021

@dudhia @smileMchen @weiwangncar
Folks,
I went through the mods, and they appear to be very clean, organized, and in final form. There is a new option that is available for RRTMG, so the README is updated accordingly. This is probably one of the options that Ming would be interested in looking at, with similar figs as above: mean T and HGT diffs with original overlap, and then the four combinations of new overlap options 4 and 5, each with the two separate decorrelation length settings. This is ready for the physics group to review.

Wei,
As a reminder, place a description to these options in the User Guide. This seems like a lot of work to go unpublished. Is there a journal article that could be included in the release notes? I am OK with these mods as is, just requesting a possible tweak in the comments.

@weiwangncar
Copy link
Collaborator Author

@davegill Will try to get some information from the developers.

@davegill
Copy link
Contributor

davegill commented Dec 7, 2021

@weiwangncar
Wei,
No need to go back to the developers - just add in that the new overlap schemes have an optional decorrelation length option - either fixed at 2500 m, or a function of lat. I just wanted that ALSO mentioned in the commit message, nothing big.

@davegill davegill changed the base branch from master to develop December 8, 2021 16:14
@davegill
Copy link
Contributor

davegill commented Dec 9, 2021

@weiwangncar
This PR passes the new reg tests:

Here are the restart tests (0=pass):

0 = STATUS test_002m em em_real 34 basic
0 = STATUS test_002m em em_real 34 dfi
0 = STATUS test_002m em em_real 34 diff_opt_2
0 = STATUS test_002m em em_real 34 km_opt_1
0 = STATUS test_002m em em_real 34 km_opt_2
0 = STATUS test_002m em em_real 34 km_opt_3
0 = STATUS test_002m em em_real 34 nest_starts_later
0 = STATUS test_002m em em_real 34 nwp_diag
0 = STATUS test_002m em em_real 34 w_damping

Here is the DA builds:

SUCCESS BUILD WRFPlus
56884 -rwxr-xr-x 1 wrfuser wrf 58247080 Dec  8 16:35 main/wrfplus.exe
WRF DA Build WRFDA-4DVar
Wed Dec  8 16:35:49 UTC 2021
Wed Dec  8 16:40:03 UTC 2021
SUCCESS BUILD WRFDA-4DVar
   36 -rwxr-xr-x 1 wrfuser wrf    34416 Dec  8 16:36 var/build/da_advance_time.exe
  108 -rwxr-xr-x 1 wrfuser wrf   110112 Dec  8 16:36 var/build/da_bias_airmass.exe
   96 -rwxr-xr-x 1 wrfuser wrf    97400 Dec  8 16:36 var/build/da_bias_scan.exe
   44 -rwxr-xr-x 1 wrfuser wrf    43600 Dec  8 16:36 var/build/da_bias_sele.exe
   52 -rwxr-xr-x 1 wrfuser wrf    52208 Dec  8 16:36 var/build/da_bias_verif.exe
   88 -rwxr-xr-x 1 wrfuser wrf    89192 Dec  8 16:36 var/build/da_rad_diags.exe
  124 -rwxr-xr-x 1 wrfuser wrf   125888 Dec  8 16:36 var/build/da_tune_obs_desroziers.exe
  272 -rwxr-xr-x 1 wrfuser wrf   276256 Dec  8 16:36 var/build/da_tune_obs_hollingsworth1.exe
  180 -rwxr-xr-x 1 wrfuser wrf   182016 Dec  8 16:36 var/build/da_tune_obs_hollingsworth2.exe
  164 -rwxr-xr-x 1 wrfuser wrf   164216 Dec  8 16:36 var/build/da_update_bc.exe
  120 -rwxr-xr-x 1 wrfuser wrf   119256 Dec  8 16:36 var/build/da_update_bc_ad.exe
  200 -rwxr-xr-x 1 wrfuser wrf   201504 Dec  8 16:36 var/build/da_verif_grid.exe
   88 -rwxr-xr-x 1 wrfuser wrf    88560 Dec  8 16:36 var/build/da_verif_obs.exe
77356 -rwxr-xr-x 1 wrfuser wrf 79211320 Dec  8 16:39 var/build/da_wrfvar.exe
  684 -rwxr-xr-x 1 wrfuser wrf   696448 Dec  8 16:36 var/build/gen_be_addmean.exe
  672 -rwxr-xr-x 1 wrfuser wrf   687960 Dec  8 16:36 var/build/gen_be_cov2d.exe
  688 -rwxr-xr-x 1 wrfuser wrf   700496 Dec  8 16:36 var/build/gen_be_cov2d3d_contrib.exe
  672 -rwxr-xr-x 1 wrfuser wrf   687960 Dec  8 16:36 var/build/gen_be_cov3d.exe
  688 -rwxr-xr-x 1 wrfuser wrf   700496 Dec  8 16:36 var/build/gen_be_cov3d2d_contrib.exe
  684 -rwxr-xr-x 1 wrfuser wrf   696408 Dec  8 16:36 var/build/gen_be_cov3d3d_bin3d_contrib.exe
  692 -rwxr-xr-x 1 wrfuser wrf   704592 Dec  8 16:36 var/build/gen_be_cov3d3d_contrib.exe
  664 -rwxr-xr-x 1 wrfuser wrf   679768 Dec  8 16:36 var/build/gen_be_diags.exe
  680 -rwxr-xr-x 1 wrfuser wrf   696320 Dec  8 16:36 var/build/gen_be_diags_read.exe
  676 -rwxr-xr-x 1 wrfuser wrf   692160 Dec  8 16:36 var/build/gen_be_ensmean.exe
  692 -rwxr-xr-x 1 wrfuser wrf   708552 Dec  8 16:36 var/build/gen_be_ensrf.exe
  752 -rwxr-xr-x 1 wrfuser wrf   769880 Dec  8 16:36 var/build/gen_be_ep1.exe
  824 -rwxr-xr-x 1 wrfuser wrf   843720 Dec  8 16:36 var/build/gen_be_ep2.exe
  736 -rwxr-xr-x 1 wrfuser wrf   749808 Dec  8 16:37 var/build/gen_be_etkf.exe
  680 -rwxr-xr-x 1 wrfuser wrf   696152 Dec  8 16:36 var/build/gen_be_hist.exe
  784 -rwxr-xr-x 1 wrfuser wrf   799280 Dec  8 16:36 var/build/gen_be_stage0_gsi.exe
  804 -rwxr-xr-x 1 wrfuser wrf   819496 Dec  8 16:36 var/build/gen_be_stage0_wrf.exe
  720 -rwxr-xr-x 1 wrfuser wrf   737120 Dec  8 16:36 var/build/gen_be_stage1.exe
  716 -rwxr-xr-x 1 wrfuser wrf   733024 Dec  8 16:36 var/build/gen_be_stage1_1dvar.exe
  720 -rwxr-xr-x 1 wrfuser wrf   733232 Dec  8 16:36 var/build/gen_be_stage1_gsi.exe
  696 -rwxr-xr-x 1 wrfuser wrf   712544 Dec  8 16:36 var/build/gen_be_stage2.exe
  716 -rwxr-xr-x 1 wrfuser wrf   733080 Dec  8 16:36 var/build/gen_be_stage2_1dvar.exe
  168 -rwxr-xr-x 1 wrfuser wrf   171856 Dec  8 16:36 var/build/gen_be_stage2_gsi.exe
  684 -rwxr-xr-x 1 wrfuser wrf   700256 Dec  8 16:36 var/build/gen_be_stage2a.exe
  692 -rwxr-xr-x 1 wrfuser wrf   708448 Dec  8 16:36 var/build/gen_be_stage3.exe
  672 -rwxr-xr-x 1 wrfuser wrf   687968 Dec  8 16:36 var/build/gen_be_stage4_global.exe
  708 -rwxr-xr-x 1 wrfuser wrf   721008 Dec  8 16:36 var/build/gen_be_stage4_regional.exe
  668 -rwxr-xr-x 1 wrfuser wrf   683928 Dec  8 16:36 var/build/gen_be_vertloc.exe
  804 -rwxr-xr-x 1 wrfuser wrf   819224 Dec  8 16:37 var/build/gen_mbe_stage2.exe
Wed Dec  8 16:40:03 UTC 2021
SUCCESSFUL BUILD
0 -rw-r--r-- 1 wrfuser wrf 0 Dec  8 16:35 SUCCESS_BUILD_WRF_d01_wrfplus_18_WRFPlus
0 -rw-r--r-- 1 wrfuser wrf 0 Dec  8 16:40 SUCCESS_BUILD_WRF_d01_all_wrfvar_18_WRFDA-4DVar
Wed Dec  8 16:40:03 UTC 2021

@davegill davegill removed the Feature label Dec 9, 2021
@davegill
Copy link
Contributor

@dudhia @weiwangncar @smileMchen
Folks,
This PR seems ready for a review.

@weiwangncar
Copy link
Collaborator Author

The updated code has passed all Jenkins tests:

Test Type              | Expected  | Received |  Failed
= = = = = = = = = = = = = = = = = = = = = = = =  = = = =
Number of Tests        : 23           24
Number of Builds       : 60           58
Number of Simulations  : 159           155        0
Number of Comparisons  : 96           93        0

Failed Simulations are: 
None
Which comparisons are not bit-for-bit: 
None

@davegill
Copy link
Contributor

@dudhia @smileMchen

Jimy,
Nice catch on the code missing on the SW part.

Jimy and Ming,
This is now ready for review, as suggestions have been addressed.

@dudhia
Copy link
Collaborator

dudhia commented Dec 16, 2021

Is cldovrlp = 5 a namelist option?

@weiwangncar
Copy link
Collaborator Author

weiwangncar commented Dec 16, 2021

@dudhia Yes, it is. You can use it to set all five options for cloud overlapping. For ARW, the default is 2 from registry.

@dudhia
Copy link
Collaborator

dudhia commented Dec 16, 2021 via email

@weiwangncar
Copy link
Collaborator Author

@dudhia Release note updated. This namelist isn't introduced in the PR. It's ben in the code since option 4 was added.

@dudhia
Copy link
Collaborator

dudhia commented Dec 16, 2021 via email

@davegill davegill self-requested a review December 16, 2021 22:46
@weiwangncar weiwangncar merged commit 85d43e1 into wrf-model:develop Dec 17, 2021
vlakshmanan-scala pushed a commit to scala-computing/WRF that referenced this pull request Apr 4, 2024
TYPE: enhancement

KEYWORDS: RRTMG, cloud overlap option

SOURCE: Mike Iacono, John Henderson, AER

DESCRIPTION OF CHANGES:

This PR adds a new cloud overlap option for RRTMG, exponential-random (cldovrlp=5). Also added is an optional decorrelation length option - either fixed at 2500 m (idcor=0), or a function of latitude (idcor=1). The decorrelation option works for cloud overlapping option 4 (exponential) too. Tests from one 15 km, east Asian domain in summer, suggests that compared to the default option, maximum-random, this exponential-random produces a slightly cooler atmosphere.

For more information, please see this [report](https://dtcenter.org/sites/default/files/visitor-projects/Iacono_Henderson_DTC_VP_FinalReport_2019-2020.pdf).

LIST OF MODIFIED FILES: 
M       Registry/Registry.EM_COMMON
M       dyn_em/module_first_rk_step_part1.F
M       phys/module_ra_rrtmg_lw.F
M       phys/module_ra_rrtmg_sw.F
M       phys/module_radiation_driver.F
M       run/README.namelist
M       share/output_wrf.F

TESTS CONDUCTED: 
1. Tested in a 15 km east Asian domain. 
2. The Jenkins tests are all passing.

RELEASE NOTE: This PR adds a new cloud overlap option for RRTMG, exponential-random set by namelist cldovrlp=5. A decorrelation length option, idcor, is also added. It is set either fixed at 2500 m (idcor=0), or a function of latitude (idcor=1). The decorrelation option works for cloud overlapping option 4 (exponential) too.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants