! Generated by TAPENADE (INRIA, Tropics team) ! Tapenade 3.6 (r4343) - 10 Feb 2012 10:52 ! !WRF:MODEL_LAYER:INITIALIZATION ! ! This MODULE holds the routines which are used to perform model start-up operations ! for the individual domains. This is the stage after inputting wrfinput and before ! calling 'integrate'. ! This MODULE CONTAINS the following routines: MODULE G_MODULE_PHYSICS_INIT USE module_state_description USE module_model_constants USE module_configure, ONLY : grid_config_rec_type USE module_physics_init, ONLY : landuse_init CONTAINS ! Differentiation of phy_init in forward (tangent) mode: ! variations of useful results: rublten raincv rqccuten rthcuten ! rainnc rainc rqvblten w0avg snownc rvblten ! graupelnc rthblten graupelncv snowncv rqvcuten ! rainncv ! with respect to varying inputs: rublten raincv rqccuten rthcuten ! rainnc rainc rqvblten w0avg snownc rvblten ! graupelnc rthblten graupelncv snowncv rqvcuten ! rainncv ! RW status of diff variables: rublten:in-out raincv:in-out rqccuten:in-out ! rthcuten:in-out rainnc:in-out rainc:in-out rqvblten:in-out ! w0avg:in-out snownc:in-out rvblten:in-out ! graupelnc:in-out rthblten:in-out graupelncv:in-out ! snowncv:in-out rqvcuten:in-out rainncv:in-out !================================================================= !EMBCK new ! temporary ! num_roof_layers,num_wall_layers, & ! num_road_layers, & ! WA 12/21/09 ! WA 9/27/10 ! WA 2/22/11 SUBROUTINE G_PHY_INIT(id, config_flags, dt, restart, zfull, zhalf, & & p_top, tsk, radt, bldt, cudt, mpdt, rucuten, rvcuten, rthcuten& & , rthcutend, rqvcuten, rqvcutend, rqrcuten, rqccuten, rqccutend, & & rqscuten, rqicuten, rushten, rvshten, rthshten, rqvshten, rqrshten, & & rqcshten, rqsshten, rqishten, rqgshten, rublten, rubltend, rvblten, & & rvbltend, rthblten, rthbltend, rqvblten, rqvbltend, rqcblten, & & rqiblten, rthraten, rthratenlw, rthratensw, stepbl, stepra, stepcu, & & w0avg, w0avgd, rainnc, rainncd, rainc, raincd, raincv, raincvd, & & rainncv, rainncvd, snownc, snowncd, snowncv, snowncvd, graupelnc, & & graupelncd, graupelncv, graupelncvd, nca, swrad_scat, cldefi, lowlyr& & , mass_flux, rthften, rqvften, cldfra, cldfra_old, glw, gsw, emiss, & & embck, lu_index, landuse_isice, landuse_lucats, landuse_luseas, & & landuse_isn, lu_state, xlat, xlong, albedo, albbck, gmt, julyr, & & julday, levsiz, n_ozmixm, n_aerosolc, paerlev, alevsiz, no_src_types, tmn, xland, znt, z0, & & ust, mol, pblh, tke_pbl, exch_h, thc, snowc, mavail, hfx, qfx, & & rainbl, tslb, zs, dzs, num_soil_layers, warm_rain, adv_moist_cond, & & apr_gr, apr_w, apr_mc, apr_st, apr_as, apr_capma, apr_capme, & & apr_capmi, xice, xicem, vegfra, snow, canwat, smstav, smstot, & & sfcrunoff, udrunoff, grdflx, acsnow, acsnom, ivgtyp, isltyp, sfcevp& & , smois, sh2o, snowh, smfr3d, snoalb, dx, dy, f_ice_phy, f_rain_phy& & , f_rimef_phy, mp_restart_state, tbpvs_state, tbpvs0_state, & & allowed_to_read, moved, start_of_simulation, lagday, ids, ide, jds, & & jde, kds, kde, ims, ime, jms, jme, kms, kme, its, ite, jts, jte, kts& & , kte, num_urban_layers, NUM_URBAN_HI, raincv_a, raincv_b, gd_cloud, gd_cloud2, & & gd_cloud_a, gd_cloud2_a, gd_cloud_b, gd_cloud2_b, ozmixm, pin, & & m_ps_1, m_ps_2, m_hybi, aerosolc_1, aerosolc_2, rundgdten, rvndgdten& & , rthndgdten, rphndgdten, rqvndgdten, rmundgdten, fgdt, stepfg, & & cugd_tten, cugd_ttens, cugd_qvten, cugd_qvtens, cugd_qcten, isnowxy& & , zsnsoxy, tsnoxy, snicexy, snliqxy, tvxy, tgxy, canicexy, canliqxy& & , eahxy, tahxy, cmxy, chxy, fwetxy, sneqvoxy, alboldxy, qsnowxy, & & wslakexy, zwtxy, waxy, wtxy, lfmassxy, rtmassxy, stmassxy, woodxy, & & stblcpxy, fastcpxy, xsaixy, t2mvxy, t2mbxy, chstarxy, dzr, dzb, dzg& & , tr_urb2d, tb_urb2d, tg_urb2d, tc_urb2d, qc_urb2d, xxxr_urb2d, & & xxxb_urb2d, xxxg_urb2d, xxxc_urb2d, trl_urb3d, tbl_urb3d, tgl_urb3d& & , sh_urb2d, lh_urb2d, g_urb2d, rn_urb2d, ts_urb2d, frc_urb2d, & & utype_urb2d, trb_urb4d, tw1_urb4d, tw2_urb4d, tgb_urb4d, tlev_urb3d& & , qlev_urb3d, tw1lev_urb3d, tw2lev_urb3d, tglev_urb3d, tflev_urb3d, & & sf_ac_urb3d, lf_ac_urb3d, cm_ac_urb3d, sfvent_urb3d, lfvent_urb3d, & & sfwin1_urb3d, sfwin2_urb3d, sfw1_urb3d, sfw2_urb3d, sfr_urb3d, & & sfg_urb3d, a_u_bep, a_v_bep, a_t_bep, a_q_bep, a_e_bep, b_u_bep, & & b_v_bep, b_t_bep, b_q_bep, b_e_bep, dlg_bep, dl_u_bep, sf_bep, & & vl_bep, tml, t0ml, hml, h0ml, huml, hvml, tmoml, itimestep, fdob, & & t00, p00, tlp, tyr, tyra, tdly, tlag, nyear, nday, tmn_update, achfx& & , aclhf, acgrdflx, te_temf, cf3d_temf, wm_temf, massflux_edkf, & & entr_edkf, detr_edkf, thl_up, thv_up, rt_up, rv_up, rc_up, u_up, & & v_up, frac_up) USE module_domain USE module_wrf_error IMPLICIT NONE !----------------------------------------------------------------- !----------------------------------------------------------------- TYPE(GRID_CONFIG_REC_TYPE) :: config_flags INTEGER, INTENT(IN) :: id INTEGER, INTENT(IN), OPTIONAL :: tmn_update LOGICAL, INTENT(OUT) :: warm_rain, adv_moist_cond ! LOGICAL , INTENT (IN) :: FNDSOILW, FNDSNOWH LOGICAL, PARAMETER :: fndsoilw=.true., fndsnowh=.true. INTEGER, INTENT(IN) :: ids, ide, jds, jde, kds, kde, ims, ime, jms, & & jme, kms, kme, its, ite, jts, jte, kts, kte INTEGER, INTENT(IN) :: num_soil_layers INTEGER, INTENT(IN) :: lagday INTEGER, INTENT(OUT), OPTIONAL :: nyear REAL, INTENT(OUT), OPTIONAL :: nday LOGICAL, INTENT(IN) :: start_of_simulation REAL, INTENT(IN) :: dt, p_top, dx, dy LOGICAL, INTENT(IN) :: restart REAL, INTENT(IN) :: radt, bldt, cudt, mpdt REAL, INTENT(IN) :: swrad_scat REAL, DIMENSION(kms:kme), INTENT(IN) :: zfull, zhalf REAL, DIMENSION(ims:ime, jms:jme), INTENT(IN) :: tsk, xlat, xlong INTEGER, INTENT(IN) :: levsiz, n_ozmixm INTEGER, INTENT(IN) :: paerlev, n_aerosolc INTEGER, INTENT(IN ) :: alevsiz, no_src_types REAL, DIMENSION(ims:ime, levsiz, jms:jme, n_ozmixm), OPTIONAL, & & INTENT(INOUT) :: ozmixm REAL, DIMENSION(levsiz), OPTIONAL, INTENT(INOUT) :: pin REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: m_ps_1& & , m_ps_2 REAL, DIMENSION(paerlev), OPTIONAL, INTENT(INOUT) :: m_hybi REAL, DIMENSION(ims:ime, paerlev, jms:jme, n_aerosolc), OPTIONAL, & & INTENT(INOUT) :: aerosolc_1, aerosolc_2 REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), INTENT(INOUT) ::& & smois, sh2o, tslb REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), INTENT(OUT) :: & & smfr3d REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: snow, snowc, & & snowh, canwat, smstav, smstot, sfcrunoff, udrunoff, sfcevp, grdflx, & & acsnow, xice, xicem, vegfra, acsnom REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: achfx& & , aclhf, acgrdflx INTEGER, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: ivgtyp, & & isltyp ! rad REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rthraten& & , rthratenlw, rthratensw, cldfra REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) ::& & cldfra_old !EMBCK new REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: gsw, albedo, & & albbck, glw, emiss, embck REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: snoalb REAL, INTENT(IN) :: gmt INTEGER, INTENT(OUT) :: stepra, stepbl, stepcu INTEGER, INTENT(IN) :: julyr, julday ! cps REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rucuten, & & rvcuten, rthcuten, rqvcuten, rqrcuten, rqccuten, rqscuten, rqicuten& & , rushten, rvshten, rthshten, rqvshten, rqrshten, rqcshten, rqsshten& & , rqishten, rqgshten REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rthcutend& & , rqvcutend, rqccutend REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(OUT) :: w0avg REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(OUT) :: w0avgd REAL, DIMENSION(ims:ime, jms:jme), INTENT(OUT) :: mass_flux, apr_gr& & , apr_w, apr_mc, apr_st, apr_as, apr_capma, apr_capme, apr_capmi REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rthften, & & rqvften REAL, DIMENSION(ims:ime, jms:jme), INTENT(OUT) :: rainnc, rainc, & & raincv, rainncv, snownc, snowncv, graupelnc, graupelncv REAL, DIMENSION(ims:ime, jms:jme), INTENT(OUT) :: rainncd, raincd, & & raincvd, rainncvd, snowncd, snowncvd, graupelncd, graupelncvd REAL, DIMENSION(ims:ime, jms:jme), INTENT(OUT) :: cldefi, nca INTEGER, DIMENSION(ims:ime, jms:jme), INTENT(OUT) :: lowlyr !pbl ! soil layer REAL, DIMENSION(num_soil_layers), INTENT(INOUT) :: zs, dzs REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rublten, & & rvblten, rthblten, rqvblten, rqcblten, rqiblten, exch_h, tke_pbl REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rubltend& & , rvbltend, rthbltend, rqvbltend REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(OUT), OPTIONAL ::& & massflux_edkf, entr_edkf, detr_edkf, thl_up, thv_up, rt_up, rv_up, & & rc_up, u_up, v_up, frac_up REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(OUT) ::& & cugd_tten, cugd_ttens, cugd_qvten, cugd_qvtens, cugd_qcten REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: xland, znt, z0, & & ust, mol, lu_index, pblh, thc, mavail, hfx, qfx, rainbl INTEGER, INTENT(INOUT) :: landuse_isice, landuse_lucats INTEGER, INTENT(INOUT) :: landuse_luseas, landuse_isn REAL, DIMENSION(:), INTENT(INOUT) :: lu_state REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: tmn REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT), OPTIONAL :: tyr REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT), OPTIONAL :: tyra REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT), OPTIONAL :: tdly REAL, DIMENSION(ims:ime, lagday, jms:jme), INTENT(INOUT), OPTIONAL & & :: tlag REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: gd_cloud, gd_cloud2, gd_cloud_a, gd_cloud2_a, gd_cloud_b, & & gd_cloud2_b REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: raincv_a, & & raincv_b !Noah-MP INTEGER, DIMENSION(ims:ime, jms:jme), OPTIONAL :: isnowxy REAL, DIMENSION(ims:ime, -2:num_soil_layers, jms:jme), OPTIONAL :: & & zsnsoxy REAL, DIMENSION(ims:ime, -2:0, jms:jme), OPTIONAL :: tsnoxy REAL, DIMENSION(ims:ime, -2:0, jms:jme), OPTIONAL :: snicexy REAL, DIMENSION(ims:ime, -2:0, jms:jme), OPTIONAL :: snliqxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: tvxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: tgxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: canicexy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: canliqxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: eahxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: tahxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: cmxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: chxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: fwetxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: sneqvoxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: alboldxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: qsnowxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: wslakexy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: zwtxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: waxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: wtxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: lfmassxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: rtmassxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: stmassxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: woodxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: stblcpxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: fastcpxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: xsaixy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: t2mvxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: t2mbxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: chstarxy !mp REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(OUT) :: f_ice_phy& & , f_rain_phy, f_rimef_phy REAL, DIMENSION(:), INTENT(INOUT) :: mp_restart_state, tbpvs_state, & & tbpvs0_state LOGICAL, INTENT(IN) :: allowed_to_read, moved ! ocean mixed layer REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: tml, & & t0ml, hml, h0ml, huml, hvml, tmoml !fdda REAL, OPTIONAL, INTENT(IN) :: fgdt INTEGER, OPTIONAL, INTENT(OUT) :: stepfg REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(OUT) ::& & rundgdten, rvndgdten, rthndgdten, rphndgdten, rqvndgdten REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(OUT) :: & & rmundgdten !URBAN ! REAL, DIMENSION(1:num_roof_layers), INTENT(INOUT) :: DZR !urban ! REAL, DIMENSION(1:num_wall_layers), INTENT(INOUT) :: DZB !urban ! REAL, DIMENSION(1:num_road_layers), INTENT(INOUT) :: DZG !urban !urban REAL, DIMENSION(num_soil_layers), OPTIONAL, INTENT(INOUT) :: dzr !urban REAL, DIMENSION(num_soil_layers), OPTIONAL, INTENT(INOUT) :: dzb !urban REAL, DIMENSION(num_soil_layers), OPTIONAL, INTENT(INOUT) :: dzg !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & tr_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & tb_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & tg_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & tc_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & qc_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & xxxr_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & xxxb_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & xxxg_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & xxxc_urb2d ! REAL, DIMENSION(ims:ime, 1:num_roof_layers, jms:jme), INTENT(INOUT) :: TRL_URB3D !urban ! REAL, DIMENSION(ims:ime, 1:num_wall_layers, jms:jme), INTENT(INOUT) :: TBL_URB3D !urban ! REAL, DIMENSION(ims:ime, 1:num_road_layers, jms:jme), INTENT(INOUT) :: TGL_URB3D !urban !urban REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: trl_urb3d !urban REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tbl_urb3d !urban REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tgl_urb3d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & sh_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & lh_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & g_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & rn_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & ts_urb2d !urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & frc_urb2d !urban INTEGER, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & utype_urb2d INTEGER, INTENT(IN) :: num_urban_layers INTEGER, INTENT(IN) :: num_urban_hi ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: trb_urb4d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tw1_urb4d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tw2_urb4d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tgb_urb4d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tlev_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: qlev_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tw1lev_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tw2lev_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tglev_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tflev_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & lf_ac_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & sf_ac_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & cm_ac_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & sfvent_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & lfvent_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfwin1_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfwin2_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfg_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfr_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfw1_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfw2_urb3d REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_u_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_v_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_t_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_q_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_e_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_u_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_v_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_t_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_q_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_e_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: vl_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: dlg_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: sf_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: dl_u_bep !obs fdda INTEGER, OPTIONAL, INTENT(IN) :: itimestep TYPE(FDOB_TYPE), OPTIONAL, INTENT(INOUT) :: fdob ! for obs-nudging base-state calcn REAL, OPTIONAL, INTENT(IN) :: p00, t00, tlp ! WA 12/21/09 REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(OUT) ::& & te_temf, cf3d_temf ! WA 2/22/11 REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(OUT) :: wm_temf ! Local data REAL :: alblnd, zzlnd, zzwtr, thinld, xmava, cen_lat, pptop REAL, DIMENSION(kms:kme) :: sfull, shalf REAL :: obs_twindo_cg, obs_twindo CHARACTER(len=256) :: mminlu_loc CHARACTER(len=512) :: message INTEGER :: iswater INTEGER :: isice INTEGER :: isurban INTEGER :: sf_urban_physics INTEGER :: sf_ocean_physics REAL :: oml_hml0 LOGICAL :: usemonalb LOGICAL :: rdmaxalb INTEGER :: mfshconv INTEGER :: i, j, k, itf, jtf, ktf, n INTEGER :: myproc INTRINSIC MIN0 INTRINSIC TRIM INTRINSIC PRESENT !----------------------------------------------------------------- sf_urban_physics = config_flags%sf_urban_physics usemonalb = config_flags%usemonalb rdmaxalb = config_flags%rdmaxalb mfshconv = config_flags%mfshconv obs_twindo_cg = model_config_rec%obs_twindo(1) obs_twindo = config_flags%obs_twindo oml_hml0 = config_flags%oml_hml0 sf_ocean_physics=config_flags%sf_ocean_physics !-- should be from the namelist sfull = 0. shalf = 0. CALL WRF_DEBUG(100, 'top of phy_init') WRITE(wrf_err_message, *) 'phy_init: start_of_simulation = ', & & start_of_simulation CALL WRF_DEBUG(100, TRIM(wrf_err_message)) IF (ite .GT. ide - 1) THEN itf = ide - 1 ELSE itf = ite END IF IF (jte .GT. jde - 1) THEN jtf = jde - 1 ELSE jtf = jte END IF IF (kte .GT. kde - 1) THEN ktf = kde - 1 ELSE ktf = kte END IF zzlnd = 0.1 zzwtr = 0.0001 thinld = 0.04 alblnd = 0.2 xmava = 0.3 CALL NL_GET_CEN_LAT(id, cen_lat) CALL WRF_DEBUG(100, & & 'calling nl_get_iswater, nl_get_isice, nl_get_mminlu_loc') CALL NL_GET_ISWATER(id, iswater) CALL NL_GET_ISICE(id, isice) CALL NL_GET_ISURBAN(id, isurban) CALL NL_GET_MMINLU(1, mminlu_loc) CALL WRF_DEBUG(100, & & 'after nl_get_iswater, nl_get_isice, nl_get_mminlu_loc') !-- temporary fix by ww landuse_isice = isice IF (.NOT.restart) THEN !-- initialize common variables IF (.NOT.moved) THEN DO j=jts,jtf DO i=its,itf xland(i, j) = 1. gsw(i, j) = 0. glw(i, j) = 0. !-- initialize ust to a small value ust(i, j) = 0.0001 mol(i, j) = 0.0 pblh(i, j) = 0.0 hfx(i, j) = 0. qfx(i, j) = 0. rainbl(i, j) = 0. rainncvd(i, j) = 0.0 rainncv(i, j) = 0. snowncvd(i, j) = 0.0 snowncv(i, j) = 0. graupelncvd(i, j) = 0.0 graupelncv(i, j) = 0. acsnow(i, j) = 0. !wig, 17-May-2006: Added for idealized chem. runs DO k=kms,kme exch_h(i, k, j) = 0. END DO END DO END DO END IF ! IF (PRESENT(tmn_update)) THEN IF (tmn_update .EQ. 1) THEN nyear = 1 nday = 0. DO j=jts,jtf DO i=its,itf tyr(i, j) = tmn(i, j) tyra(i, j) = 0.0 tdly(i, j) = 0.0 DO n=1,lagday tlag(i, n, j) = tmn(i, j) END DO END DO END DO END IF END IF ! ! DO j=jts,jtf DO i=its,itf IF (xland(i, j) .LT. 1.5) THEN IF (mminlu_loc .EQ. ' ') albbck(i, j) = alblnd embck(i, j) = 0.85 albedo(i, j) = albbck(i, j) emiss(i, j) = embck(i, j) thc(i, j) = thinld znt(i, j) = zzlnd z0(i, j) = zzlnd mavail(i, j) = xmava ELSE IF (mminlu_loc .EQ. ' ') albbck(i, j) = 0.08 albedo(i, j) = albbck(i, j) embck(i, j) = 0.98 emiss(i, j) = embck(i, j) thc(i, j) = thinld znt(i, j) = zzwtr z0(i, j) = zzwtr mavail(i, j) = 1.0 END IF END DO END DO IF (config_flags%cu_diag .EQ. 1) THEN DO j=jts,jtf DO k=kts,ktf DO i=its,itf gd_cloud(i, k, j) = 0. gd_cloud2(i, k, j) = 0. gd_cloud_a(i, k, j) = 0. gd_cloud2_a(i, k, j) = 0. END DO END DO END DO END IF DO j=jts,jtf DO i=its,itf raincv_a(i, j) = 0. raincv_b(i, j) = 0. END DO END DO CALL WRF_DEBUG(200, & & 'module_start: phy_init: Before call to landuse_init') IF (mminlu_loc .NE. ' ') CALL LANDUSE_INIT(lu_index, snowc, & & albedo, albbck, snoalb, & & mavail, emiss, embck, znt& & , z0, thc, xland, xice, & & xicem, julday, cen_lat, & & iswater, TRIM(mminlu_loc)& & , landuse_isice, & & landuse_lucats, & & landuse_luseas, & & landuse_isn, config_flags%& & fractional_seaice, & & lu_state, allowed_to_read& & , usemonalb, ids, ide, jds& & , jde, kds, kde, ims, ime& & , jms, jme, kms, kme, its& & , ite, jts, jte, kts, kte) !-- initialize surface properties END IF !-- initialize physics !-- ra: radiation !-- bl: pbl !-- cu: cumulus !-- mp: microphysics CALL WRF_DEBUG(200, 'module_start: phy_init: Before call to bl_init'& & ) ! WA !Optional urban !Optional urban !Optional urban !Optional urban !Optional urban !Optional urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional oml CALL BL_INIT_D(stepbl, bldt, dt, rublten, rubltend, rvblten, & & rvbltend, rthblten, rthbltend, rqvblten, rqvbltend, & & rqcblten, rqiblten, tsk, tmn, config_flags, restart, & & ust, lowlyr, tslb, zs, dzs, num_soil_layers, tke_pbl, & & mfshconv, massflux_edkf, entr_edkf, detr_edkf, thl_up, & & thv_up, rt_up, rv_up, rc_up, u_up, v_up, frac_up, exch_h, & & vegfra, snow, snowc, canwat, smstav, smstot, sfcrunoff, & & udrunoff, acsnow, acsnom, ivgtyp, isltyp, isurban, smois, & & smfr3d, mavail, snowh, sh2o, snoalb, fndsoilw, fndsnowh, & & rdmaxalb, znt, xland, xice, sfcevp, grdflx, TRIM(mminlu_loc& & ), isnowxy, zsnsoxy, tsnoxy, snicexy, snliqxy, tvxy, tgxy, & & canicexy, canliqxy, eahxy, tahxy, cmxy, chxy, fwetxy, & & sneqvoxy, alboldxy, qsnowxy, wslakexy, zwtxy, waxy, wtxy, & & lfmassxy, rtmassxy, stmassxy, woodxy, stblcpxy, fastcpxy, & & xsaixy, t2mvxy, t2mbxy, chstarxy, allowed_to_read, & & start_of_simulation, te_temf, cf3d_temf, wm_temf, dzr, dzb& & , dzg, tr_urb2d, tb_urb2d, tg_urb2d, tc_urb2d, qc_urb2d, & & xxxr_urb2d, xxxb_urb2d, xxxg_urb2d, xxxc_urb2d, trl_urb3d, & & tbl_urb3d, tgl_urb3d, sh_urb2d, lh_urb2d, g_urb2d, rn_urb2d& & , ts_urb2d, frc_urb2d, utype_urb2d, sf_urban_physics, & & num_urban_layers, trb_urb4d, tw1_urb4d, tw2_urb4d, & & tgb_urb4d, tlev_urb3d, qlev_urb3d, tw1lev_urb3d, & & tw2lev_urb3d, tglev_urb3d, tflev_urb3d, sf_ac_urb3d, & & lf_ac_urb3d, cm_ac_urb3d, sfvent_urb3d, lfvent_urb3d, & & sfwin1_urb3d, sfwin2_urb3d, sfw1_urb3d, sfw2_urb3d, & & sfr_urb3d, sfg_urb3d, a_u_bep, a_v_bep, a_t_bep, a_q_bep, & & a_e_bep, b_u_bep, b_v_bep, b_t_bep, b_q_bep, b_e_bep, & & dlg_bep, dl_u_bep, sf_bep, vl_bep, ids, ide, jds, jde, kds& & , kde, ims, ime, jms, jme, kms, kme, its, ite, jts, jte, & & kts, kte, achfx, aclhf, acgrdflx, oml_hml0, sf_ocean_physics, tml, & & t0ml, hml, h0ml, huml, hvml, tmoml) !Optional oml CALL WRF_DEBUG(200, 'module_start: phy_init: Before call to cu_init'& & ) CALL CU_INIT_D(dx, stepcu, cudt, dt, rucuten, rvcuten, rthcuten, & & rthcutend, rqvcuten, rqvcutend, rqrcuten, rqccuten, & & rqccutend, rqscuten, rqicuten, nca, rainc, raincd, raincv, & & raincvd, w0avg, w0avgd, config_flags, restart, cldefi, & & lowlyr, mass_flux, rthften, rqvften, apr_gr, apr_w, apr_mc& & , apr_st, apr_as, apr_capma, apr_capme, apr_capmi, & & cugd_tten, cugd_ttens, cugd_qvten, cugd_qvtens, cugd_qcten& & , allowed_to_read, start_of_simulation, ids, ide, jds, jde& & , kds, kde, ims, ime, jms, jme, kms, kme, its, ite, jts, & & jte, kts, kte) CALL WRF_DEBUG(200, 'module_start: phy_init: Before call to mp_init'& & ) CALL MP_INIT_D(rainnc, rainncd, snownc, snowncd, graupelnc, & & graupelncd, config_flags, restart, warm_rain, & & adv_moist_cond, mpdt, dt, dx, dy, lowlyr, f_ice_phy, & & f_rain_phy, f_rimef_phy, mp_restart_state, tbpvs_state, & & tbpvs0_state, allowed_to_read, start_of_simulation, ids, & & ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, its& & , ite, jts, jte, kts, kte) END SUBROUTINE G_PHY_INIT !================================================================= ! Differentiation of bl_init in forward (tangent) mode: ! variations of useful results: rublten rqvblten rvblten ! rthblten ! with respect to varying inputs: rublten rqvblten rvblten ! rthblten !===================================================================== ! WA ! num_roof_layers,num_wall_layers,num_road_layers,& !Optional urban !Optional urban !Optional urban !Optional urban !Optional urban !Optional urban !Optional urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional multi-layer urban !Optional oml SUBROUTINE BL_INIT_D(stepbl, bldt, dt, rublten, rubltend, rvblten, & & rvbltend, rthblten, rthbltend, rqvblten, rqvbltend, rqcblten, & & rqiblten, tsk, tmn, config_flags, restart, ust, lowlyr, tslb, & & zs, dzs, num_soil_layers, tke_pbl, mfshconv, massflux_edkf, & & entr_edkf, detr_edkf, thl_up, thv_up, rt_up, rv_up, rc_up, u_up, & & v_up, frac_up, exch_h, vegfra, snow, snowc, canwat, smstav, smstot, & & sfcrunoff, udrunoff, acsnow, acsnom, ivgtyp, isltyp, isurban, smois& & , smfr3d, mavail, snowh, sh2o, snoalb, fndsoilw, fndsnowh, rdmaxalb& & , znt, xland, xice, sfcevp, grdflx, mminlu, isnowxy, zsnsoxy, tsnoxy& & , snicexy, snliqxy, tvxy, tgxy, canicexy, canliqxy, eahxy, tahxy, & & cmxy, chxy, fwetxy, sneqvoxy, alboldxy, qsnowxy, wslakexy, zwtxy, & & waxy, wtxy, lfmassxy, rtmassxy, stmassxy, woodxy, stblcpxy, fastcpxy& & , xsaixy, t2mvxy, t2mbxy, chstarxy, allowed_to_read, & & start_of_simulation, te_temf, cf3d_temf, wm_temf, dzr, dzb, dzg, & & tr_urb2d, tb_urb2d, tg_urb2d, tc_urb2d, qc_urb2d, xxxr_urb2d, & & xxxb_urb2d, xxxg_urb2d, xxxc_urb2d, trl_urb3d, tbl_urb3d, tgl_urb3d& & , sh_urb2d, lh_urb2d, g_urb2d, rn_urb2d, ts_urb2d, frc_urb2d, & & utype_urb2d, sf_urban_physics, num_urban_layers, trb_urb4d, & & tw1_urb4d, tw2_urb4d, tgb_urb4d, tlev_urb3d, qlev_urb3d, & & tw1lev_urb3d, tw2lev_urb3d, tglev_urb3d, tflev_urb3d, sf_ac_urb3d, & & lf_ac_urb3d, cm_ac_urb3d, sfvent_urb3d, lfvent_urb3d, sfwin1_urb3d, & & sfwin2_urb3d, sfw1_urb3d, sfw2_urb3d, sfr_urb3d, sfg_urb3d, a_u_bep& & , a_v_bep, a_t_bep, a_q_bep, a_e_bep, b_u_bep, b_v_bep, b_t_bep, & & b_q_bep, b_e_bep, dlg_bep, dl_u_bep, sf_bep, vl_bep, ids, ide, jds, & & jde, kds, kde, ims, ime, jms, jme, kms, kme, its, ite, jts, jte, kts& & , kte, achfx, aclhf, acgrdflx, oml_hml0, sf_ocean_physics, tml, t0ml, hml, & & h0ml, huml, hvml, tmoml) USE module_sf_sfclay USE g_module_bl_surface_drag, ONLY : surface_drag_init_d IMPLICIT NONE !Optional oml !-------------------------------------------------------------------- !-------------------------------------------------------------------- TYPE(GRID_CONFIG_REC_TYPE) :: config_flags LOGICAL, INTENT(IN) :: restart LOGICAL, INTENT(IN) :: fndsoilw, fndsnowh LOGICAL, INTENT(IN) :: rdmaxalb INTEGER, INTENT(IN) :: ids, ide, jds, jde, kds, kde, ims, ime, jms, & & jme, kms, kme, its, ite, jts, jte, kts, kte INTEGER, INTENT(IN) :: num_soil_layers INTEGER, INTENT(IN) :: sf_urban_physics REAL, INTENT(IN) :: dt, bldt INTEGER, INTENT(INOUT) :: stepbl REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), INTENT(OUT) :: & & smfr3d REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), INTENT(INOUT) ::& & smois, sh2o, tslb REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: snow, snowh, & & snowc, snoalb, canwat, mavail, smstav, smstot, sfcrunoff, udrunoff, & & acsnow, vegfra, acsnom, sfcevp, grdflx, ust, znt, xland, xice INTEGER, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: ivgtyp, & & isltyp, lowlyr REAL, DIMENSION(num_soil_layers), INTENT(INOUT) :: zs, dzs REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rublten, & & rvblten, exch_h, rthblten, rqvblten, rqcblten, rqiblten, tke_pbl REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rubltend& & , rvbltend, rthbltend, rqvbltend REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(OUT), OPTIONAL ::& & massflux_edkf, entr_edkf, detr_edkf, thl_up, thv_up, rt_up, rv_up, & & rc_up, u_up, v_up, frac_up ! WRF JP INTEGER, INTENT(IN) :: mfshconv REAL, DIMENSION(ims:ime, jms:jme), INTENT(IN) :: tsk REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: tmn CHARACTER(len=*), INTENT(IN) :: mminlu LOGICAL, INTENT(IN) :: allowed_to_read INTEGER, INTENT(IN) :: isurban INTEGER :: isn, isfc INTEGER :: k !WA REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(OUT) ::& & te_temf, cf3d_temf REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(OUT) :: wm_temf !Noah-MP INTEGER, DIMENSION(ims:ime, jms:jme), OPTIONAL :: isnowxy REAL, DIMENSION(ims:ime, -2:num_soil_layers, jms:jme), OPTIONAL :: & & zsnsoxy REAL, DIMENSION(ims:ime, -2:0, jms:jme), OPTIONAL :: tsnoxy REAL, DIMENSION(ims:ime, -2:0, jms:jme), OPTIONAL :: snicexy REAL, DIMENSION(ims:ime, -2:0, jms:jme), OPTIONAL :: snliqxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: tvxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: tgxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: canicexy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: canliqxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: eahxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: tahxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: cmxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: chxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: fwetxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: sneqvoxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: alboldxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: qsnowxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: wslakexy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: zwtxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: waxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: wtxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: lfmassxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: rtmassxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: stmassxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: woodxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: stblcpxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: fastcpxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: xsaixy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: t2mvxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: t2mbxy REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL :: chstarxy !URBAN ! REAL, DIMENSION(1:num_roof_layers), INTENT(INOUT) :: DZR !Optional urban ! REAL, DIMENSION(1:num_wall_layers), INTENT(INOUT) :: DZB !Optional urban ! REAL, DIMENSION(1:num_road_layers), INTENT(INOUT) :: DZG !Optional urban !Optional urban REAL, DIMENSION(num_soil_layers), OPTIONAL, INTENT(INOUT) :: dzr !Optional urban REAL, DIMENSION(num_soil_layers), OPTIONAL, INTENT(INOUT) :: dzb !Optional urban REAL, DIMENSION(num_soil_layers), OPTIONAL, INTENT(INOUT) :: dzg !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & tr_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & tb_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & tg_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & tc_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & qc_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & xxxr_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & xxxb_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & xxxg_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & xxxc_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & sh_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & lh_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & g_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & rn_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & ts_urb2d !Optional urban REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & frc_urb2d !Optional urban INTEGER, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & utype_urb2d ! REAL, DIMENSION( ims:ime, 1:num_roof_layers, jms:jme ), INTENT(INOUT) :: TRL_URB3D !Optional urban ! REAL, DIMENSION( ims:ime, 1:num_wall_layers, jms:jme ), INTENT(INOUT) :: TBL_URB3D !Optional urban ! REAL, DIMENSION( ims:ime, 1:num_road_layers, jms:jme ), INTENT(INOUT) :: TGL_URB3D !Optional urban !Optional urban REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: trl_urb3d !Optional urban REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tbl_urb3d !Optional urban REAL, DIMENSION(ims:ime, num_soil_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tgl_urb3d INTEGER, INTENT(IN) :: num_urban_layers !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: trb_urb4d !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tw1_urb4d !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tw2_urb4d !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tgb_urb4d !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tlev_urb3d !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: qlev_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tw1lev_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tw2lev_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tglev_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: tflev_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & lf_ac_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & sf_ac_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & cm_ac_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & sfvent_urb3d !multi-layer UCM REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: & & lfvent_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfwin1_urb3d ! multi-layer UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfwin2_urb3d !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfw1_urb3d !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfw2_urb3d !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfr_urb3d !Optional UCM REAL, DIMENSION(ims:ime, num_urban_layers, jms:jme), OPTIONAL, & & INTENT(INOUT) :: sfg_urb3d REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_u_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_v_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_t_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_q_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: a_e_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_u_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_v_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_t_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_q_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: b_e_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: vl_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: dlg_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: sf_bep REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: dl_u_bep REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: achfx& & , aclhf, acgrdflx ! Optional OML variables REAL, DIMENSION(ims:ime, jms:jme), OPTIONAL, INTENT(INOUT) :: tml, & & t0ml, hml, h0ml, huml, hvml, tmoml INTEGER, OPTIONAL, INTENT(IN) :: sf_ocean_physics REAL, OPTIONAL, INTENT(IN) :: oml_hml0 LOGICAL, INTENT(IN) :: start_of_simulation INTEGER :: i, j !local mynn INTEGER :: mynn_closure_level INTRINSIC MAX INTRINSIC PRESENT INTRINSIC NINT !-- calculate pbl time step stepbl = NINT(bldt*60./dt) IF (stepbl .LT. 1) THEN stepbl = 1 ELSE stepbl = stepbl END IF !-- initialization IF (PRESENT(achfx)) THEN IF (.NOT.restart) THEN DO j=jts,jte DO i=its,ite achfx(i, j) = 0. aclhf(i, j) = 0. acgrdflx(i, j) = 0. sfcevp(i, j) = 0. END DO END DO END IF END IF !-- initialize surface layer scheme SELECT CASE (config_flags%sf_sfclay_physics) CASE (sfclayscheme) CALL SFCLAYINIT(allowed_to_read) isfc = 1 CASE DEFAULT !mynn !-- initialize pbl scheme SELECT CASE (config_flags%bl_pbl_physics) CASE (surfdragscheme) CALL SURFACE_DRAG_INIT_D(rublten, rubltend, rvblten, rvbltend, & & rthblten, rthbltend, rqvblten, rqvbltend, & & restart, ids, ide, jds, jde, kds, kde, ims, & & ime, jms, jme, kms, kme, its, ite, jts, jte, & & kts, kte) END SELECT END SELECT END SUBROUTINE BL_INIT_D !===================================================================== ! Differentiation of cu_init in forward (tangent) mode: ! variations of useful results: raincv rqccuten rthcuten rainc ! w0avg rqvcuten ! with respect to varying inputs: raincv rqccuten rthcuten rainc ! w0avg rqvcuten !================================================================== SUBROUTINE CU_INIT_D(dx, stepcu, cudt, dt, rucuten, rvcuten, rthcuten, & & rthcutend, rqvcuten, rqvcutend, rqrcuten, rqccuten, rqccutend, & & rqscuten, rqicuten, nca, rainc, raincd, raincv, raincvd, w0avg, & & w0avgd, config_flags, restart, cldefi, lowlyr, mass_flux, rthften, & & rqvften, apr_gr, apr_w, apr_mc, apr_st, apr_as, apr_capma, apr_capme& & , apr_capmi, cugd_tten, cugd_ttens, cugd_qvten, cugd_qvtens, & & cugd_qcten, allowed_to_read, start_of_simulation, ids, ide, jds, jde& & , kds, kde, ims, ime, jms, jme, kms, kme, its, ite, jts, jte, kts, & & kte) USE g_module_cu_du, ONLY : ducuinit_d IMPLICIT NONE !------------------------------------------------------------------ !------------------------------------------------------------------ TYPE(GRID_CONFIG_REC_TYPE) :: config_flags LOGICAL, INTENT(IN) :: restart INTEGER, INTENT(IN) :: ids, ide, jds, jde, kds, kde, ims, ime, jms, & & jme, kms, kme, its, ite, jts, jte, kts, kte REAL, INTENT(IN) :: dt, cudt, dx LOGICAL, INTENT(IN) :: start_of_simulation LOGICAL, INTENT(IN) :: allowed_to_read INTEGER, INTENT(INOUT) :: stepcu REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rucuten& & , rvcuten, rthcuten, rqvcuten, rqccuten, rqrcuten, rqicuten, & & rqscuten REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: & & rthcutend, rqvcutend, rqccutend REAL, DIMENSION(ims:ime, kms:kme, jms:jme), OPTIONAL, INTENT(INOUT) & & :: cugd_tten, cugd_ttens, cugd_qvten, cugd_qvtens, cugd_qcten REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(OUT) :: w0avg REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(OUT) :: w0avgd REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: rthften, & & rqvften REAL, DIMENSION(ims:ime, jms:jme), INTENT(OUT) :: rainc, raincv REAL, DIMENSION(ims:ime, jms:jme), INTENT(OUT) :: raincd, raincvd REAL, DIMENSION(ims:ime, jms:jme), INTENT(OUT) :: cldefi REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: nca REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: mass_flux, & & apr_gr, apr_w, apr_mc, apr_st, apr_as, apr_capma, apr_capme, & & apr_capmi INTEGER, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: lowlyr ! LOCAL VAR INTEGER :: i, j, itf, jtf INTRINSIC MIN0 INTRINSIC MAX INTRINSIC NINT IF (ite .GT. ide - 1) THEN itf = ide - 1 ELSE itf = ite END IF IF (jte .GT. jde - 1) THEN jtf = jde - 1 ELSE jtf = jte END IF ! stepcu = NINT(cudt*60./dt) IF (stepcu .LT. 1) THEN stepcu = 1 ELSE stepcu = stepcu END IF !-- initialization IF (start_of_simulation) THEN DO j=jts,jtf DO i=its,itf raincd(i, j) = 0.0 rainc(i, j) = 0. raincvd(i, j) = 0.0 raincv(i, j) = 0. END DO END DO END IF !-- deep convection and hybrid deep-shallow convection schemes SELECT CASE (config_flags%cu_physics) CASE (ducuscheme) CALL DUCUINIT_D(rthcuten, rthcutend, rqvcuten, rqvcutend, rqccuten& & , rqccutend, rqrcuten, rqicuten, rqscuten, nca, w0avg, & & w0avgd, p_qc, p_qr, svp1, svp2, svp3, svpt0, & & param_first_scalar, restart, allowed_to_read, ids, ide, & & jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, its, & & ite, jts, jte, kts, kte) END SELECT END SUBROUTINE CU_INIT_D !================================================================== ! Differentiation of mp_init in forward (tangent) mode: ! variations of useful results: rainnc snownc graupelnc ! with respect to varying inputs: rainnc snownc graupelnc !================================================================== ! for eta mp ! for eta mp ! eta mp SUBROUTINE MP_INIT_D(rainnc, rainncd, snownc, snowncd, graupelnc, & & graupelncd, config_flags, restart, warm_rain, adv_moist_cond, mpdt, & & dt, dx, dy, lowlyr, f_ice_phy, f_rain_phy, f_rimef_phy, & & mp_restart_state, tbpvs_state, tbpvs0_state, allowed_to_read, & & start_of_simulation, ids, ide, jds, jde, kds, kde, ims, ime, jms, & & jme, kms, kme, its, ite, jts, jte, kts, kte) IMPLICIT NONE !------------------------------------------------------------------ ! Arguments TYPE(GRID_CONFIG_REC_TYPE) :: config_flags LOGICAL, INTENT(IN) :: restart LOGICAL, INTENT(OUT) :: warm_rain, adv_moist_cond REAL, INTENT(IN) :: mpdt, dt, dx, dy LOGICAL, INTENT(IN) :: start_of_simulation INTEGER, INTENT(IN) :: ids, ide, jds, jde, kds, kde, ims, ime, jms, & & jme, kms, kme, its, ite, jts, jte, kts, kte INTEGER, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: lowlyr REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: rainnc, snownc, & & graupelnc REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT) :: rainncd, snowncd& & , graupelncd REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: & & f_ice_phy, f_rain_phy, f_rimef_phy REAL, DIMENSION(:), INTENT(INOUT) :: mp_restart_state, tbpvs_state, & & tbpvs0_state LOGICAL, INTENT(IN) :: allowed_to_read ! Local INTEGER :: i, j, itf, jtf INTRINSIC MIN0 warm_rain = .false. adv_moist_cond = .true. IF (ite .GT. ide - 1) THEN itf = ide - 1 ELSE itf = ite END IF IF (jte .GT. jde - 1) THEN jtf = jde - 1 ELSE jtf = jte END IF IF (start_of_simulation) THEN DO j=jts,jtf DO i=its,itf rainncd(i, j) = 0.0 rainnc(i, j) = 0. snowncd(i, j) = 0.0 snownc(i, j) = 0. graupelncd(i, j) = 0.0 graupelnc(i, j) = 0. END DO END DO END IF SELECT CASE (config_flags%mp_physics) END SELECT END SUBROUTINE MP_INIT_D !================================================================== SUBROUTINE Z2SIGMA(zf, zh, sf, sh, p_top, pptop, config_flags, & & allowed_to_read, kds, kde, kms, kme, kts, kte) IMPLICIT NONE ! Arguments INTEGER, INTENT(IN) :: kds, kde, kms, kme, kts, kte REAL, DIMENSION(kms:kme), INTENT(IN) :: zf, zh REAL, DIMENSION(kms:kme), INTENT(OUT) :: sf, sh REAL, INTENT(IN) :: p_top REAL, INTENT(OUT) :: pptop TYPE(GRID_CONFIG_REC_TYPE) :: config_flags LOGICAL, INTENT(IN) :: allowed_to_read ! Local REAL :: r, g, ts, gamma, ps, ztrop, tstrat, ptrop, z, t, p, ztop, & & ptop INTEGER :: k REAL :: pwx1 REAL :: pwy1 REAL :: pwr1 REAL :: arg1 INTRINSIC EXP IF (zf(kde/2) .GT. 1.0) THEN ! Height levels assumed (zeta coordinate) ! Convert to sigma using standard atmosphere for pressure-height relation ! constants for standard atmosphere definition r = 287.05 g = 9.80665 ts = 288.15 gamma = -(6.5/1000.) ps = 1013.25 ztrop = 11000. tstrat = ts + gamma*ztrop pwx1 = tstrat/ts pwy1 = -(g/(gamma*r)) pwr1 = pwx1**pwy1 ptrop = ps*pwr1 DO k=kde,kds,-1 ! full levels z = zf(k) IF (z .LE. ztrop) THEN t = ts + gamma*z pwx1 = t/ts pwy1 = -(g/(gamma*r)) pwr1 = pwx1**pwy1 p = ps*pwr1 ELSE t = tstrat arg1 = -(g*(z-ztrop)/(r*tstrat)) p = ptrop*EXP(arg1) END IF IF (k .EQ. kde) THEN ztop = zf(k) ptop = p END IF sf(k) = (p-ptop)/(ps-ptop) ! half levels IF (k .NE. kds) THEN z = 0.5*(zf(k)+zf(k-1)) IF (z .LE. ztrop) THEN t = ts + gamma*z pwx1 = t/ts pwy1 = -(g/(gamma*r)) pwr1 = pwx1**pwy1 p = ps*pwr1 ELSE t = tstrat arg1 = -(g*(z-ztrop)/(r*tstrat)) p = ptrop*EXP(arg1) END IF sh(k-1) = (p-ptop)/(ps-ptop) END IF END DO pptop = ptop/10. ELSE ! Levels are already sigma/eta DO k=kde,kds,-1 ! sf(k)=zf(kde-k+kds) ! if(k .ne. kde)sh(k)=zh(kde-1-k+kds) sf(k) = zf(k) IF (k .NE. kde) sh(k) = zh(k) END DO pptop = p_top/1000. END IF END SUBROUTINE Z2SIGMA END MODULE G_MODULE_PHYSICS_INIT