diff --git a/phys/MYNN-SFC b/phys/MYNN-SFC index 9dd8e52ead..b66187f54a 160000 --- a/phys/MYNN-SFC +++ b/phys/MYNN-SFC @@ -1 +1 @@ -Subproject commit 9dd8e52ead0bc67e4da0cd41fea31956a3af9d54 +Subproject commit b66187f54a0c169e36592552634867474e603ac0 diff --git a/phys/module_surface_driver.F b/phys/module_surface_driver.F index 50cbc820be..22f96580fb 100644 --- a/phys/module_surface_driver.F +++ b/phys/module_surface_driver.F @@ -112,9 +112,8 @@ SUBROUTINE surface_driver( & & ,zi3d, watsat3d, csol3d, tkmg3d & !lake & ,tkdry3d, tksatu3d, LakeModel, lake_min_elev & !lake #if (EM_CORE==1) - ! & ,lakemask, lakeflag & !lake - & ,lakemask & !lake - , restart_flag & ! restart_flag + & ,lakemask & !lake + , restart_flag & ! restart_flag #endif ! KIM TOFD ,kim_tofd & ! kim tofd @@ -2477,6 +2476,7 @@ SUBROUTINE surface_driver( & its=i_start(ij),ite=i_end(ij),jts=j_start(ij), & jte=j_end(ij), kts=kts,kte=kte, & ustm=ustm,ck=ck,cka=cka,cd=cd,cda=cda, & + xice=xice,xice_threshold=xice_threshold, & sf_mynn_sfcflux_land=sf_mynn_sfcflux_land, & sf_mynn_sfcflux_water=sf_mynn_sfcflux_water, & isfflx=isfflx,shalwater_z0=shalwater_z0, & @@ -5374,8 +5374,9 @@ SUBROUTINE mynn_seaice_wrapper(U3D,V3D,T3D,QV3D, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte, & - ustm,ck,cka,cd,cda,sf_mynn_sfcflux_land, & - sf_mynn_sfcflux_water,flag_lsm,restart, & + ustm,ck,cka,cd,cda, & + sf_mynn_sfcflux_land,sf_mynn_sfcflux_water, & + flag_lsm,restart, & shalwater_z0,water_depth,lakemask, & errmsg,errflg ) @@ -5565,6 +5566,7 @@ SUBROUTINE mynn_seaice_wrapper(U3D,V3D,T3D,QV3D, & ims=ims,ime=ime,jms=jms,jme=jme,kms=kms,kme=kme,& its=its,ite=ite,jts=jts,jte=jte,kts=kts,kte=kte,& ustm=ustm,ck=ck,cka=cka,cd=cd,cda=cda, & + xice=xice,xice_threshold=xice_threshold, & sf_mynn_sfcflux_land=sf_mynn_sfcflux_land, & sf_mynn_sfcflux_water=sf_mynn_sfcflux_water, & isfflx=isfflx,shalwater_z0=shalwater_z0, & @@ -5632,6 +5634,7 @@ SUBROUTINE mynn_seaice_wrapper(U3D,V3D,T3D,QV3D, & ims=ims,ime=ime, jms=jms,jme=jme, kms=kms,kme=kme, & its=its,ite=ite, jts=jts,jte=jte, kts=kts,kte=kte, & ustm=ustm_sea,ck=ck_sea,cka=cka_sea,cd=cd_sea,cda=cda_sea, & + xice=xice,xice_threshold=xice_threshold, & sf_mynn_sfcflux_land=sf_mynn_sfcflux_land, & sf_mynn_sfcflux_water=sf_mynn_sfcflux_water,isfflx=isfflx, & flag_lsm=flag_lsm,restart=restart, & diff --git a/run/README.namelist b/run/README.namelist index 28d6b33e45..931cf092a6 100644 --- a/run/README.namelist +++ b/run/README.namelist @@ -775,6 +775,7 @@ use_rap_aero_icbc = .false. ! Set to .true. to ingest real-t (replacing bl_mynn_tkebudget in prior versions) bl_mynn_closure = 2.5: level 2.5 2.6: level 2.6 (default) + 2.7: level 2.7 (TTE) 3.0: level 3.0 bl_mynn_tkeadvect (max_dom) = .false., ! default off; = .true. do MYNN tke advection icloud_bl option to couple the subgrid-scale clouds from the PBL scheme (MYNN only) @@ -793,7 +794,7 @@ use_rap_aero_icbc = .false. ! Set to .true. to ingest real-t 0: original (Sommeria and Deardorf 1977); 1: Kuwano et al 2010, similar to option 0, but uses resolved scale gradients as opposed to higher order moments ; - 2: from Chaboureau and Bechtold (2002, JAS, with mods, default) + 2: heavily modified Chaboureau and Bechtold (2002, JAS, default) bl_mynn_edmf (max_dom) = 1 ! activate mass-flux option in MYNN, 0: mass-flux option off Related (hidden) options: bl_mynn_edmf_mom (max_dom) = 1 ! - activates momentum transport in MYNN mass-flux scheme @@ -802,6 +803,16 @@ use_rap_aero_icbc = .false. ! Set to .true. to ingest real-t bl_mynn_edmf_tke (max_dom) = 0, ! default; = 1 activates TKE transport in MYNN mass-flux scheme (assuming bl_mynn_edmf > 0) 0 no TKE transport (default);1: activate TKE transport + bl_mynn_output (max_dom) = 0, ! = 0: do not output extra arrays (default) + 1 ! = 1: allocate and output extra 3D arrays from MYNN-EDMF mass-flux component + bl_mynn_mixaerosols (max_dom) = 1, ! = 0: do not mix qnwfa, qnifa, or qnbca + ! = 1: mix aerosols (default) + bl_mynn_mixnumcon (max_dom) = 0, ! = 0: do not mix number concentrations (qnc, qni) (default) + ! = 1: mix number concentrations + bl_mynn_mixscalars (max_dom) = 0 ! = 0: off, 1: activate mixing of scalars (must set up 4D scalar array manually in WRF) + bl_mynn_ess (max_dom) ! = 0: use Chaboureau and Bechtold (2002) buoyancy-flux relationships to estimate effect static stability + ! = 1: use modified OGorman (2011, JAS) effective static stability + bl_mynn_mixqt (max_dom) = 0 ! mixing moisture species separately, 1: mixing total water (do not use) bl_mynn_edmf_dd (max_dom) = 0, ! default = 0; 1 activates downdraft component in MYNN bl_mynn_output (max_dom) = 0, ! default = 0; do not output extra arrays ! =1:allocate and output extra 3D arrays from MYNN PBL for mass flux updrafts @@ -1053,8 +1064,12 @@ use_rap_aero_icbc = .false. ! Set to .true. to ingest real-t seaice_thickness_default = 3.0 ! default value of seaice thickness for seaice_thickness_opt=0 tice2tsk_if2cold = .false. ! set Tice to Tsk to avoid unrealistically low sea ice temperatures iz0tlnd = 0 ! thermal roughness length for sfclay (0 = old, 1 = veg dependent Chen-Zhang Czil, - 2 = Zilitinkevitch (czil=0.1)) - for mynn sfc (0=Zilitinkevitch (def),1=Chen-Zhang,2=mod Yang,3=const zt) + ! 2 = Zilitinkevitch (czil=0.1)). + ! Note for mynn sfc, this iz0tlnd option has been replaced by sf_mynn_sfcflux_water and sf_mynn_sfcflux_land + sf_mynn_sfcflux_water = 1 ! MYNN-SFC option for the calculation of z0, zt and zq over water + ! 0:COARE3.0, 1:COARE3.5 (default), 2:Davis/COARE3.5 ,3:Davis/Garratt, 4:Taylor-Yelland + sf_mynn_sfcflux_land = 0 ! MYNN-SFC option for the calculation of zt and zq over land + ! 0:constant Czil (default), 1:variable Czil, 2:Yang (2008) mp_tend_lim = 10., ! limit on temp tendency from mp latent heating from radar data assimilation prec_acc_dt (max_dom) = 0., ! number of minutes in precipitation bucket - will add three new 2d output fields: prec_acc_c, prec_acc_nc and snow_acc_nc