MODULE module_microphysics_driver
CONTAINS
SUBROUTINE microphysics_driver( &
th, rho, pi_phy, p &
,ht, dz8w, p8w, dt,dx,dy &
,mp_physics, spec_zone &
,specified, channel_switch &
,warm_rain &
,t8w &
,chem_opt, progn &
,cldfra, cldfra_old, exch_h, nsource &
,qlsink, precr, preci, precs, precg &
,xland,snowh,itimestep &
,f_ice_phy,f_rain_phy,f_rimef_phy &
,lowlyr,sr, id &
,ids,ide, jds,jde, kds,kde &
,ims,ime, jms,jme, kms,kme &
,ips,ipe, jps,jpe, kps,kpe &
,i_start,i_end,j_start,j_end,kts,kte &
,num_tiles, naer &
,dlf,dlf2,t_phy,p_hyd,p8w_hyd,tke_pbl,z_at_w,qfx &
,rliq,turbtype3d,smaw3d,wsedl3d,cldfra_old_mp &
,cldfra_mp,cldfra_mp_all,lradius,iradius &
,cldfrai,cldfral,cldfra_conv &
,alt &
,accum_mode,aitken_mode,coarse_mode &
,icwmrsh3d,icwmrdp3d,shfrc3d,cmfmc3d,cmfmc2_3d &
,config_flags,fnm,fnp,rh_old_mp,lcd_old_mp &
,qv_curr,qc_curr,qr_curr,qi_curr,qs_curr,qg_curr &
,qic_curr,qip_curr,qid_curr &
,qnic_curr,qnip_curr,qnid_curr &
,qndrop_curr,qni_curr,qh_curr,qnh_curr &
,qzr_curr,qzi_curr,qzs_curr,qzg_curr,qzh_curr &
,qns_curr,qnr_curr,qng_curr,qnn_curr,qnc_curr &
,qnwfa_curr,qnifa_curr &
,f_qnwfa,f_qnifa &
,qvolg_curr,qvolh_curr &
,qir_curr,qib_curr &
,effr_curr,ice_effr_curr,tot_effr_curr &
,qic_effr_curr,qip_effr_curr,qid_effr_curr &
,f_qv,f_qc,f_qr,f_qi,f_qs,f_qg,f_qndrop,f_qni &
,f_qns,f_qnr,f_qng,f_qnc,f_qnn,f_qh,f_qnh &
, f_qzr,f_qzi,f_qzs,f_qzg,f_qzh &
,f_qvolg,f_qvolh &
,f_qic,f_qip,f_qid &
,f_qnic,f_qnip,f_qnid &
,f_qir,f_qib &
,f_effr,f_ice_effr,f_tot_effr &
,f_qic_effr,f_qip_effr,f_qid_effr &
,cu_used &
,qrcuten, qscuten, qicuten, qccuten &
,qt_curr,f_qt &
,mp_restart_state,tbpvs_state,tbpvs0_state &
,hail,ice2 &
,u,v,w,z &
,rainnc, rainncv &
,snownc, snowncv &
,hailnc, hailncv &
,graupelnc, graupelncv &
,qnwfa2d &
,refl_10cm &
,vmi3d &
,di3d &
,rhopo3d &
,ri_curr &
,diagflag, do_radar_ref &
,re_cloud, re_ice, re_snow &
,has_reqc, has_reqi, has_reqs &
,ccn_conc &
,scalar,num_scalar &
,kext_ql,kext_qs,kext_qg &
,kext_qh,kext_qa &
,kext_qic,kext_qid,kext_qip &
,kext_ft_qic,kext_ft_qid,kext_ft_qip &
,kext_ft_qs,kext_ft_qg &
,height,tempc &
,TH_OLD &
,QV_OLD &
,xlat,xlong,ivgtyp &
,qrimef_curr,f_qrimef &
)
USE module_state_description, ONLY : &
KESSLERSCHEME, LINSCHEME, SBU_YLINSCHEME, WSM3SCHEME, WSM5SCHEME &
,WSM6SCHEME, ETAMPNEW, FER_MP_HIRES, etamp_HWRF,THOMPSON, THOMPSONAERO, MORR_TWO_MOMENT &
,GSFCGCESCHEME, WDM5SCHEME, WDM6SCHEME, NSSL_2MOM, NSSL_2MOMCCN, NSSL_2MOMG &
,NSSL_1MOM,NSSL_1MOMLFO, FER_MP_HIRES_ADVECT &
,MILBRANDT2MOM, P3_1CATEGORY, P3_1CATEGORY_NC
USE module_dm, ONLY : &
local_communicator, mytask, wrf_dm_min_real, wrf_dm_max_real
USE module_model_constants
USE module_wrf_error
USE module_configure, only: grid_config_rec_type
USE module_mp_kessler
USE module_mp_lin
USE module_mp_sbu_ylin
USE module_mp_wsm3
USE module_mp_wsm5
USE module_mp_wsm6
USE module_mp_etanew
USE module_mp_fer_hires
USE module_mp_thompson
USE module_mp_full_sbm
USE module_mp_fast_sbm
USE module_mp_gsfcgce
USE module_mp_morr_two_moment
USE module_mp_p3
USE module_mp_wdm5
USE module_mp_wdm6
USE module_mp_milbrandt2mom
USE module_mp_nssl_2mom
USE module_mp_HWRF
USE module_mixactivate, only: prescribe_aerosol_mixactivate
USE module_utility, ONLY: WRFU_Clock, WRFU_Alarm
USE module_domain, ONLY : HISTORY_ALARM, Is_alarm_tstep
IMPLICIT NONE
INTEGER,parameter :: iunit=6
INTEGER :: mpi_error_code=1
TYPE(grid_config_rec_type), INTENT(IN ) , OPTIONAL :: config_flags
INTEGER, INTENT(IN ) :: mp_physics
LOGICAL, INTENT(IN ) :: specified
INTEGER, OPTIONAL, INTENT(IN ) :: chem_opt, progn
INTEGER, OPTIONAL, INTENT(IN ) :: hail, ice2
INTEGER, INTENT(IN ) :: ids,ide, jds,jde, kds,kde
INTEGER, INTENT(IN ) :: ims,ime, jms,jme, kms,kme,num_scalar
INTEGER, OPTIONAL, INTENT(IN ) :: ips,ipe, jps,jpe, kps,kpe
INTEGER, INTENT(IN ) :: kts,kte
INTEGER, INTENT(IN ) :: itimestep,num_tiles,spec_zone
INTEGER, DIMENSION(num_tiles), INTENT(IN) :: &
& i_start,i_end,j_start,j_end
LOGICAL, INTENT(IN ) :: warm_rain
REAL, DIMENSION( ims:ime , kms:kme , jms:jme ), &
INTENT(INOUT) :: th
REAL, DIMENSION( ims:ime , kms:kme , jms:jme ), &
INTENT(IN ) :: &
rho, &
dz8w, &
p8w, &
pi_phy, &
p
REAL, DIMENSION( ims:ime , kms:kme , jms:jme ),INTENT(INOUT), OPTIONAL :: th_old,qv_old
REAL,DIMENSION(ims:ime,kms:kme,jms:jme,num_scalar),INTENT(INOUT), OPTIONAL :: scalar
INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(IN), OPTIONAL:: IVGTYP
REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN), OPTIONAL :: XLAT, XLONG
REAL,INTENT(IN), OPTIONAL ::accum_mode,aitken_mode,coarse_mode
REAL , DIMENSION( kms:kme ) , &
INTENT(IN ) , OPTIONAL :: fnm, &
fnp
REAL, DIMENSION( ims:ime, jms:jme ), &
INTENT(IN), OPTIONAL :: &
qfx, &
rliq
REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
INTENT(IN), OPTIONAL :: &
dlf, &
dlf2, &
t_phy, &
p_hyd, &
p8w_hyd, &
z_at_w, &
tke_pbl, &
turbtype3d, &
smaw3d, &
alt, &
icwmrsh3d, &
icwmrdp3d, &
shfrc3d, &
cmfmc3d, &
cmfmc2_3d
REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
INTENT(INOUT) , OPTIONAL :: &
cldfra_old_mp, &
rh_old_mp, &
lcd_old_mp
REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
INTENT(INOUT) , OPTIONAL:: &
wsedl3d, &
cldfra_mp, &
cldfra_mp_all, &
cldfrai, &
cldfral, &
lradius, &
iradius, &
cldfra_conv
REAL, INTENT(INOUT), DIMENSION(ims:ime, kms:kme, jms:jme ) :: &
F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY
REAL, OPTIONAL, INTENT(OUT), DIMENSION(ims:ime, kms:kme, jms:jme ) :: &
qlsink, &
precr, &
preci, &
precs, &
precg
REAL , DIMENSION( ims:ime , jms:jme ) , INTENT(IN) :: XLAND
REAL , DIMENSION( ims:ime , jms:jme ) , INTENT(IN), OPTIONAL :: SNOWH
REAL , DIMENSION( ims:ime , jms:jme ) , INTENT(OUT) :: SR
REAL, INTENT(IN ) :: dt,dx,dy
INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT) :: LOWLYR
REAL, OPTIONAL, DIMENSION( ims:ime , kms:kme, jms:jme ) , INTENT(OUT) :: refl_10cm
REAL, OPTIONAL, DIMENSION( ims:ime , kms:kme, jms:jme ) , INTENT(OUT) :: vmi3d,di3d,rhopo3d
LOGICAL, OPTIONAL, INTENT(IN ) :: channel_switch
REAL, OPTIONAL, INTENT(INOUT ) :: naer
REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) , OPTIONAL :: qnwfa2d
REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
OPTIONAL, &
INTENT(INOUT ) :: &
u,v,w, z, t8w &
,cldfra, cldfra_old, exch_h &
,qv_curr,qc_curr,qr_curr,qi_curr,qs_curr,qg_curr &
,qt_curr,qndrop_curr,qni_curr,qh_curr,qnh_curr &
,qns_curr,qnr_curr,qng_curr,qnn_curr,qnc_curr &
,qic_curr,qip_curr,qid_curr &
,qnic_curr,qnip_curr,qnid_curr &
,qzr_curr,qzi_curr,qzs_curr,qzg_curr,qzh_curr &
,qir_curr,qib_curr &
,effr_curr,ice_effr_curr,tot_effr_curr &
,qic_effr_curr,qip_effr_curr,qid_effr_curr &
,kext_ql,kext_qs,kext_qg &
,kext_qh,kext_qa &
,kext_qic,kext_qip,kext_qid,tempc,height &
,kext_ft_qic,kext_ft_qip,kext_ft_qid &
,kext_ft_qs,kext_ft_qg &
,qnwfa_curr,qnifa_curr &
,qvolg_curr,qvolh_curr, qrimef_curr
REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
OPTIONAL, &
INTENT(IN) :: qrcuten, qscuten, qicuten, qccuten
INTEGER, INTENT(IN), optional :: cu_used
REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
OPTIONAL, &
INTENT(INOUT) :: ri_curr
REAL, DIMENSION(ims:ime, kms:kme, jms:jme ), &
OPTIONAL, &
INTENT(OUT ) :: &
nsource
REAL, DIMENSION( ims:ime , jms:jme ), &
INTENT(INOUT), &
OPTIONAL :: &
RAINNC &
,RAINNCV &
,SNOWNC &
,SNOWNCV &
,GRAUPELNC &
,GRAUPELNCV &
,HAILNC &
,HAILNCV
INTEGER,OPTIONAL,INTENT(IN ) :: id
REAL , DIMENSION( ims:ime , jms:jme ) , OPTIONAL , &
INTENT(IN) :: ht
REAL, DIMENSION (:), OPTIONAL, INTENT(INOUT) :: mp_restart_state &
,tbpvs_state,tbpvs0_state
LOGICAL, OPTIONAL :: f_qv,f_qc,f_qr,f_qi,f_qs,f_qg,f_qndrop,f_qni,f_qt &
,f_qns,f_qnr,f_qng,f_qnn,f_qnc,f_qh,f_qnh,f_qzr &
,f_effr,f_ice_effr,f_tot_effr &
,f_qic_effr,f_qip_effr,f_qid_effr &
,f_qic,f_qip,f_qid &
,f_qnic,f_qnip,f_qnid &
,f_qzi,f_qzs,f_qzg,f_qzh,f_qvolg,f_qvolh &
,f_qrimef &
,f_qir,f_qib &
,f_qnwfa, f_qnifa
LOGICAL, OPTIONAL, INTENT(IN) :: diagflag
REAL, INTENT(IN) :: ccn_conc
INTEGER, OPTIONAL, INTENT(IN) :: do_radar_ref
REAL, DIMENSION(ims:ime,kms:kme,jms:jme), INTENT(INOUT) :: &
re_cloud, re_ice, re_snow
INTEGER, INTENT(IN):: has_reqc, has_reqi, has_reqs
INTEGER :: i,j,k,its,ite,jts,jte,ij,sz,n
LOGICAL :: channel
LOGICAL :: nssl_progn = .false.
REAL :: z0, z1, z2, w1, w2
integer, parameter :: ntot = 50
real :: wmin, wmax
integer :: ierr
channel = .FALSE.
IF ( PRESENT ( channel_switch ) ) channel = channel_switch
if (mp_physics .eq. 0) return
IF( specified ) THEN
sz = spec_zone
ELSE
sz = 0
ENDIF
IF ( .false. ) THEN
wmax = maxval( w(ips:ipe,kps:kpe,jps:jpe) )
wmin = minval( w(ips:ipe,kps:kpe,jps:jpe) )
wmax = wrf_dm_max_real ( wmax )
wmin = wrf_dm_min_real ( wmin )
WRITE( wrf_err_message , * ) 'microphysics_driver: GLOBAL w max/min = ', wmax, wmin
CALL wrf_message ( wrf_err_message )
ENDIF
!$OMP PARALLEL DO &
!$OMP PRIVATE ( ij, its, ite, jts, jte, i,j,k,n )
DO ij = 1 , num_tiles
IF (channel) THEN
its = max(i_start(ij),ids)
ite = min(i_end(ij),ide-1)
ELSE
its = max(i_start(ij),ids+sz)
ite = min(i_end(ij),ide-1-sz)
ENDIF
jts = max(j_start(ij),jds+sz)
jte = min(j_end(ij),jde-1-sz)
IF( PRESENT(qlsink) ) qlsink(its:ite,kts:kte,jts:jte) = 0.
IF( PRESENT(precr ) ) precr(its:ite,kts:kte,jts:jte) = 0.
IF( PRESENT(preci ) ) preci(its:ite,kts:kte,jts:jte) = 0.
IF( PRESENT(precs ) ) precs(its:ite,kts:kte,jts:jte) = 0.
IF( PRESENT(precg ) ) precg(its:ite,kts:kte,jts:jte) = 0.
IF( PRESENT(chem_opt) .AND. PRESENT(progn) ) THEN
IF ( mp_physics==NSSL_2MOMCCN .or. mp_physics==NSSL_2MOM .or. mp_physics==NSSL_2MOMG ) THEN
IF ( progn > 0 ) THEN
IF ( .not. (chem_opt == 0 .or. chem_opt == 401) ) nssl_progn = .true.
ELSE
nssl_progn = .false.
ENDIF
ENDIF
IF( (chem_opt==0 .OR. chem_opt==401) .AND. progn==1 .AND. (mp_physics==LINSCHEME .OR. mp_physics==MORR_TWO_MOMENT)) THEN
IF( PRESENT( QNDROP_CURR ) ) THEN
CALL wrf_debug ( 100 , 'microphysics_driver: calling prescribe_aerosol_mixactivate' )
call prescribe_aerosol_mixactivate ( &
id, itimestep, dt, naer, &
ccn_conc, chem_opt, &
rho, th, pi_phy, w, cldfra, cldfra_old, &
z, dz8w, p8w, t8w, exch_h, &
qv_curr, qc_curr, qi_curr, qndrop_curr, &
nsource, &
ids,ide, jds,jde, kds,kde, &
ims,ime, jms,jme, kms,kme, &
its,ite, jts,jte, kts,kte, &
F_QC=f_qc, F_QI=f_qi )
END IF
ELSEIF ( (chem_opt==0 .OR. chem_opt==401) .AND. progn==1 .AND. (mp_physics==NSSL_2MOMCCN .or. &
mp_physics==NSSL_2MOM .or. mp_physics==NSSL_2MOMG)) THEN
ELSEIF ( progn==1 .AND. mp_physics/=LINSCHEME .AND. mp_physics/=MORR_TWO_MOMENT &
.AND. mp_physics/=NSSL_2MOM .AND. mp_physics/=NSSL_2MOMCCN .AND. mp_physics/=NSSL_2MOMG ) THEN
call wrf_error_fatal3("<stdin>",585,&
"SETTINGS ERROR: Prognostic cloud droplet number can only be used with the mp_physics=LINSCHEME or MORRISON or NSSL_2MOM.")
END IF
END IF
micro_select: SELECT CASE(mp_physics)
CASE (KESSLERSCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling kessler' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) .AND. &
PRESENT( Z )) THEN
CALL kessler( &
T=th &
,QV=qv_curr &
,QC=qc_curr &
,QR=qr_curr &
,RHO=rho, PII=pi_phy,DT_IN=dt, Z=z, XLV=xlv, CP=cp &
,EP2=ep_2,SVP1=svp1,SVP2=svp2 &
,SVP3=svp3,SVPT0=svpt0,RHOWATER=rhowater &
,DZ8W=dz8w &
,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_fatal3("<stdin>",613,&
'arguments not present for calling kessler' )
ENDIF
CASE (THOMPSONAERO)
CALL wrf_debug ( 100 , 'microphysics_driver: calling thompson' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. PRESENT ( QI_CURR ) .AND. &
PRESENT( QS_CURR ) .AND. PRESENT ( QG_CURR ) .AND. &
PRESENT( QNR_CURR) .AND. PRESENT ( QNI_CURR) .AND. &
PRESENT( QNC_CURR) .AND. PRESENT ( QNWFA_CURR) .AND. &
PRESENT( QNIFA_CURR).AND.PRESENT ( QNWFA2D) .AND. &
PRESENT( SNOWNC) .AND. PRESENT ( SNOWNCV) .AND. &
PRESENT( GRAUPELNC).AND. PRESENT ( GRAUPELNCV) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) ) THEN
CALL mp_gt_driver( &
QV=qv_curr, &
QC=qc_curr, &
QR=qr_curr, &
QI=qi_curr, &
QS=qs_curr, &
QG=qg_curr, &
NI=qni_curr, &
NR=qnr_curr, &
NC=qnc_curr, &
NWFA=qnwfa_curr, &
NIFA=qnifa_curr, &
NWFA2D=qnwfa2d, &
TH=th, &
PII=pi_phy, &
P=p, &
W=w, &
DZ=dz8w, &
DT_IN=dt, &
ITIMESTEP=itimestep, &
RAINNC=RAINNC, &
RAINNCV=RAINNCV, &
SNOWNC=SNOWNC, &
SNOWNCV=SNOWNCV, &
GRAUPELNC=GRAUPELNC, &
GRAUPELNCV=GRAUPELNCV, &
SR=SR, &
REFL_10CM=refl_10cm, &
diagflag=diagflag, &
do_radar_ref=do_radar_ref, &
re_cloud=re_cloud, &
re_ice=re_ice, &
re_snow=re_snow, &
has_reqc=has_reqc, &
has_reqi=has_reqi, &
has_reqs=has_reqs, &
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_fatal3("<stdin>",669,&
'arguments not present for calling thompson_et_al' )
ENDIF
CASE (THOMPSON)
CALL wrf_debug ( 100 , 'microphysics_driver: calling thompson' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. PRESENT ( QI_CURR ) .AND. &
PRESENT( QS_CURR ) .AND. PRESENT ( QG_CURR ) .AND. &
PRESENT( QNR_CURR) .AND. PRESENT ( QNI_CURR) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) ) THEN
CALL mp_gt_driver( &
QV=qv_curr, &
QC=qc_curr, &
QR=qr_curr, &
QI=qi_curr, &
QS=qs_curr, &
QG=qg_curr, &
NI=qni_curr, &
NR=qnr_curr, &
TH=th, &
PII=pi_phy, &
P=p, &
W=w, &
DZ=dz8w, &
DT_IN=dt, &
ITIMESTEP=itimestep, &
RAINNC=RAINNC, &
RAINNCV=RAINNCV, &
SNOWNC=SNOWNC, &
SNOWNCV=SNOWNCV, &
GRAUPELNC=GRAUPELNC, &
GRAUPELNCV=GRAUPELNCV, &
SR=SR, &
REFL_10CM=refl_10cm, &
diagflag=diagflag, &
do_radar_ref=do_radar_ref, &
re_cloud=re_cloud, &
re_ice=re_ice, &
re_snow=re_snow, &
has_reqc=has_reqc, &
has_reqi=has_reqi, &
has_reqs=has_reqs, &
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_fatal3("<stdin>",718,&
'arguments not present for calling thompson_et_al' )
ENDIF
CASE (MORR_TWO_MOMENT)
CALL wrf_debug(100, 'microphysics_driver: calling morrison two moment')
IF (PRESENT (QV_CURR) .AND. PRESENT (QC_CURR) .AND. &
PRESENT (QR_CURR) .AND. PRESENT (QI_CURR) .AND. &
PRESENT (QS_CURR) .AND. PRESENT (QG_CURR) .AND. &
PRESENT (QR_CURR) .AND. PRESENT (QI_CURR) .AND. &
PRESENT (QNS_CURR) .AND. PRESENT (QNI_CURR).AND. &
PRESENT (QNR_CURR) .AND. PRESENT (QNG_CURR).AND. &
PRESENT (QSCUTEN).AND. &
PRESENT (QRCUTEN) .AND. PRESENT (QICUTEN).AND. &
PRESENT (RAINNC ) .AND. PRESENT (RAINNCV) .AND. &
PRESENT ( W ) ) THEN
CALL mp_morr_two_moment( &
ITIMESTEP=itimestep, &
TH=th, &
QV=qv_curr, &
QC=qc_curr, &
QR=qr_curr, &
QI=qi_curr, &
QS=qs_curr, &
QG=qg_curr, &
NI=qni_curr, &
NS=qns_curr, &
NR=qnr_curr, &
NG=qng_curr, &
RHO=rho, &
PII=pi_phy, &
P=p, &
DT_IN=dt, &
DZ=dz8w, &
HT=ht, &
W=w &
,RAINNC=RAINNC &
,RAINNCV=RAINNCV &
,SNOWNC=SNOWNC &
,SNOWNCV=SNOWNCV &
,GRAUPELNC=GRAUPELNC &
,GRAUPELNCV=GRAUPELNCV &
,SR=SR &
,REFL_10CM=refl_10cm &
,diagflag=diagflag &
,do_radar_ref=do_radar_ref &
,qrcuten=qrcuten &
,qscuten=qscuten &
,qicuten=qicuten &
,F_QNDROP=f_qndrop &
,QNDROP=qndrop_curr &
,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 &
,QLSINK=qlsink &
,PRECR=precr,PRECI=preci,PRECS=precs,PRECG=precg &
)
ELSE
Call wrf_error_fatal3("<stdin>",778,&
'arguments not present for calling morrison two moment')
ENDIF
CASE (MILBRANDT2MOM)
CALL wrf_debug(100, 'microphysics_driver: calling milbrandt2mom')
IF (PRESENT (QV_CURR) .AND. &
PRESENT (QC_CURR) .AND. PRESENT (QNC_CURR) .AND. &
PRESENT (QR_CURR) .AND. PRESENT (QNR_CURR) .AND. &
PRESENT (QI_CURR) .AND. PRESENT (QNI_CURR) .AND. &
PRESENT (QS_CURR) .AND. PRESENT (QNS_CURR) .AND. &
PRESENT (QG_CURR) .AND. PRESENT (QNG_CURR) .AND. &
PRESENT (QH_CURR) .AND. PRESENT (QNH_CURR) .AND. &
PRESENT (RAINNC ) .AND. PRESENT (RAINNCV) .AND. &
PRESENT (SNOWNC ) .AND. PRESENT (SNOWNCV) .AND. &
PRESENT (HAILNC ) .AND. PRESENT (HAILNCV) .AND. &
PRESENT (GRAUPELNC).AND.PRESENT (GRAUPELNCV).AND. &
PRESENT ( W ) ) THEN
CALL mp_milbrandt2mom_driver( &
ITIMESTEP=itimestep, &
p8w=p8w, &
TH=th, &
QV=qv_curr, &
QC=qc_curr, &
QR=qr_curr, &
QI=qi_curr, &
QS=qs_curr, &
QG=qg_curr, &
QH=qh_curr, &
NC=qnc_curr, &
NR=qnr_curr, &
NI=qni_curr, &
NS=qns_curr, &
NG=qng_curr, &
NH=qnh_curr, &
PII=pi_phy, &
P=p, &
DT_IN=dt, &
DZ=dz8w, &
W=w, &
RAINNC = RAINNC, &
RAINNCV = RAINNCV, &
SNOWNC = SNOWNC, &
SNOWNCV = SNOWNCV, &
HAILNC = HAILNC, &
HAILNCV = HAILNCV, &
GRPLNC = GRAUPELNC, &
GRPLNCV = GRAUPELNCV, &
SR=SR, &
Zet = refl_10cm, &
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_fatal3("<stdin>",837,&
'arguments not present for calling milbrandt2mom')
ENDIF
CASE (NSSL_1MOM)
CALL wrf_debug(100, 'microphysics_driver: calling nssl1mom')
IF (PRESENT (QV_CURR) .AND. &
PRESENT (QC_CURR) .AND. &
PRESENT (QR_CURR) .AND. &
PRESENT (QI_CURR) .AND. &
PRESENT (QS_CURR) .AND. &
PRESENT (QG_CURR) .AND. &
PRESENT (QH_CURR) .AND. &
PRESENT (RAINNC ) .AND. PRESENT (RAINNCV) .AND. &
PRESENT ( W ) .AND. &
PRESENT (QVOLG_CURR) ) THEN
CALL nssl_2mom_driver( &
ITIMESTEP=itimestep, &
TH=th, &
QV=qv_curr, &
QC=qc_curr, &
QR=qr_curr, &
QI=qi_curr, &
QS=qs_curr, &
QH=qg_curr, &
QHL=qh_curr, &
VHW=qvolg_curr, &
PII=pi_phy, &
P=p, &
W=w, &
DZ=dz8w, &
DTP=dt, &
DN=rho, &
RAINNC = RAINNC, &
RAINNCV = RAINNCV, &
SNOWNC = SNOWNC, &
SNOWNCV = SNOWNCV, &
HAILNC = HAILNC, &
HAILNCV = HAILNCV, &
GRPLNC = GRAUPELNC, &
GRPLNCV = GRAUPELNCV, &
SR=SR, &
dbz = refl_10cm, &
diagflag = diagflag, &
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_fatal3("<stdin>",943,&
'arguments not present for calling nssl_1mom')
ENDIF
CASE (NSSL_1MOMLFO)
CALL wrf_debug(100, 'microphysics_driver: calling nssl1mom')
IF (PRESENT (QV_CURR) .AND. &
PRESENT (QC_CURR) .AND. &
PRESENT (QR_CURR) .AND. &
PRESENT (QI_CURR) .AND. &
PRESENT (QS_CURR) .AND. &
PRESENT (QG_CURR) .AND. &
PRESENT (RAINNC ) .AND. PRESENT (RAINNCV) .AND. &
PRESENT ( W ) ) THEN
CALL nssl_2mom_driver( &
ITIMESTEP=itimestep, &
TH=th, &
QV=qv_curr, &
QC=qc_curr, &
QR=qr_curr, &
QI=qi_curr, &
QS=qs_curr, &
QH=qg_curr, &
PII=pi_phy, &
P=p, &
W=w, &
DZ=dz8w, &
DTP=dt, &
DN=rho, &
RAINNC = RAINNC, &
RAINNCV = RAINNCV, &
SNOWNC = SNOWNC, &
SNOWNCV = SNOWNCV, &
GRPLNC = GRAUPELNC, &
GRPLNCV = GRAUPELNCV, &
SR=SR, &
dbz = refl_10cm, &
diagflag = diagflag, &
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_fatal3("<stdin>",989,&
'arguments not present for calling nssl_1momlfo')
ENDIF
CASE (NSSL_2MOM)
CALL wrf_debug(100, 'microphysics_driver: calling nssl2mom')
IF (PRESENT (QV_CURR) .AND. &
PRESENT (QC_CURR) .AND. PRESENT (QNdrop_CURR) .AND. &
PRESENT (QR_CURR) .AND. PRESENT (QNR_CURR) .AND. &
PRESENT (QI_CURR) .AND. PRESENT (QNI_CURR) .AND. &
PRESENT (QS_CURR) .AND. PRESENT (QNS_CURR) .AND. &
PRESENT (QG_CURR) .AND. PRESENT (QNG_CURR) .AND. &
PRESENT (QH_CURR) .AND. PRESENT (QNH_CURR) .AND. &
PRESENT (RAINNC ) .AND. PRESENT (RAINNCV) .AND. &
PRESENT ( W ) .AND. &
PRESENT (QVOLG_CURR) .AND. F_QVOLG .AND. &
PRESENT (QVOLH_CURR) .AND. F_QVOLH ) THEN
CALL nssl_2mom_driver( &
ITIMESTEP=itimestep, &
TH=th, &
QV=qv_curr, &
QC=qc_curr, &
QR=qr_curr, &
QI=qi_curr, &
QS=qs_curr, &
QH=qg_curr, &
QHL=qh_curr, &
CCW=qndrop_curr, &
CRW=qnr_curr, &
CCI=qni_curr, &
CSW=qns_curr, &
CHW=qng_curr, &
CHL=qnh_curr, &
VHW=qvolg_curr, &
VHL=qvolh_curr, &
PII=pi_phy, &
P=p, &
W=w, &
DZ=dz8w, &
DTP=dt, &
DN=rho, &
RAINNC = RAINNC, &
RAINNCV = RAINNCV, &
SNOWNC = SNOWNC, &
SNOWNCV = SNOWNCV, &
HAILNC = HAILNC, &
HAILNCV = HAILNCV, &
GRPLNC = GRAUPELNC, &
GRPLNCV = GRAUPELNCV, &
SR=SR, &
dbz = refl_10cm, &
nssl_progn=nssl_progn, &
diagflag = diagflag, &
cu_used=cu_used, &
qrcuten=qrcuten, &
qscuten=qscuten, &
qicuten=qicuten, &
qccuten=qccuten, &
re_cloud=re_cloud, &
re_ice=re_ice, &
re_snow=re_snow, &
has_reqc=has_reqc, &
has_reqi=has_reqi, &
has_reqs=has_reqs, &
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_fatal3("<stdin>",1062,&
'arguments not present for calling nssl_2mom')
ENDIF
CASE (NSSL_2MOMG)
CALL wrf_debug(100, 'microphysics_driver: calling nssl2mom')
IF (PRESENT (QV_CURR) .AND. &
PRESENT (QC_CURR) .AND. PRESENT (QNdrop_CURR) .AND. &
PRESENT (QR_CURR) .AND. PRESENT (QNR_CURR) .AND. &
PRESENT (QI_CURR) .AND. PRESENT (QNI_CURR) .AND. &
PRESENT (QS_CURR) .AND. PRESENT (QNS_CURR) .AND. &
PRESENT (QG_CURR) .AND. PRESENT (QNG_CURR) .AND. &
PRESENT (RAINNC ) .AND. PRESENT (RAINNCV) .AND. &
PRESENT ( W ) .AND. &
PRESENT (QVOLG_CURR) .AND. F_QVOLG ) THEN
CALL nssl_2mom_driver( &
ITIMESTEP=itimestep, &
TH=th, &
QV=qv_curr, &
QC=qc_curr, &
QR=qr_curr, &
QI=qi_curr, &
QS=qs_curr, &
QH=qg_curr, &
CCW=qndrop_curr, &
CRW=qnr_curr, &
CCI=qni_curr, &
CSW=qns_curr, &
CHW=qng_curr, &
VHW=qvolg_curr, &
PII=pi_phy, &
P=p, &
W=w, &
DZ=dz8w, &
DTP=dt, &
DN=rho, &
RAINNC = RAINNC, &
RAINNCV = RAINNCV, &
SNOWNC = SNOWNC, &
SNOWNCV = SNOWNCV, &
HAILNC = HAILNC, &
HAILNCV = HAILNCV, &
GRPLNC = GRAUPELNC, &
GRPLNCV = GRAUPELNCV, &
SR=SR, &
dbz = refl_10cm, &
nssl_progn=nssl_progn, &
diagflag = diagflag, &
cu_used=cu_used, &
qrcuten=qrcuten, &
qscuten=qscuten, &
qicuten=qicuten, &
qccuten=qccuten, &
re_cloud=re_cloud, &
re_ice=re_ice, &
re_snow=re_snow, &
has_reqc=has_reqc, &
has_reqi=has_reqi, &
has_reqs=has_reqs, &
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_fatal3("<stdin>",1130,&
'arguments not present for calling nssl_2momg')
ENDIF
CASE (NSSL_2MOMCCN)
CALL wrf_debug(100, 'microphysics_driver: calling nssl_2momccn')
IF (PRESENT (QV_CURR) .AND. &
PRESENT (QC_CURR) .AND. PRESENT (QNDROP_CURR) .AND. &
PRESENT (QR_CURR) .AND. PRESENT (QNR_CURR) .AND. &
PRESENT (QI_CURR) .AND. PRESENT (QNI_CURR) .AND. &
PRESENT (QS_CURR) .AND. PRESENT (QNS_CURR) .AND. &
PRESENT (QG_CURR) .AND. PRESENT (QNG_CURR) .AND. &
PRESENT (QH_CURR) .AND. PRESENT (QNH_CURR) .AND. &
PRESENT (RAINNC ) .AND. PRESENT (RAINNCV) .AND. &
PRESENT ( W ) .AND. &
PRESENT (QVOLG_CURR) .AND. F_QVOLG .AND. &
PRESENT (QVOLH_CURR) .AND. F_QVOLH .AND. &
PRESENT( QNN_CURR ) ) THEN
CALL nssl_2mom_driver( &
ITIMESTEP=itimestep, &
TH=th, &
QV=qv_curr, &
QC=qc_curr, &
QR=qr_curr, &
QI=qi_curr, &
QS=qs_curr, &
QH=qg_curr, &
QHL=qh_curr, &
CCW=qndrop_curr, &
CRW=qnr_curr, &
CCI=qni_curr, &
CSW=qns_curr, &
CHW=qng_curr, &
CHL=qnh_curr, &
VHW=qvolg_curr, &
VHL=qvolh_curr, &
cn=qnn_curr, &
PII=pi_phy, &
P=p, &
W=w, &
DZ=dz8w, &
DTP=dt, &
DN=rho, &
RAINNC = RAINNC, &
RAINNCV = RAINNCV, &
SNOWNC = SNOWNC, &
SNOWNCV = SNOWNCV, &
HAILNC = HAILNC, &
HAILNCV = HAILNCV, &
GRPLNC = GRAUPELNC, &
GRPLNCV = GRAUPELNCV, &
SR=SR, &
dbz = refl_10cm, &
nssl_progn=nssl_progn, &
diagflag = diagflag, &
cu_used=cu_used, &
qrcuten=qrcuten, &
qscuten=qscuten, &
qicuten=qicuten, &
qccuten=qccuten, &
re_cloud=re_cloud, &
re_ice=re_ice, &
re_snow=re_snow, &
has_reqc=has_reqc, &
has_reqi=has_reqi, &
has_reqs=has_reqs, &
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_fatal3("<stdin>",1204,&
'arguments not present for calling nssl_2momccn')
ENDIF
CASE (GSFCGCESCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling GSFCGCE' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. PRESENT ( QI_CURR ) .AND. &
PRESENT( QS_CURR ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) .AND. &
PRESENT( HAIL ) .AND. PRESENT ( ICE2 ) .AND. &
PRESENT( Z ) .AND. PRESENT ( W ) ) THEN
CALL gsfcgce( &
TH=th &
,QV=qv_curr &
,QL=qc_curr &
,QR=qr_curr &
,QI=qi_curr &
,QS=qs_curr &
,RHO=rho, PII=pi_phy, P=p, DT_IN=dt, Z=z &
,HT=ht, DZ8W=dz8w, GRAV=G &
,RHOWATER=rhowater, RHOSNOW=rhosnow &
,ITIMESTEP=itimestep &
,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 &
,RAINNC=rainnc, RAINNCV=rainncv &
,SNOWNC=snownc, SNOWNCV=snowncv ,SR=sr &
,GRAUPELNC=graupelnc ,GRAUPELNCV=graupelncv &
,REFL_10CM=refl_10cm &
,diagflag=diagflag &
,do_radar_ref=do_radar_ref &
,F_QG=f_qg &
,QG=qg_curr &
,IHAIL=hail, ICE2=ice2 &
)
ELSE
CALL wrf_error_fatal3("<stdin>",1248,&
'arguments not present for calling GSFCGCE' )
ENDIF
CASE (LINSCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling lin_et_al' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. PRESENT ( QI_CURR ) .AND. &
PRESENT( QS_CURR ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) .AND. &
PRESENT( Z ) ) THEN
CALL lin_et_al( &
TH=th &
,QV=qv_curr &
,QL=qc_curr &
,QR=qr_curr &
,QI=qi_curr &
,QS=qs_curr &
,QLSINK=qlsink &
,RHO=rho, PII=pi_phy, P=p, DT_IN=dt, Z=z &
,HT=ht, DZ8W=dz8w, GRAV=G, CP=cp &
,RAIR=r_d, RVAPOR=R_v &
,XLS=xls, XLV=xlv, XLF=xlf &
,RHOWATER=rhowater, RHOSNOW=rhosnow &
,EP2=ep_2,SVP1=svp1,SVP2=svp2 &
,SVP3=svp3,SVPT0=svpt0 &
,RAINNC=rainnc, RAINNCV=rainncv &
,SNOWNC=snownc, SNOWNCV=snowncv &
,GRAUPELNC=graupelnc, GRAUPELNCV=graupelncv, SR=sr &
,REFL_10CM=refl_10cm &
,diagflag=diagflag &
,do_radar_ref=do_radar_ref &
,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 &
,PRECR=precr,PRECI=preci,PRECS=precs,PRECG=precg &
,F_QG=f_qg, F_QNDROP=f_qndrop &
,QG=qg_curr &
,QNDROP=qndrop_curr &
)
ELSE
CALL wrf_error_fatal3("<stdin>",1289,&
'arguments not present for calling lin_et_al' )
ENDIF
CASE (SBU_YLINSCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling sbu_ylin' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. PRESENT ( QI_CURR ) .AND. &
PRESENT( QS_CURR ) .AND. &
PRESENT( RI_CURR ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) .AND. &
PRESENT( Z ) ) THEN
CALL sbu_ylin( &
TH=th &
,QV=qv_curr &
,QL=qc_curr &
,QR=qr_curr &
,QI=qi_curr &
,QS=qs_curr &
,RI3D=ri_curr &
,RHO=rho, PII=pi_phy, P=p, DT_IN=dt, Z=z &
,HT=ht, DZ8W=dz8w &
,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_fatal3("<stdin>",1331,&
'arguments not present for calling sbu_ylin' )
ENDIF
CASE (WSM3SCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling wsm3' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) .AND. &
PRESENT( W ) ) THEN
CALL wsm3( &
TH=th &
,Q=qv_curr &
,QCI=qc_curr &
,QRS=qr_curr &
,W=w,DEN=rho,PII=pi_phy,P=p,DELZ=dz8w &
,DELT=dt,G=g,CPD=cp,CPV=cpv &
,RD=r_d,RV=r_v,T0C=svpt0 &
,EP1=ep_1, EP2=ep_2, QMIN=epsilon &
,XLS=xls, XLV0=xlv, XLF0=xlf &
,DEN0=rhoair0, DENR=rhowater &
,CLIQ=cliq,CICE=cice,PSAT=psat &
,RAIN=rainnc ,RAINNCV=rainncv &
,SNOW=snownc ,SNOWNCV=snowncv &
,SR=sr &
,has_reqc=has_reqc &
,has_reqi=has_reqi &
,has_reqs=has_reqs &
,re_cloud=re_cloud &
,re_ice=re_ice &
,re_snow=re_snow &
,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_fatal3("<stdin>",1368,&
'arguments not present for calling wsm3' )
ENDIF
CASE (WSM5SCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling wsm5' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. PRESENT ( QI_CURR ) .AND. &
PRESENT( QS_CURR ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) ) THEN
CALL wsm5( &
TH=th &
,Q=qv_curr &
,QC=qc_curr &
,QR=qr_curr &
,QI=qi_curr &
,QS=qs_curr &
,DEN=rho,PII=pi_phy,P=p,DELZ=dz8w &
,DELT=dt,G=g,CPD=cp,CPV=cpv &
,RD=r_d,RV=r_v,T0C=svpt0 &
,EP1=ep_1, EP2=ep_2, QMIN=epsilon &
,XLS=xls, XLV0=xlv, XLF0=xlf &
,DEN0=rhoair0, DENR=rhowater &
,CLIQ=cliq,CICE=cice,PSAT=psat &
,RAIN=rainnc ,RAINNCV=rainncv &
,SNOW=snownc ,SNOWNCV=snowncv &
,SR=sr &
,REFL_10CM=refl_10cm &
,diagflag=diagflag &
,do_radar_ref=do_radar_ref &
,has_reqc=has_reqc &
,has_reqi=has_reqi &
,has_reqs=has_reqs &
,re_cloud=re_cloud &
,re_ice=re_ice &
,re_snow=re_snow &
,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_fatal3("<stdin>",1409,&
'arguments not present for calling wsm5' )
ENDIF
CASE (WSM6SCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling wsm6' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. PRESENT ( QI_CURR ) .AND. &
PRESENT( QS_CURR ) .AND. PRESENT ( QG_CURR ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) ) THEN
CALL wsm6( &
TH=th &
,Q=qv_curr &
,QC=qc_curr &
,QR=qr_curr &
,QI=qi_curr &
,QS=qs_curr &
,QG=qg_curr &
,DEN=rho,PII=pi_phy,P=p,DELZ=dz8w &
,DELT=dt,G=g,CPD=cp,CPV=cpv &
,RD=r_d,RV=r_v,T0C=svpt0 &
,EP1=ep_1, EP2=ep_2, QMIN=epsilon &
,XLS=xls, XLV0=xlv, XLF0=xlf &
,DEN0=rhoair0, DENR=rhowater &
,CLIQ=cliq,CICE=cice,PSAT=psat &
,RAIN=rainnc ,RAINNCV=rainncv &
,SNOW=snownc ,SNOWNCV=snowncv &
,SR=sr &
,REFL_10CM=refl_10cm &
,diagflag=diagflag &
,do_radar_ref=do_radar_ref &
,GRAUPEL=graupelnc ,GRAUPELNCV=graupelncv &
,has_reqc=has_reqc &
,has_reqi=has_reqi &
,has_reqs=has_reqs &
,re_cloud=re_cloud &
,re_ice=re_ice &
,re_snow=re_snow &
,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_fatal3("<stdin>",1452,&
'arguments not present for calling wsm6' )
ENDIF
CASE (WDM5SCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling wdm5' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. PRESENT ( QI_CURR ) .AND. &
PRESENT( QS_CURR ) .AND. PRESENT( QNN_CURR ) .AND. &
PRESENT ( QNC_CURR ) .AND. PRESENT( QNR_CURR ).AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) ) THEN
CALL wdm5( &
TH=th &
,Q=qv_curr &
,QC=qc_curr &
,QR=qr_curr &
,QI=qi_curr &
,QS=qs_curr &
,NN=qnn_curr &
,NC=qnc_curr &
,NR=qnr_curr &
,DEN=rho,PII=pi_phy,P=p,DELZ=dz8w &
,DELT=dt,G=g,CPD=cp,CPV=cpv,CCN0=ccn_conc &
,RD=r_d,RV=r_v,T0C=svpt0 &
,EP1=ep_1, EP2=ep_2, QMIN=epsilon &
,XLS=xls, XLV0=xlv, XLF0=xlf &
,DEN0=rhoair0, DENR=rhowater &
,CLIQ=cliq,CICE=cice,PSAT=psat &
,RAIN=rainnc ,RAINNCV=rainncv &
,SNOW=snownc ,SNOWNCV=snowncv &
,SR=sr &
,REFL_10CM=refl_10cm &
,diagflag=diagflag &
,do_radar_ref=do_radar_ref &
,has_reqc=has_reqc &
,has_reqi=has_reqi &
,has_reqs=has_reqs &
,re_cloud=re_cloud &
,re_ice=re_ice &
,re_snow=re_snow &
,ITIMESTEP=itimestep &
,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_fatal3("<stdin>",1498,&
'arguments not present for calling wdm5')
ENDIF
CASE (WDM6SCHEME)
CALL wrf_debug ( 100 , 'microphysics_driver: calling wdm6' )
IF ( PRESENT( QV_CURR ) .AND. PRESENT ( QC_CURR ) .AND. &
PRESENT( QR_CURR ) .AND. PRESENT ( QI_CURR ) .AND. &
PRESENT( QS_CURR ) .AND. PRESENT ( QG_CURR ) .AND. &
PRESENT( QNN_CURR ) .AND. PRESENT ( QNC_CURR ) .AND. &
PRESENT( QNR_CURR ).AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) ) THEN
CALL wdm6( &
TH=th &
,Q=qv_curr &
,QC=qc_curr &
,QR=qr_curr &
,QI=qi_curr &
,QS=qs_curr &
,QG=qg_curr &
,NN=qnn_curr &
,NC=qnc_curr &
,NR=qnr_curr &
,DEN=rho,PII=pi_phy,P=p,DELZ=dz8w &
,DELT=dt,G=g,CPD=cp,CPV=cpv,CCN0=ccn_conc &
,RD=r_d,RV=r_v,T0C=svpt0 &
,EP1=ep_1, EP2=ep_2, QMIN=epsilon &
,XLS=xls, XLV0=xlv, XLF0=xlf &
,DEN0=rhoair0, DENR=rhowater &
,CLIQ=cliq,CICE=cice,PSAT=psat &
,RAIN=rainnc ,RAINNCV=rainncv &
,SNOW=snownc ,SNOWNCV=snowncv &
,SR=sr &
,REFL_10CM=refl_10cm &
,diagflag=diagflag &
,do_radar_ref=do_radar_ref &
,GRAUPEL=graupelnc ,GRAUPELNCV=graupelncv &
,ITIMESTEP=itimestep &
,has_reqc=has_reqc &
,has_reqi=has_reqi &
,has_reqs=has_reqs &
,re_cloud=re_cloud &
,re_ice=re_ice &
,re_snow=re_snow &
,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_fatal3("<stdin>",1547,&
'arguments not present for calling wdm6')
ENDIF
CASE (ETAMP_HWRF)
CALL wrf_debug ( 100 , 'microphysics_driver: calling etampnew_HWRF')
IF ( PRESENT( qv_curr ) .AND. PRESENT( qt_curr ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) .AND. &
PRESENT( mp_restart_state ) .AND. &
PRESENT( tbpvs_state ) .AND. &
PRESENT( tbpvs0_state ) ) THEN
CALL ETAMP_NEW_HWRF( &
ITIMESTEP=itimestep,DT=dt,DX=dx,DY=dy, GID=id &
,RAINNC=rainnc,RAINNCV=rainncv &
,DZ8W=dz8w,RHO_PHY=rho,P_PHY=p,PI_PHY=pi_phy,TH_PHY=th &
,QV=qv_curr &
,QT=qt_curr &
,LOWLYR=LOWLYR,SR=SR &
,F_ICE_PHY=F_ICE_PHY,F_RAIN_PHY=F_RAIN_PHY &
,F_RIMEF_PHY=F_RIMEF_PHY &
,QC=qc_curr,QR=Qr_curr,QI=Qi_curr &
,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_fatal3("<stdin>",1574,&
'arguments not present for calling etampnew' )
ENDIF
CASE (ETAMPNEW)
CALL wrf_debug ( 100 , 'microphysics_driver: calling etampnew')
IF ( PRESENT( qv_curr ) .AND. PRESENT( qt_curr ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) .AND. &
PRESENT( mp_restart_state ) .AND. &
PRESENT( tbpvs_state ) .AND. &
PRESENT( tbpvs0_state ) ) THEN
CALL ETAMP_NEW( &
ITIMESTEP=itimestep,DT=dt,DX=dx,DY=dy &
,DZ8W=dz8w,RHO_PHY=rho,P_PHY=p,PI_PHY=pi_phy,TH_PHY=th &
,QV=qv_curr &
,QC=qc_curr &
,QS=qs_curr &
,QR=qr_curr &
,QT=qt_curr &
,LOWLYR=LOWLYR,SR=SR &
,F_ICE_PHY=F_ICE_PHY,F_RAIN_PHY=F_RAIN_PHY &
,F_RIMEF_PHY=F_RIMEF_PHY &
,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 &
,MP_RESTART_STATE=mp_restart_state &
,TBPVS_STATE=tbpvs_state,TBPVS0_STATE=tbpvs0_state &
)
ELSE
CALL wrf_error_fatal3("<stdin>",1604,&
'arguments not present for calling etampnew' )
ENDIF
CASE (FER_MP_HIRES)
CALL wrf_debug ( 100 , 'microphysics_driver: calling etampnew')
IF ( PRESENT( qv_curr ) .AND. PRESENT( qt_curr ) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) .AND. &
PRESENT( mp_restart_state ) .AND. &
PRESENT( tbpvs_state ) .AND. &
PRESENT( tbpvs0_state ) ) THEN
CALL FER_HIRES( &
ITIMESTEP=itimestep,DT=dt,DX=dx,DY=dy, GID=id &
,RAINNC=rainnc,RAINNCV=rainncv &
,DZ8W=dz8w,RHO_PHY=rho,P_PHY=p,PI_PHY=pi_phy,TH_PHY=th &
,QV=qv_curr &
,QT=qt_curr &
,LOWLYR=LOWLYR,SR=SR &
,F_ICE_PHY=F_ICE_PHY,F_RAIN_PHY=F_RAIN_PHY &
,F_RIMEF_PHY=F_RIMEF_PHY &
,QC=qc_curr,QR=Qr_curr,QI=Qi_curr &
,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_fatal3("<stdin>",1637,&
'arguments not present for calling fer_hires' )
ENDIF
CASE (FER_MP_HIRES_ADVECT)
CALL wrf_debug ( 100 , 'microphysics_driver: calling etampnew')
IF ( PRESENT( qv_curr ) .AND. PRESENT( qi_curr ) .AND. &
PRESENT( qc_curr ) .and. PRESENT(qrimef_curr) .AND. &
PRESENT( RAINNC ) .AND. PRESENT ( RAINNCV ) .AND. &
PRESENT( mp_restart_state ) .AND. &
PRESENT( tbpvs_state ) .AND. &
PRESENT( tbpvs0_state ) ) THEN
CALL FER_HIRES_ADVECT( &
ITIMESTEP=itimestep,DT=dt,DX=dx,DY=dy, GID=id &
,RAINNC=rainnc,RAINNCV=rainncv &
,DZ8W=dz8w,RHO_PHY=rho,P_PHY=p,PI_PHY=pi_phy,TH_PHY=th &
,QV=qv_curr &
,LOWLYR=LOWLYR,SR=SR &
,QC=qc_curr,QR=Qr_curr,QI=Qi_curr,QRIMEF=qrimef_curr &
,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_fatal3("<stdin>",1670,&
'arguments not present for calling fer_hires' )
ENDIF
CASE DEFAULT
WRITE( wrf_err_message , * ) 'The microphysics option does not exist: mp_physics = ', mp_physics
CALL wrf_error_fatal3("<stdin>",1678,&
wrf_err_message )
END SELECT micro_select
ENDDO
!$OMP END PARALLEL DO
CALL wrf_debug ( 200 , 'microphysics_driver: returning from' )
RETURN
END SUBROUTINE microphysics_driver
END MODULE module_microphysics_driver