Skip to content

Adding new Ristic microphysics scheme #1308

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

Closed
Nikolar97 wants to merge 42 commits into from
Closed

Adding new Ristic microphysics scheme #1308

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
42 commits
Select commit Hold shift + click to select a range
00f89fc
Adding the new microphysics
Nikolar97 Oct 29, 2020
307c9da
Adding descriptions of custom fields
Nikolar97 Mar 3, 2021
599c6ce
Adding descriptions to custom fields registry
Nikolar97 Mar 3, 2021
e9deaba
improve MAD-WRF section (#1519)
dudhia May 18, 2021
1602bd1
Add got_var_sso to history file for DA cycling application (#1517)
weiwangncar May 18, 2021
e4632b0
Add mod_clddet_geoir.o as compilation dependency to da_radiance.o for…
honnorat Jul 2, 2021
a80a449
Correct units for several NTU microphysics variables (#1532)
weiwangncar Sep 15, 2021
e774092
Add radiation driver dependency on eclipse module (#1544)
weiwangncar Sep 15, 2021
a42b7c6
Default values in URBPARM_LCZ.TBL made more generic (#1520)
matthiasdemuzere Sep 15, 2021
cde45fd
Dimension fix for perturbed bc (#1525)
weiwangncar Sep 16, 2021
7aff4cd
Fix one word: required (#1529)
kwadrat Sep 27, 2021
a4faba8
Fix spelling of "pressure" (#1523)
kwadrat Sep 27, 2021
8401e54
Fix moving nest + Intel tools/registry segfault (#1555)
davegill Sep 27, 2021
30b7ad6
Add CAMtr_volume_mixing_ratio for SSPs (#1553)
jesusff Sep 27, 2021
d55c6cf
revert print and when to turn on w_damping to pre 4.3 code (#1563)
weiwangncar Oct 1, 2021
f4c1c06
NSSL_2mom Snow Reflectivity Correction (#1550)
MicroTed Oct 4, 2021
a743208
Bug fix longwave radiation for bep and bem urban options (#1538)
andreazonato Oct 4, 2021
82e13fb
Increase total number of time periods available (mostly for climate s…
davegill Oct 18, 2021
7d26c10
improvement of the photovoltaic panels scheme (#1551)
andreazonato Oct 18, 2021
84ac4f2
Bug fix for nested aerosol-aware Thompson Scheme QNWFA2D (#1542)
smileMchen Oct 18, 2021
9ed6b94
Updated P3 code with bug fixes and cleanup (#1534)
hmorrison100 Oct 18, 2021
535683f
initializing smstav and a few other local variables before used (#1568)
weiwangncar Oct 22, 2021
7ef7c13
correction for arguments to module_stoch for rand_perturb=1 (#1566)
dudhia Oct 27, 2021
40a4582
Prepare for the release of WRFv4.3.1 (#1556)
smileMchen Oct 27, 2021
4ab0bb5
Merge remote-tracking branch 'origin/release-v4.3.1'
smileMchen Oct 28, 2021
7c6fd57
Small physics syntax errors (#1577)
davegill Nov 29, 2021
ca220f4
bug fix: initialization of 1000-h fuel moisture content, computation …
tmgiannaros Nov 29, 2021
997951e
Fix misspelling in comment: pressure (#1588)
kwadrat Nov 29, 2021
eeb3f2b
smooth_cg_topo=T requires first-guess soil elevation data (#1580)
davegill Nov 29, 2021
4ff1e3d
configure script missing quotes in two echo statments (#1576)
davegill Nov 29, 2021
eacf0fc
Remove photolysis rates from default history stream
jordanschnell Nov 29, 2021
28c33ad
Changes the description of "dustload5" to be load in size bin 5, not …
jordanschnell Nov 29, 2021
30dbdc8
Correct variable name for CLM packaging (#1589)
weiwangncar Nov 30, 2021
238a7d2
Another physics syntax error - similar to #1577 (#1594)
davegill Dec 9, 2021
8f8984e
fix WRFDA CRTM Cray Fortran build (#1595)
jamiebresch Dec 9, 2021
fc99e4e
Add new stanza for ARM aarch64/openmpi (#1557)
yip111 Dec 13, 2021
712e9cc
PBL + LES restrictions (#1596)
davegill Dec 14, 2021
8a928d8
Add time-varying greenhouse gases to RRTMG shortwave options (#1598)
weiwangncar Dec 14, 2021
e2ea09e
Update top-level README and inc/version_decl files to v4.3.2 in prepa…
kkeene44 Dec 14, 2021
57f8f55
Finalize WRFV4.3.2 by merging release-v4.3.2 branch onto master.
kkeene44 Dec 14, 2021
cafc166
Merge branch 'master' of https://github.com/wrf-model/WRF into wrf-mo…
Nikolar97 Dec 25, 2021
0e70827
Merge branch 'wrf-model-master'
Nikolar97 Dec 25, 2021
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
8 changes: 7 additions & 1 deletion Registry/Registry.EM_COMMON
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1189,7 +1189,12 @@ state real F_ICE_PHY ikj misc 1 - rhdu "F
state real F_RAIN_PHY ikj misc 1 - rhdu "F_RAIN_PHY" "FRACTION OF RAIN " ""
state real F_RIMEF_PHY ikj misc 1 - rhdu "F_RIMEF_PHY" "MASS RATIO OF RIMED ICE " ""
state real qndropsource ikj misc 1 - r "qndropsource" "Droplet number source" "#/kg/s"

# ivanr microphysics variables
state real cld ikj misc 1 - r "CLD" "w2id " ""
state real tw2 ikj misc 1 - r "TW2" "w2id " ""
state real pw2 ikj misc 1 - r "PW2" "w2id " ""
state real qw2 ikj misc 1 - r "QW2" "w2id " ""
state real f_snow ikj misc 1 - r "F_SNOW" "w2id " ""

@davegill davegill Feb 9, 2021

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@Nikolar97
Please include a description of these fields with units. An example is the line directly above.
Since these are "r" restart variables, have you tested that you do indeed get reproducible results with a restart?

@Nikolar97 Nikolar97 Mar 3, 2021

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@davegill
Added the descriptions to the fields in the EM_COMMON and NMM registry. The variables should be r, they were used in a similar way to other r fields in other custom schemes already in the model.

# cyl 3DPWP variables
state real OM_TMP i{nocnl}j misc 1 Z i012rhdu=(copy_fcnm)f=(c2f_interp:grid_id) "OM_TMP" "temperature" "k"
state real OM_S i{nocnl}j misc 1 Z i012rhdu=(copy_fcnm)f=(c2f_interp:grid_id) "OM_S" "salinity" "degree"
Expand Down Expand Up @@ -2824,6 +2829,7 @@ package p3_1category mp_physics==50 - moist:qv,qc
package p3_1category_nc mp_physics==51 - moist:qv,qc,qr,qi;scalar:qnc,qni,qnr,qir,qib;state:re_cloud,re_ice,vmi3d,rhopo3d,di3d,refl_10cm,th_old,qv_old
package p3_2category mp_physics==52 - moist:qv,qc,qr,qi,qi2;scalar:qnc,qni,qnr,qir,qib,qni2,qir2,qib2;state:re_cloud,re_ice,vmi3d,rhopo3d,di3d,vmi3d_2,rhopo3d_2,di3d_2,refl_10cm,th_old,qv_old
package morr_tm_aero mp_physics==40 - moist:qv,qc,qr,qi,qs,qg;scalar:qnc,qni,qns,qnr,qng;state:rqrcuten,rqscuten,rqicuten,EFCG,EFIG,EFSG,WACT,CCN1_GS,CCN2_GS,CCN3_GS,CCN4_GS,CCN5_GS,CCN6_GS,CCN7_GS,re_cloud,re_ice,re_snow,mskf_refl_10cm
package ivanr_micro mp_physics==150 - moist:qv,qc,qr,qs;scalar:qt;state:cld,tw2,pw2,qw2,f_snow

@davegill davegill Feb 9, 2021

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@dudhia @weiwangncar
Is there a better MP scheme value than 150?

@Nikolar97 Nikolar97 Mar 3, 2021

Copy link
Copy Markdown
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@davegill
We just used a random one that does not overlap with the others, whatever value is free can be used if 150 is too big/not following a rule.

package jensen_ishmael mp_physics==55 - moist:qv,qc,qr,qi,qi2,qi3;scalar:qnr,qni,qvoli,qaoli,qni2,qvoli2,qaoli2,qni3,qvoli3,qaoli3;state:re_cloud,re_ice,vmi3d,rhopo3d,di3d,phii3d,itype,vmi3d_2,rhopo3d_2,di3d_2,phii3d_2,itype_2,vmi3d_3,rhopo3d_3,di3d_3,phii3d_3,itype_3,refl_10cm
package etampnew mp_physics==95 - moist:qv,qc,qr,qs;scalar:qt;state:f_ice_phy,f_rain_phy,f_rimef_phy
package gsfcgcescheme mp_physics==97 - moist:qv,qc,qr,qi,qs,qg
Expand Down
6 changes: 6 additions & 0 deletions Registry/Registry.NMM
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1057,6 +1057,11 @@ state REAL CPRATE ij misc 1 - rh023
state real F_ICE_PHY ikj misc 1 - - "F_ICE_PHY" "FRACTION OF ICE" ""
state real F_RAIN_PHY ikj misc 1 - - "F_RAIN_PHY" "FRACTION OF RAIN " ""
state real F_RIMEF_PHY ikj misc 1 - - "F_RIMEF_PHY" "MASS RATIO OF RIMED ICE " ""
state real cld ikj misc 1 - r "CLD" "w2id " ""
state real tw2 ikj misc 1 - r "TW2" "w2id " ""
state real pw2 ikj misc 1 - r "PW2" "w2id " ""
state real qw2 ikj misc 1 - r "QW2" "w2id " ""
state real f_snow ikj misc 1 - r "F_SNOW" "w2id " ""
state real MASS_FLUX ij misc 1 - r "MASS_FLUX" "DOWNDRAFT MASS FLUX FOR IN GRELL CUMULUS SCHEME" "mb/hour"
state real apr_gr ij misc 1 - r "APR_GR" "PRECIP FROM CLOSURE OLD_GRELL " "mm/hour"
state real apr_w ij misc 1 - r "APR_W" "PRECIP FROM CLOSURE W " "mm/hour"
Expand Down Expand Up @@ -1775,6 +1780,7 @@ package wdm7scheme mp_physics==26 - moist:qv,qc,q
package p3_1category mp_physics==50 - moist:qv
package p3_1category_nc mp_physics==51 - moist:qv
package p3_2category mp_physics==52 - moist:qv
package ivanr_micro mp_physics==150 - moist:qv,qc,qr,qs;state:cld,tw2,pw2,qw2,f_snow
package jensen_ishmael mp_physics==86 - moist:qv
package etamp_hwrf mp_physics==85 - moist:qv,qc,qr,qi,qs;state:f_ice,f_rain,f_rimef
package etampnew mp_physics==95 - moist:qv,qc,qr,qs;state:f_ice,f_rain,f_rimef
Expand Down
5 changes: 5 additions & 0 deletions dyn_em/solve_em.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3783,6 +3783,11 @@ SUBROUTINE solve_em ( grid , config_flags &
! YLIN
! RI_CURR INPUT
& , RI_CURR=grid%rimi &
& , cld=grid%cld & ! w2id
& , tw2=grid%tw2 & ! w2id
& , pw2=grid%pw2 & ! w2id
& , qw2=grid%qw2 & ! w2id
& , f_snow=grid%f_snow & ! w2id
& , re_cloud=grid%re_cloud, re_ice=grid%re_ice, re_snow=grid%re_snow & ! G. Thompson
& , has_reqc=grid%has_reqc, has_reqi=grid%has_reqi, has_reqs=grid%has_reqs & ! G. Thompson
& , qnwfa2d=grid%qnwfa2d, qnifa2d=grid%qnifa2d & ! G. Thompson
Expand Down
5 changes: 5 additions & 0 deletions dyn_nmm/module_PHYSICS_CALLS.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3423,6 +3423,11 @@ SUBROUTINE GSMDRIVE(NTSD,DT,NPHS,N_MOIST &
& ,ID=grid%id &
& ,num_scalar=1 & !mchen temporary
& ,refl_10cm=grid%refl_10cm & ! to calc. radar reflectivity
& ,cld=grid%cld & ! w2id
& ,tw2=grid%tw2 & ! w2id
& ,pw2=grid%pw2 & ! w2id
& ,qw2=grid%qw2 & ! w2id
& ,f_snow=grid%f_snow & ! w2id
& ,re_cloud=grid%re_cloud & ! G. Thompson
& ,re_ice=grid%re_ice & ! G. Thompson
& ,re_snow=grid%re_snow & ! G. Thompson
Expand Down
1 change: 1 addition & 0 deletions dyn_nmm/solve_nmm.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -356,6 +356,7 @@ SUBROUTINE SOLVE_NMM(GRID,CONFIG_FLAGS &
!
IF (CONFIG_FLAGS%MP_PHYSICS == ETAMPNEW .OR. &
& CONFIG_FLAGS%MP_PHYSICS == FER_MP_HIRES .OR. &
& CONFIG_FLAGS%MP_PHYSICS == IVANR_MICRO .OR. & ! w2id
& CONFIG_FLAGS%MP_PHYSICS == ETAMP_HWRF ) THEN
!
ETAMP_PHYSICS=.TRUE.
Expand Down
3 changes: 2 additions & 1 deletion phys/Makefile
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -84,6 +84,7 @@ MODULES = \
module_mp_wsm6.o \
module_mp_wsm7.o \
module_mp_etanew.o \
module_mp_ivanr_micro.o \
module_mp_fer_hires.o \
module_mp_HWRF.o \
module_mp_thompson.o \
Expand Down Expand Up @@ -192,7 +193,7 @@ MODULES = \
module_wind_fitch.o \
module_sf_lake.o \
module_diagnostics_driver.o \
module_irrigation.o
module_irrigation.o

FIRE_MODULES = \
module_fr_fire_driver.o \
Expand Down
35 changes: 33 additions & 2 deletions phys/module_microphysics_driver.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -123,6 +123,7 @@ SUBROUTINE microphysics_driver( &
,ri_curr &
,diagflag, do_radar_ref &
,ke_diag &
,cld, tw2, pw2, qw2, f_snow & ! w2id
,re_cloud, re_ice, re_snow & ! G. Thompson
,has_reqc, has_reqi, has_reqs & ! G. Thompson
,ccn_conc & ! RAS
Expand Down Expand Up @@ -157,7 +158,7 @@ SUBROUTINE microphysics_driver( &
,NSSL_1MOM,NSSL_1MOMLFO, FER_MP_HIRES_ADVECT &
,WSM7SCHEME, WDM7SCHEME &
,NUWRF4ICESCHEME &
,MILBRANDT2MOM, P3_1CATEGORY, P3_1CATEGORY_NC, JENSEN_ISHMAEL !, P3_2CATEGORY ,MILBRANDT3MOM
,MILBRANDT2MOM, P3_1CATEGORY, P3_1CATEGORY_NC, JENSEN_ISHMAEL,IVANR_MICRO !, P3_2CATEGORY ,MILBRANDT3MOM
#else
USE module_state_description, ONLY : &
KESSLERSCHEME, LINSCHEME, SBU_YLINSCHEME, WSM3SCHEME, WSM5SCHEME &
Expand All @@ -166,7 +167,7 @@ SUBROUTINE microphysics_driver( &
,NSSL_1MOM,NSSL_1MOMLFO, FER_MP_HIRES_ADVECT & ! ,NSSL_3MOM &
,WSM7SCHEME, WDM7SCHEME &
,NUWRF4ICESCHEME &
,MILBRANDT2MOM , CAMMGMPSCHEME,FULL_KHAIN_LYNN, P3_1CATEGORY, P3_1CATEGORY_NC, P3_2CATEGORY, MORR_TM_AERO, JENSEN_ISHMAEL,SPRINKLER !,MILBRANDT3MOM
,MILBRANDT2MOM , CAMMGMPSCHEME,FULL_KHAIN_LYNN, P3_1CATEGORY, P3_1CATEGORY_NC, P3_2CATEGORY, MORR_TM_AERO, JENSEN_ISHMAEL,SPRINKLER,IVANR_MICRO !,MILBRANDT3MOM
#if ( WRFPLUS == 1 )
USE module_state_description, ONLY : LSCONDSCHEME, MKESSLERSCHEME
#endif
Expand Down Expand Up @@ -212,6 +213,7 @@ SUBROUTINE microphysics_driver( &
USE module_mp_wsm6
USE module_mp_wsm7
USE module_mp_etanew
USE module_mp_ivanr_micro
USE module_mp_fer_hires
USE module_mp_thompson
USE module_mp_full_sbm
Expand Down Expand Up @@ -696,6 +698,8 @@ SUBROUTINE microphysics_driver( &
INTEGER, OPTIONAL, INTENT(IN) :: ke_diag ! tells reflectivity calculation whether to do full depth or only k=1
REAL, INTENT(IN) :: ccn_conc ! RAS
INTEGER, OPTIONAL, INTENT(IN) :: do_radar_ref
REAL, DIMENSION(ims:ime,kms:kme,jms:jme), INTENT(INOUT) :: & ! w2id
cld, tw2, pw2, qw2, f_snow
REAL, DIMENSION(ims:ime,kms:kme,jms:jme), INTENT(INOUT) :: & ! G. Thompson
re_cloud, re_ice, re_snow
INTEGER, INTENT(IN):: has_reqc, has_reqi, has_reqs
Expand Down Expand Up @@ -934,7 +938,34 @@ SUBROUTINE microphysics_driver( &
#endif

micro_select: SELECT CASE(mp_physics)
CASE (IVANR_MICRO) !-- Operational 4-km High-Resolution Window (HRW)version IVANR_MICRO
CALL wrf_debug ( 100 , 'microphysics_driver: calling ivanr_micro')

IF ( PRESENT( qv_curr ) .AND. PRESENT( qt_curr ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) ) THEN
CALL IVANRM_NEW( &
ITIMESTEP=itimestep,DT=dt &
,DZ8W=dz8w,RHO_PHY=rho,P_PHY=p,PI_PHY=pi_phy,TH_PHY=th &
,QV=qv_curr &
,QT=qt_curr &
,SR=SR &
,QC=qc_curr &
,QS=qs_curr &
,QR=qr_curr &
,xland=xland &
,CLD=cld &
,T0=tw2 &
,P0=pw2 &
,Q0=qw2 &
,F_SNOW=f_snow &
,RAINNC=rainnc,RAINNCV=rainncv &
,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde &
,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme &
,ITS=its,ITE=ite, JTS=jts,JTE=jte, KTS=kts,KTE=kte &
)
ELSE
CALL wrf_error_fatal ( 'arguments not present for calling ivanr_micro' )
ENDIF
CASE (KESSLERSCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling kessler' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT( QC_CURR ) .AND. &
Expand Down
Loading