PageRenderTime 74ms CodeModel.GetById 17ms RepoModel.GetById 0ms app.codeStats 1ms

/wrfv2_fire/phys/module_physics_init.F

http://github.com/jbeezley/wrf-fire
FORTRAN Legacy | 2890 lines | 2027 code | 331 blank | 532 comment | 29 complexity | d714e1d2e2743694f30eb848ea228815 MD5 | raw file
Possible License(s): AGPL-1.0 
    

  1. !WRF:MODEL_LAYER:INITIALIZATION
    2. 

  2. !
    3. 

  3. 

  4. !  This MODULE holds the routines which are used to perform model start-up operations
    5. 

  5. !  for the individual domains.  This is the stage after inputting wrfinput and before
    6. 

  6. !  calling 'integrate'.
    7. 

  7. 

  8. !  This MODULE CONTAINS the following routines:
    9. 

  9. 

  10. 

  11. MODULE module_physics_init
    12. 

  12. 

  13.    USE module_state_description
    14. 

  14.    USE module_model_constants
    15. 

  15.    USE module_configure, ONLY : grid_config_rec_type
    16. 

  16. !  USE module_ssib_veg , ONLY : init_module_ssib_veg !fds (SSiB constants)
    17. 

  17. 

  18. CONTAINS
    19. 

  19. 

  20. 

  21. !=================================================================
    22. 

  22.    SUBROUTINE phy_init ( id, config_flags, DT, restart, zfull, zhalf,     &
    23. 

  23.                          p_top, TSK,RADT,BLDT,CUDT,MPDT,         &
    24. 

  24.                          RUCUTEN, RVCUTEN,                       &
    25. 

  25.                          RTHCUTEN, RQVCUTEN, RQRCUTEN,           &
    26. 

  26.                          RQCCUTEN, RQSCUTEN, RQICUTEN,           &
    27. 

  27.                          RUSHTEN,  RVSHTEN,  RTHSHTEN,           &
    28. 

  28.                          RQVSHTEN, RQRSHTEN, RQCSHTEN,           &
    29. 

  29.                          RQSSHTEN, RQISHTEN, RQGSHTEN,           &
    30. 

  30.                          RUBLTEN,RVBLTEN,RTHBLTEN,               &
    31. 

  31.                          RQVBLTEN,RQCBLTEN,RQIBLTEN,             &
    32. 

  32.                          RTHRATEN,RTHRATENLW,RTHRATENSW,         &
    33. 

  33.                          STEPBL,STEPRA,STEPCU,                   &
    34. 

  34.                          W0AVG, RAINNC, RAINC, RAINCV, RAINNCV,  &
    35. 

  35.                          SNOWNC, SNOWNCV, GRAUPELNC, GRAUPELNCV, &
    36. 

  36.                          NCA,swrad_scat,                         &
    37. 

  37.                          CLDEFI,LOWLYR,                          &
    38. 

  38.                          MASS_FLUX,                              &
    39. 

  39.                          RTHFTEN, RQVFTEN,                       &
    40. 

  40.                          CLDFRA,CLDFRA_OLD,GLW,GSW,EMISS,EMBCK,  & !EMBCK new
    41. 

  41.                          LU_INDEX,                               &
    42. 

  42.                          landuse_ISICE, landuse_LUCATS,          &
    43. 

  43.                          landuse_LUSEAS, landuse_ISN,            &
    44. 

  44.                          lu_state,                               &
    45. 

  45.                          XLAT,XLONG,ALBEDO,ALBBCK,GMT,JULYR,JULDAY,&
    46. 

  46.                          levsiz, n_ozmixm, n_aerosolc, paerlev,  &
    47. 

  47.                          TMN,XLAND,ZNT,Z0,UST,MOL,PBLH,TKE_PBL,  &
    48. 

  48.                          EXCH_H,THC,SNOWC,MAVAIL,HFX,QFX,RAINBL, &
    49. 

  49.                          TSLB,ZS,DZS,num_soil_layers,warm_rain,  &
    50. 

  50.                          adv_moist_cond,                         &
    51. 

  51.                          APR_GR,APR_W,APR_MC,APR_ST,APR_AS,      &
    52. 

  52.                          APR_CAPMA,APR_CAPME,APR_CAPMI,          &
    53. 

  53.                          XICE,XICEM,VEGFRA,SNOW,CANWAT,SMSTAV,   &
    54. 

  54.                          SMSTOT, SFCRUNOFF,UDRUNOFF,GRDFLX,ACSNOW,&
    55. 

  55.                          ACSNOM,IVGTYP,ISLTYP, SFCEVP, SMOIS,    &
    56. 

  56.                          SH2O, SNOWH, SMFR3D,                    &  ! temporary
    57. 

  57.                          SNOALB,                                 &
    58. 

  58.                          DX,DY,F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY, &
    59. 

  59.                          mp_restart_state,tbpvs_state,tbpvs0_state,&
    60. 

  60.                          allowed_to_read, moved, start_of_simulation,&
    61. 

  61.                          LAGDAY,                                 &
    62. 

  62.                          ids, ide, jds, jde, kds, kde,           &
    63. 

  63.                          ims, ime, jms, jme, kms, kme,           &
    64. 

  64.                          its, ite, jts, jte, kts, kte,           &
    65. 

  65.                          NUM_URBAN_LAYERS,                       &
    66. 

  66.                          raincv_a,raincv_b,                      &
    67. 

  67.                          gd_cloud,gd_cloud2,                     &    ! Optional
    68. 

  68.                          gd_cloud_a,gd_cloud2_a,                 &    ! Optional
    69. 

  69.                          gd_cloud_b,gd_cloud2_b,                 &    ! Optional
    70. 

  70.                          ozmixm,pin,                             &    ! Optional
    71. 

  71.                          m_ps_1,m_ps_2,m_hybi,aerosolc_1,aerosolc_2,& ! Optional
    72. 

  72.                          RUNDGDTEN,RVNDGDTEN,RTHNDGDTEN,         &    ! Optional
    73. 

  73.                          RPHNDGDTEN,RQVNDGDTEN,RMUNDGDTEN,       &    ! Optional
    74. 

  74.                          FGDT,STEPFG,                            &    ! Optional
    75. 

  75.                          cugd_tten,cugd_ttens,cugd_qvten,        &    ! Optional
    76. 

  76.                          cugd_qvtens,cugd_qcten,                 &    ! Optional
    77. 

  77.                          ISNOWXY, ZSNSOXY, TSNOXY,                        & ! Optional Noah-MP
    78. 

  78.                          SNICEXY, SNLIQXY, TVXY, TGXY, CANICEXY,          & ! Optional Noah-MP
    79. 

  79.                          CANLIQXY, EAHXY, TAHXY, CMXY,                    & ! Optional Noah-MP
    80. 

  80.                          CHXY, FWETXY, SNEQVOXY, ALBOLDXY, QSNOWXY,       & ! Optional Noah-MP
    81. 

  81.                          WSLAKEXY, ZWTXY, WAXY, WTXY, LFMASSXY, RTMASSXY, & ! Optional Noah-MP
    82. 

  82.                          STMASSXY, WOODXY, STBLCPXY, FASTCPXY,            & ! Optional Noah-MP
    83. 

  83.                          XSAIXY,                                          & ! Optional Noah-MP
    84. 

  84.                          T2MVXY, T2MBXY, CHSTARXY ,                       & ! Optional Noah-MP
    85. 

  85. !                        num_roof_layers,num_wall_layers,        & !Optional urban
    86. 

  86. !                        num_road_layers,                        & !Optional urban
    87. 

  87.                          DZR, DZB, DZG,                          & !Optional urban
    88. 

  88.                          TR_URB2D,TB_URB2D,TG_URB2D,TC_URB2D,    & !Optional urban
    89. 

  89.                          QC_URB2D, XXXR_URB2D,XXXB_URB2D,        & !Optional urban
    90. 

  90.                          XXXG_URB2D, XXXC_URB2D,                 & !Optional urban
    91. 

  91.                          TRL_URB3D, TBL_URB3D, TGL_URB3D,        & !Optional urban
    92. 

  92.                          SH_URB2D, LH_URB2D, G_URB2D, RN_URB2D,  & !Optional urban
    93. 

  93.                          TS_URB2D, FRC_URB2D, UTYPE_URB2D,       & !Optional urban
    94. 

  94.                          TRB_URB4D,TW1_URB4D,TW2_URB4D,          & !Optional multi-layer urban
    95. 

  95.                          TGB_URB4D,TLEV_URB3D,QLEV_URB3D,        & !Optional multi-layer urban
    96. 

  96.                          TW1LEV_URB3D,TW2LEV_URB3D,              & !Optional multi-layer urban
    97. 

  97.                          TGLEV_URB3D,TFLEV_URB3D,                & !Optional multi-layer urban
    98. 

  98.                          SF_AC_URB3D,LF_AC_URB3D,CM_AC_URB3D,    & !Optional multi-layer urban
    99. 

  99.                          SFVENT_URB3D,LFVENT_URB3D,              & !Optional multi-layer urban
    100. 

  100.                          SFWIN1_URB3D,SFWIN2_URB3D,              & !Optional multi-layer urban
    101. 

  101.                          SFW1_URB3D,SFW2_URB3D,                  & !Optional multi-layer urban
    102. 

  102.                          SFR_URB3D,SFG_URB3D,                    & !Optional multi-layer urban
    103. 

  103.                          A_U_BEP,A_V_BEP,A_T_BEP,A_Q_BEP,        & !Optional multi-layer urban
    104. 

  104.                          A_E_BEP,B_U_BEP,B_V_BEP,                & !Optional multi-layer urban
    105. 

  105.                          B_T_BEP,B_Q_BEP,B_E_BEP,DLG_BEP,        & !Optional multi-layer urban
    106. 

  106.                          DL_U_BEP,SF_BEP,VL_BEP,                 & !Optional multi-layer urban
    107. 

  107.                          TML,T0ML,HML,H0ML,HUML,HVML,TMOML,      & !Optional oml
    108. 

  108.                          itimestep,                              & !Optional obs fdda
    109. 

  109. #if ( EM_CORE == 1 )
    110. 

  110.                          fdob,                                   & !Optional obs fdda
    111. 

  111. #endif
    112. 

  112.                          t00, p00, tlp,                          & !for obs-nudging
    113. 

  113.                          TYR,TYRA,TDLY,TLAG,NYEAR,NDAY,tmn_update,   &
    114. 

  114.                          ACHFX,ACLHF,ACGRDFLX                    &
    115. 

  115.    ! OPTIONAL 
    116. 

  116.                          ,te_temf          & ! WA 12/21/09
    117. 

  117.                          ,cf3d_temf        & ! WA 9/27/10
    118. 

  118.                          ,wm_temf          & ! WA 2/22/11
    119. 

  119.                          ,massflux_EDKF, entr_EDKF, detr_EDKF  & ! Optional for qnse
    120. 

  120.                          ,thl_up, thv_up, rt_up                & ! Optional for qnse
    121. 

  121.                          ,rv_up, rc_up, u_up, v_up, frac_up    & ! Optional for qnse
    122. 

  122.                          )
    123. 

  123. 

  124. !-----------------------------------------------------------------
    125. 

  125.    USE module_domain
    126. 

  126.    USE module_wrf_error
    127. 

  127.    USE module_wind_generic
    128. 

  128.    USE module_wind_fitch
    129. 

  129.    IMPLICIT NONE
    130. 

  130. !-----------------------------------------------------------------
    131. 

  131.    TYPE (grid_config_rec_type)              :: config_flags
    132. 

  132. 

  133.    INTEGER , INTENT(IN)        :: id
    134. 

  134.    INTEGER , INTENT(IN) ,OPTIONAL       :: tmn_update
    135. 

  135.    LOGICAL , INTENT(OUT)       :: warm_rain,adv_moist_cond
    136. 

  136. !   LOGICAL , INTENT (IN)       :: FNDSOILW, FNDSNOWH
    137. 

  137.    LOGICAL, PARAMETER          :: FNDSOILW=.true., FNDSNOWH=.true.
    138. 

  138.    INTEGER , INTENT(IN)        :: ids, ide, jds, jde, kds, kde,  &
    139. 

  139.                                   ims, ime, jms, jme, kms, kme,  &
    140. 

  140.                                   its, ite, jts, jte, kts, kte
    141. 

  141. 

  142.    INTEGER , INTENT(IN)        :: num_soil_layers
    143. 

  143.    INTEGER , INTENT(IN)        :: lagday
    144. 

  144.    INTEGER , INTENT(OUT) ,OPTIONAL      :: nyear
    145. 

  145.    REAL    , INTENT(OUT) ,OPTIONAL      :: nday
    146. 

  146. 

  147.    LOGICAL,  INTENT(IN)        :: start_of_simulation
    148. 

  148.    REAL,     INTENT(IN)        :: DT, p_top, DX, DY
    149. 

  149.    LOGICAL,  INTENT(IN)        :: restart
    150. 

  150.    REAL,     INTENT(IN)        :: RADT,BLDT,CUDT,MPDT
    151. 

  151.    REAL,     INTENT(IN)        :: swrad_scat
    152. 

  152. 

  153.    REAL,     DIMENSION( kms:kme ) , INTENT(IN) :: zfull, zhalf
    154. 

  154.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(IN) :: TSK, XLAT, XLONG
    155. 

  155. 

  156.    INTEGER,      INTENT(IN   )    ::   levsiz, n_ozmixm
    157. 

  157.    INTEGER,      INTENT(IN   )    ::   paerlev, n_aerosolc
    158. 

  158. 

  159.    REAL,  DIMENSION( ims:ime, levsiz, jms:jme, n_ozmixm ), OPTIONAL, &
    160. 

  160.           INTENT(INOUT) ::                                  OZMIXM
    161. 

  161. 

  162.    REAL,  DIMENSION(levsiz), OPTIONAL, INTENT(INOUT)  ::        PIN
    163. 

  163. 

  164.    REAL,  DIMENSION(ims:ime,jms:jme), OPTIONAL, INTENT(INOUT)  :: m_ps_1,m_ps_2
    165. 

  165.    REAL,  DIMENSION(paerlev), OPTIONAL,INTENT(INOUT)  ::          m_hybi
    166. 

  166.    REAL,  DIMENSION( ims:ime, paerlev, jms:jme, n_aerosolc ), OPTIONAL, &
    167. 

  167.           INTENT(INOUT) ::                    aerosolc_1, aerosolc_2
    168. 

  168. 

  169.    REAL,     DIMENSION( ims:ime , 1:num_soil_layers , jms:jme ),&
    170. 

  170.                  INTENT(INOUT) :: SMOIS, SH2O,TSLB
    171. 

  171.    REAL,     DIMENSION( ims:ime , 1:num_soil_layers , jms:jme ), INTENT(OUT) :: SMFR3D
    172. 

  172. 

  173.    REAL,    DIMENSION( ims:ime, jms:jme )                     , &
    174. 

  174.             INTENT(INOUT)    ::                           SNOW, &
    175. 

  175.                                                          SNOWC, &
    176. 

  176.                                                          SNOWH, &
    177. 

  177.                                                         CANWAT, &
    178. 

  178.                                                         SMSTAV, &
    179. 

  179.                                                         SMSTOT, &
    180. 

  180.                                                      SFCRUNOFF, &
    181. 

  181.                                                       UDRUNOFF, &
    182. 

  182.                                                         SFCEVP, &
    183. 

  183.                                                         GRDFLX, &
    184. 

  184.                                                         ACSNOW, &
    185. 

  185.                                                           XICE, &
    186. 

  186.                                                          XICEM, &
    187. 

  187.                                                         VEGFRA, &
    188. 

  188.                                                         ACSNOM
    189. 

  189. 

  190.    REAL,    DIMENSION( ims:ime, jms:jme )                     , &
    191. 

  191.             OPTIONAL, INTENT(INOUT)    ::                ACHFX, &
    192. 

  192.                                                          ACLHF, &
    193. 

  193.                                                       ACGRDFLX
    194. 

  194. 

  195.    INTEGER, DIMENSION( ims:ime, jms:jme )                     , &
    196. 

  196.             INTENT(INOUT)    ::                         IVGTYP, &
    197. 

  197.                                                         ISLTYP
    198. 

  198. 

  199. ! rad
    200. 

  200. 

  201.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) ::    &
    202. 

  202.              RTHRATEN, RTHRATENLW, RTHRATENSW, CLDFRA
    203. 

  203. 

  204.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , OPTIONAL, INTENT(OUT) :: &
    205. 

  205.              CLDFRA_OLD
    206. 

  206. 

  207.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT) ::         &
    208. 

  208.              GSW,ALBEDO,ALBBCK,GLW,EMISS,EMBCK                          !EMBCK new
    209. 

  209.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT) ::   SNOALB
    210. 

  210. 

  211. 

  212.    REAL,     INTENT(IN) :: GMT
    213. 

  213. 

  214.    INTEGER , INTENT(OUT) :: STEPRA, STEPBL, STEPCU
    215. 

  215.    INTEGER , INTENT(IN) :: JULYR, JULDAY
    216. 

  216. 

  217. ! cps
    218. 

  218. 

  219.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) :: &
    220. 

  220.              RUCUTEN, RVCUTEN, RTHCUTEN, RQVCUTEN, RQRCUTEN, RQCCUTEN, &
    221. 

  221.              RQSCUTEN, RQICUTEN,                                       &
    222. 

  222.              RUSHTEN, RVSHTEN, RTHSHTEN, RQVSHTEN, RQRSHTEN, RQCSHTEN, &
    223. 

  223.              RQSSHTEN, RQISHTEN, RQGSHTEN
    224. 

  224. 

  225.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) :: W0AVG
    226. 

  226. 

  227.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(OUT) :: MASS_FLUX,   &
    228. 

  228.                       APR_GR,APR_W,APR_MC,APR_ST,APR_AS,          &
    229. 

  229.                       APR_CAPMA,APR_CAPME,APR_CAPMI
    230. 

  230. 

  231.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) ::    &
    232. 

  232.              RTHFTEN, RQVFTEN
    233. 

  233. 

  234.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(OUT) ::           &
    235. 

  235.              RAINNC, RAINC, RAINCV, RAINNCV,  &
    236. 

  236.              SNOWNC, SNOWNCV, GRAUPELNC, GRAUPELNCV
    237. 

  237. 

  238.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(OUT) :: CLDEFI, NCA
    239. 

  239. 

  240.    INTEGER,  DIMENSION( ims:ime , jms:jme ) , INTENT(OUT) :: LOWLYR
    241. 

  241. 

  242. !pbl
    243. 

  243. 

  244.    ! soil layer
    245. 

  245. 

  246. 

  247.    REAL,     DIMENSION(1:num_soil_layers),      INTENT(INOUT) :: ZS,DZS
    248. 

  248. 

  249.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) ::    &
    250. 

  250.              RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN,RQCBLTEN,RQIBLTEN,EXCH_H,TKE_PBL
    251. 

  251. 

  252.   REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT), OPTIONAL :: &
    253. 

  253.                                             massflux_EDKF, entr_EDKF, detr_EDKF & 
    254. 

  254.                                                    ,thl_up, thv_up, rt_up       &
    255. 

  255.                                                    ,rv_up, rc_up, u_up, v_up    &
    256. 

  256.                                                    ,frac_up
    257. 

  257. 

  258. 

  259.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , OPTIONAL, INTENT(OUT) ::    &
    260. 

  260.              cugd_tten,cugd_ttens,cugd_qvten,                &
    261. 

  261.              cugd_qvtens,cugd_qcten
    262. 

  262.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT) ::         &
    263. 

  263.              XLAND,ZNT,Z0,UST,MOL,LU_INDEX,                         &
    264. 

  264.              PBLH,THC,MAVAIL,HFX,QFX,RAINBL
    265. 

  265.    INTEGER , INTENT(INOUT)  :: landuse_ISICE, landuse_LUCATS
    266. 

  266.    INTEGER , INTENT(INOUT)  :: landuse_LUSEAS, landuse_ISN
    267. 

  267.    REAL    , INTENT(INOUT)  , DIMENSION( : ) :: lu_state
    268. 

  268. 

  269.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT) :: TMN
    270. 

  270.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT),OPTIONAL :: TYR
    271. 

  271.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT),OPTIONAL :: TYRA
    272. 

  272.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT),OPTIONAL :: TDLY
    273. 

  273.    REAL,     DIMENSION( ims:ime , 1:lagday , jms:jme ) , INTENT(INOUT),OPTIONAL :: TLAG
    274. 

  274. 

  275. 

  276.    REAL,  DIMENSION( ims:ime , kms:kme , jms:jme )         ,         &
    277. 

  277.           OPTIONAL,                                                  &
    278. 

  278.           INTENT(INOUT   ) ::                                        &
    279. 

  279.                                gd_cloud, gd_cloud2,                  &
    280. 

  280.                                gd_cloud_a, gd_cloud2_a,              &
    281. 

  281.                                gd_cloud_b, gd_cloud2_b
    282. 

  282. 

  283.    REAL,  DIMENSION( ims:ime ,  jms:jme )         ,         &
    284. 

  284.           INTENT(INOUT   ) ::                                        &
    285. 

  285.                                raincv_a,raincv_b
    286. 

  286. 

  287. !Noah-MP
    288. 

  288. 

  289. 

  290.    INTEGER, OPTIONAL, DIMENSION(ims:ime,jms:jme) :: ISNOWXY
    291. 

  291.    REAL,    OPTIONAL, DIMENSION(ims:ime,-2:num_soil_layers, jms:jme) :: ZSNSOXY
    292. 

  292.    REAL,    OPTIONAL, DIMENSION(ims:ime,-2:0, jms:jme) :: TSNOXY
    293. 

  293.    REAL,    OPTIONAL, DIMENSION(ims:ime,-2:0, jms:jme) :: SNICEXY
    294. 

  294.    REAL,    OPTIONAL, DIMENSION(ims:ime,-2:0, jms:jme) :: SNLIQXY
    295. 

  295.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: TVXY
    296. 

  296.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: TGXY
    297. 

  297.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CANICEXY
    298. 

  298.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CANLIQXY
    299. 

  299.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: EAHXY
    300. 

  300.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: TAHXY
    301. 

  301.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CMXY
    302. 

  302.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CHXY
    303. 

  303.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: FWETXY
    304. 

  304.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: SNEQVOXY
    305. 

  305.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: ALBOLDXY
    306. 

  306.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: QSNOWXY
    307. 

  307.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: WSLAKEXY
    308. 

  308.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: ZWTXY
    309. 

  309.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: WAXY
    310. 

  310.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: WTXY
    311. 

  311.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: LFMASSXY
    312. 

  312.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: RTMASSXY
    313. 

  313.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: STMASSXY
    314. 

  314.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: WOODXY
    315. 

  315.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: STBLCPXY
    316. 

  316.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: FASTCPXY
    317. 

  317.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: XSAIXY
    318. 

  318.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: T2MVXY
    319. 

  319.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: T2MBXY
    320. 

  320.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CHSTARXY
    321. 

  321. 

  322. !mp
    323. 

  323.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) ::   &
    324. 

  324.              F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY
    325. 

  325.    REAL, DIMENSION(:), INTENT(INOUT)   :: mp_restart_state,tbpvs_state,tbpvs0_state
    326. 

  326.    LOGICAL,  INTENT(IN)  :: allowed_to_read, moved
    327. 

  327. 

  328. ! ocean mixed layer
    329. 

  329.    REAL,     DIMENSION( ims:ime , jms:jme ) , OPTIONAL, INTENT(INOUT) ::    &
    330. 

  330.              TML,T0ML,HML,H0ML,HUML,HVML,TMOML
    331. 

  331. 

  332. !fdda
    333. 

  333.    REAL,     OPTIONAL, INTENT(IN) :: FGDT
    334. 

  334.    INTEGER , OPTIONAL, INTENT(OUT) :: STEPFG
    335. 

  335.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , OPTIONAL, INTENT(OUT) ::    &
    336. 

  336.              RUNDGDTEN, RVNDGDTEN, RTHNDGDTEN, RPHNDGDTEN, RQVNDGDTEN
    337. 

  337.    REAL,     DIMENSION( ims:ime , jms:jme ) , OPTIONAL, INTENT(OUT) ::    &
    338. 

  338.              RMUNDGDTEN
    339. 

  339. 

  340. !URBAN
    341. 

  341. !   REAL, DIMENSION(1:num_roof_layers), INTENT(INOUT) :: DZR   !urban
    342. 

  342. !   REAL, DIMENSION(1:num_wall_layers), INTENT(INOUT) :: DZB   !urban
    343. 

  343. !   REAL, DIMENSION(1:num_road_layers), INTENT(INOUT) :: DZG   !urban
    344. 

  344.    REAL, OPTIONAL, DIMENSION(1:num_soil_layers), INTENT(INOUT) :: DZR    !urban
    345. 

  345.    REAL, OPTIONAL, DIMENSION(1:num_soil_layers), INTENT(INOUT) :: DZB    !urban
    346. 

  346.    REAL, OPTIONAL, DIMENSION(1:num_soil_layers), INTENT(INOUT) :: DZG    !urban
    347. 

  347. 

  348.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TR_URB2D !urban
    349. 

  349.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TB_URB2D !urban
    350. 

  350.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TG_URB2D !urban
    351. 

  351.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TC_URB2D !urban
    352. 

  352.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: QC_URB2D !urban
    353. 

  353.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: XXXR_URB2D !urban
    354. 

  354.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: XXXB_URB2D !urban
    355. 

  355.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: XXXG_URB2D !urban
    356. 

  356.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: XXXC_URB2D !urban
    357. 

  357. 

  358. !   REAL, DIMENSION(ims:ime, 1:num_roof_layers, jms:jme), INTENT(INOUT) :: TRL_URB3D !urban
    359. 

  359. !   REAL, DIMENSION(ims:ime, 1:num_wall_layers, jms:jme), INTENT(INOUT) :: TBL_URB3D !urban
    360. 

  360. !   REAL, DIMENSION(ims:ime, 1:num_road_layers, jms:jme), INTENT(INOUT) :: TGL_URB3D !urban
    361. 

  361.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_soil_layers, jms:jme), INTENT(INOUT) :: TRL_URB3D  !urban
    362. 

  362.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_soil_layers, jms:jme), INTENT(INOUT) :: TBL_URB3D  !urban
    363. 

  363.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_soil_layers, jms:jme), INTENT(INOUT) :: TGL_URB3D  !urban
    364. 

  364. 

  365.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: SH_URB2D !urban
    366. 

  366.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: LH_URB2D !urban
    367. 

  367.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: G_URB2D !urban
    368. 

  368.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: RN_URB2D !urban
    369. 

  369.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TS_URB2D !urban
    370. 

  370.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: FRC_URB2D !urban
    371. 

  371.    INTEGER, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: UTYPE_URB2D !urban
    372. 

  372. 

  373.    INTEGER , INTENT(IN)        :: num_urban_layers
    374. 

  374.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: TRB_URB4D ! multi-layer UCM
    375. 

  375.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: TW1_URB4D ! multi-layer UCM
    376. 

  376.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: TW2_URB4D ! multi-layer UCM
    377. 

  377.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: TGB_URB4D ! multi-layer UCM
    378. 

  378.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: TLEV_URB3D ! multi-layer UCM
    379. 

  379.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: QLEV_URB3D ! multi-layer UCM
    380. 

  380.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TW1LEV_URB3D ! multi-layer UCM
    381. 

  381.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TW2LEV_URB3D ! multi-layer UCM
    382. 

  382.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TGLEV_URB3D ! multi-layer UCM
    383. 

  383.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TFLEV_URB3D ! multi-layer UCM
    384. 

  384.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: LF_AC_URB3D !multi-layer UCM
    385. 

  385.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: SF_AC_URB3D !multi-layer UCM
    386. 

  386.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: CM_AC_URB3D !multi-layer UCM
    387. 

  387.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: SFVENT_URB3D !multi-layer UCM
    388. 

  388.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: LFVENT_URB3D !multi-layer UCM
    389. 

  389.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: SFWIN1_URB3D ! multi-layer UCM
    390. 

  390.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: SFWIN2_URB3D ! multi-layer UCM
    391. 

  391.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: SFG_URB3D ! multi-layer UCM
    392. 

  392.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: SFR_URB3D ! multi-layer UCM
    393. 

  393.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: SFW1_URB3D ! multi-layer UCM
    394. 

  394.    REAL, OPTIONAL, DIMENSION(ims:ime, 1:num_urban_layers, jms:jme), INTENT(INOUT) :: SFW2_URB3D ! multi-layer UCM
    395. 

  395.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_U_BEP
    396. 

  396.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_V_BEP
    397. 

  397.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_T_BEP
    398. 

  398.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_Q_BEP
    399. 

  399.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_E_BEP
    400. 

  400.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_U_BEP
    401. 

  401.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_V_BEP
    402. 

  402.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_T_BEP
    403. 

  403.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_Q_BEP
    404. 

  404.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_E_BEP
    405. 

  405.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: VL_BEP
    406. 

  406.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: DLG_BEP
    407. 

  407.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme,jms:jme), INTENT(INOUT) :: SF_BEP
    408. 

  408.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: DL_U_BEP
    409. 

  409. 

  410. !obs fdda
    411. 

  411.    INTEGER, OPTIONAL, INTENT(IN) :: itimestep
    412. 

  412. #if ( EM_CORE == 1 )
    413. 

  413.    TYPE(fdob_type), OPTIONAL, INTENT(INOUT) :: fdob
    414. 

  414. #endif
    415. 

  415.    REAL, OPTIONAL, INTENT(IN) :: p00, t00, tlp   ! for obs-nudging base-state calcn
    416. 

  416. 

  417. ! WA 12/21/09
    418. 

  418.    REAL,OPTIONAL, DIMENSION( ims:ime , kms:kme , jms:jme ) , &
    419. 

  419.           INTENT(OUT) ::    te_temf, cf3d_temf    
    420. 

  420. ! WA 2/22/11
    421. 

  421.    REAL,OPTIONAL, DIMENSION( ims:ime , jms:jme ) , &
    422. 

  422.           INTENT(OUT) ::    wm_temf
    423. 

  423. 

  424. ! Local data
    425. 

  425. 

  426.    REAL    :: ALBLND,ZZLND,ZZWTR,THINLD,XMAVA,CEN_LAT,pptop
    427. 

  427.    REAL,     DIMENSION( kms:kme )  :: sfull, shalf
    428. 

  428.    REAL :: obs_twindo_cg, obs_twindo
    429. 

  429. 

  430.    CHARACTER*256 :: MMINLU_loc
    431. 

  431.    CHARACTER*80 :: message
    432. 

  432.    INTEGER :: ISWATER
    433. 

  433.    INTEGER :: ISICE
    434. 

  434.    INTEGER :: ISURBAN
    435. 

  435.    INTEGER :: sf_urban_physics
    436. 

  436.    INTEGER :: omlcall
    437. 

  437.    REAL    :: oml_hml0
    438. 

  438.    LOGICAL :: usemonalb
    439. 

  439.    LOGICAL :: rdmaxalb
    440. 

  440.    INTEGER :: mfshconv
    441. 

  441. 

  442.    INTEGER :: i, j, k, itf, jtf, ktf, n
    443. 

  443. integer myproc
    444. 

  444. 

  445. !-----------------------------------------------------------------
    446. 

  446. 

  447.    sf_urban_physics=config_flags%sf_urban_physics
    448. 

  448.    usemonalb=config_flags%usemonalb
    449. 

  449.    rdmaxalb=config_flags%rdmaxalb
    450. 

  450.    mfshconv=config_flags%mfshconv
    451. 

  451. #if ( EM_CORE == 1 )
    452. 

  452.    obs_twindo_cg=model_config_rec%obs_twindo(1)
    453. 

  453.    obs_twindo=config_flags%obs_twindo
    454. 

  454.    oml_hml0=config_flags%oml_hml0
    455. 

  455.    omlcall=config_flags%omlcall
    456. 

  456. #endif
    457. 

  457. 

  458. !-- should be from the namelist
    459. 

  459. 

  460.    sfull = 0.
    461. 

  461.    shalf = 0.
    462. 

  462. 

  463.    CALL wrf_debug(100,'top of phy_init')
    464. 

  464. 

  465.    WRITE(wrf_err_message,*) 'phy_init:  start_of_simulation = ',start_of_simulation
    466. 

  466.    CALL wrf_debug ( 100, TRIM(wrf_err_message) )
    467. 

  467. 

  468.    itf=min0(ite,ide-1)
    469. 

  469.    jtf=min0(jte,jde-1)
    470. 

  470.    ktf=min0(kte,kde-1)
    471. 

  471. 

  472.    ZZLND=0.1
    473. 

  473.    ZZWTR=0.0001
    474. 

  474.    THINLD=0.04
    475. 

  475.    ALBLND=0.2
    476. 

  476.    XMAVA=0.3
    477. 

  477. 

  478. #if (NMM_CORE == 1)
    479. 

  479.    if (.not.usemonalb) CALL wrf_error_fatal('usemonalb should always be true for NMM')
    480. 

  480. #endif
    481. 

  481. 

  482.    CALL nl_get_cen_lat(id,cen_lat)
    483. 

  483.    CALL wrf_debug(100,'calling nl_get_iswater, nl_get_isice, nl_get_mminlu_loc')
    484. 

  484.    CALL nl_get_iswater(id,iswater)
    485. 

  485.    CALL nl_get_isice(id,isice)
    486. 

  486.    CALL nl_get_isurban(id,isurban)
    487. 

  487.    CALL nl_get_mminlu( 1, mminlu_loc )
    488. 

  488.    CALL wrf_debug(100,'after nl_get_iswater, nl_get_isice, nl_get_mminlu_loc')
    489. 

  489. !-- temporary fix by ww
    490. 

  490.    landuse_ISICE = isice
    491. 

  491. 

  492. ! Added for Wind Turbine parameterization code -- This will only read in an optional 
    493. 

  493. ! configuration file with information that will be used by inividual turbine init routines 
    494. 

  494. ! as each domain is initialized.
    495. 

  495.   IF ( id .EQ. 1 ) THEN
    496. 

  496.     CALL init_module_wind_generic
    497. 

  497.     CALL init_module_wind_fitch
    498. 

  498.   ENDIF
    499. 

  499. !
    500. 

  500. 

  501.   IF(.not.restart)THEN
    502. 

  502. !-- initialize common variables
    503. 

  503. 

  504.    IF ( .NOT. moved ) THEN
    505. 

  505.    DO j=jts,jtf
    506. 

  506.    DO i=its,itf
    507. 

  507.       XLAND(i,j)=1.
    508. 

  508.       GSW(i,j)=0.
    509. 

  509.       GLW(i,j)=0.
    510. 

  510. !-- initialize ust to a small value
    511. 

  511.       UST(i,j)=0.0001
    512. 

  512.       MOL(i,j)=0.0
    513. 

  513.       PBLH(i,j)=0.0
    514. 

  514.       HFX(i,j)=0.
    515. 

  515.       QFX(i,j)=0.
    516. 

  516.       RAINBL(i,j)=0.
    517. 

  517.       RAINNCV(i,j)=0.
    518. 

  518.       SNOWNCV(i,j)=0.
    519. 

  519.       GRAUPELNCV(i,j)=0.
    520. 

  520.       ACSNOW(i,j)=0.
    521. 

  521.       DO k=kms,kme  !wig, 17-May-2006: Added for idealized chem. runs
    522. 

  522.          EXCH_H(i,k,j) = 0.
    523. 

  523.       END DO
    524. 

  524.    ENDDO
    525. 

  525.    ENDDO
    526. 

  526.    ENDIF
    527. 

  527. !
    528. 

  528.    IF(PRESENT(TMN_UPDATE))THEN
    529. 

  529.    if(tmn_update.eq.1) then
    530. 

  530.    nyear=1
    531. 

  531.    nday=0.
    532. 

  532.    DO j=jts,jtf
    533. 

  533.    DO i=its,itf
    534. 

  534.       TYR(i,j)=TMN(i,j)
    535. 

  535.       TYRA(i,j)=0.0
    536. 

  536.       TDLY(i,j)=0.0
    537. 

  537.     DO n=1,lagday
    538. 

  538.       TLAG(i,n,j)=TMN(i,j)
    539. 

  539.     ENDDO
    540. 

  540.    ENDDO
    541. 

  541.    ENDDO
    542. 

  542.    endif
    543. 

  543.    ENDIF
    544. 

  544. !
    545. 

  545. 

  546. !
    547. 

  547.    DO j=jts,jtf
    548. 

  548.    DO i=its,itf
    549. 

  549.      IF(XLAND(i,j) .LT. 1.5)THEN
    550. 

  550.        IF(mminlu_loc .EQ. '    ') ALBBCK(i,j)=ALBLND
    551. 

  551.        EMBCK(i,j)=0.85
    552. 

  552.        ALBEDO(i,j)=ALBBCK(i,j)
    553. 

  553.        EMISS(i,j)=EMBCK(i,j)
    554. 

  554.        THC(i,j)=THINLD
    555. 

  555.        ZNT(i,j)=ZZLND
    556. 

  556. #if  ! ( NMM_CORE == 1 )
    557. 

  557.        Z0(i,j)=ZZLND
    558. 

  558. #endif
    559. 

  559.        MAVAIL(i,j)=XMAVA
    560. 

  560.      ELSE
    561. 

  561.        IF(mminlu_loc .EQ. '    ') ALBBCK(i,j)=0.08
    562. 

  562.        ALBEDO(i,j)=ALBBCK(i,j)
    563. 

  563.        EMBCK(i,j)=0.98
    564. 

  564.        EMISS(i,j)=EMBCK(i,j)
    565. 

  565.        THC(i,j)=THINLD
    566. 

  566.        ZNT(i,j)=ZZWTR
    567. 

  567. #if  ! ( NMM_CORE == 1 )
    568. 

  568.        Z0(i,j)=ZZWTR
    569. 

  569. #endif
    570. 

  570.        MAVAIL(i,j)=1.0
    571. 

  571.      ENDIF
    572. 

  572. 

  573.    ENDDO
    574. 

  574.    ENDDO
    575. 

  575. 

  576.     if (config_flags%cu_diag == 1 )then
    577. 

  577.     do j=jts,jtf
    578. 

  578.        do k=kts,ktf
    579. 

  579.           do i=its,itf
    580. 

  580.              gd_cloud(i,k,j) = 0.
    581. 

  581.              gd_cloud2(i,k,j) = 0.
    582. 

  582.              gd_cloud_a(i,k,j) = 0.
    583. 

  583.              gd_cloud2_a(i,k,j) = 0.
    584. 

  584.              gd_cloud_b(i,k,j) = 0.
    585. 

  585.              gd_cloud2_b(i,k,j) = 0.
    586. 

  586.           end do
    587. 

  587.        end do
    588. 

  588.     end do
    589. 

  589.     endif
    590. 

  590. 

  591.     do j=jts,jtf
    592. 

  592.          do i=its,itf
    593. 

  593.            raincv_a(i,j)=0.
    594. 

  594.            raincv_b(i,j)=0.
    595. 

  595.       end do
    596. 

  596.     end do
    597. 

  597. 

  598.    CALL wrf_debug ( 200 , 'module_start: phy_init: Before call to landuse_init' )
    599. 

  599. 

  600.    IF(mminlu_loc .ne. '    ')THEN
    601. 

  601. !-- initialize surface properties
    602. 

  602. 

  603.      CALL landuse_init(lu_index, snowc, albedo, albbck, snoalb, mavail, emiss, embck, &
    604. 

  604.                 znt, Z0, thc, xland, xice, xicem, julday, cen_lat, iswater, &
    605. 

  605.                 TRIM ( mminlu_loc ) ,                               &
    606. 

  606.                 landuse_ISICE, landuse_LUCATS,                      &
    607. 

  607.                 landuse_LUSEAS, landuse_ISN,                        &
    608. 

  608.                 config_flags%fractional_seaice,                      &
    609. 

  609.                 lu_state,                                           &
    610. 

  610.                 allowed_to_read , usemonalb ,                       &
    611. 

  611.                 ids, ide, jds, jde, kds, kde,                       &
    612. 

  612.                 ims, ime, jms, jme, kms, kme,                       &
    613. 

  613.                 its, ite, jts, jte, kts, kte                       )
    614. 

  614.    ENDIF
    615. 

  615. 

  616.   ENDIF
    617. 

  617. 

  618. !-- convert zfull and zhalf to sigma values for ra_init (Eta CO2 needs these)
    619. 

  619. !-- zfull/zhalf may be either zeta or eta
    620. 

  620. !-- what is done here depends on coordinate (check this code if adding new coordinates)
    621. 

  621.    CALL z2sigma(zfull,zhalf,sfull,shalf,p_top,pptop,config_flags, &
    622. 

  622.                 allowed_to_read,                                  &
    623. 

  623.                 kds,kde,kms,kme,kts,kte)
    624. 

  624. 

  625. !-- initialize physics
    626. 

  626. !-- ra: radiation
    627. 

  627. !-- bl: pbl
    628. 

  628. !-- cu: cumulus
    629. 

  629. !-- mp: microphysics
    630. 

  630. 

  631.    CALL wrf_debug ( 200 , 'module_start: phy_init: Before call to ra_init' )
    632. 

  632. 

  633.    CALL ra_init(id=id,STEPRA=STEPRA,RADT=RADT,DT=DT,RTHRATEN=RTHRATEN,RTHRATENLW=RTHRATENLW,             &
    634. 

  634.                 RTHRATENSW=RTHRATENSW,CLDFRA=CLDFRA,EMISS=EMISS,cen_lat=cen_lat,JULYR=JULYR,JULDAY=JULDAY,GMT=GMT,    &
    635. 

  635.                 levsiz=levsiz,XLAT=XLAT,n_ozmixm=n_ozmixm,                           &
    636. 

  636.                 cldfra_old=cldfra_old,                                     & ! Optional
    637. 

  637.                 ozmixm=ozmixm,pin=pin,                                     & ! Optional
    638. 

  638.                 m_ps_1=m_ps_1,m_ps_2=m_ps_2,m_hybi=m_hybi,aerosolc_1=aerosolc_1,aerosolc_2=aerosolc_2,     & ! Optional
    639. 

  639.                 paerlev=paerlev,n_aerosolc=n_aerosolc,                             &
    640. 

  640.                 sfull=sfull,shalf=shalf,pptop=pptop,swrad_scat=swrad_scat,p_top=p_top,       &
    641. 

  641.                 config_flags=config_flags,restart=restart,                           &
    642. 

  642.                 allowed_to_read=allowed_to_read, start_of_simulation=start_of_simulation,           &
    643. 

  643.                 ids=ids, ide=ide, jds=jds, jde=jde, kds=kds, kde=kde,                   &
    644. 

  644.                 ims=ims, ime=ime, jms=jms, jme=jme, kms=kms, kme=kme,                   &
    645. 

  645.                 its=its, ite=ite, jts=jts, jte=jte, kts=kts, kte=kte                    )
    646. 

  646. 

  647.    CALL wrf_debug ( 200 , 'module_start: phy_init: Before call to bl_init' )
    648. 

  648.    CALL bl_init(STEPBL,BLDT,DT,RUBLTEN,RVBLTEN,RTHBLTEN,        &
    649. 

  649.                 RQVBLTEN,RQCBLTEN,RQIBLTEN,TSK,TMN,             &
    650. 

  650.                 config_flags,restart,UST,LOWLYR,TSLB,ZS,DZS,    &
    651. 

  651.                 num_soil_layers,TKE_PBL,mfshconv,               &
    652. 

  652.                 massflux_EDKF, entr_EDKF, detr_EDKF, & 
    653. 

  653.                 thl_up, thv_up, rt_up,       &
    654. 

  654.                 rv_up, rc_up, u_up, v_up,    &
    655. 

  655.                 frac_up, &
    656. 

  656.                 EXCH_H,VEGFRA,                                  &
    657. 

  657.                 SNOW,SNOWC, CANWAT,SMSTAV,                      &
    658. 

  658.                 SMSTOT, SFCRUNOFF,UDRUNOFF,ACSNOW,ACSNOM,       &
    659. 

  659.                 IVGTYP,ISLTYP,ISURBAN,SMOIS,SMFR3D,MAVAIL,      &
    660. 

  660.                 SNOWH,SH2O,SNOALB,FNDSOILW,FNDSNOWH,RDMAXALB,   &
    661. 

  661. #if (NMM_CORE == 1)
    662. 

  662.                 Z0,XLAND,XICE,                                  &
    663. 

  663. #else
    664. 

  664.                 ZNT,XLAND,XICE,                                 &
    665. 

  665. #endif
    666. 

  666.                 SFCEVP,GRDFLX,                                  &
    667. 

  667.                 TRIM (MMINLU_LOC),                              &
    668. 

  668.                 ISNOWXY, ZSNSOXY, TSNOXY,                       &
    669. 

  669.                 SNICEXY, SNLIQXY, TVXY, TGXY, CANICEXY,         &
    670. 

  670.                 CANLIQXY, EAHXY, TAHXY, CMXY,                   &
    671. 

  671.                 CHXY, FWETXY, SNEQVOXY, ALBOLDXY, QSNOWXY,      &
    672. 

  672.                 WSLAKEXY, ZWTXY, WAXY, WTXY, LFMASSXY, RTMASSXY,&
    673. 

  673.                 STMASSXY, WOODXY, STBLCPXY, FASTCPXY,           &
    674. 

  674.                 XSAIXY,                                         &
    675. 

  675.                 T2MVXY,T2MBXY,CHSTARXY ,                        &
    676. 

  676.                 allowed_to_read ,                               &
    677. 

  677.                 start_of_simulation ,                           &
    678. 

  678.                 te_temf,cf3d_temf,wm_temf,                      & ! WA
    679. 

  679.                 DZR, DZB, DZG,                                  & !Optional urban
    680. 

  680.                 TR_URB2D,TB_URB2D,TG_URB2D,TC_URB2D,QC_URB2D,   & !Optional urban
    681. 

  681.                 XXXR_URB2D,XXXB_URB2D,XXXG_URB2D,XXXC_URB2D,    & !Optional urban
    682. 

  682.                 TRL_URB3D, TBL_URB3D, TGL_URB3D,                & !Optional urban
    683. 

  683.                 SH_URB2D, LH_URB2D, G_URB2D, RN_URB2D,          & !Optional urban
    684. 

  684.                 TS_URB2D, FRC_URB2D, UTYPE_URB2D,               & 
    685. 

  685.                 SF_URBAN_PHYSICS,                               & !Optional urban
    686. 

  686.                 NUM_URBAN_LAYERS,                               & !Optional multi-layer urban
    687. 

  687.                 TRB_URB4D,TW1_URB4D,TW2_URB4D,                  & !Optional multi-layer urban
    688. 

  688.                 TGB_URB4D,TLEV_URB3D,QLEV_URB3D,                & !Optional multi-layer urban
    689. 

  689.                 TW1LEV_URB3D,TW2LEV_URB3D,                      & !Optional multi-layer urban
    690. 

  690.                 TGLEV_URB3D,TFLEV_URB3D,                        & !Optional multi-layer urban
    691. 

  691.                 SF_AC_URB3D,LF_AC_URB3D,CM_AC_URB3D,            & !Optional multi-layer urban
    692. 

  692.                 SFVENT_URB3D,LFVENT_URB3D,                      & !Optional multi-layer urban
    693. 

  693.                 SFWIN1_URB3D,SFWIN2_URB3D,                      & !Optional multi-layer urban
    694. 

  694.                 SFW1_URB3D,SFW2_URB3D,                          & !Optional multi-layer urban
    695. 

  695.                 SFR_URB3D,SFG_URB3D,                            & !Optional multi-layer urban
    696. 

  696.                 A_U_BEP,A_V_BEP,A_T_BEP,A_Q_BEP,                & !Optional multi-layer urban
    697. 

  697.                 A_E_BEP,B_U_BEP,B_V_BEP,                        & !Optional multi-layer urban
    698. 

  698.                 B_T_BEP,B_Q_BEP,B_E_BEP,DLG_BEP,                & !Optional multi-layer urban
    699. 

  699.                 DL_U_BEP,SF_BEP,VL_BEP,                         & !Optional multi-layer urban
    700. 

  700.                 ids, ide, jds, jde, kds, kde,                   &
    701. 

  701.                 ims, ime, jms, jme, kms, kme,                   &
    702. 

  702.                 its, ite, jts, jte, kts, kte,                   &
    703. 

  703.                 ACHFX,ACLHF,ACGRDFLX,                           &
    704. 

  704.                 oml_hml0, omlcall,                              & !Optional oml
    705. 

  705.                 TML,T0ML,HML,H0ML,HUML,HVML,TMOML               ) !Optional oml
    706. 

  706. 

  707.    CALL wrf_debug ( 200 , 'module_start: phy_init: Before call to cu_init' )
    708. 

  708. 

  709.    CALL cu_init(STEPCU,CUDT,DT,RUCUTEN,RVCUTEN,RTHCUTEN,        &
    710. 

  710.                 RQVCUTEN,RQRCUTEN,RQCCUTEN,RQSCUTEN,RQICUTEN,   &
    711. 

  711.                 NCA,RAINC,RAINCV,W0AVG,config_flags,restart,    &
    712. 

  712.                 CLDEFI,LOWLYR,MASS_FLUX,                        &
    713. 

  713.                 RTHFTEN, RQVFTEN,                               &
    714. 

  714.                 APR_GR,APR_W,APR_MC,APR_ST,APR_AS,              &
    715. 

  715.                 APR_CAPMA,APR_CAPME,APR_CAPMI,                  &
    716. 

  716.                 cugd_tten,cugd_ttens,cugd_qvten,                &
    717. 

  717.                 cugd_qvtens,cugd_qcten,                         &
    718. 

  718.                 allowed_to_read, start_of_simulation,           &
    719. 

  719.                 ids, ide, jds, jde, kds, kde,                   &
    720. 

  720.                 ims, ime, jms, jme, kms, kme,                   &
    721. 

  721.                 its, ite, jts, jte, kts, kte                    )
    722. 

  722. 

  723.    CALL wrf_debug ( 200 , 'module_start: phy_init: Before call to shcu_init' )
    724. 

  724. 

  725.    CALL shcu_init(STEPCU,CUDT,DT,RUSHTEN,RVSHTEN,RTHSHTEN,      &
    726. 

  726.                 RQVSHTEN,RQRSHTEN,RQCSHTEN,                     &
    727. 

  727.                 RQSSHTEN,RQISHTEN,RQGSHTEN,                     &
    728. 

  728.                 NCA,RAINC,RAINCV,config_flags,restart,          &
    729. 

  729.                 allowed_to_read, start_of_simulation,           &
    730. 

  730.                 ids, ide, jds, jde, kds, kde,                   &
    731. 

  731.                 ims, ime, jms, jme, kms, kme,                   &
    732. 

  732.                 its, ite, jts, jte, kts, kte                    )
    733. 

  733. 

  734.    CALL wrf_debug ( 200 , 'module_start: phy_init: Before call to mp_init' )
    735. 

  735. 

  736.    CALL mp_init(RAINNC,SNOWNC,GRAUPELNC,config_flags,restart,warm_rain,          &
    737. 

  737.                 adv_moist_cond,                                 &
    738. 

  738.                 MPDT, DT, DX, DY, LOWLYR,                       &
    739. 

  739.                 F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY,               &
    740. 

  740.                 mp_restart_state,tbpvs_state,tbpvs0_state,      &
    741. 

  741.                 allowed_to_read, start_of_simulation,           &
    742. 

  742.                 ids, ide, jds, jde, kds, kde,                   &
    743. 

  743.                 ims, ime, jms, jme, kms, kme,                   &
    744. 

  744.                 its, ite, jts, jte, kts, kte                    )
    745. 

  745. 

  746. #if  ( EM_CORE == 1 )
    747. 

  747.    CALL wrf_debug ( 200 , 'module_start: phy_init: Before call to fg_init' )
    748. 

  748. 

  749.    CALL fg_init(STEPFG,FGDT,DT,id,RUNDGDTEN,RVNDGDTEN,          &
    750. 

  750.                 RTHNDGDTEN,RPHNDGDTEN,RQVNDGDTEN,RMUNDGDTEN,    &
    751. 

  751.                 config_flags,restart,                           &
    752. 

  752.                 allowed_to_read ,                               &
    753. 

  753.                 ids, ide, jds, jde, kds, kde,                   &
    754. 

  754.                 ims, ime, jms, jme, kms, kme,                   &
    755. 

  755.                 its, ite, jts, jte, kts, kte                    )
    756. 

  756. 

  757.    CALL wrf_debug ( 200 , 'module_start: phy_init: Before call to fdob_init' )
    758. 

  758. 

  759.    CALL fdob_init(model_config_rec%obs_nudge_opt,               &
    760. 

  760.                   model_config_rec%max_dom,                     &
    761. 

  761.                   id,                                           &
    762. 

  762.                   model_config_rec%parent_id,                   &
    763. 

  763.                   model_config_rec%obs_idynin,                  &
    764. 

  764.                   model_config_rec%obs_dtramp,                  &
    765. 

  765.                   model_config_rec%fdda_end,                    &
    766. 

  766.                   model_config_rec%restart,                     &
    767. 

  767.                   obs_twindo_cg, obs_twindo,                    &
    768. 

  768.                   itimestep,                                    &
    769. 

  769.                   model_config_rec%obs_no_pbl_nudge_uv,         &
    770. 

  770.                   model_config_rec%obs_no_pbl_nudge_t,          &
    771. 

  771.                   model_config_rec%obs_no_pbl_nudge_q,          &
    772. 

  772.                   model_config_rec%obs_sfc_scheme_horiz,        &
    773. 

  773.                   model_config_rec%obs_sfc_scheme_vert,         &
    774. 

  774.                   model_config_rec%obs_max_sndng_gap,           &
    775. 

  775.                   model_config_rec%obs_sfcfact,                 &
    776. 

  776.                   model_config_rec%obs_sfcfacr,                 &
    777. 

  777.                   model_config_rec%obs_dpsmx,                   &
    778. 

  778.                   model_config_rec%obs_nudge_wind,              &
    779. 

  779.                   model_config_rec%obs_nudge_temp,              &
    780. 

  780.                   model_config_rec%obs_nudge_mois,              &
    781. 

  781.                   model_config_rec%obs_nudgezfullr1_uv,         &
    782. 

  782.                   model_config_rec%obs_nudgezrampr1_uv,         &
    783. 

  783.                   model_config_rec%obs_nudgezfullr2_uv,         &
    784. 

  784.                   model_config_rec%obs_nudgezrampr2_uv,         &
    785. 

  785.                   model_config_rec%obs_nudgezfullr4_uv,         &
    786. 

  786.                   model_config_rec%obs_nudgezrampr4_uv,         &
    787. 

  787.                   model_config_rec%obs_nudgezfullr1_t,          &
    788. 

  788.                   model_config_rec%obs_nudgezrampr1_t,          &
    789. 

  789.                   model_config_rec%obs_nudgezfullr2_t,          &
    790. 

  790.                   model_config_rec%obs_nudgezrampr2_t,          &
    791. 

  791.                   model_config_rec%obs_nudgezfullr4_t,          &
    792. 

  792.                   model_config_rec%obs_nudgezrampr4_t,          &
    793. 

  793.                   model_config_rec%obs_nudgezfullr1_q,          &
    794. 

  794.                   model_config_rec%obs_nudgezrampr1_q,          &
    795. 

  795.                   model_config_rec%obs_nudgezfullr2_q,          &
    796. 

  796.                   model_config_rec%obs_nudgezrampr2_q,          &
    797. 

  797.                   model_config_rec%obs_nudgezfullr4_q,          &
    798. 

  798.                   model_config_rec%obs_nudgezrampr4_q,          &
    799. 

  799.                   model_config_rec%obs_nudgezfullmin,           &
    800. 

  800.                   model_config_rec%obs_nudgezrampmin,           &
    801. 

  801.                   model_config_rec%obs_nudgezmax,               &
    802. 

  802.                   xlat,                                         &
    803. 

  803.                   xlong,                                        &
    804. 

  804.                   model_config_rec%start_year(id),              &
    805. 

  805.                   model_config_rec%start_month(id),             &
    806. 

  806.                   model_config_rec%start_day(id),               &
    807. 

  807.                   model_config_rec%start_hour(id),              &
    808. 

  808.                   model_config_rec%start_minute(id),            &
    809. 

  809.                   model_config_rec%start_second(id),            &
    810. 

  810.                   p00, t00, tlp,                                &
    811. 

  811.                   zhalf, p_top,                                 &
    812. 

  812.                   fdob,                                         &
    813. 

  813.                   model_config_rec%obs_ipf_init,                &
    814. 

  814.                   ids, ide, jds, jde, kds, kde,                 &
    815. 

  815.                   ims, ime, jms, jme, kms, kme,                 &
    816. 

  816.                   its, ite, jts, jte, kts, kte                  )
    817. 

  817. 

  818. #endif
    819. 

  819. 

  820.    END SUBROUTINE phy_init
    821. 

  821. 

  822. !=====================================================================
    823. 

  823.    SUBROUTINE landuse_init(lu_index, snowc, albedo, albbck, snoalb, mavail, emiss, embck, &
    824. 

  824.                 znt,Z0,thc,xland, xice, xicem, julday, cen_lat, iswater, mminlu,  &
    825. 

  825.                 ISICE, LUCATS, LUSEAS, ISN,                         &
    826. 

  826.                 fractional_seaice,                                  &
    827. 

  827.                 lu_state,                                           &
    828. 

  828.                 allowed_to_read , usemonalb ,                       &
    829. 

  829.                 ids, ide, jds, jde, kds, kde,                       &
    830. 

  830.                 ims, ime, jms, jme, kms, kme,                       &
    831. 

  831.                 its, ite, jts, jte, kts, kte                       )
    832. 

  832. 

  833.    USE module_wrf_error
    834. 

  834.    IMPLICIT NONE
    835. 

  835. 

  836. !---------------------------------------------------------------------
    837. 

  837.    INTEGER , INTENT(IN)           :: ids, ide, jds, jde, kds, kde,   &
    838. 

  838.                                      ims, ime, jms, jme, kms, kme,   &
    839. 

  839.                                      its, ite, jts, jte, kts, kte
    840. 

  840. 

  841.    INTEGER , INTENT(IN)           :: iswater, julday
    842. 

  842.    REAL    , INTENT(IN)           :: cen_lat
    843. 

  843.    CHARACTER(LEN=*), INTENT(IN)        :: mminlu
    844. 

  844.    LOGICAL,  INTENT(IN)           :: allowed_to_read , usemonalb
    845. 

  845.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(IN   ) :: lu_index, snowc, xice, snoalb 
    846. 

  846.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(OUT  ) :: albedo, albbck, mavail, emiss, &
    847. 

  847.                                                                embck,                         &
    848. 

  848.                                                                znt, Z0, thc, xland, xicem
    849. 

  849.    INTEGER , INTENT(INOUT)  :: ISICE, LUCATS, LUSEAS, ISN, fractional_seaice
    850. 

  850.    REAL    , INTENT(INOUT)  , DIMENSION( : ) :: lu_state
    851. 

  851. 

  852.    REAL :: xice_threshold
    853. 

  853. !---------------------------------------------------------------------
    854. 

  854. ! Local
    855. 

  855.    CHARACTER*256 LUTYPE
    856. 

  856.    CHARACTER*80 :: message
    857. 

  857.    INTEGER  :: landuse_unit, LS, LC, LI, LUN, NSN
    858. 

  858.    INTEGER  :: i, j, itf, jtf, is, cats, seas, curs
    859. 

  859.    INTEGER , PARAMETER :: OPEN_OK = 0
    860. 

  860.    INTEGER :: ierr
    861. 

  861.    INTEGER , PARAMETER :: max_cats = 100 , max_seas = 12
    862. 

  862.    REAL    , DIMENSION( max_cats, max_seas ) :: ALBD, SLMO, SFEM, SFZ0, THERIN, SFHC
    863. 

  863.    REAL    , DIMENSION( max_cats )     :: SCFX
    864. 

  864. ! save these fields in case nest moves or has to be reinitialized
    865. 

  865. ! and this routine is called with allowed_to_read set to false
    866. 

  866. ! note that by saving these, we're locking in the same landuse for
    867. 

  867. ! the duration of a run; possible implications for long climate runs
    868. 

  868.    LOGICAL :: found_lu, end_of_file
    869. 

  869.    LOGICAL, EXTERNAL :: wrf_dm_on_monitor
    870. 

  870. 

  871. !---------------------------------------------------------------------
    872. 

  872. 

  873.    CALL wrf_debug( 100 , 'top of landuse_init' )
    874. 

  874. 

  875.    NSN=-1  ! set this to suppress uninitalized data messages from tools
    876. 

  876. 

  877.   if ( fractional_seaice == 0 ) then
    878. 

  878.      xice_threshold = 0.5
    879. 

  879.   else if ( fractional_seaice == 1 ) then
    880. 

  880.      xice_threshold = 0.02
    881. 

  881.   endif
    882. 

  882. 

  883. ! recover LU variables from state
    884. 

  884.    IF ( 6*(max_cats*max_seas)+1*max_cats .GT. 7501 ) THEN
    885. 

  885.       WRITE(message,*)'landuse_init: lu_state overflow. Make Registry dimspec p > ',6*(max_cats*max_seas)+1*max_cats
    886. 

  886.    ENDIF
    887. 

  887.    curs = 1
    888. 

  888.    DO cats = 1, max_cats
    889. 

  889.      SCFX(cats) =           lu_state(curs)         ; curs = curs + 1
    890. 

  890.      DO seas = 1, max_seas
    891. 

  891.        ALBD(cats,seas) =    lu_state(curs)         ; curs = curs + 1
    892. 

  892.        SLMO(cats,seas) =    lu_state(curs)         ; curs = curs + 1
    893. 

  893.        SFEM(cats,seas) =    lu_state(curs)         ; curs = curs + 1
    894. 

  894.        SFZ0(cats,seas) =    lu_state(curs)         ; curs = curs + 1
    895. 

  895.        SFHC(cats,seas) =    lu_state(curs)         ; curs = curs + 1
    896. 

  896.        THERIN(cats,seas) =  lu_state(curs)         ; curs = curs + 1
    897. 

  897.      ENDDO
    898. 

  898.    ENDDO
    899. 

  899. 

  900. ! Determine season (summer=1, winter=2)
    901. 

  901.    ISN=1
    902. 

  902.    IF(JULDAY.LT.105.OR.JULDAY.GT.288)ISN=2
    903. 

  903.    IF(CEN_LAT.LT.0.0)ISN=3-ISN
    904. 

  904. 

  905.    FOUND_LU = .TRUE.
    906. 

  906.    IF ( allowed_to_read ) THEN
    907. 

  907.       landuse_unit = 29
    908. 

  908.       IF ( wrf_dm_on_monitor() ) THEN
    909. 

  909.         OPEN(landuse_unit, FILE='LANDUSE.TBL',FORM='FORMATTED',STATUS='OLD',IOSTAT=ierr)
    910. 

  910.         IF ( ierr .NE. OPEN_OK ) THEN
    911. 

  911.           WRITE(message,FMT='(A)') &
    912. 

  912.           'module_physics_init.F: LANDUSE_INIT: open failure for LANDUSE.TBL'
    913. 

  913.           CALL wrf_error_fatal ( message )
    914. 

  914.         END IF
    915. 

  915.       ENDIF
    916. 

  916. 

  917. ! Read info from file LANDUSE.TBL
    918. 

  918. !     IF(MMINLU.EQ.'OLD ')THEN
    919. 

  919. !       ISWATER=7
    920. 

  920. !       ISICE=11
    921. 

  921. !     ELSE IF(MMINLU.EQ.'USGS')THEN
    922. 

  922. !       ISWATER=16
    923. 

  923. !       ISICE=24
    924. 

  924. !     ELSE IF(MMINLU.EQ.'SiB ')THEN
    925. 

  925. !       ISWATER=15
    926. 

  926. !       ISICE=16
    927. 

  927. !     ELSE IF(MMINLU.EQ.'LW12')THEN
    928. 

  928. !       ISWATER=15
    929. 

  929. !       ISICE=3
    930. 

  930. !     ELSE IF (MMINLU .EQ. 'MODIFIED_IGBP_MODIS_NOAH') THEN
    931. 

  931. !       ISICE = 15
    932. 

  932. !     ELSE
    933. 

  933. !        call wrf_error_fatal ("INPUT LandUse not found: "//TRIM(MMINLU))
    934. 

  934. !     ENDIF
    935. 

  935.       call wrf_message ( 'INPUT LandUse = "' // TRIM(MMINLU) // '"' )
    936. 

  936.       FOUND_LU = .FALSE.
    937. 

  937.       end_of_file = .FALSE.
    938. 

  938. !!! BEGINNING OF 1999 LOOP
    939. 

  939.  1999 CONTINUE
    940. 

  940.       IF ( wrf_dm_on_monitor() ) THEN
    941. 

  941.         READ (landuse_unit,*,END=2002)LUTYPE
    942. 

  942.         GOTO 2003
    943. 

  943.  2002   CONTINUE
    944. 

  944.         CALL wrf_message( 'INPUT FILE FOR LANDUSE REACHED END OF FILE' )
    945. 

  945.         end_of_file = .TRUE.
    946. 

  946.  2003   CONTINUE
    947. 

  947.         IF ( .NOT. end_of_file ) READ (landuse_unit,*)LUCATS,LUSEAS
    948. 

  948.         FOUND_LU = LUTYPE.EQ.MMINLU
    949. 

  949.       ENDIF
    950. 

  950.       CALL wrf_dm_bcast_bytes (end_of_file, LWORDSIZE )
    951. 

  951.       IF ( .NOT. end_of_file ) THEN
    952. 

  952.         CALL wrf_dm_bcast_string(lutype, 256)
    953. 

  953.         CALL wrf_dm_bcast_bytes (lucats,  IWORDSIZE )
    954. 

  954.         CALL wrf_dm_bcast_bytes (luseas,  IWORDSIZE )
    955. 

  955.         CALL wrf_dm_bcast_bytes (found_lu,  LWORDSIZE )
    956. 

  956.         IF(FOUND_LU)THEN
    957. 

  957.           LUN=LUCATS
    958. 

  958.           NSN=LUSEAS
    959. 

  959.             IF(LUTYPE.NE.'SSIB') THEN !this is not really true for ssib lsm (fds)
    960. 

  960.             PRINT *, 'LANDUSE TYPE = "' // TRIM (LUTYPE) // '" FOUND',        &
    961. 

  961.                    LUCATS,' CATEGORIES',LUSEAS,' SEASONS',     &
    962. 

  962.                    ' WATER CATEGORY = ',ISWATER,               &
    963. 

  963.                    ' SNOW CATEGORY = ',ISICE
    964. 

  964.             ENDIF
    965. 

  965.         ENDIF
    966. 

  966.         DO ls=1,luseas
    967. 

  967.           if ( wrf_dm_on_monitor() ) then
    968. 

  968.             READ (landuse_unit,*)
    969. 

  969.           endif
    970. 

  970.           DO LC=1,LUCATS
    971. 

  971.             IF(found_lu)THEN
    972. 

  972.               IF ( wrf_dm_on_monitor() ) THEN
    973. 

  973.                 READ (landuse_unit,*)LI,ALBD(LC,LS),SLMO(LC,LS),SFEM(LC,LS),        &
    974. 

  974.                            SFZ0(LC,LS),THERIN(LC,LS),SCFX(LC),SFHC(LC,LS)
    975. 

  975.               ENDIF
    976. 

  976.               CALL wrf_dm_bcast_bytes (LI,  IWORDSIZE )
    977. 

  977.               IF(LC.NE.LI)CALL wrf_error_fatal ( 'module_start: MISSING LANDUSE UNIT ' )
    978. 

  978.             ELSE
    979. 

  979.               IF ( wrf_dm_on_monitor() ) THEN
    980. 

  980.                 READ (landuse_unit,*)
    981. 

  981.               ENDIF
    982. 

  982.             ENDIF
    983. 

  983.           ENDDO
    984. 

  984.         ENDDO
    985. 

  985.         IF(NSN.EQ.1.AND.FOUND_LU) THEN
    986. 

  986.            ISN = 1
    987. 

  987.         END IF
    988. 

  988.         CALL wrf_dm_bcast_bytes (albd,   max_cats * max_seas * RWORDSIZE )
    989. 

  989.         CALL wrf_dm_bcast_bytes (slmo,   max_cats * max_seas * RWORDSIZE )
    990. 

  990.         CALL wrf_dm_bcast_bytes (sfem,   max_cats * max_seas * RWORDSIZE )
    991. 

  991.         CALL wrf_dm_bcast_bytes (sfz0,   max_cats * max_seas * RWORDSIZE )
    992. 

  992.         CALL wrf_dm_bcast_bytes (therin, max_cats * max_seas * RWORDSIZE )
    993. 

  993.         CALL wrf_dm_bcast_bytes (sfhc,   max_cats * max_seas * RWORDSIZE )
    994. 

  994.         CALL wrf_dm_bcast_bytes (scfx,   max_cats *            RWORDSIZE )
    995. 

  995.       ENDIF
    996. 

  996. 

  997.       IF(.NOT. found_lu .AND. .NOT. end_of_file ) GOTO 1999
    998. 

  998. !!! END OF 1999 LOOP
    999. 

  999. 

  1000.       IF(.NOT. found_lu .OR. end_of_file )THEN
    1001. 

  1001.         CALL wrf_message ( 'LANDUSE IN INPUT FILE DOES NOT MATCH LUTABLE: TABLE NOT USED' )
    1002. 

  1002.       ENDIF
    1003. 

  1003.     ENDIF  ! allowed_to_read
    1004. 

  1004. 

  1005.     IF(FOUND_LU)THEN
    1006. 

  1006. ! Set arrays according to lu_index
    1007. 

  1007.       itf = min0(ite, ide-1)
    1008. 

  1008.       jtf = min0(jte, jde-1)
    1009. 

  1009.       IF(usemonalb)CALL wrf_message ( 'Climatological albedo is used instead of table values' )
    1010. 

  1010.       DO j = jts, jtf
    1011. 

  1011.         DO i = its, itf
    1012. 

  1012.           IS=nint(lu_index(i,j))
    1013. 

  1013.           ! only do this check on read-in data
    1014. 

  1014.           IF(allowed_to_read)THEN
    1015. 

  1015.              IF(IS.LT.0.OR.IS.GT.LUN)THEN
    1016. 

  1016.                WRITE ( wrf_err_message , * ) 'ERROR: LANDUSE OUTSIDE RANGE =',IS,' AT ',I,J,' LUN= ',LUN
    1017. 

  1017.                CALL wrf_error_fatal ( TRIM ( wrf_err_message ) )
    1018. 

  1018.              ENDIF
    1019. 

  1019.           ENDIF
    1020. 

  1020. !   SET NO-DATA POINTS (IS=0) TO WATER
    1021. 

  1021.           IF(IS.EQ.0)THEN
    1022. 

  1022.             IS=ISWATER
    1023. 

  1023.           ENDIF
    1024. 

  1024.           IF(.NOT.usemonalb)ALBBCK(I,J)=ALBD(IS,ISN)/100.
    1025. 

  1025.           ALBEDO(I,J)=ALBBCK(I,J)
    1026. 

  1026.           IF(SNOWC(I,J) .GT. 0.5) THEN 
    1027. 

  1027.              IF (usemonalb) THEN
    1028. 

  1028.                  ALBEDO(I,J)=SNOALB(I,J)
    1029. 

  1029.              ELSE
    1030. 

  1030.                  ALBEDO(I,J)=ALBBCK(I,J)*(1.+SCFX(IS))
    1031. 

  1031.              ENDIF
    1032. 

  1032.           ENDIF
    1033. 

  1033.           THC(I,J)=THERIN(IS,ISN)/100.
    1034. 

  1034.           Z0(I,J)=SFZ0(IS,ISN)/100.
    1035. 

  1035.           ZNT(I,J)=Z0(I,J)
    1036. 

  1036.           EMBCK(I,J)=SFEM(IS,ISN)
    1037. 

  1037.           EMISS(I,J)=EMBCK(I,J)
    1038. 

  1038.           MAVAIL(I,J)=SLMO(IS,ISN)
    1039. 

  1039.           IF(IS.NE.ISWATER)THEN
    1040. 

  1040.             XLAND(I,J)=1.0
    1041. 

  1041.           ELSE
    1042. 

  1042.             XLAND(I,J)=2.0
    1043. 

  1043.           ENDIF
    1044. 

  1044. !    SET SEA-ICE POINTS TO LAND WITH ICE/SNOW SURFACE PROPERTIES
    1045. 

  1045.           XICEM(I,J)=XICE(I,J)
    1046. 

  1046.           IF(XICE(I,J).GE.xice_threshold)THEN
    1047. 

  1047.             XLAND(I,J)=1.0
    1048. 

  1048.             ALBBCK(I,J)=ALBD(ISICE,ISN)/100.
    1049. 

  1049.             EMBCK(I,J)=SFEM(ISICE,ISN)
    1050. 

  1050.             IF (FRACTIONAL_SEAICE == 1) THEN
    1051. 

  1051.                ! The 0.08 value is the albedo over open water.
    1052. 

  1052.                ! The 0.98 value is the emissivity over open water.
    1053. 

  1053.                ALBEDO(I,J) = ( XICE(I,J) * ALBBCK(I,J) ) + ( (1.0-XICE(I,J)) * 0.08 )
    1054. 

  1054.                EMISS(I,J)  = ( XICE(I,J) * EMBCK(I,J)  ) + ( (1.0-XICE(I,J)) * 0.98 )
    1055. 

  1055.             ELSE
    1056. 

  1056.                ALBEDO(I,J)=ALBBCK(I,J)
    1057. 

  1057.                EMISS(I,J)=EMBCK(I,J)
    1058. 

  1058.             ENDIF
    1059. 

  1059.             THC(I,J)=THERIN(ISICE,ISN)/100.
    1060. 

  1060.             Z0(I,J)=SFZ0(ISICE,ISN)/100.
    1061. 

  1061.             ZNT(I,J)=Z0(I,J)
    1062. 

  1062.             MAVAIL(I,J)=SLMO(ISICE,ISN)
    1063. 

  1063.           ENDIF
    1064. 

  1064.         ENDDO
    1065. 

  1065.       ENDDO
    1066. 

  1066.     ENDIF
    1067. 

  1067.     if ( wrf_dm_on_monitor() .and. allowed_to_read ) then
    1068. 

  1068.       CLOSE (landuse_unit)
    1069. 

  1069.     endif
    1070. 

  1070.     CALL wrf_debug( 100 , 'returning from of landuse_init' )
    1071. 

  1071. 

  1072. ! restore LU variables from state
    1073. 

  1073.     curs = 1
    1074. 

  1074.     DO cats = 1, max_cats
    1075. 

  1075.       lu_state(curs) = SCFX(cats)                 ; curs = curs + 1
    1076. 

  1076.       DO seas = 1, max_seas
    1077. 

  1077.         lu_state(curs) = ALBD(cats,seas)          ; curs = curs + 1
    1078. 

  1078.         lu_state(curs) = SLMO(cats,seas)          ; curs = curs + 1
    1079. 

  1079.         lu_state(curs) = SFEM(cats,seas)          ; curs = curs + 1
    1080. 

  1080.         lu_state(curs) = SFZ0(cats,seas)          ; curs = curs + 1
    1081. 

  1081.         lu_state(curs) = SFHC(cats,seas)          ; curs = curs + 1
    1082. 

  1082.         lu_state(curs) = THERIN(cats,seas)        ; curs = curs + 1
    1083. 

  1083.       ENDDO
    1084. 

  1084.     ENDDO
    1085. 

  1085. !
    1086. 

  1086. !-- fds (07/2010)
    1087. 

  1087. !-- SSIB's 12-category vegetation parameters are defined in module_ssib_veg
    1088. 

  1088. !-- Data in LANDUSE.TBL is only used temporarily until the first
    1089. 

  1089. !-- call to SSIB, and also to set water/ice points properties
    1090. 

  1090. !    IF(MMINLU.EQ.'SSIB')THEN
    1091. 

  1091. !       CALL init_module_ssib_veg
    1092. 

  1092. !       CALL wrf_message ( 'READING SSIB VEGETATION PARAMETERS' )
    1093. 

  1093. !    ENDIF
    1094. 

  1094. !
    1095. 

  1095. 

  1096.    END SUBROUTINE landuse_init
    1097. 

  1097. 

  1098. !=====================================================================
    1099. 

  1099.    SUBROUTINE ra_init(id,STEPRA,RADT,DT,RTHRATEN,RTHRATENLW,       &
    1100. 

  1100.                       RTHRATENSW,CLDFRA,EMISS,cen_lat,JULYR,JULDAY,GMT,    &
    1101. 

  1101.                       levsiz,XLAT,n_ozmixm,                           &
    1102. 

  1102.                       cldfra_old,                                     & ! Optional
    1103. 

  1103.                       ozmixm,pin,                                     & ! Optional
    1104. 

  1104.                       m_ps_1,m_ps_2,m_hybi,aerosolc_1,aerosolc_2,     & ! Optional
    1105. 

  1105.                       paerlev,n_aerosolc,                             &
    1106. 

  1106.                       sfull,shalf,pptop,swrad_scat,p_top,            &
    1107. 

  1107.                       config_flags,restart,                          &
    1108. 

  1108.                       allowed_to_read, start_of_simulation,          &
    1109. 

  1109.                       ids, ide, jds, jde, kds, kde,                  &
    1110. 

  1110.                       ims, ime, jms, jme, kms, kme,                  &
    1111. 

  1111.                       its, ite, jts, jte, kts, kte                   )
    1112. 

  1112. !---------------------------------------------------------------------
    1113. 

  1113.    USE module_ra_rrtm      , ONLY : rrtminit
    1114. 

  1114.    USE module_ra_rrtmg_lw  , ONLY : rrtmg_lwinit
    1115. 

  1115.    USE module_ra_rrtmg_sw  , ONLY : rrtmg_swinit
    1116. 

  1116.    USE module_ra_cam       , ONLY : camradinit
    1117. 

  1117.    USE module_ra_sw        , ONLY : swinit
    1118. 

  1118.    USE module_ra_gsfcsw    , ONLY : gsfc_swinit
    1119. 

  1119.    USE module_ra_gfdleta   , ONLY : gfdletainit
    1120. 

  1120. #if(NMM_CORE==1)
    1121. 

  1121.    USE module_ra_hwrf      , ONLY : hwrfrainit
    1122. 

  1122. #endif
    1123. 

  1123.    USE module_ra_hs        , ONLY : hsinit
    1124. 

  1124.    USE module_domain
    1125. 

  1125. !---------------------------------------------------------------------
    1126. 

  1126.    IMPLICIT NONE
    1127. 

  1127. !---------------------------------------------------------------------
    1128. 

  1128.    INTEGER,  INTENT(IN)           :: id
    1129. 

  1129.    TYPE (grid_config_rec_type)    :: config_flags
    1130. 

  1130.    LOGICAL , INTENT(IN)           :: restart
    1131. 

  1131.    LOGICAL,  INTENT(IN)           :: allowed_to_read
    1132. 

  1132. 

  1133.    INTEGER , INTENT(IN)           :: ids, ide, jds, jde, kds, kde,   &
    1134. 

  1134.                                      ims, ime, jms, jme, kms, kme,   &
    1135. 

  1135.                                      its, ite, jts, jte, kts, kte
    1136. 

  1136. 

  1137.    INTEGER , INTENT(IN)           :: JULDAY,JULYR
    1138. 

  1138.    REAL ,    INTENT(IN)           :: DT, RADT, cen_lat, GMT, pptop,  &
    1139. 

  1139.                                      swrad_scat, p_top
    1140. 

  1140.    LOGICAL,  INTENT(IN)           :: start_of_simulation
    1141. 

  1141. 

  1142.    INTEGER,      INTENT(IN   )    ::   levsiz, n_ozmixm
    1143. 

  1143.    INTEGER,      INTENT(IN   )    ::   paerlev, n_aerosolc
    1144. 

  1144. 

  1145.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(IN) ::  XLAT
    1146. 

  1146. 

  1147.    REAL,  DIMENSION( ims:ime, levsiz, jms:jme, n_ozmixm ), OPTIONAL,      &
    1148. 

  1148.           INTENT(INOUT) ::                                  OZMIXM
    1149. 

  1149. 

  1150.    REAL,  DIMENSION(ims:ime,jms:jme), OPTIONAL, INTENT(INOUT)  :: m_ps_1,m_ps_2
    1151. 

  1151.    REAL,  DIMENSION(paerlev), OPTIONAL, INTENT(INOUT)  ::         m_hybi
    1152. 

  1152.    REAL,  DIMENSION( ims:ime, paerlev, jms:jme, n_aerosolc ), OPTIONAL,     &
    1153. 

  1153.           INTENT(INOUT) ::                      aerosolc_1, aerosolc_2
    1154. 

  1154. 

  1155.    REAL,  DIMENSION(levsiz), OPTIONAL, INTENT(INOUT)  ::          PIN
    1156. 

  1156. 

  1157.    INTEGER , INTENT(INOUT)        :: STEPRA
    1158. 

  1158.    INTEGER :: isn
    1159. 

  1159. 

  1160.    REAL , DIMENSION( kms:kme ) , INTENT(IN) :: sfull, shalf
    1161. 

  1161.    REAL , DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) ::           &
    1162. 

  1162.                                                            RTHRATEN, &
    1163. 

  1163.                                                          RTHRATENLW, &
    1164. 

  1164.                                                          RTHRATENSW, &
    1165. 

  1165.                                                              CLDFRA
    1166. 

  1166. 

  1167.    REAL , DIMENSION( ims:ime , kms:kme , jms:jme ) , OPTIONAL, INTENT(OUT) :: &
    1168. 

  1168.                                                          CLDFRA_OLD
    1169. 

  1169. 

  1170.    REAL , DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT) :: EMISS
    1171. 

  1171.    LOGICAL :: etalw = .false.
    1172. 

  1172.    LOGICAL :: hwrflw= .false.
    1173. 

  1173.    LOGICAL :: camlw = .false.
    1174. 

  1174. !   LOGICAL :: etamp = .false.
    1175. 

  1175.    LOGICAL :: acswalloc = .false.
    1176. 

  1176.    LOGICAL :: aclwalloc = .false.
    1177. 

  1177.    integer :: month,iday
    1178. 

  1178.    INTEGER :: i, j, k, itf, jtf, ktf
    1179. 

  1179. !---------------------------------------------------------------------
    1180. 

  1180. 

  1181.    jtf=min0(jte,jde-1)
    1182. 

  1182.    ktf=min0(kte,kde-1)
    1183. 

  1183.    itf=min0(ite,ide-1)
    1184. 

  1184. 

  1185. !---------------------------------------------------------------------
    1186. 

  1186. 

  1187. !-- calculate radiation time step
    1188. 

  1188. 

  1189.     STEPRA = nint(RADT*60./DT)
    1190. 

  1190.     STEPRA = max(STEPRA,1)
    1191. 

  1191. 

  1192. !-- initialization
    1193. 

  1193. 

  1194.    IF(start_of_simulation)THEN
    1195. 

  1195.      DO j=jts,jtf
    1196. 

  1196.      DO k=kts,ktf
    1197. 

  1197.      DO i=its,itf
    1198. 

  1198.         RTHRATEN(i,k,j)=0.
    1199. 

  1199.         RTHRATENLW(i,k,j)=0.
    1200. 

  1200.         RTHRATENSW(i,k,j)=0.
    1201. 

  1201.         CLDFRA(i,k,j)=0.
    1202. 

  1202.      ENDDO
    1203. 

  1203.      ENDDO
    1204. 

  1204.      ENDDO
    1205. 

  1205. 

  1206.      if( present(cldfra_old) ) then
    1207. 

  1207.         DO j=jts,jtf
    1208. 

  1208.         DO k=kts,ktf
    1209. 

  1209.         DO i=its,itf
    1210. 

  1210.            cldfra_old(i,k,j) = 0.
    1211. 

  1211.         ENDDO
    1212. 

  1212.         ENDDO
    1213. 

  1213.         ENDDO
    1214. 

  1214.      end if
    1215. 

  1215.    ENDIF
    1216. 

  1216. 

  1217. !-- find out which microphysics option is used first
    1218. 

  1218. 

  1219. !   mp_select: SELECT CASE(config_flags%mp_physics)
    1220. 

  1220. !
    1221. 

  1221. !        CASE (ETAMPNEW)
    1222. 

  1222. !             etamp = .true.
    1223. 

  1223. !
    1224. 

  1224. !   END SELECT mp_select
    1225. 

  1225. 

  1226. !-- chose long wave radiation scheme
    1227. 

  1227. 

  1228.    lwrad_select: SELECT CASE(config_flags%ra_lw_physics)
    1229. 

  1229. 

  1230.         CASE (RRTMSCHEME)
    1231. 

  1231.              CALL rrtminit(                                 &
    1232. 

  1232.                            p_top, allowed_to_read ,         &
    1233. 

  1233.                            ids, ide, jds, jde, kds, kde,    &
    1234. 

  1234.                            ims, ime, jms, jme, kms, kme,    &
    1235. 

  1235.                            its, ite, jts, jte, kts, kte     )
    1236. 

  1236. 

  1237.         CASE (CAMLWSCHEME)
    1238. 

  1238. #ifdef MAC_KLUDGE
    1239. 

  1239.              CALL wrf_error_fatal ( 'CAM radiation scheme not supported under the chosen build configuration' )
    1240. 

  1240. #endif
    1241. 

  1241.              IF ( PRESENT( OZMIXM ) .AND. PRESENT( PIN ) .AND. &
    1242. 

  1242.                   PRESENT(M_PS_1) .AND. PRESENT(M_PS_2) .AND.  &
    1243. 

  1243.                   PRESENT(M_HYBI) .AND. PRESENT(AEROSOLC_1)    &
    1244. 

  1244.                   .AND. PRESENT(AEROSOLC_2)) THEN
    1245. 

  1245.              CALL camradinit(                                  &
    1246. 

  1246.                          R_D,R_V,CP,G,STBOLT,EP_2,shalf,pptop, &
    1247. 

  1247.                          ozmixm,pin,levsiz,XLAT,n_ozmixm,      &
    1248. 

  1248.                          m_ps_1,m_ps_2,m_hybi,aerosolc_1,aerosolc_2,&
    1249. 

  1249.                          paerlev, n_aerosolc,              &
    1250. 

  1250.                          ids, ide, jds, jde, kds, kde,     &
    1251. 

  1251.                          ims, ime, jms, jme, kms, kme,     &
    1252. 

  1252.                          its, ite, jts, jte, kts, kte      )
    1253. 

  1253.              ELSE
    1254. 

  1254.                 CALL wrf_error_fatal ( 'arguments not present for calling cam radiation' )
    1255. 

  1255.              ENDIF
    1256. 

  1256. 

  1257.              camlw = .true.
    1258. 

  1258.              aclwalloc = .true.
    1259. 

  1259. 

  1260.         CASE (RRTMG_LWSCHEME)
    1261. 

  1261.              CALL rrtmg_lwinit(                             &
    1262. 

  1262.                            p_top, allowed_to_read ,         &
    1263. 

  1263.                            ids, ide, jds, jde, kds, kde,    &
    1264. 

  1264.                            ims, ime, jms, jme, kms, kme,    &
    1265. 

  1265.                            its, ite, jts, jte, kts, kte     )
    1266. 

  1266. 

  1267.              aclwalloc = .true.
    1268. 

  1268. 

  1269.         CASE (GFDLLWSCHEME)
    1270. 

  1270.              CALL nl_get_start_month(id,month)
    1271. 

  1271.              CALL nl_get_start_day(id,iday)
    1272. 

  1272.              CALL gfdletainit(emiss,sfull,shalf,pptop,      &
    1273. 

  1273.                               julyr,month,iday,gmt,         &
    1274. 

  1274.                               config_flags,allowed_to_read, &
    1275. 

  1275.                               ids, ide, jds, jde, kds, kde, &
    1276. 

  1276.                               ims, ime, jms, jme, kms, kme, &
    1277. 

  1277.                               its, ite, jts, jte, kts, kte  )
    1278. 

  1278.              etalw = .true.
    1279. 

  1279. #if(NMM_CORE==1)
    1280. 

  1280.         CASE (HWRFLWSCHEME)
    1281. 

  1281.              CALL nl_get_start_month(id,month)
    1282. 

  1282.              CALL nl_get_start_day(id,iday)
    1283. 

  1283. ! test this with standard jul-day calls
    1284. 

  1284. !     CALL nl_get_start_year(id,start_year)
    1285. 

  1285. !     CALL nl_get_start_month(id,start_month)
    1286. 

  1286. !     CALL nl_get_start_day(id,start_day)
    1287. 

  1287. !     CALL nl_get_start_hour(id,start_hour)
    1288. 

  1288. !     CALL nl_get_start_minute(id,start_minute)
    1289. 

  1289. !     CALL nl_get_start_second(id,start_second)
    1290. 

  1290. !     CALL jdn_sec(day_in_sec,start_year,start_month,start_day,0,0,0)
    1291. 

  1291. !     CALL jdn_sec(day_in_sec_ref,start_year,1,1,0,0,0)
    1292. 

  1292. !     julyr_start=start_year
    1293. 

  1293. !     julday_start=(day_in_sec-day_in_sec_ref)/(3600.*24.)+1
    1294. 

  1294. !     gmt_start=start_hour+real(start_minute)/60.+real(start_second)/3600.
    1295. 

  1295.               CALL hwrfrainit(sfull,shalf,pptop,JULYR,MONTH,IDAY,GMT,&
    1296. 

  1296. !             CALL hwrfrainit(sfull,shalf,pptop,JULYR_start,MONTH,IDAY,GMT_start,&
    1297. 

  1297.                             config_flags,allowed_to_read ,                         &
    1298. 

  1298.                            kds, kde, kms, kme, kts, kte     )
    1299. 

  1299.              hwrflw = .true.
    1300. 

  1300. #endif
    1301. 

  1301.         CASE (HELDSUAREZ)
    1302. 

  1302.              CALL hsinit(RTHRATEN,restart,             &
    1303. 

  1303.                          ids, ide, jds, jde, kds, kde, &
    1304. 

  1304.                          ims, ime, jms, jme, kms, kme, &
    1305. 

  1305.                          its, ite, jts, jte, kts, kte )
    1306. 

  1306. 

  1307.         CASE (FLGLWSCHEME)
    1308. 

  1308. 

  1309.         CASE DEFAULT
    1310. 

  1310. 

  1311.    END SELECT lwrad_select
    1312. 

  1312. !-- initialize short wave radiation scheme
    1313. 

  1313. 

  1314.    swrad_select: SELECT CASE(config_flags%ra_sw_physics)
    1315. 

  1315. 

  1316.         CASE (SWRADSCHEME)
    1317. 

  1317.              CALL swinit(                                  &
    1318. 

  1318.                          swrad_scat,                       &
    1319. 

  1319.                          allowed_to_read ,                 &
    1320. 

  1320.                          ids, ide, jds, jde, kds, kde,     &
    1321. 

  1321.                          ims, ime, jms, jme, kms, kme,     &
    1322. 

  1322.                          its, ite, jts, jte, kts, kte      )
    1323. 

  1323. 

  1324.         CASE (CAMSWSCHEME)
    1325. 

  1325. #ifdef MAC_KLUDGE
    1326. 

  1326.              CALL wrf_error_fatal ( 'CAM radiation scheme not supported under the chosen build configuration' )
    1327. 

  1327. #endif
    1328. 

  1328.              IF(.not.camlw)THEN
    1329. 

  1329.              CALL camradinit(                              &
    1330. 

  1330.                          R_D,R_V,CP,G,STBOLT,EP_2,shalf,pptop,               &
    1331. 

  1331.                          ozmixm,pin,levsiz,XLAT,n_ozmixm,     &
    1332. 

  1332.                          m_ps_1,m_ps_2,m_hybi,aerosolc_1,aerosolc_2,&
    1333. 

  1333.                          paerlev, n_aerosolc,              &
    1334. 

  1334.                          ids, ide, jds, jde, kds, kde,     &
    1335. 

  1335.                          ims, ime, jms, jme, kms, kme,     &
    1336. 

  1336.                          its, ite, jts, jte, kts, kte      )
    1337. 

  1337.              ENDIF
    1338. 

  1338.              acswalloc = .true.
    1339. 

  1339. 

  1340.         CASE (GSFCSWSCHEME)
    1341. 

  1341.              CALL gsfc_swinit(cen_lat, allowed_to_read )
    1342. 

  1342. 

  1343.         CASE (RRTMG_SWSCHEME)
    1344. 

  1344.              CALL rrtmg_swinit(                             &
    1345. 

  1345.                            allowed_to_read ,                &
    1346. 

  1346.                            ids, ide, jds, jde, kds, kde,    &
    1347. 

  1347.                            ims, ime, jms, jme, kms, kme,    &
    1348. 

  1348.                            its, ite, jts, jte, kts, kte     )
    1349. 

  1349. 

  1350.              acswalloc = .true.
    1351. 

  1351. 

  1352.         CASE (GFDLSWSCHEME)
    1353. 

  1353.              IF(.not.etalw)THEN
    1354. 

  1354.              CALL nl_get_start_month(id,month)
    1355. 

  1355.              CALL nl_get_start_day(id,iday)
    1356. 

  1356.              CALL gfdletainit(emiss,sfull,shalf,pptop,      &
    1357. 

  1357.                               julyr,month,iday,gmt,         &
    1358. 

  1358.                               config_flags,allowed_to_read, &
    1359. 

  1359.                               ids, ide, jds, jde, kds, kde, &
    1360. 

  1360.                               ims, ime, jms, jme, kms, kme, &
    1361. 

  1361.                               its, ite, jts, jte, kts, kte  )
    1362. 

  1362.              ENDIF
    1363. 

  1363. #if(NMM_CORE==1)
    1364. 

  1364.         CASE (HWRFSWSCHEME)
    1365. 

  1365.              IF(.not.hwrflw)THEN
    1366. 

  1366.              CALL nl_get_start_month(id,month)
    1367. 

  1367.              CALL nl_get_start_day(id,iday)
    1368. 

  1368.              CALL hwrfrainit(sfull,shalf,pptop,JULYR,MONTH,IDAY,GMT,&
    1369. 

  1369.                             config_flags,allowed_to_read,                          &
    1370. 

  1370.                            kds, kde, kms, kme, kts, kte     )
    1371. 

  1371.              ENDIF
    1372. 

  1372. #endif
    1373. 

  1373.         CASE (FLGSWSCHEME)
    1374. 

  1374. 

  1375.         CASE DEFAULT
    1376. 

  1376. 

  1377.    END SELECT swrad_select
    1378. 

  1378. 

  1379. #if ( EM_CORE == 1 ) 
    1380. 

  1380.    ! test for conditionally allocated arrays when using bucket_J
    1381. 

  1381. 

  1382.    IF(config_flags%bucket_J .gt. 0.0)THEN
    1383. 

  1383.      IF(.not. (acswalloc .and. aclwalloc))THEN
    1384. 

  1384.            CALL wrf_error_fatal ( 'Need CAM or RRTMG radiation for bucket_J option')
    1385. 

  1385.      ENDIF
    1386. 

  1386.    ENDIF
    1387. 

  1387. #endif
    1388. 

  1388. 

  1389.    END SUBROUTINE ra_init
    1390. 

  1390. 

  1391.    SUBROUTINE bl_init(STEPBL,BLDT,DT,RUBLTEN,RVBLTEN,RTHBLTEN,  &
    1392. 

  1392.                 RQVBLTEN,RQCBLTEN,RQIBLTEN,TSK,TMN,             &
    1393. 

  1393.                 config_flags,restart,UST,LOWLYR,TSLB,ZS,DZS,    &
    1394. 

  1394.                 num_soil_layers,TKE_PBL,mfshconv,               &
    1395. 

  1395.                 massflux_EDKF, entr_EDKF, detr_EDKF, & 
    1396. 

  1396.                 thl_up, thv_up, rt_up,       &
    1397. 

  1397.                 rv_up, rc_up, u_up, v_up,    &
    1398. 

  1398.                 frac_up, &
    1399. 

  1399.                 EXCH_H,VEGFRA,                                  &
    1400. 

  1400.                 SNOW,SNOWC, CANWAT,SMSTAV,                      &
    1401. 

  1401.                 SMSTOT, SFCRUNOFF,UDRUNOFF,ACSNOW,ACSNOM,       &
    1402. 

  1402.                 IVGTYP,ISLTYP,ISURBAN,SMOIS,SMFR3D,mavail,      &
    1403. 

  1403.                 SNOWH,SH2O,SNOALB,FNDSOILW,FNDSNOWH,RDMAXALB,   &
    1404. 

  1404. #if  ( NMM_CORE == 1 )
    1405. 

  1405.                 Z0,XLAND,XICE,                                  &
    1406. 

  1406. #else
    1407. 

  1407.                 ZNT,XLAND,XICE,                                 &
    1408. 

  1408. #endif
    1409. 

  1409.                 SFCEVP,GRDFLX,                                  &
    1410. 

  1410.                 MMINLU,                                         &
    1411. 

  1411.                 ISNOWXY, ZSNSOXY, TSNOXY,                       &
    1412. 

  1412.                 SNICEXY, SNLIQXY, TVXY, TGXY, CANICEXY,         &
    1413. 

  1413.                 CANLIQXY, EAHXY, TAHXY, CMXY,                   &
    1414. 

  1414.                 CHXY, FWETXY, SNEQVOXY, ALBOLDXY, QSNOWXY,      &
    1415. 

  1415.                 WSLAKEXY, ZWTXY, WAXY, WTXY, LFMASSXY, RTMASSXY,&
    1416. 

  1416.                 STMASSXY, WOODXY, STBLCPXY, FASTCPXY,           &
    1417. 

  1417.                 XSAIXY,                                         &
    1418. 

  1418.                 T2MVXY, T2MBXY ,CHSTARXY,                       &
    1419. 

  1419.                 allowed_to_read,                                &
    1420. 

  1420.                 start_of_simulation,                            &
    1421. 

  1421.                 te_temf,cf3d_temf,wm_temf,                      & ! WA
    1422. 

  1422. !                num_roof_layers,num_wall_layers,num_road_layers,& !Optional urban
    1423. 

  1423.                 DZR, DZB, DZG,                                  & !Optional urban
    1424. 

  1424.                 TR_URB2D,TB_URB2D,TG_URB2D,TC_URB2D,QC_URB2D,   & !Optional urban
    1425. 

  1425.                 XXXR_URB2D,XXXB_URB2D,XXXG_URB2D,XXXC_URB2D,    & !Optional urban
    1426. 

  1426.                 TRL_URB3D, TBL_URB3D, TGL_URB3D,                & !Optional urban
    1427. 

  1427.                 SH_URB2D,LH_URB2D,G_URB2D,RN_URB2D,             & !Optional urban
    1428. 

  1428.                 TS_URB2D, FRC_URB2D, UTYPE_URB2D,               &
    1429. 

  1429.                 SF_URBAN_PHYSICS,                               & !Optional urban
    1430. 

  1430.                 NUM_URBAN_LAYERS,                               & !Optional multi-layer urban
    1431. 

  1431.                 TRB_URB4D,TW1_URB4D,TW2_URB4D,                  & !Optional multi-layer urban
    1432. 

  1432.                 TGB_URB4D,TLEV_URB3D,QLEV_URB3D,                & !Optional multi-layer urban
    1433. 

  1433.                 TW1LEV_URB3D,TW2LEV_URB3D,                      & !Optional multi-layer urban
    1434. 

  1434.                 TGLEV_URB3D,TFLEV_URB3D,                        & !Optional multi-layer urban
    1435. 

  1435.                 SF_AC_URB3D,LF_AC_URB3D,CM_AC_URB3D,            & !Optional multi-layer urban
    1436. 

  1436.                 SFVENT_URB3D,LFVENT_URB3D,                      & !Optional multi-layer urban
    1437. 

  1437.                 SFWIN1_URB3D,SFWIN2_URB3D,                      & !Optional multi-layer urban
    1438. 

  1438.                 SFW1_URB3D,SFW2_URB3D,                          & !Optional multi-layer urban
    1439. 

  1439.                 SFR_URB3D,SFG_URB3D,                            & !Optional multi-layer urban
    1440. 

  1440.                 A_U_BEP,A_V_BEP,A_T_BEP,A_Q_BEP,                & !Optional multi-layer urban
    1441. 

  1441.                 A_E_BEP,B_U_BEP,B_V_BEP,                        & !Optional multi-layer urban
    1442. 

  1442.                 B_T_BEP,B_Q_BEP,B_E_BEP,DLG_BEP,                & !Optional multi-layer urban
    1443. 

  1443.                 DL_U_BEP,SF_BEP,VL_BEP,                         & !Optional multi-layer urban
    1444. 

  1444.                 ids, ide, jds, jde, kds, kde,                   &
    1445. 

  1445.                 ims, ime, jms, jme, kms, kme,                   &
    1446. 

  1446.                 its, ite, jts, jte, kts, kte,                   &
    1447. 

  1447.                 ACHFX,ACLHF,ACGRDFLX,                           &
    1448. 

  1448.                 oml_hml0, omlcall,                              & !Optional oml
    1449. 

  1449.                 TML,T0ML,HML,H0ML,HUML,HVML,TMOML               ) !Optional oml
    1450. 

  1450. !--------------------------------------------------------------------
    1451. 

  1451.    USE module_sf_sfclay
    1452. 

  1452.    USE module_sf_sfclayrev
    1453. 

  1453.    USE module_sf_slab
    1454. 

  1454.    USE module_sf_pxsfclay
    1455. 

  1455.    USE module_bl_ysu
    1456. 

  1456.    USE module_bl_mrf
    1457. 

  1457.    USE module_bl_gfs
    1458. 

  1458.    USE module_bl_gfs2011, only : gfs2011init
    1459. 

  1459.    USE module_bl_acm
    1460. 

  1460.    USE module_sf_myjsfc
    1461. 

  1461.    USE module_sf_qnsesfc
    1462. 

  1462.    USE module_sf_noahdrv
    1463. 

  1463.    USE module_sf_noahmpdrv
    1464. 

  1464.    USE module_sf_urban
    1465. 

  1465.    USE module_sf_bep                                  !BEP
    1466. 

  1466.    USE module_sf_bep_bem
    1467. 

  1467.    USE module_sf_ruclsm
    1468. 

  1468.    USE module_sf_pxlsm
    1469. 

  1469.    USE module_sf_oml
    1470. 

  1470.    USE module_bl_myjpbl
    1471. 

  1471.    USE module_bl_myjurb
    1472. 

  1472.    USE module_bl_boulac
    1473. 

  1473.    USE module_bl_camuwpbl_driver, ONLY :  camuwpblinit
    1474. 

  1474.    USE module_bl_qnsepbl
    1475. 

  1475.    USE module_bl_qnsepbl09
    1476. 

  1476.    USE module_bl_mfshconvpbl
    1477. 

  1477. #if ( EM_CORE == 1 )
    1478. 

  1478.    USE module_bl_mynn
    1479. 

  1479.    USE module_bl_temf
    1480. 

  1480.    USE module_sf_temfsfclay
    1481. 

  1481.    USE module_sf_mynn
    1482. 

  1482. #endif
    1483. 

  1483. 

  1484. #if (NMM_CORE == 1)
    1485. 

  1485.    USE module_sf_gfdl
    1486. 

  1486. #endif
    1487. 

  1487. !--------------------------------------------------------------------
    1488. 

  1488.    IMPLICIT NONE
    1489. 

  1489. !--------------------------------------------------------------------
    1490. 

  1490.    TYPE (grid_config_rec_type) ::     config_flags
    1491. 

  1491.    LOGICAL , INTENT(IN)        :: restart
    1492. 

  1492.    LOGICAL, INTENT(IN)         ::   FNDSOILW, FNDSNOWH
    1493. 

  1493.    LOGICAL, INTENT(IN)         ::   RDMAXALB
    1494. 

  1494. 

  1495.    INTEGER , INTENT(IN)        ::     ids, ide, jds, jde, kds, kde, &
    1496. 

  1496.                                       ims, ime, jms, jme, kms, kme, &
    1497. 

  1497.                                       its, ite, jts, jte, kts, kte
    1498. 

  1498.    INTEGER , INTENT(IN)        ::     num_soil_layers
    1499. 

  1499.    INTEGER , INTENT(IN)        ::     SF_URBAN_PHYSICS 
    1500. 

  1500. 

  1501.    REAL ,    INTENT(IN)        ::     DT, BLDT
    1502. 

  1502.    INTEGER , INTENT(INOUT)     ::     STEPBL
    1503. 

  1503. 

  1504.    REAL,     DIMENSION( ims:ime , 1:num_soil_layers , jms:jme ),    &
    1505. 

  1505.              INTENT(OUT) :: SMFR3D
    1506. 

  1506. 

  1507.    REAL,     DIMENSION( ims:ime , 1:num_soil_layers , jms:jme ),&
    1508. 

  1508.                    INTENT(INOUT) :: SMOIS,SH2O,TSLB
    1509. 

  1509. 

  1510.    REAL,    DIMENSION( ims:ime, jms:jme )                     , &
    1511. 

  1511.             INTENT(INOUT)    ::                           SNOW, &
    1512. 

  1512.                                                          SNOWH, &
    1513. 

  1513.                                                          SNOWC, &
    1514. 

  1514.                                                         SNOALB, &
    1515. 

  1515.                                                         CANWAT, &
    1516. 

  1516.                                                         MAVAIL, &
    1517. 

  1517.                                                         SMSTAV, &
    1518. 

  1518.                                                         SMSTOT, &
    1519. 

  1519.                                                      SFCRUNOFF, &
    1520. 

  1520.                                                       UDRUNOFF, &
    1521. 

  1521.                                                         ACSNOW, &
    1522. 

  1522.                                                         VEGFRA, &
    1523. 

  1523.                                                         ACSNOM, &
    1524. 

  1524.                                                         SFCEVP, &
    1525. 

  1525.                                                         GRDFLX, &
    1526. 

  1526.                                                            UST, &
    1527. 

  1527. #if ( NMM_CORE == 1 )
    1528. 

  1528.                                                             Z0, &
    1529. 

  1529. #else
    1530. 

  1530.                                                            ZNT, &
    1531. 

  1531. #endif
    1532. 

  1532.                                                          XLAND, &
    1533. 

  1533.                                                          XICE
    1534. 

  1534. 

  1535.    INTEGER, DIMENSION( ims:ime, jms:jme )                     , &
    1536. 

  1536.             INTENT(INOUT)    ::                         IVGTYP, &
    1537. 

  1537.                                                         ISLTYP, &
    1538. 

  1538.                                                         LOWLYR
    1539. 

  1539. 

  1540. 

  1541.    REAL,     DIMENSION(1:num_soil_layers), INTENT(INOUT)  ::  ZS,DZS
    1542. 

  1542. 

  1543.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) ::       &
    1544. 

  1544.                                                            RUBLTEN, &
    1545. 

  1545.                                                            RVBLTEN, &
    1546. 

  1546.                                                           EXCH_H,   &
    1547. 

  1547.                                                           RTHBLTEN, &
    1548. 

  1548.                                                           RQVBLTEN, &
    1549. 

  1549.                                                           RQCBLTEN, &
    1550. 

  1550.                                                           RQIBLTEN, &
    1551. 

  1551.                                                           TKE_PBL
    1552. 

  1552.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT), OPTIONAL :: &
    1553. 

  1553.                               massflux_EDKF, entr_EDKF, detr_EDKF & 
    1554. 

  1554.                                      ,thl_up, thv_up, rt_up       &
    1555. 

  1555.                                      ,rv_up, rc_up, u_up, v_up    &
    1556. 

  1556.                                      ,frac_up
    1557. 

  1557. 

  1558.    INTEGER, INTENT(IN)           ::  mfshconv ! WRF JP
    1559. 

  1559. 

  1560.    REAL,  DIMENSION( ims:ime , jms:jme ) , INTENT(IN) ::     TSK
    1561. 

  1561.    REAL,  DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT) ::  TMN
    1562. 

  1562.    CHARACTER(LEN=*), INTENT(IN)   :: MMINLU
    1563. 

  1563.    LOGICAL,  INTENT(IN)           :: allowed_to_read
    1564. 

  1564.    INTEGER,  INTENT(IN)           :: ISURBAN
    1565. 

  1565.    INTEGER :: isn, isfc
    1566. 

  1566.    INTEGER :: k
    1567. 

  1567. 

  1568.    REAL, OPTIONAL, DIMENSION( ims:ime , kms:kme , jms:jme ) , &
    1569. 

  1569.             INTENT(OUT) :: te_temf, cf3d_temf !WA 
    1570. 

  1570.    REAL, OPTIONAL, DIMENSION( ims:ime , jms:jme ) , &
    1571. 

  1571.             INTENT(OUT) :: wm_temf
    1572. 

  1572. 

  1573. !Noah-MP
    1574. 

  1574. 

  1575.    INTEGER, OPTIONAL, DIMENSION(ims:ime,jms:jme) :: ISNOWXY
    1576. 

  1576.    REAL,    OPTIONAL, DIMENSION(ims:ime,-2:num_soil_layers, jms:jme) :: ZSNSOXY
    1577. 

  1577.    REAL,    OPTIONAL, DIMENSION(ims:ime,-2:0, jms:jme) :: TSNOXY
    1578. 

  1578.    REAL,    OPTIONAL, DIMENSION(ims:ime,-2:0, jms:jme) :: SNICEXY
    1579. 

  1579.    REAL,    OPTIONAL, DIMENSION(ims:ime,-2:0, jms:jme) :: SNLIQXY
    1580. 

  1580.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: TVXY
    1581. 

  1581.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: TGXY
    1582. 

  1582.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CANICEXY
    1583. 

  1583.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CANLIQXY
    1584. 

  1584.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: EAHXY
    1585. 

  1585.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: TAHXY
    1586. 

  1586.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CMXY
    1587. 

  1587.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CHXY
    1588. 

  1588.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: FWETXY
    1589. 

  1589.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: SNEQVOXY
    1590. 

  1590.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: ALBOLDXY
    1591. 

  1591.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: QSNOWXY
    1592. 

  1592.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: WSLAKEXY
    1593. 

  1593.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: ZWTXY
    1594. 

  1594.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: WAXY
    1595. 

  1595.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: WTXY
    1596. 

  1596.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: LFMASSXY
    1597. 

  1597.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: RTMASSXY
    1598. 

  1598.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: STMASSXY
    1599. 

  1599.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: WOODXY
    1600. 

  1600.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: STBLCPXY
    1601. 

  1601.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: FASTCPXY
    1602. 

  1602.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: XSAIXY
    1603. 

  1603.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: T2MVXY
    1604. 

  1604.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: T2MBXY
    1605. 

  1605.    REAL,    OPTIONAL, DIMENSION(ims:ime,jms:jme) :: CHSTARXY
    1606. 

  1606. 

  1607. !URBAN
    1608. 

  1608. !   REAL, DIMENSION(1:num_roof_layers), INTENT(INOUT) :: DZR  !Optional urban
    1609. 

  1609. !   REAL, DIMENSION(1:num_wall_layers), INTENT(INOUT) :: DZB  !Optional urban
    1610. 

  1610. !   REAL, DIMENSION(1:num_road_layers), INTENT(INOUT) :: DZG  !Optional urban
    1611. 

  1611.     REAL, OPTIONAL, DIMENSION(1:num_soil_layers), INTENT(INOUT) :: DZR  !Optional urban
    1612. 

  1612.     REAL, OPTIONAL, DIMENSION(1:num_soil_layers), INTENT(INOUT) :: DZB  !Optional urban
    1613. 

  1613.     REAL, OPTIONAL, DIMENSION(1:num_soil_layers), INTENT(INOUT) :: DZG  !Optional urban
    1614. 

  1614.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TR_URB2D !Optional urban
    1615. 

  1615.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TB_URB2D !Optional urban
    1616. 

  1616.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TG_URB2D !Optional urban
    1617. 

  1617.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TC_URB2D !Optional urban
    1618. 

  1618.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: QC_URB2D !Optional urban
    1619. 

  1619.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: XXXR_URB2D !Optional urban
    1620. 

  1620.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: XXXB_URB2D !Optional urban
    1621. 

  1621.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: XXXG_URB2D !Optional urban
    1622. 

  1622.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: XXXC_URB2D !Optional urban
    1623. 

  1623.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: SH_URB2D !Optional urban
    1624. 

  1624.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: LH_URB2D !Optional urban
    1625. 

  1625.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: G_URB2D !Optional urban
    1626. 

  1626.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: RN_URB2D !Optional urban
    1627. 

  1627.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: TS_URB2D !Optional urban
    1628. 

  1628.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: FRC_URB2D !Optional urban
    1629. 

  1629.     INTEGER, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: UTYPE_URB2D !Optional urban
    1630. 

  1630. !    REAL, DIMENSION( ims:ime, 1:num_roof_layers, jms:jme ), INTENT(INOUT) :: TRL_URB3D !Optional urban
    1631. 

  1631. !    REAL, DIMENSION( ims:ime, 1:num_wall_layers, jms:jme ), INTENT(INOUT) :: TBL_URB3D !Optional urban
    1632. 

  1632. !    REAL, DIMENSION( ims:ime, 1:num_road_layers, jms:jme ), INTENT(INOUT) :: TGL_URB3D !Optional urban
    1633. 

  1633.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_soil_layers, jms:jme ), INTENT(INOUT) :: TRL_URB3D !Optional urban
    1634. 

  1634.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_soil_layers, jms:jme ), INTENT(INOUT) :: TBL_URB3D !Optional urban
    1635. 

  1635.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_soil_layers, jms:jme ), INTENT(INOUT) :: TGL_URB3D !Optional urban
    1636. 

  1636. 

  1637.     INTEGER , INTENT(IN)        ::     num_urban_layers
    1638. 

  1638.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TRB_URB4D !Optional UCM
    1639. 

  1639.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TW1_URB4D !Optional UCM
    1640. 

  1640.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TW2_URB4D !Optional UCM
    1641. 

  1641.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TGB_URB4D !Optional UCM
    1642. 

  1642.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TLEV_URB3D !Optional UCM
    1643. 

  1643.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: QLEV_URB3D !Optional UCM
    1644. 

  1644.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TW1LEV_URB3D ! multi-layer UCM
    1645. 

  1645.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TW2LEV_URB3D ! multi-layer UCM
    1646. 

  1646.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TGLEV_URB3D ! multi-layer UCM
    1647. 

  1647.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: TFLEV_URB3D ! multi-layer UCM
    1648. 

  1648.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: LF_AC_URB3D  !multi-layer UCM
    1649. 

  1649.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: SF_AC_URB3D  !multi-layer UCM
    1650. 

  1650.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: CM_AC_URB3D  !multi-layer UCM
    1651. 

  1651.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: SFVENT_URB3D !multi-layer UCM
    1652. 

  1652.     REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: LFVENT_URB3D !multi-layer UCM
    1653. 

  1653.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: SFWIN1_URB3D ! multi-layer UCM
    1654. 

  1654.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: SFWIN2_URB3D ! multi-layer UCM
    1655. 

  1655.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: SFW1_URB3D !Optional UCM
    1656. 

  1656.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: SFW2_URB3D !Optional UCM
    1657. 

  1657.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: SFR_URB3D !Optional UCM
    1658. 

  1658.     REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: SFG_URB3D !Optional UCM
    1659. 

  1659.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_U_BEP 
    1660. 

  1660.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_V_BEP 
    1661. 

  1661.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_T_BEP 
    1662. 

  1662.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_Q_BEP 
    1663. 

  1663.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: A_E_BEP 
    1664. 

  1664.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_U_BEP 
    1665. 

  1665.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_V_BEP 
    1666. 

  1666.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_T_BEP 
    1667. 

  1667.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_Q_BEP 
    1668. 

  1668.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: B_E_BEP 
    1669. 

  1669.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: VL_BEP
    1670. 

  1670.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: DLG_BEP 
    1671. 

  1671.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme,jms:jme),INTENT(INOUT) :: SF_BEP
    1672. 

  1672.    REAL, OPTIONAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) :: DL_U_BEP
    1673. 

  1673. 

  1674.    REAL,  DIMENSION( ims:ime , jms:jme ) , OPTIONAL, INTENT(INOUT) ::    &
    1675. 

  1675.                                         ACHFX,ACLHF,ACGRDFLX
    1676. 

  1676. ! Optional OML variables
    1677. 

  1677.    REAL,  DIMENSION( ims:ime , jms:jme ) , OPTIONAL, INTENT(INOUT) ::    &
    1678. 

  1678.                                         TML,T0ML,HML,H0ML,HUML,HVML,TMOML
    1679. 

  1679.    INTEGER,  OPTIONAL,  INTENT(IN) :: omlcall
    1680. 

  1680.    REAL,  OPTIONAL,  INTENT(IN) :: oml_hml0
    1681. 

  1681.    LOGICAL,  INTENT(IN) :: start_of_simulation
    1682. 

  1682.    INTEGER :: i,j
    1683. 

  1683. 

  1684. 

  1685. #if ( EM_CORE == 1 )
    1686. 

  1686. !local mynn 
    1687. 

  1687.    INTEGER :: mynn_closure_level
    1688. 

  1688. #endif
    1689. 

  1689. 

  1690. !-- calculate pbl time step
    1691. 

  1691. 

  1692.    STEPBL = nint(BLDT*60./DT)
    1693. 

  1693.    STEPBL = max(STEPBL,1)
    1694. 

  1694. !-- initialization
    1695. 

  1695. 

  1696.    IF(PRESENT(ACHFX))THEN
    1697. 

  1697.    IF(.not.restart)THEN
    1698. 

  1698.      DO j=jts,jte
    1699. 

  1699.      DO i=its,ite
    1700. 

  1700.         ACHFX(i,j)=0.
    1701. 

  1701.         ACLHF(i,j)=0.
    1702. 

  1702.         ACGRDFLX(i,j)=0.
    1703. 

  1703.         SFCEVP(i,j)=0.
    1704. 

  1704.      ENDDO
    1705. 

  1705.      ENDDO
    1706. 

  1706.    ENDIF
    1707. 

  1707.    ENDIF
    1708. 

  1708. 

  1709. !-- initialize surface layer scheme
    1710. 

  1710. 

  1711.    sfclay_select: SELECT CASE(config_flags%sf_sfclay_physics)
    1712. 

  1712. 

  1713.       CASE (SFCLAYSCHEME)
    1714. 

  1714.            CALL sfclayinit( allowed_to_read )
    1715. 

  1715.            isfc = 1
    1716. 

  1716.       CASE (SFCLAYREVSCHEME)
    1717. 

  1717. !           CALL sfclayinit( allowed_to_read )
    1718. 

  1718.            isfc = 1
    1719. 

  1719.       CASE (PXSFCSCHEME)
    1720. 

  1720.            CALL pxsfclayinit( allowed_to_read )
    1721. 

  1721.            isfc = 7
    1722. 

  1722.       CASE (MYJSFCSCHEME)
    1723. 

  1723.            CALL myjsfcinit(LOWLYR,UST,                         &
    1724. 

  1724. #if ( NMM_CORE == 1 )
    1725. 

  1725.                                       Z0,                      &
    1726. 

  1726. #else
    1727. 

  1727.                                       ZNT,                     &
    1728. 

  1728. #endif
    1729. 

  1729.                                           XLAND,XICE,          &
    1730. 

  1730.                          IVGTYP,restart,                       &
    1731. 

  1731.                          allowed_to_read ,                     &
    1732. 

  1732.                          ids, ide, jds, jde, kds, kde,         &
    1733. 

  1733.                          ims, ime, jms, jme, kms, kme,         &
    1734. 

  1734.                          its, ite, jts, jte, kts, kte          )
    1735. 

  1735.            isfc = 2
    1736. 

  1736. 

  1737.       CASE (QNSESFCSCHEME)
    1738. 

  1738.            CALL qnsesfcinit(LOWLYR,UST,                         &
    1739. 

  1739. #if ( NMM_CORE == 1 )
    1740. 

  1740.                                       Z0,                      &
    1741. 

  1741. #else
    1742. 

  1742.                                       ZNT,                     &
    1743. 

  1743. #endif
    1744. 

  1744.                                           XLAND,XICE,          &
    1745. 

  1745.                          IVGTYP,restart,                       &
    1746. 

  1746.                          allowed_to_read ,                     &
    1747. 

  1747.                          ids, ide, jds, jde, kds, kde,         &
    1748. 

  1748.                          ims, ime, jms, jme, kms, kme,         &
    1749. 

  1749.                          its, ite, jts, jte, kts, kte          )
    1750. 

  1750.            isfc = 4
    1751. 

  1751. 

  1752.       CASE (GFSSFCSCHEME)
    1753. 

  1753.            CALL myjsfcinit(LOWLYR,UST,                         &
    1754. 

  1754. #if ( NMM_CORE == 1 )
    1755. 

  1755.                                       Z0,                      &
    1756. 

  1756. #else
    1757. 

  1757.                                       ZNT,                     &
    1758. 

  1758. #endif
    1759. 

  1759.                                           XLAND,XICE,          &
    1760. 

  1760.                          IVGTYP,restart,                       &
    1761. 

  1761.                          allowed_to_read ,                     &
    1762. 

  1762.                          ids, ide, jds, jde, kds, kde,         &
    1763. 

  1763.                          ims, ime, jms, jme, kms, kme,         &
    1764. 

  1764.                          its, ite, jts, jte, kts, kte          )
    1765. 

  1765.            isfc = 2
    1766. 

  1766. #if (NMM_CORE==1)
    1767. 

  1767.       CASE (GFDLSFCSCHEME)
    1768. 

  1768.            CALL myjsfcinit(LOWLYR,UST,                         &
    1769. 

  1769.                                       Z0,                      &
    1770. 

  1770.                                           XLAND,XICE,          &
    1771. 

  1771.                          IVGTYP,restart,                       &
    1772. 

  1772.                          allowed_to_read ,                     &
    1773. 

  1773.                          ids, ide, jds, jde, kds, kde,         &
    1774. 

  1774.                          ims, ime, jms, jme, kms, kme,         &
    1775. 

  1775.                          its, ite, jts, jte, kts, kte          )
    1776. 

  1776.            isfc = 2
    1777. 

  1777. #endif
    1778. 

  1778. 

  1779. #if ( EM_CORE == 1 )
    1780. 

  1780. !mynn 
    1781. 

  1781. 

  1782.         CASE (MYNNSFCSCHEME)
    1783. 

  1783. 

  1784.            CALL mynn_sf_init_driver(allowed_to_read)
    1785. 

  1785.            isfc=5
    1786. 

  1786. !          isfc=3
    1787. 

  1787. 

  1788.       CASE (TEMFSFCSCHEME)
    1789. 

  1789.            CALL wrf_debug( 100, 'calling temfsfclayinit' )
    1790. 

  1790.            CALL temfsfclayinit( restart, allowed_to_read ,        &
    1791. 

  1791.                              wm_temf,                          &
    1792. 

  1792.                          ids, ide, jds, jde, kds, kde,         &
    1793. 

  1793.                          ims, ime, jms, jme, kms, kme,         &
    1794. 

  1794.                          its, ite, jts, jte, kts, kte          )
    1795. 

  1795. #endif
    1796. 

  1796. 

  1797.       CASE DEFAULT
    1798. 

  1798. 

  1799.    END SELECT sfclay_select
    1800. 

  1800. 

  1801. 

  1802. !-- initialize surface scheme
    1803. 

  1803. 

  1804.    sfc_select: SELECT CASE(config_flags%sf_surface_physics)
    1805. 

  1805. 

  1806.       CASE (SLABSCHEME)
    1807. 

  1807. 

  1808.            CALL slabinit(TSK,TMN,                              &
    1809. 

  1809.                          TSLB,ZS,DZS,num_soil_layers,          &
    1810. 

  1810.                          allowed_to_read ,start_of_simulation ,&
    1811. 

  1811.                          ids, ide, jds, jde, kds, kde,         &
    1812. 

  1812.                          ims, ime, jms, jme, kms, kme,         &
    1813. 

  1813.                          its, ite, jts, jte, kts, kte          )
    1814. 

  1814. 

  1815. #if (NMM_CORE == 1)
    1816. 

  1816.      CASE (GFDLSLAB)
    1817. 

  1817.            CALL hwrfsfcinit(isn,XICE,VEGFRA,SNOW,SNOWC, CANWAT,SMSTAV, &
    1818. 

  1818.                      SMSTOT, SFCRUNOFF,UDRUNOFF,GRDFLX,ACSNOW,     &
    1819. 

  1819.                      ACSNOM,IVGTYP,ISLTYP,TSLB,SMOIS,DZS,SFCEVP,   &
    1820. 

  1820.                      TMN,                                          &
    1821. 

  1821.                      num_soil_layers,                              &
    1822. 

  1822.                      allowed_to_read ,                             &
    1823. 

  1823.                      ids,ide, jds,jde, kds,kde,                    &
    1824. 

  1824.                      ims,ime, jms,jme, kms,kme,                    &
    1825. 

  1825.                      its,ite, jts,jte, kts,kte                     )
    1826. 

  1826. 

  1827. #endif
    1828. 

  1828.       CASE (LSMSCHEME)
    1829. 

  1829.           CALL LSMINIT(VEGFRA,SNOW,SNOWC,SNOWH,CANWAT,SMSTAV,  &
    1830. 

  1830.                      SMSTOT, SFCRUNOFF,UDRUNOFF,ACSNOW,        &
    1831. 

  1831.                      ACSNOM,IVGTYP,ISLTYP,TSLB,SMOIS,SH2O,ZS,DZS, &
    1832. 

  1832.                      MMINLU,                                   &
    1833. 

  1833.                      SNOALB, FNDSOILW, FNDSNOWH, RDMAXALB,     &
    1834. 

  1834.                      num_soil_layers, restart,                 &
    1835. 

  1835.                      allowed_to_read ,                         &
    1836. 

  1836.                      ids,ide, jds,jde, kds,kde,                &
    1837. 

  1837.                      ims,ime, jms,jme, kms,kme,                &
    1838. 

  1838.                      its,ite, jts,jte, kts,kte                 )
    1839. 

  1839. 

  1840. !URBAN
    1841. 

  1841.           IF ((SF_URBAN_PHYSICS.eq.1).OR.(SF_URBAN_PHYSICS.EQ.2).OR.(SF_URBAN_PHYSICS.EQ.3)) THEN
    1842. 

  1842. 

  1843.              IF ( PRESENT( FRC_URB2D ) .AND. PRESENT( UTYPE_URB2D )) THEN
    1844. 

  1844.                 
    1845. 

  1845.                 CALL urban_param_init(DZR,DZB,DZG,num_soil_layers,                   & !urban
    1846. 

  1846.                                 sf_urban_physics)
    1847. 

  1847. !                                num_roof_layers,num_wall_layers,road_soil_layers)   !urban
    1848. 

  1848.                                
    1849. 

  1849.                 
    1850. 

  1850.                 CALL urban_var_init(ISURBAN,TSK,TSLB,TMN,IVGTYP,                     & !urban
    1851. 

  1851.                               ims,ime,jms,jme,kms,kme,num_soil_layers,               & !urban
    1852. 

  1852. !                              num_roof_layers,num_wall_layers,num_road_layers, & !urban
    1853. 

  1853.                               restart,sf_urban_physics,                        & !urban
    1854. 

  1854.                               XXXR_URB2D,XXXB_URB2D,XXXG_URB2D,XXXC_URB2D,     & !urban
    1855. 

  1855.                               TR_URB2D,TB_URB2D,TG_URB2D,TC_URB2D,QC_URB2D,    & !urban
    1856. 

  1856.                               TRL_URB3D,TBL_URB3D,TGL_URB3D,                   & !urban
    1857. 

  1857.                               SH_URB2D,LH_URB2D,G_URB2D,RN_URB2D, TS_URB2D,    & !urban
    1858. 

  1858.                               num_urban_layers,                                & !urban
    1859. 

  1859.                               TRB_URB4D,TW1_URB4D,TW2_URB4D,TGB_URB4D,         & !urban
    1860. 

  1860.                               TLEV_URB3D,QLEV_URB3D,                           & !urban
    1861. 

  1861.                               TW1LEV_URB3D,TW2LEV_URB3D,                       & !urban
    1862. 

  1862.                               TGLEV_URB3D,TFLEV_URB3D,                         & !urban
    1863. 

  1863.                               SF_AC_URB3D,LF_AC_URB3D,CM_AC_URB3D,             & !urban
    1864. 

  1864.                               SFVENT_URB3D,LFVENT_URB3D,                       & !urban
    1865. 

  1865.                               SFWIN1_URB3D,SFWIN2_URB3D,                       & !urban
    1866. 

  1866.                               SFW1_URB3D,SFW2_URB3D,SFR_URB3D,SFG_URB3D,       & !urban
    1867. 

  1867.                               A_U_BEP,A_V_BEP,A_T_BEP,A_Q_BEP,                 & !multi-layer urban
    1868. 

  1868.                               A_E_BEP,B_U_BEP,B_V_BEP,                         & !multi-layer urban
    1869. 

  1869.                               B_T_BEP,B_Q_BEP,B_E_BEP,DLG_BEP,                 & !multi-layer urban
    1870. 

  1870.                               DL_U_BEP,SF_BEP,VL_BEP,                          & !multi-layer urban
    1871. 

  1871.                               FRC_URB2D, UTYPE_URB2D)                            !urban
    1872. 

  1872.              ELSE
    1873. 

  1873.                 CALL wrf_error_fatal ( 'arguments not present for calling urban model' )
    1874. 

  1874.              ENDIF
    1875. 

  1875.           ENDIF
    1876. 

  1876. !
    1877. 

  1877.       CASE (NOAHMPSCHEME)
    1878. 

  1878.           CALL NOAHMP_INIT(MMINLU, SNOW,SNOWH,CANWAT,ISLTYP,                                &
    1879. 

  1879.                      TSLB,SMOIS,SH2O,DZS, FNDSOILW, FNDSNOWH,                               &
    1880. 

  1880.                      TSK,isnowxy  ,tvxy     ,tgxy     ,canicexy ,                           &
    1881. 

  1881.                      canliqxy ,eahxy    ,tahxy    ,cmxy     ,chxy     ,                     &
    1882. 

  1882.                      fwetxy   ,sneqvoxy ,alboldxy ,qsnowxy  ,wslakexy ,zwtxy    ,waxy     , &
    1883. 

  1883.                      wtxy     ,tsnoxy   ,zsnsoxy  ,snicexy  ,snliqxy  ,lfmassxy ,rtmassxy , &
    1884. 

  1884.                      stmassxy ,woodxy   ,stblcpxy ,fastcpxy ,xsaixy   ,                     &
    1885. 

  1885.                      t2mvxy   ,t2mbxy   ,chstarxy ,            &
    1886. 

  1886.                      num_soil_layers, restart,                 &
    1887. 

  1887.                      allowed_to_read ,                         &
    1888. 

  1888.                      ids,ide, jds,jde, kds,kde,                &
    1889. 

  1889.                      ims,ime, jms,jme, kms,kme,                &
    1890. 

  1890.                      its,ite, jts,jte, kts,kte                 )
    1891. 

  1891. 

  1892.       CASE (RUCLSMSCHEME)
    1893. 

  1893. !          if(isfc .ne. 2)CALL wrf_error_fatal &
    1894. 

  1894. !           ( 'module_physics_init: use myjsfc and myjpbl scheme for this lsm option' )
    1895. 

  1895.            CALL ruclsminit( SH2O,SMFR3D,TSLB,SMOIS,ISLTYP,IVGTYP,MMINLU,XICE,  &
    1896. 

  1896.                      mavail,num_soil_layers, config_flags%iswater,      &
    1897. 

  1897. #if (NMM_CORE == 1)
    1898. 

  1898.                      config_flags%isice, z0, restart,                  &
    1899. 

  1899. #else
    1900. 

  1900.                      config_flags%isice, znt, restart,                  &
    1901. 

  1901. #endif
    1902. 

  1902.                      allowed_to_read ,                             &
    1903. 

  1903.                      ids,ide, jds,jde, kds,kde,                    &
    1904. 

  1904.                      ims,ime, jms,jme, kms,kme,                    &
    1905. 

  1905.                      its,ite, jts,jte, kts,kte                     )
    1906. 

  1906. 

  1907.       CASE (PXLSMSCHEME)
    1908. 

  1908.           if(config_flags%num_land_cat .ne. 24) CALL wrf_error_fatal &
    1909. 

  1909.             ( 'module_physics_init: 24 cat USGS must be used with PX LSM option' )
    1910. 

  1910.           CALL LSMINIT(VEGFRA,SNOW,SNOWC,SNOWH,CANWAT,SMSTAV,  &
    1911. 

  1911.                      SMSTOT, SFCRUNOFF,UDRUNOFF,ACSNOW,        &
    1912. 

  1912.                      ACSNOM,IVGTYP,ISLTYP,TSLB,SMOIS,SH2O,ZS,DZS, &
    1913. 

  1913.                      MMINLU,                                   &
    1914. 

  1914.                      SNOALB, FNDSOILW, FNDSNOWH, RDMAXALB,     &
    1915. 

  1915.                      num_soil_layers, restart,                 &
    1916. 

  1916.                      allowed_to_read ,                         &
    1917. 

  1917.                      ids,ide, jds,jde, kds,kde,                &
    1918. 

  1918.                      ims,ime, jms,jme, kms,kme,                &
    1919. 

  1919.                      its,ite, jts,jte, kts,kte                 )
    1920. 

  1920. 

  1921. !--------------fds (06/2010)-----------------------------------
    1922. 

  1922.      CASE  (SSIBSCHEME)
    1923. 

  1923. !SSiB only works with sfclay and YSU schemes. Check this here!
    1924. 

  1924.           if(isfc .ne. 1)CALL wrf_error_fatal &
    1925. 

  1925.            ( 'module_physics_init: use sfclay scheme with SSiB' )
    1926. 

  1926.           if(config_flags%bl_pbl_physics .ne. 1)CALL wrf_error_fatal &
    1927. 

  1927.            ( 'module_physics_init: use ysu scheme with SSiB' )
    1928. 

  1928. ! Add radiation scheme 4 (RRTMg) for SSiB, By Zhenxin 2011-06-20 ************************
    1929. 

  1929.           if(config_flags%ra_lw_physics .eq. 2 .or. config_flags%ra_lw_physics .gt. 4)CALL wrf_error_fatal &
    1930. 

  1930.            ( 'module_physics_init: SSiB only works with rrtm, cam scheme or rrtmg scheme (lw_phys=1,3,4)' )
    1931. 

  1931.           if(config_flags%ra_sw_physics .eq. 2 .or. config_flags%ra_sw_physics .gt. 4)CALL wrf_error_fatal &
    1932. 

  1932.            ( 'module_physics_init: SSiB only works with rrtm, cam scheme or rrtmg scheme (sw_phys=1,3,4)' )
    1933. 

  1933. ! End of Adding radiation scheme 4 (RRTMg) for SSiB, By Zhenxin 2011-06-20 **************
    1934. 

  1934. !--------------------------------------------------------------
    1935. 

  1935. 

  1936.       CASE DEFAULT
    1937. 

  1937. 

  1938.    END SELECT sfc_select
    1939. 

  1939. 

  1940.    IF(PRESENT(OMLCALL))THEN
    1941. 

  1941.      IF (omlcall .EQ. 1) THEN
    1942. 

  1942.         CALL omlinit(oml_hml0, tsk,                           &
    1943. 

  1943.                      tml,t0ml,hml,h0ml,huml,hvml,tmoml,       &
    1944. 

  1944.                      allowed_to_read, start_of_simulation,    &
    1945. 

  1945.                      ids,ide, jds,jde, kds,kde,               &
    1946. 

  1946.                      ims,ime, jms,jme, kms,kme,               &
    1947. 

  1947.                      its,ite, jts,jte, kts,kte                )
    1948. 

  1948.      ENDIF
    1949. 

  1949.    ENDIF
    1950. 

  1950. 

  1951. !-- initialize pbl scheme
    1952. 

  1952. 

  1953.    pbl_select: SELECT CASE(config_flags%bl_pbl_physics)
    1954. 

  1954. 

  1955.       CASE (YSUSCHEME)
    1956. 

  1956.            if(isfc .ne. 1)CALL wrf_error_fatal &
    1957. 

  1957.             ( 'module_physics_init: use sfclay scheme for this pbl option' )
    1958. 

  1958.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    1959. 

  1959.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    1960. 

  1960.            CALL ysuinit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN,    &
    1961. 

  1961.                         RQCBLTEN,RQIBLTEN,P_QI,               &
    1962. 

  1962.                         PARAM_FIRST_SCALAR,                   &
    1963. 

  1963.                         restart,                              &
    1964. 

  1964.                         allowed_to_read ,                     &
    1965. 

  1965.                         ids, ide, jds, jde, kds, kde,         &
    1966. 

  1966.                         ims, ime, jms, jme, kms, kme,         &
    1967. 

  1967.                         its, ite, jts, jte, kts, kte          )
    1968. 

  1968.       CASE (MRFSCHEME)
    1969. 

  1969.            if(isfc .ne. 1)CALL wrf_error_fatal &
    1970. 

  1970.             ( 'module_physics_init: use sfclay scheme for this pbl option' )
    1971. 

  1971.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    1972. 

  1972.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    1973. 

  1973.            CALL mrfinit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN,    &
    1974. 

  1974.                         RQCBLTEN,RQIBLTEN,P_QI,               &
    1975. 

  1975.                         PARAM_FIRST_SCALAR,                   &
    1976. 

  1976.                         restart,                              &
    1977. 

  1977.                         allowed_to_read ,                     &
    1978. 

  1978.                         ids, ide, jds, jde, kds, kde,         &
    1979. 

  1979.                         ims, ime, jms, jme, kms, kme,         &
    1980. 

  1980.                         its, ite, jts, jte, kts, kte          )
    1981. 

  1981.       CASE (ACMPBLSCHEME)
    1982. 

  1982.            if(isfc .ne. 1 .and. isfc .ne. 7)CALL wrf_error_fatal &
    1983. 

  1983.             ( 'module_physics_init: use sfclay or pxsfc scheme for this pbl option' )
    1984. 

  1984.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    1985. 

  1985.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    1986. 

  1986.            CALL acminit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN,    &
    1987. 

  1987.                         RQCBLTEN,RQIBLTEN,P_QI,               &
    1988. 

  1988.                         PARAM_FIRST_SCALAR,                   &
    1989. 

  1989.                         restart,                              &
    1990. 

  1990.                         allowed_to_read ,                     &
    1991. 

  1991.                         ids, ide, jds, jde, kds, kde,         &
    1992. 

  1992.                         ims, ime, jms, jme, kms, kme,         &
    1993. 

  1993.                         its, ite, jts, jte, kts, kte          )
    1994. 

  1994.       CASE (GFSSCHEME)
    1995. 

  1995.            if(isfc .ne. 2)CALL wrf_error_fatal &
    1996. 

  1996.             ( 'module_physics_init: use myjsfc scheme for this pbl option' )
    1997. 

  1997.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    1998. 

  1998.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    1999. 

  1999.            CALL gfsinit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN,    &
    2000. 

  2000.                         RQCBLTEN,RQIBLTEN,P_QI,               &
    2001. 

  2001.                         PARAM_FIRST_SCALAR,                   &
    2002. 

  2002.                         restart,                              &
    2003. 

  2003.                         allowed_to_read ,                     &
    2004. 

  2004.                         ids, ide, jds, jde, kds, kde,         &
    2005. 

  2005.                         ims, ime, jms, jme, kms, kme,         &
    2006. 

  2006.                         its, ite, jts, jte, kts, kte          )
    2007. 

  2007. #if (NMM_CORE == 1)
    2008. 

  2008.       CASE (GFS2011SCHEME)
    2009. 

  2009.            if(isfc .ne. 2)CALL wrf_error_fatal &
    2010. 

  2010.             ( 'module_physics_init: use myjsfc scheme for this pbl option' )
    2011. 

  2011.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    2012. 

  2012.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    2013. 

  2013.            CALL gfs2011init(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN,&
    2014. 

  2014.                         RQCBLTEN,RQIBLTEN,P_QI,               &
    2015. 

  2015.                         PARAM_FIRST_SCALAR,                   &
    2016. 

  2016.                         restart,                              &
    2017. 

  2017.                         allowed_to_read ,                     &
    2018. 

  2018.                         ids, ide, jds, jde, kds, kde,         &
    2019. 

  2019.                         ims, ime, jms, jme, kms, kme,         &
    2020. 

  2020.                         its, ite, jts, jte, kts, kte          )
    2021. 

  2021. #endif
    2022. 

  2022.       CASE (MYJPBLSCHEME)
    2023. 

  2023.            if(isfc .ne. 2)CALL wrf_error_fatal &
    2024. 

  2024.             ( 'module_physics_init: use myjsfc scheme for this pbl option' )
    2025. 

  2025.           IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) THEN
    2026. 

  2026.            CALL myjurbinit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN, &
    2027. 

  2027.                         TKE_PBL,EXCH_H,restart,               &
    2028. 

  2028.                         allowed_to_read ,                     &
    2029. 

  2029.                         ids, ide, jds, jde, kds, kde,         &
    2030. 

  2030.                         ims, ime, jms, jme, kms, kme,         &
    2031. 

  2031.                         its, ite, jts, jte, kts, kte          )
    2032. 

  2032.           ELSE
    2033. 

  2033. 

  2034.            CALL myjpblinit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN, &
    2035. 

  2035.                         TKE_PBL,EXCH_H,restart,               &
    2036. 

  2036.                         allowed_to_read ,                     &
    2037. 

  2037.                         ids, ide, jds, jde, kds, kde,         &
    2038. 

  2038.                         ims, ime, jms, jme, kms, kme,         &
    2039. 

  2039.                         its, ite, jts, jte, kts, kte          )
    2040. 

  2040.           END IF
    2041. 

  2041.       CASE (QNSEPBLSCHEME)
    2042. 

  2042.            if(isfc .ne. 4)CALL wrf_error_fatal &
    2043. 

  2043.             ( 'module_physics_init: use qnsesfc scheme for this pbl option' )
    2044. 

  2044.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    2045. 

  2045.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    2046. 

  2046.            CALL qnsepblinit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN, &
    2047. 

  2047.                         TKE_PBL,EXCH_H,restart,               &
    2048. 

  2048.                         allowed_to_read ,                     &
    2049. 

  2049.                         ids, ide, jds, jde, kds, kde,         &
    2050. 

  2050.                         ims, ime, jms, jme, kms, kme,         &
    2051. 

  2051.                         its, ite, jts, jte, kts, kte          )
    2052. 

  2052.            
    2053. 

  2053. !          IF ( PRESENT (mfshconv) ) THEN
    2054. 

  2054.               if (mfshconv.EQ.1) &
    2055. 

  2055.               CALL mfshconvpblinit( massflux_EDKF, entr_EDKF, detr_EDKF & 
    2056. 

  2056.                                     ,thl_up, thv_up, rt_up              &
    2057. 

  2057.                                     ,rv_up, rc_up, u_up, v_up           &
    2058. 

  2058.                                     ,frac_up, restart,                  &
    2059. 

  2059.                                     allowed_to_read ,                   &
    2060. 

  2060.                                     ids, ide, jds, jde, kds, kde,       &
    2061. 

  2061.                                     ims, ime, jms, jme, kms, kme,       &
    2062. 

  2062.                                     its, ite, jts, jte, kts, kte   )  
    2063. 

  2063. !          ENDIF
    2064. 

  2064. 

  2065.       CASE (QNSEPBL09SCHEME)
    2066. 

  2066.            if(isfc .ne. 4)CALL wrf_error_fatal &
    2067. 

  2067.             ( 'module_physics_init: use qnsesfc scheme for this pbl option' )
    2068. 

  2068.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    2069. 

  2069.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    2070. 

  2070.            CALL qnsepblinit09(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN, &
    2071. 

  2071.                         TKE_PBL,EXCH_H,restart,               &
    2072. 

  2072.                         allowed_to_read ,                     &
    2073. 

  2073.                         ids, ide, jds, jde, kds, kde,         &
    2074. 

  2074.                         ims, ime, jms, jme, kms, kme,         &
    2075. 

  2075.                         its, ite, jts, jte, kts, kte          )
    2076. 

  2076. 

  2077. #if (NMM_CORE != 1)
    2078. 

  2078.       CASE (BOULACSCHEME)
    2079. 

  2079.            if(isfc .ne. 1 .and. isfc .ne. 2)CALL wrf_error_fatal &
    2080. 

  2080.             ( 'module_physics_init: use sfclay or myjsfc scheme for this pbl option' )
    2081. 

  2081.            CALL boulacinit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN, &
    2082. 

  2082.                          TKE_PBL,EXCH_H,restart,               &
    2083. 

  2083.                          allowed_to_read ,                     &
    2084. 

  2084.                          ids, ide, jds, jde, kds, kde,         &
    2085. 

  2085.                          ims, ime, jms, jme, kms, kme,         &
    2086. 

  2086.                          its, ite, jts, jte, kts, kte          )
    2087. 

  2087.         CASE (CAMUWPBLSCHEME)
    2088. 

  2088.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    2089. 

  2089.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    2090. 

  2090.            CALL camuwpblinit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN, &
    2091. 

  2091.                 restart,TKE_PBL,config_flags%grid_id,           &
    2092. 

  2092.                 ids, ide, jds, jde, kds, kde,         &
    2093. 

  2093.                 ims, ime, jms, jme, kms, kme,         &
    2094. 

  2094.                 its, ite, jts, jte, kts, kte          )
    2095. 

  2095. #endif
    2096. 

  2096. 

  2097. #if ( EM_CORE == 1 )
    2098. 

  2098. 

  2099. !mynn 
    2100. 

  2100.            
    2101. 

  2101.         CASE (MYNNPBLSCHEME2, MYNNPBLSCHEME3)
    2102. 

  2102.            IF(isfc .NE. 5 .AND. isfc .NE. 1 .AND. isfc .NE. 2) CALL wrf_error_fatal &
    2103. 

  2103.                 ( 'module_physics_init: use mynnsfc or sfclay or myjsfc scheme for this pbl option')
    2104. 

  2104.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    2105. 

  2105.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    2106. 

  2106.            
    2107. 

  2107.            SELECT CASE(config_flags%bl_pbl_physics)
    2108. 

  2108. 

  2109.              CASE(MYNNPBLSCHEME2)
    2110. 

  2110.                 mynn_closure_level=2
    2111. 

  2111. 

  2112.              CASE(MYNNPBLSCHEME3)
    2113. 

  2113.                 mynn_closure_level=3
    2114. 

  2114. 

  2115.              CASE DEFAULT
    2116. 

  2116. 

  2117.            END SELECT
    2118. 

  2118. 

  2119.            CALL mynn_bl_init_driver(&
    2120. 

  2120.                 &RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN,RQCBLTEN&
    2121. 

  2121.                 &,restart,allowed_to_read,mynn_closure_level &
    2122. 

  2122.                 &,IDS,IDE,JDS,JDE,KDS,KDE                    &
    2123. 

  2123.                 &,IMS,IME,JMS,JME,KMS,KME                    &
    2124. 

  2124.                 &,ITS,ITE,JTS,JTE,KTS,KTE)
    2125. 

  2125. 

  2126.       CASE (TEMFPBLSCHEME)
    2127. 

  2127.            ! if(isfc .ne. 0)CALL wrf_error_fatal &
    2128. 

  2128.            !  ( 'module_physics_init: use sfclay scheme = 0 for this pbl option' )
    2129. 

  2129.            IF ((SF_URBAN_PHYSICS.eq.2).OR.(SF_URBAN_PHYSICS.EQ.3)) CALL wrf_error_fatal &
    2130. 

  2130.             ( 'module_physics_init: use myj (option 2) or boulac (option 8) with BEP/BEM urban scheme' )
    2131. 

  2131.          IF ( PRESENT( te_temf ) .AND. PRESENT( cf3d_temf )) THEN
    2132. 

  2132.            CALL temfinit(RUBLTEN,RVBLTEN,RTHBLTEN,RQVBLTEN,    &
    2133. 

  2133.                         RQCBLTEN,RQIBLTEN,P_QI,               &
    2134. 

  2134.                         PARAM_FIRST_SCALAR,                   &
    2135. 

  2135.                         restart,                              &
    2136. 

  2136.                         allowed_to_read ,                     &
    2137. 

  2137.                         te_temf,cf3d_temf,                    & ! WA
    2138. 

  2138.                         ids, ide, jds, jde, kds, kde,         &
    2139. 

  2139.                         ims, ime, jms, jme, kms, kme,         &
    2140. 

  2140.                         its, ite, jts, jte, kts, kte          )
    2141. 

  2141.              ELSE
    2142. 

  2142.                 CALL wrf_error_fatal ( 'arguments not present for calling TEMF scheme' )
    2143. 

  2143.          ENDIF
    2144. 

  2144. 

  2145. #endif
    2146. 

  2146. 

  2147.       CASE DEFAULT
    2148. 

  2148. 

  2149.    END SELECT pbl_select
    2150. 

  2150. 

  2151. 

  2152.    END SUBROUTINE bl_init
    2153. 

  2153. 

  2154. !==================================================================
    2155. 

  2155.    SUBROUTINE cu_init(STEPCU,CUDT,DT,RUCUTEN,RVCUTEN,RTHCUTEN,     &
    2156. 

  2156.                       RQVCUTEN,RQRCUTEN,RQCCUTEN,RQSCUTEN,RQICUTEN,&
    2157. 

  2157.                       NCA,RAINC,RAINCV,W0AVG,config_flags,restart, &
    2158. 

  2158.                       CLDEFI,LOWLYR,MASS_FLUX,                     &
    2159. 

  2159.                       RTHFTEN, RQVFTEN,                            &
    2160. 

  2160.                       APR_GR,APR_W,APR_MC,APR_ST,APR_AS,           &
    2161. 

  2161.                       APR_CAPMA,APR_CAPME,APR_CAPMI,               &
    2162. 

  2162.                       cugd_tten,cugd_ttens,cugd_qvten,             &
    2163. 

  2163.                       cugd_qvtens,cugd_qcten,                      &
    2164. 

  2164.                       allowed_to_read, start_of_simulation,        &
    2165. 

  2165.                       ids, ide, jds, jde, kds, kde,                &
    2166. 

  2166.                       ims, ime, jms, jme, kms, kme,                &
    2167. 

  2167.                       its, ite, jts, jte, kts, kte                 )
    2168. 

  2168. !------------------------------------------------------------------
    2169. 

  2169.    USE module_cu_kf
    2170. 

  2170.    USE module_cu_kfeta
    2171. 

  2171.    USE MODULE_CU_BMJ
    2172. 

  2172.    USE module_cu_gd,  ONLY : GDINIT
    2173. 

  2173.    USE module_cu_g3,  ONLY : G3INIT
    2174. 

  2174.    USE module_cu_sas
    2175. 

  2175.    USE module_cu_osas
    2176. 

  2176.    USE module_cu_camzm_driver, ONLY : zm_conv_init
    2177. 

  2177.    USE module_cu_nsas
    2178. 

  2178.    USE module_cu_tiedtke
    2179. 

  2179. !------------------------------------------------------------------
    2180. 

  2180.    IMPLICIT NONE
    2181. 

  2181. !------------------------------------------------------------------
    2182. 

  2182.    TYPE (grid_config_rec_type) ::     config_flags
    2183. 

  2183.    LOGICAL , INTENT(IN)        :: restart
    2184. 

  2184. 

  2185. 

  2186.    INTEGER , INTENT(IN)        :: ids, ide, jds, jde, kds, kde,   &
    2187. 

  2187.                                   ims, ime, jms, jme, kms, kme,   &
    2188. 

  2188.                                   its, ite, jts, jte, kts, kte
    2189. 

  2189. 

  2190.    REAL ,    INTENT(IN)        :: DT, CUDT
    2191. 

  2191.    LOGICAL , INTENT(IN)        :: start_of_simulation
    2192. 

  2192.    LOGICAL , INTENT(IN)        :: allowed_to_read
    2193. 

  2193.    INTEGER , INTENT(INOUT)     :: STEPCU
    2194. 

  2194. 

  2195.    REAL ,   DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(INOUT) ::    &
    2196. 

  2196.             RUCUTEN, RVCUTEN, RTHCUTEN, &
    2197. 

  2197.             RQVCUTEN, RQCCUTEN, RQRCUTEN, RQICUTEN, RQSCUTEN
    2198. 

  2198.    REAL ,   DIMENSION( ims:ime , kms:kme , jms:jme ) , OPTIONAL, INTENT(INOUT) ::    &
    2199. 

  2199.                         cugd_tten,cugd_ttens,cugd_qvten,            &
    2200. 

  2200.                         cugd_qvtens,cugd_qcten
    2201. 

  2201. 

  2202.    REAL ,   DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) :: W0AVG
    2203. 

  2203. 

  2204.    REAL,    DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) ::    &
    2205. 

  2205.             RTHFTEN, RQVFTEN
    2206. 

  2206. 

  2207.    REAL ,   DIMENSION( ims:ime , jms:jme ), INTENT(OUT):: RAINC, RAINCV
    2208. 

  2208. 

  2209.    REAL ,   DIMENSION( ims:ime , jms:jme ), INTENT(OUT):: CLDEFI
    2210. 

  2210. 

  2211.    REAL ,   DIMENSION( ims:ime , jms:jme ), INTENT(INOUT):: NCA
    2212. 

  2212. 

  2213.    REAL ,   DIMENSION( ims:ime , jms:jme ), INTENT(INOUT):: MASS_FLUX,   &
    2214. 

  2214.                                    APR_GR,APR_W,APR_MC,APR_ST,APR_AS,    &
    2215. 

  2215.                                    APR_CAPMA,APR_CAPME,APR_CAPMI
    2216. 

  2216. 

  2217.    INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT):: LOWLYR
    2218. 

  2218. 

  2219. ! LOCAL VAR
    2220. 

  2220. 

  2221.   INTEGER :: i,j,itf,jtf
    2222. 

  2222. 

  2223. !--------------------------------------------------------------------
    2224. 

  2224. 

  2225. !-- calculate cumulus parameterization time step
    2226. 

  2226. 

  2227.    itf=min0(ite,ide-1)
    2228. 

  2228.    jtf=min0(jte,jde-1)
    2229. 

  2229. !
    2230. 

  2230.    STEPCU = nint(CUDT*60./DT)
    2231. 

  2231.    STEPCU = max(STEPCU,1)
    2232. 

  2232. 

  2233. !-- initialization
    2234. 

  2234. 

  2235.    IF(start_of_simulation)THEN
    2236. 

  2236.      DO j=jts,jtf
    2237. 

  2237.      DO i=its,itf
    2238. 

  2238.         RAINC(i,j)=0.
    2239. 

  2239.         RAINCV(i,j)=0.
    2240. 

  2240.      ENDDO
    2241. 

  2241.      ENDDO
    2242. 

  2242.    ENDIF
    2243. 

  2243. 

  2244. !-- deep convection and hybrid deep-shallow convection schemes
    2245. 

  2245.    cps_select: SELECT CASE(config_flags%cu_physics)
    2246. 

  2246. 

  2247.      CASE (KFSCHEME)
    2248. 

  2248.           CALL kfinit(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQRCUTEN,        &
    2249. 

  2249.                       RQICUTEN,RQSCUTEN,NCA,W0AVG,P_QI,P_QS,      &
    2250. 

  2250.                       PARAM_FIRST_SCALAR,restart,                 &
    2251. 

  2251.                       allowed_to_read ,                           &
    2252. 

  2252.                       ids, ide, jds, jde, kds, kde,               &
    2253. 

  2253.                       ims, ime, jms, jme, kms, kme,               &
    2254. 

  2254.                       its, ite, jts, jte, kts, kte                )
    2255. 

  2255. 

  2256.      CASE (BMJSCHEME)
    2257. 

  2257.           CALL bmjinit(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQRCUTEN,       &
    2258. 

  2258.                       CLDEFI,LOWLYR,cp,r_d,restart,               &
    2259. 

  2259.                       allowed_to_read ,                           &
    2260. 

  2260.                       ids, ide, jds, jde, kds, kde,               &
    2261. 

  2261.                       ims, ime, jms, jme, kms, kme,               &
    2262. 

  2262.                       its, ite, jts, jte, kts, kte                )
    2263. 

  2263. 

  2264.      CASE (KFETASCHEME)
    2265. 

  2265.           CALL kf_eta_init(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQRCUTEN,   &
    2266. 

  2266.                       RQICUTEN,RQSCUTEN,NCA,W0AVG,P_QI,P_QS,      &
    2267. 

  2267.                       SVP1,SVP2,SVP3,SVPT0,                       &
    2268. 

  2268.                       PARAM_FIRST_SCALAR,restart,                 &
    2269. 

  2269.                       allowed_to_read ,                           &
    2270. 

  2270.                       ids, ide, jds, jde, kds, kde,               &
    2271. 

  2271.                       ims, ime, jms, jme, kms, kme,               &
    2272. 

  2272.                       its, ite, jts, jte, kts, kte                )
    2273. 

  2273.      CASE (GDSCHEME)
    2274. 

  2274.           CALL gdinit(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQICUTEN,        &
    2275. 

  2275.                       MASS_FLUX,cp,restart,                       &
    2276. 

  2276.                       P_QC,P_QI,PARAM_FIRST_SCALAR,               &
    2277. 

  2277.                       RTHFTEN, RQVFTEN,                           &
    2278. 

  2278.                       APR_GR,APR_W,APR_MC,APR_ST,APR_AS,          &
    2279. 

  2279.                       APR_CAPMA,APR_CAPME,APR_CAPMI,              &
    2280. 

  2280.                       allowed_to_read ,                           &
    2281. 

  2281.                       ids, ide, jds, jde, kds, kde,               &
    2282. 

  2282.                       ims, ime, jms, jme, kms, kme,               &
    2283. 

  2283.                       its, ite, jts, jte, kts, kte                )
    2284. 

  2284.      CASE (NSASSCHEME)
    2285. 

  2285.          CALL nsasinit(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQICUTEN,       &
    2286. 

  2286.                       RUCUTEN,RVCUTEN,                            & 
    2287. 

  2287.                       restart,P_QC,P_QI,PARAM_FIRST_SCALAR,       &
    2288. 

  2288.                       allowed_to_read ,                           &
    2289. 

  2289.                       ids, ide, jds, jde, kds, kde,               &
    2290. 

  2290.                       ims, ime, jms, jme, kms, kme,               &
    2291. 

  2291.                       its, ite, jts, jte, kts, kte                )
    2292. 

  2292. 

  2293. #if ( EM_CORE == 1 )
    2294. 

  2294.      CASE (G3SCHEME)
    2295. 

  2295.           CALL g3init(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQICUTEN,        &
    2296. 

  2296.                       MASS_FLUX,cp,restart,                       &
    2297. 

  2297.                       P_QC,P_QI,PARAM_FIRST_SCALAR,               &
    2298. 

  2298.                       RTHFTEN, RQVFTEN,                           &
    2299. 

  2299.                       APR_GR,APR_W,APR_MC,APR_ST,APR_AS,          &
    2300. 

  2300.                       APR_CAPMA,APR_CAPME,APR_CAPMI,              &
    2301. 

  2301.                       cugd_tten,cugd_ttens,cugd_qvten,            &
    2302. 

  2302.                       cugd_qvtens,cugd_qcten,                     &
    2303. 

  2303.                       allowed_to_read ,                           &
    2304. 

  2304.                       ids, ide, jds, jde, kds, kde,               &
    2305. 

  2305.                       ims, ime, jms, jme, kms, kme,               &
    2306. 

  2306.                       its, ite, jts, jte, kts, kte                )
    2307. 

  2307. #endif
    2308. 

  2308.      CASE (SASSCHEME)
    2309. 

  2309.           CALL sasinit(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQICUTEN,       &
    2310. 

  2310.                       RUCUTEN,RVCUTEN,                            &   ! gopal's doing for SAS
    2311. 

  2311.                       restart,P_QC,P_QI,PARAM_FIRST_SCALAR,       &
    2312. 

  2312.                       allowed_to_read ,                           &
    2313. 

  2313.                       ids, ide, jds, jde, kds, kde,               &
    2314. 

  2314.                       ims, ime, jms, jme, kms, kme,               &
    2315. 

  2315.                       its, ite, jts, jte, kts, kte                )
    2316. 

  2316. 

  2317.      CASE (OSASSCHEME)
    2318. 

  2318.           CALL osasinit(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQICUTEN,      &
    2319. 

  2319.                       RUCUTEN,RVCUTEN,                            &   ! gopal's doing for SAS
    2320. 

  2320.                       restart,P_QC,P_QI,PARAM_FIRST_SCALAR,       &
    2321. 

  2321.                       allowed_to_read ,                           &
    2322. 

  2322.                       ids, ide, jds, jde, kds, kde,               &
    2323. 

  2323.                       ims, ime, jms, jme, kms, kme,               &
    2324. 

  2324.                       its, ite, jts, jte, kts, kte                )
    2325. 

  2325. 

  2326.      CASE (CAMZMSCHEME)
    2327. 

  2327.           CALL zm_conv_init(rucuten, rvcuten, rthcuten, rqvcuten,       &
    2328. 

  2328.                       rqccuten, rqicuten,                               &
    2329. 

  2329.                       p_qc, p_qi, param_first_scalar,                   &
    2330. 

  2330.                       restart,                                          &
    2331. 

  2331.                       ids, ide, jds, jde, kds, kde,                     &
    2332. 

  2332.                       ims, ime, jms, jme, kms, kme,                     &
    2333. 

  2333.                       its, ite, jts, jte, kts, kte                      )
    2334. 

  2334. 

  2335. ! Tiedtke Scheme - ZCX&YQW
    2336. 

  2336.       CASE (TIEDTKESCHEME)
    2337. 

  2337.           CALL tiedtkeinit(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQICUTEN,   &
    2338. 

  2338.                       RUCUTEN,RVCUTEN,                            & 
    2339. 

  2339.                       restart,P_QC,P_QI,PARAM_FIRST_SCALAR,       &
    2340. 

  2340.                       allowed_to_read ,                           &
    2341. 

  2341.                       ids, ide, jds, jde, kds, kde,               &
    2342. 

  2342.                       ims, ime, jms, jme, kms, kme,               &
    2343. 

  2343.                       its, ite, jts, jte, kts, kte                )
    2344. 

  2344. 

  2345.      CASE DEFAULT
    2346. 

  2346. 

  2347.    END SELECT cps_select
    2348. 

  2348. 

  2349.    END SUBROUTINE cu_init
    2350. 

  2350. 

  2351. !==================================================================
    2352. 

  2352.    SUBROUTINE shcu_init(STEPCU,CUDT,DT,RUSHTEN,RVSHTEN,RTHSHTEN,   &
    2353. 

  2353.                       RQVSHTEN,RQRSHTEN,RQCSHTEN,                  &
    2354. 

  2354.                       RQSSHTEN,RQISHTEN,RQGSHTEN,                  &
    2355. 

  2355.                       NCA,RAINC,RAINCV,config_flags,restart,       &
    2356. 

  2356.                       allowed_to_read, start_of_simulation,        &
    2357. 

  2357.                       ids, ide, jds, jde, kds, kde,                &
    2358. 

  2358.                       ims, ime, jms, jme, kms, kme,                &
    2359. 

  2359.                       its, ite, jts, jte, kts, kte                 )
    2360. 

  2360. !------------------------------------------------------------------
    2361. 

  2361.    USE uwshcu,        ONLY: init_uwshcu
    2362. 

  2362.    USE physconst,     ONLY: cpair, gravit, latice, latvap, mwdry, mwh2o, &
    2363. 

  2363.                             rair, zvir
    2364. 

  2364.    USE shr_kind_mod,  ONLY: r8 => shr_kind_r8
    2365. 

  2365. 

  2366. !------------------------------------------------------------------
    2367. 

  2367.    IMPLICIT NONE
    2368. 

  2368. !------------------------------------------------------------------
    2369. 

  2369.    TYPE (grid_config_rec_type) :: config_flags
    2370. 

  2370.    LOGICAL , INTENT(IN)        :: restart
    2371. 

  2371. 

  2372. 

  2373.    INTEGER , INTENT(IN)        :: ids, ide, jds, jde, kds, kde,   &
    2374. 

  2374.                                   ims, ime, jms, jme, kms, kme,   &
    2375. 

  2375.                                   its, ite, jts, jte, kts, kte
    2376. 

  2376. 

  2377.    REAL ,    INTENT(IN)        :: DT, CUDT
    2378. 

  2378.    LOGICAL , INTENT(IN)        :: start_of_simulation
    2379. 

  2379.    LOGICAL , INTENT(IN)        :: allowed_to_read
    2380. 

  2380.    INTEGER , INTENT(INOUT)     :: STEPCU
    2381. 

  2381. 

  2382.    REAL ,   DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(INOUT) ::    &
    2383. 

  2383.             RUSHTEN, RVSHTEN, RTHSHTEN, &
    2384. 

  2384.             RQVSHTEN, RQCSHTEN, RQRSHTEN, RQISHTEN, RQSSHTEN, RQGSHTEN
    2385. 

  2385. 

  2386.    REAL ,   DIMENSION( ims:ime , jms:jme ), INTENT(OUT):: RAINC, RAINCV
    2387. 

  2387. 

  2388.    REAL ,   DIMENSION( ims:ime , jms:jme ), INTENT(INOUT):: NCA
    2389. 

  2389. 

  2390. ! LOCAL VAR
    2391. 

  2391. 

  2392.   INTEGER :: i,j,itf,jtf
    2393. 

  2393. 

  2394. !--------------------------------------------------------------------
    2395. 

  2395. 

  2396. ! Some of this stuff is redundant with deep convection, but redo it
    2397. 

  2397. ! in case deep is turned off...
    2398. 

  2398. 

  2399. !-- calculate cumulus parameterization time step
    2400. 

  2400. 

  2401.    itf=min0(ite,ide-1)
    2402. 

  2402.    jtf=min0(jte,jde-1)
    2403. 

  2403. !
    2404. 

  2404.    STEPCU = nint(CUDT*60./DT)
    2405. 

  2405.    STEPCU = max(STEPCU,1)
    2406. 

  2406. 

  2407. !-- initialization
    2408. 

  2408. 

  2409.    IF(start_of_simulation)THEN
    2410. 

  2410.      DO j=jts,jtf
    2411. 

  2411.      DO i=its,itf
    2412. 

  2412.         RAINC(i,j)=0.
    2413. 

  2413.         RAINCV(i,j)=0.
    2414. 

  2414.      ENDDO
    2415. 

  2415.      ENDDO
    2416. 

  2416.    ENDIF
    2417. 

  2417. 

  2418. !-- independent shallow convection schemes
    2419. 

  2419.    shcu_select: SELECT CASE(config_flags%shcu_physics)
    2420. 

  2420. 

  2421.    CASE (CAMUWSHCUSCHEME)
    2422. 

  2422.       CALL init_uwshcu(r8,latvap,cpair,latice,zvir,rair,gravit,      &
    2423. 

  2423.            mwh2o/mwdry,                                              &
    2424. 

  2424.            rushten, rvshten, rthshten, rqvshten,                     &
    2425. 

  2425.            rqcshten, rqrshten, rqishten, rqsshten, rqgshten,         &
    2426. 

  2426.            p_qc, p_qr, p_qi, p_qs, p_qg,                             &
    2427. 

  2427.            config_flags%bl_pbl_physics, param_first_scalar, restart, &
    2428. 

  2428.            config_flags%grid_id,                                     &
    2429. 

  2429.            ids, ide, jds, jde, kds, kde,                             &
    2430. 

  2430.            ims, ime, jms, jme, kms, kme,                             &
    2431. 

  2431.            its, ite, jts, jte, kts, kte                              )
    2432. 

  2432. 

  2433.    CASE DEFAULT
    2434. 

  2434. 

  2435.    END SELECT shcu_select
    2436. 

  2436. 

  2437.    END SUBROUTINE shcu_init
    2438. 

  2438. 

  2439. !==================================================================
    2440. 

  2440.    SUBROUTINE mp_init(RAINNC,SNOWNC,GRAUPELNC,config_flags,restart,warm_rain,      &
    2441. 

  2441.                       adv_moist_cond,                             &
    2442. 

  2442.                       MPDT, DT, DX, DY, LOWLYR,                   & ! for eta mp
    2443. 

  2443.                       F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY,           & ! for eta mp
    2444. 

  2444.                       mp_restart_state,tbpvs_state,tbpvs0_state,   & ! eta mp
    2445. 

  2445.                       allowed_to_read, start_of_simulation,       &
    2446. 

  2446.                       ids, ide, jds, jde, kds, kde,               &
    2447. 

  2447.                       ims, ime, jms, jme, kms, kme,               &
    2448. 

  2448.                       its, ite, jts, jte, kts, kte                )
    2449. 

  2449. !------------------------------------------------------------------
    2450. 

  2450.    USE module_mp_wsm3
    2451. 

  2451.    USE module_mp_wsm5
    2452. 

  2452.    USE module_mp_wsm6
    2453. 

  2453.    USE module_mp_etanew
    2454. 

  2454.    USE module_mp_etaold
    2455. 

  2455. #if (NMM_CORE == 1)
    2456. 

  2456.    USE module_mp_HWRF
    2457. 

  2457. #endif
    2458. 

  2458.    USE module_mp_thompson
    2459. 

  2459.    USE module_mp_morr_two_moment
    2460. 

  2460.    USE module_mp_milbrandt2mom
    2461. 

  2461. !  USE module_mp_milbrandt3mom
    2462. 

  2462.    USE module_mp_wdm5
    2463. 

  2463.    USE module_mp_wdm6
    2464. 

  2464. #if (EM_CORE==1)
    2465. 

  2465.    USE module_mp_nssl_2mom
    2466. 

  2466. #endif
    2467. 

  2467. !------------------------------------------------------------------
    2468. 

  2468.    IMPLICIT NONE
    2469. 

  2469. !------------------------------------------------------------------
    2470. 

  2470. ! Arguments
    2471. 

  2471.    TYPE (grid_config_rec_type) ::     config_flags
    2472. 

  2472.    LOGICAL , INTENT(IN)        :: restart
    2473. 

  2473.    LOGICAL , INTENT(OUT)       :: warm_rain,adv_moist_cond
    2474. 

  2474.    REAL    , INTENT(IN)        :: MPDT, DT, DX, DY
    2475. 

  2475.    LOGICAL , INTENT(IN)        :: start_of_simulation
    2476. 

  2476. 

  2477.    INTEGER , INTENT(IN)        :: ids, ide, jds, jde, kds, kde,   &
    2478. 

  2478.                                   ims, ime, jms, jme, kms, kme,   &
    2479. 

  2479.                                   its, ite, jts, jte, kts, kte
    2480. 

  2480. 

  2481.    INTEGER , DIMENSION( ims:ime , jms:jme ) ,INTENT(INOUT)  :: LOWLYR
    2482. 

  2482.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(INOUT) :: RAINNC,SNOWNC,GRAUPELNC
    2483. 

  2483.    REAL,     DIMENSION( ims:ime , kms:kme, jms:jme ) , INTENT(INOUT) :: &
    2484. 

  2484.                                   F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY
    2485. 

  2485.    REAL , DIMENSION(:) ,INTENT(INOUT)  :: mp_restart_state,tbpvs_state,tbpvs0_state
    2486. 

  2486.    LOGICAL , INTENT(IN)  :: allowed_to_read
    2487. 

  2487. 

  2488. ! Local
    2489. 

  2489.    INTEGER :: i, j, itf, jtf
    2490. 

  2490. 

  2491.    warm_rain = .false.
    2492. 

  2492.    adv_moist_cond = .true.
    2493. 

  2493.    itf=min0(ite,ide-1)
    2494. 

  2494.    jtf=min0(jte,jde-1)
    2495. 

  2495. 

  2496.    IF(start_of_simulation)THEN
    2497. 

  2497.      DO j=jts,jtf
    2498. 

  2498.      DO i=its,itf
    2499. 

  2499.         RAINNC(i,j) = 0.
    2500. 

  2500.         SNOWNC(i,j) = 0.
    2501. 

  2501.         GRAUPELNC(i,j) = 0.
    2502. 

  2502.      ENDDO
    2503. 

  2503.      ENDDO
    2504. 

  2504.    ENDIF
    2505. 

  2505. 

  2506.    mp_select: SELECT CASE(config_flags%mp_physics)
    2507. 

  2507. 

  2508.      CASE (KESSLERSCHEME)
    2509. 

  2509.           warm_rain = .true.
    2510. 

  2510.      CASE (WSM3SCHEME)
    2511. 

  2511.           CALL wsm3init(rhoair0,rhowater,rhosnow,cliq,cpv, allowed_to_read )
    2512. 

  2512.      CASE (WSM5SCHEME)
    2513. 

  2513.           CALL wsm5init(rhoair0,rhowater,rhosnow,cliq,cpv, allowed_to_read )
    2514. 

  2514.      CASE (WSM6SCHEME)
    2515. 

  2515.           CALL wsm6init(rhoair0,rhowater,rhosnow,cliq,cpv, allowed_to_read )
    2516. 

  2516.      CASE (ETAMPNEW)
    2517. 

  2517.          adv_moist_cond = .false.
    2518. 

  2518.          CALL etanewinit (MPDT,DT,DX,DY,LOWLYR,restart,           &
    2519. 

  2519.                           F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY,       &
    2520. 

  2520.                           mp_restart_state,tbpvs_state,tbpvs0_state,&
    2521. 

  2521.                           allowed_to_read,                        &
    2522. 

  2522.                           ids, ide, jds, jde, kds, kde,           &
    2523. 

  2523.                           ims, ime, jms, jme, kms, kme,           &
    2524. 

  2524.                           its, ite, jts, jte, kts, kte            )
    2525. 

  2525. #if(NMM_CORE==1)
    2526. 

  2526.      CASE (etamp_HWRF)
    2527. 

  2527.          CALL etanewinit_HWRF (MPDT,DT,DX,DY,LOWLYR,restart,           &
    2528. 

  2528.                           F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY,       &
    2529. 

  2529.                           allowed_to_read,                        &
    2530. 

  2530.                           ids, ide, jds, jde, kds, kde,           &
    2531. 

  2531.                           ims, ime, jms, jme, kms, kme,           &
    2532. 

  2532.                           its, ite, jts, jte, kts, kte            )
    2533. 

  2533. #endif 
    2534. 

  2534.      CASE (THOMPSON)
    2535. 

  2535. ! Cycling the WRF forecast with moving nests will cause this initialization to be
    2536. 

  2536. ! called for each nest move. This is potentially very computationally expensive.
    2537. 

  2537.          IF(start_of_simulation.or.restart.or.config_flags%cycling)CALL thompson_init
    2538. 

  2538. 

  2539.      CASE (MORR_TWO_MOMENT)
    2540. 

  2540.          CALL morr_two_moment_init
    2541. 

  2541.      CASE (MILBRANDT2MOM)
    2542. 

  2542.          CALL milbrandt2mom_init
    2543. 

  2543. !      CASE (MILBRANDT3MOM)
    2544. 

  2544. !          CALL milbrandt3mom_init
    2545. 

  2545.      CASE (WDM5SCHEME)
    2546. 

  2546.           CALL wdm5init(rhoair0,rhowater,rhosnow,cliq,cpv,n_ccn0,allowed_to_read )
    2547. 

  2547.      CASE (WDM6SCHEME)
    2548. 

  2548.           CALL wdm6init(rhoair0,rhowater,rhosnow,cliq,cpv,n_ccn0,allowed_to_read )
    2549. 

  2549. #if (EM_CORE==1)
    2550. 

  2550.      CASE (NSSL_2MOM)
    2551. 

  2551.          CALL nssl_2mom_init(ims,ime, jms,jme, kms,kme)
    2552. 

  2552.      CASE (NSSL_2MOMCCN)
    2553. 

  2553.          CALL nssl_2mom_init(ims,ime, jms,jme, kms,kme)
    2554. 

  2554. #endif
    2555. 

  2555.      CASE (ETAMPOLD)
    2556. 

  2556.          adv_moist_cond = .false.
    2557. 

  2557.          CALL etaoldinit (MPDT,DT,DX,DY,LOWLYR,restart,           &
    2558. 

  2558.                           F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY,       &
    2559. 

  2559.                           mp_restart_state,tbpvs_state,tbpvs0_state,&
    2560. 

  2560.                           allowed_to_read,                        &
    2561. 

  2561.                           ids, ide, jds, jde, kds, kde,           &
    2562. 

  2562.                           ims, ime, jms, jme, kms, kme,           &
    2563. 

  2563.                           its, ite, jts, jte, kts, kte            )
    2564. 

  2564. 

  2565.      CASE DEFAULT
    2566. 

  2566. 

  2567.    END SELECT mp_select
    2568. 

  2568. 

  2569.    END SUBROUTINE mp_init
    2570. 

  2570. 

  2571. #if  ( EM_CORE == 1 )
    2572. 

  2572. !==========================================================
    2573. 

  2573.    SUBROUTINE fg_init(STEPFG,FGDT,DT,id,RUNDGDTEN,RVNDGDTEN,    &
    2574. 

  2574.                 RTHNDGDTEN,RPHNDGDTEN,RQVNDGDTEN,RMUNDGDTEN,    &
    2575. 

  2575.                 config_flags,restart,                           &
    2576. 

  2576.                 allowed_to_read ,                               &
    2577. 

  2577.                 ids, ide, jds, jde, kds, kde,                   &
    2578. 

  2578.                 ims, ime, jms, jme, kms, kme,                   &
    2579. 

  2579.                 its, ite, jts, jte, kts, kte                    )
    2580. 

  2580. 

  2581. 

  2582. !--------------------------------------------------------------------
    2583. 

  2583.    USE module_fdda_psufddagd
    2584. 

  2584.    USE module_fdda_spnudging, ONLY : fddaspnudginginit
    2585. 

  2585. !--------------------------------------------------------------------
    2586. 

  2586.    IMPLICIT NONE
    2587. 

  2587. !--------------------------------------------------------------------
    2588. 

  2588.    TYPE (grid_config_rec_type) ::     config_flags
    2589. 

  2589.    LOGICAL , INTENT(IN)        :: restart
    2590. 

  2590. 

  2591.    INTEGER , INTENT(IN)        ::     ids, ide, jds, jde, kds, kde, &
    2592. 

  2592.                                       ims, ime, jms, jme, kms, kme, &
    2593. 

  2593.                                       its, ite, jts, jte, kts, kte
    2594. 

  2594. 

  2595.    REAL ,    INTENT(IN)        ::     DT, FGDT
    2596. 

  2596.    INTEGER , INTENT(IN)        ::     id
    2597. 

  2597.    INTEGER , INTENT(INOUT)     ::     STEPFG
    2598. 

  2598.    REAL,     DIMENSION( ims:ime , kms:kme , jms:jme ) , INTENT(OUT) ::       &
    2599. 

  2599.                                                            RUNDGDTEN, &
    2600. 

  2600.                                                            RVNDGDTEN, &
    2601. 

  2601.                                                           RTHNDGDTEN, &
    2602. 

  2602.                                                           RPHNDGDTEN, &
    2603. 

  2603.                                                           RQVNDGDTEN
    2604. 

  2604.    REAL,     DIMENSION( ims:ime , jms:jme ) , INTENT(OUT) :: RMUNDGDTEN
    2605. 

  2605. 

  2606.    LOGICAL,  INTENT(IN)           :: allowed_to_read
    2607. 

  2607. !--------------------------------------------------------------------
    2608. 

  2608. 

  2609. !-- calculate pbl time step
    2610. 

  2610. 

  2611.    STEPFG = nint(FGDT*60./DT)
    2612. 

  2612.    STEPFG = max(STEPFG,1)
    2613. 

  2613. 

  2614. 

  2615. !-- initialize fdda scheme
    2616. 

  2616. 

  2617.    fdda_select: SELECT CASE(config_flags%grid_fdda)
    2618. 

  2618. 

  2619.       CASE (PSUFDDAGD)
    2620. 

  2620.            CALL fddagdinit(id,rundgdten,rvndgdten,rthndgdten,rqvndgdten,rmundgdten,&
    2621. 

  2621.                config_flags%run_hours, &
    2622. 

  2622.                config_flags%if_no_pbl_nudging_uv, &
    2623. 

  2623.                config_flags%if_no_pbl_nudging_t, &
    2624. 

  2624.                config_flags%if_no_pbl_nudging_q, &
    2625. 

  2625.                config_flags%if_zfac_uv, &
    2626. 

  2626.                config_flags%k_zfac_uv, &
    2627. 

  2627.                config_flags%if_zfac_t, &
    2628. 

  2628.                config_flags%k_zfac_t, &
    2629. 

  2629.                config_flags%if_zfac_q, &
    2630. 

  2630.                config_flags%k_zfac_q, &
    2631. 

  2631.                config_flags%guv, &
    2632. 

  2632.                config_flags%gt, config_flags%gq, &
    2633. 

  2633.                config_flags%if_ramping, config_flags%dtramp_min, &
    2634. 

  2634.                config_flags%auxinput10_end_h, &
    2635. 

  2635.                config_flags%grid_sfdda, &
    2636. 

  2636.                config_flags%guv_sfc, &
    2637. 

  2637.                config_flags%gt_sfc, &
    2638. 

  2638.                config_flags%gq_sfc, &
    2639. 

  2639.                       restart, allowed_to_read,                    &
    2640. 

  2640.                       ids, ide, jds, jde, kds, kde,                &
    2641. 

  2641.                       ims, ime, jms, jme, kms, kme,                &
    2642. 

  2642.                       its, ite, jts, jte, kts, kte                 )
    2643. 

  2643. 

  2644.       CASE (SPNUDGING)
    2645. 

  2645.            CALL fddaspnudginginit(id,rundgdten,rvndgdten,rthndgdten,rphndgdten,&
    2646. 

  2646.                config_flags%run_hours, &
    2647. 

  2647.                config_flags%if_no_pbl_nudging_uv, &
    2648. 

  2648.                config_flags%if_no_pbl_nudging_t, &
    2649. 

  2649.                config_flags%if_no_pbl_nudging_ph, &
    2650. 

  2650.                config_flags%if_zfac_uv, &
    2651. 

  2651.                config_flags%k_zfac_uv, &
    2652. 

  2652.                config_flags%dk_zfac_uv, &
    2653. 

  2653.                config_flags%if_zfac_t, &
    2654. 

  2654.                config_flags%k_zfac_t, &
    2655. 

  2655.                config_flags%dk_zfac_t, &
    2656. 

  2656.                config_flags%if_zfac_ph, &
    2657. 

  2657.                config_flags%k_zfac_ph, &
    2658. 

  2658.                config_flags%dk_zfac_ph, &
    2659. 

  2659.                config_flags%guv, &
    2660. 

  2660.                config_flags%gt, config_flags%gph, &
    2661. 

  2661.                config_flags%if_ramping, config_flags%dtramp_min, &
    2662. 

  2662.                config_flags%auxinput9_end_h, &
    2663. 

  2663.                config_flags%xwavenum,config_flags%ywavenum, &
    2664. 

  2664.                       restart, allowed_to_read,                    &
    2665. 

  2665.                       ids, ide, jds, jde, kds, kde,                &
    2666. 

  2666.                       ims, ime, jms, jme, kms, kme,                &
    2667. 

  2667.                       its, ite, jts, jte, kts, kte                 )
    2668. 

  2668. 

  2669.       CASE DEFAULT
    2670. 

  2670. 

  2671.    END SELECT fdda_select
    2672. 

  2672. 

  2673.    END SUBROUTINE fg_init
    2674. 

  2674. 

  2675. !-------------------------------------------------------------------
    2676. 

  2676.    SUBROUTINE fdob_init(obs_nudge_opt, maxdom, inest, parid,       &
    2677. 

  2677.                         idynin, dtramp, fdaend, restart,           &
    2678. 

  2678.                         obs_twindo_cg, obs_twindo, itimestep,      &
    2679. 

  2679.                         no_pbl_nudge_uv,                           &
    2680. 

  2680.                         no_pbl_nudge_t,                            &
    2681. 

  2681.                         no_pbl_nudge_q,                            &
    2682. 

  2682.                         sfc_scheme_horiz, sfc_scheme_vert,         &
    2683. 

  2683.                         maxsnd_gap,                                &
    2684. 

  2684.                         sfcfact, sfcfacr, dpsmx,                   &
    2685. 

  2685.                         nudge_wind, nudge_temp, nudge_mois,        &
    2686. 

  2686.                         nudgezfullr1_uv, nudgezrampr1_uv,          &
    2687. 

  2687.                         nudgezfullr2_uv, nudgezrampr2_uv,          &
    2688. 

  2688.                         nudgezfullr4_uv, nudgezrampr4_uv,          &
    2689. 

  2689.                         nudgezfullr1_t,  nudgezrampr1_t,           &
    2690. 

  2690.                         nudgezfullr2_t,  nudgezrampr2_t,           &
    2691. 

  2691.                         nudgezfullr4_t,  nudgezrampr4_t,           &
    2692. 

  2692.                         nudgezfullr1_q,  nudgezrampr1_q,           &
    2693. 

  2693.                         nudgezfullr2_q,  nudgezrampr2_q,           &
    2694. 

  2694.                         nudgezfullr4_q,  nudgezrampr4_q,           &
    2695. 

  2695.                         nudgezfullmin, nudgezrampmin, nudgezmax,   &
    2696. 

  2696.                         xlat, xlong,                               &
    2697. 

  2697.                         start_year, start_month, start_day,        &
    2698. 

  2698.                         start_hour, start_minute, start_second,    &
    2699. 

  2699.                         p00, t00, tlp,                             &
    2700. 

  2700.                         znu, p_top,                                &
    2701. 

  2701.                         fdob, ipf_init,                            &
    2702. 

  2702.                         ids, ide, jds, jde, kds, kde,              &
    2703. 

  2703.                         ims, ime, jms, jme, kms, kme,              &
    2704. 

  2704.                         its, ite, jts, jte, kts, kte               )
    2705. 

  2705. 

  2706. !--------------------------------------------------------------------
    2707. 

  2707.    USE module_domain
    2708. 

  2708.    USE module_fddaobs_rtfdda
    2709. 

  2709.    USE module_llxy
    2710. 

  2710. !--------------------------------------------------------------------
    2711. 

  2711.    IMPLICIT NONE
    2712. 

  2712. !--------------------------------------------------------------------
    2713. 

  2713.    INTEGER , INTENT(IN)    :: maxdom
    2714. 

  2714.    INTEGER , INTENT(IN)    :: obs_nudge_opt(maxdom)
    2715. 

  2715.    INTEGER , INTENT(IN)    :: ids,ide, jds,jde, kds,kde,           &
    2716. 

  2716.                               ims,ime, jms,jme, kms,kme,           &
    2717. 

  2717.                               its,ite, jts,jte, kts,kte
    2718. 

  2718.    INTEGER , INTENT(IN)    :: inest
    2719. 

  2719.    INTEGER , INTENT(IN)    :: parid(maxdom)
    2720. 

  2720.    INTEGER , INTENT(IN)    :: idynin          ! flag for dynamic initialization
    2721. 

  2721.    REAL    , INTENT(IN)    :: dtramp          ! time period for ramping (idynin)
    2722. 

  2722.    REAL    , INTENT(IN)    :: fdaend(maxdom)  ! nudging end time for domain (min)
    2723. 

  2723.    LOGICAL , INTENT(IN)    :: restart
    2724. 

  2724.    REAL    , INTENT(IN)    :: obs_twindo_cg   ! twindo on course grid
    2725. 

  2725.    REAL    , INTENT(IN)    :: obs_twindo
    2726. 

  2726.    INTEGER , INTENT(IN)    :: itimestep
    2727. 

  2727.    INTEGER , INTENT(IN)    :: no_pbl_nudge_uv(maxdom)  ! flags for no wind nudging in pbl
    2728. 

  2728.    INTEGER , INTENT(IN)    :: no_pbl_nudge_t(maxdom)   ! flags for no temperature nudging in pbl
    2729. 

  2729.    INTEGER , INTENT(IN)    :: no_pbl_nudge_q(maxdom)   ! flags for no moisture nudging in pbl
    2730. 

  2730.    INTEGER , INTENT(IN)    :: sfc_scheme_horiz ! horizontal spreading scheme for surf obs (wrf or orig mm5)
    2731. 

  2731.    INTEGER , INTENT(IN)    :: sfc_scheme_vert  ! vertical   spreading scheme for surf obs (orig or regime vif)
    2732. 

  2732.    REAL    , INTENT(IN)    :: maxsnd_gap       ! max allowed pressure gap in soundings for interp (centibars) 
    2733. 

  2733.    REAL    , INTENT(IN)    :: sfcfact      ! scale factor applied to time window for surface obs
    2734. 

  2734.    REAL    , INTENT(IN)    :: sfcfacr      ! scale fac applied to horiz rad of infl for sfc obs
    2735. 

  2735.    REAL    , INTENT(IN)    :: dpsmx        ! max pressure change allowed within horiz. infl. range
    2736. 

  2736.    INTEGER , INTENT(IN)    :: nudge_wind(maxdom)       ! wind-nudging flag
    2737. 

  2737.    INTEGER , INTENT(IN)    :: nudge_temp(maxdom)       ! temperature-nudging flag
    2738. 

  2738.    INTEGER , INTENT(IN)    :: nudge_mois(maxdom)       ! moisture-nudging flag
    2739. 

  2739.    REAL    , INTENT(IN)    :: nudgezfullr1_uv  ! vert infl fcn, regime=1 full-wt   hght, winds
    2740. 

  2740.    REAL    , INTENT(IN)    :: nudgezrampr1_uv  ! vert infl fcn, regime=1 ramp down hght, winds
    2741. 

  2741.    REAL    , INTENT(IN)    :: nudgezfullr2_uv  ! vert infl fcn, regime=2 full-wt   hght, winds
    2742. 

  2742.    REAL    , INTENT(IN)    :: nudgezrampr2_uv  ! vert infl fcn, regime=2 ramp down hght, winds
    2743. 

  2743.    REAL    , INTENT(IN)    :: nudgezfullr4_uv  ! vert infl fcn, regime=4 full-wt   hght, winds
    2744. 

  2744.    REAL    , INTENT(IN)    :: nudgezrampr4_uv  ! vert infl fcn, regime=4 ramp down hght, winds
    2745. 

  2745.    REAL    , INTENT(IN)    :: nudgezfullr1_t   ! vert infl fcn, regime=1 full-wt   hght, temp
    2746. 

  2746.    REAL    , INTENT(IN)    :: nudgezrampr1_t   ! vert infl fcn, regime=1 ramp down hght, temp
    2747. 

  2747.    REAL    , INTENT(IN)    :: nudgezfullr2_t   ! vert infl fcn, regime=2 full-wt   hght, temp
    2748. 

  2748.    REAL    , INTENT(IN)    :: nudgezrampr2_t   ! vert infl fcn, regime=2 ramp down hght, temp
    2749. 

  2749.    REAL    , INTENT(IN)    :: nudgezfullr4_t   ! vert infl fcn, regime=4 full-wt   hght, temp
    2750. 

  2750.    REAL    , INTENT(IN)    :: nudgezrampr4_t   ! vert infl fcn, regime=4 ramp down hght, temp
    2751. 

  2751.    REAL    , INTENT(IN)    :: nudgezfullr1_q   ! vert infl fcn, regime=1 full-wt   hght, mois
    2752. 

  2752.    REAL    , INTENT(IN)    :: nudgezrampr1_q   ! vert infl fcn, regime=1 ramp down hght, mois
    2753. 

  2753.    REAL    , INTENT(IN)    :: nudgezfullr2_q   ! vert infl fcn, regime=2 full-wt   hght, mois
    2754. 

  2754.    REAL    , INTENT(IN)    :: nudgezrampr2_q   ! vert infl fcn, regime=2 ramp down hght, mois
    2755. 

  2755.    REAL    , INTENT(IN)    :: nudgezfullr4_q   ! vert infl fcn, regime=4 full-wt   hght, mois
    2756. 

  2756.    REAL    , INTENT(IN)    :: nudgezrampr4_q   ! vert infl fcn, regime=4 ramp down hght, mois
    2757. 

  2757.    REAL    , INTENT(IN)    :: nudgezfullmin    ! min dpth thru which vert infl fcn remains 1.0 (m)
    2758. 

  2758.    REAL    , INTENT(IN)    :: nudgezrampmin    ! min dpth thru which vif decreases 1.0 to 0.0 (m)
    2759. 

  2759.    REAL    , INTENT(IN)    :: nudgezmax        ! max dpth in which vif is nonzero (m)
    2760. 

  2760.    REAL    , INTENT(IN)    :: xlat ( ims:ime, jms:jme )        ! latitudes on mass-point grid
    2761. 

  2761.    REAL    , INTENT(IN)    :: xlong( ims:ime, jms:jme )        ! longitudes on mass-point grid
    2762. 

  2762.    INTEGER , INTENT(INOUT) :: start_year
    2763. 

  2763.    INTEGER , INTENT(INOUT) :: start_month
    2764. 

  2764.    INTEGER , INTENT(INOUT) :: start_day
    2765. 

  2765.    INTEGER , INTENT(INOUT) :: start_hour
    2766. 

  2766.    INTEGER , INTENT(INOUT) :: start_minute
    2767. 

  2767.    INTEGER , INTENT(INOUT) :: start_second
    2768. 

  2768.    REAL    , INTENT(IN)    :: p00                      ! base state pressure
    2769. 

  2769.    REAL    , INTENT(IN)    :: t00                      ! base state temperature
    2770. 

  2770.    REAL    , INTENT(IN)    :: tlp                      ! base state lapse rate
    2771. 

  2771.    REAL    , INTENT(IN)    :: znu( kms:kme )           ! eta values on half (mass) levels
    2772. 

  2772.    REAL    , INTENT(IN)    :: p_top                    ! pressure at top of model
    2773. 

  2773.    TYPE(fdob_type), INTENT(INOUT)  :: fdob
    2774. 

  2774. 

  2775.    INTEGER                 :: e_sn         ! ending   north-south grid index
    2776. 

  2776.    LOGICAL                 :: ipf_init     ! print warnings detected at initialzn
    2777. 

  2777. !--------------------------------------------------------------------
    2778. 

  2778. !-- initialize fdda obs-nudging scheme
    2779. 

  2779. 

  2780.       IF ( obs_nudge_opt(inest) .eq. 0 ) RETURN
    2781. 

  2781. 

  2782.       e_sn = jde
    2783. 

  2783.       CALL fddaobs_init(obs_nudge_opt, maxdom, inest, parid,       &
    2784. 

  2784.                         idynin, dtramp, fdaend, restart,           &
    2785. 

  2785.                         obs_twindo_cg,                             &
    2786. 

  2786.                         obs_twindo, itimestep,                     &
    2787. 

  2787.                         no_pbl_nudge_uv,                           &
    2788. 

  2788.                         no_pbl_nudge_t,                            &
    2789. 

  2789.                         no_pbl_nudge_q,                            &
    2790. 

  2790.                         sfc_scheme_horiz, sfc_scheme_vert,         &
    2791. 

  2791.                         maxsnd_gap,                                &
    2792. 

  2792.                         sfcfact, sfcfacr, dpsmx,                   &
    2793. 

  2793.                         nudge_wind, nudge_temp, nudge_mois,        &
    2794. 

  2794.                         nudgezfullr1_uv, nudgezrampr1_uv,          &
    2795. 

  2795.                         nudgezfullr2_uv, nudgezrampr2_uv,          &
    2796. 

  2796.                         nudgezfullr4_uv, nudgezrampr4_uv,          &
    2797. 

  2797.                         nudgezfullr1_t,  nudgezrampr1_t,           &
    2798. 

  2798.                         nudgezfullr2_t,  nudgezrampr2_t,           &
    2799. 

  2799.                         nudgezfullr4_t,  nudgezrampr4_t,           &
    2800. 

  2800.                         nudgezfullr1_q,  nudgezrampr1_q,           &
    2801. 

  2801.                         nudgezfullr2_q,  nudgezrampr2_q,           &
    2802. 

  2802.                         nudgezfullr4_q,  nudgezrampr4_q,           &
    2803. 

  2803.                         nudgezfullmin,  nudgezrampmin, nudgezmax,  &
    2804. 

  2804.                         xlat, xlong,                               &
    2805. 

  2805.                         start_year, start_month, start_day,        &
    2806. 

  2806.                         start_hour, start_minute, start_second,    &
    2807. 

  2807.                         p00, t00, tlp,                             &
    2808. 

  2808.                         znu, p_top,                                &
    2809. 

  2809.                         fdob, ipf_init,                            &
    2810. 

  2810.                         ids,ide, jds,jde, kds,kde,                 &
    2811. 

  2811.                         ims,ime, jms,jme, kms,kme,                 &
    2812. 

  2812.                         its,ite, jts,jte, kts,kte)
    2813. 

  2813. 

  2814.    END SUBROUTINE fdob_init
    2815. 

  2815. #endif
    2816. 

  2816. 

  2817. !--------------------------------------------------------------------
    2818. 

  2818.    SUBROUTINE z2sigma(zf,zh,sf,sh,p_top,pptop,config_flags, &
    2819. 

  2819.                 allowed_to_read , &
    2820. 

  2820.                 kds,kde,kms,kme,kts,kte)
    2821. 

  2821.    IMPLICIT NONE
    2822. 

  2822. ! Arguments
    2823. 

  2823.    INTEGER, INTENT(IN) :: kds,kde,kms,kme,kts,kte
    2824. 

  2824.    REAL , DIMENSION( kms:kme ), INTENT(IN) :: zf,zh
    2825. 

  2825.    REAL , DIMENSION( kms:kme ), INTENT(OUT):: sf,sh
    2826. 

  2826.    REAL , INTENT(IN) :: p_top
    2827. 

  2827.    REAL , INTENT(OUT) :: pptop
    2828. 

  2828.    TYPE (grid_config_rec_type)              :: config_flags
    2829. 

  2829.    LOGICAL , INTENT(IN) :: allowed_to_read
    2830. 

  2830. ! Local
    2831. 

  2831.    REAL R, G, TS, GAMMA, PS, ZTROP, TSTRAT, PTROP, Z, T, P, ZTOP, PTOP
    2832. 

  2832.    INTEGER K
    2833. 

  2833. 

  2834.    IF(zf(kde/2) .GT. 1.0)THEN
    2835. 

  2835. ! Height levels assumed (zeta coordinate)
    2836. 

  2836. ! Convert to sigma using standard atmosphere for pressure-height relation
    2837. 

  2837. ! constants for standard atmosphere definition
    2838. 

  2838.       r=287.05
    2839. 

  2839.       g=9.80665
    2840. 

  2840.       ts=288.15
    2841. 

  2841.       gamma=-6.5/1000.
    2842. 

  2842.       ps=1013.25
    2843. 

  2843.       ztrop=11000.
    2844. 

  2844.       tstrat=ts+gamma*ztrop
    2845. 

  2845.       ptrop=ps*(tstrat/ts)**(-g/(gamma*r))
    2846. 

  2846. 

  2847.       do k=kde,kds,-1
    2848. 

  2848. ! full levels
    2849. 

  2849.         z=zf(k)
    2850. 

  2850.         if(z.le.ztrop)then
    2851. 

  2851.           t=ts+gamma*z
    2852. 

  2852.           p=ps*(t/ts)**(-g/(gamma*r))
    2853. 

  2853.         else
    2854. 

  2854.           t=tstrat
    2855. 

  2855.           p=ptrop*exp(-g*(z-ztrop)/(r*tstrat))
    2856. 

  2856.         endif
    2857. 

  2857.         if(k.eq.kde)then
    2858. 

  2858.           ztop=zf(k)
    2859. 

  2859.           ptop=p
    2860. 

  2860.         endif
    2861. 

  2861.         sf(k)=(p-ptop)/(ps-ptop)
    2862. 

  2862. ! half levels
    2863. 

  2863.         if(k.ne.kds)then
    2864. 

  2864.         z=0.5*(zf(k)+zf(k-1))
    2865. 

  2865.         if(z.le.ztrop)then
    2866. 

  2866.           t=ts+gamma*z
    2867. 

  2867.           p=ps*(t/ts)**(-g/(gamma*r))
    2868. 

  2868.         else
    2869. 

  2869.           t=tstrat
    2870. 

  2870.           p=ptrop*exp(-g*(z-ztrop)/(r*tstrat))
    2871. 

  2871.         endif
    2872. 

  2872.         sh(k-1)=(p-ptop)/(ps-ptop)
    2873. 

  2873.         endif
    2874. 

  2874.       enddo
    2875. 

  2875.       pptop=ptop/10.
    2876. 

  2876.    ELSE
    2877. 

  2877. !  Levels are already sigma/eta
    2878. 

  2878.       do k=kde,kds,-1
    2879. 

  2879. !        sf(k)=zf(kde-k+kds)
    2880. 

  2880. !        if(k .ne. kde)sh(k)=zh(kde-1-k+kds)
    2881. 

  2881.          sf(k)=zf(k)
    2882. 

  2882.          if(k .ne. kde)sh(k)=zh(k)
    2883. 

  2883.       enddo
    2884. 

  2884.       pptop=p_top/1000.
    2885. 

  2885. 

  2886.    ENDIF
    2887. 

  2887. 

  2888.    END SUBROUTINE z2sigma
    2889. 

  2889. 

  2890. END MODULE module_physics_init
    2891. 

  2891.