PageRenderTime 54ms CodeModel.GetById 16ms 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

Large files files are truncated, but you can click here to view the full file

    

  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.     E…
    1004.

Large files files are truncated, but you can click here to view the full file