PageRenderTime 60ms CodeModel.GetById 22ms RepoModel.GetById 1ms app.codeStats 1ms

/wrfv2_fire/phys/module_surface_driver.F

http://github.com/jbeezley/wrf-fire
FORTRAN Legacy | 4662 lines | 3098 code | 345 blank | 1219 comment | 27 complexity | 2cebd156424e3b0f3c100847e14ccdf1 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:MEDIATION_LAYER:PHYSICS
    2. 

  2. !
    3. 

  3. MODULE module_surface_driver
    4. 

  4. CONTAINS
    5. 

  5. 

  6.    SUBROUTINE surface_driver(                                         &
    7. 

  7.      &           acgrdflx,achfx,aclhf                                 &
    8. 

  8.      &          ,acsnom,acsnow,akhs,akms,albedo,br,canwat             &
    9. 

  9.      &          ,chklowq,dt,dx,dz8w,dzs,glw                           &
    10. 

  10.      &          ,grdflx,gsw,swdown,gz1oz0,hfx,ht,ifsnow,isfflx        &
    11. 

  11.      &          ,fractional_seaice,seaice_albedo_opt,tice2tsk_if2cold &
    12. 

  12.      &          ,isltyp,itimestep,julian_in,ivgtyp,lowlyr,mavail,rmol &
    13. 

  13.      &          ,num_soil_layers,p8w,pblh,pi_phy,pshltr,psih          &
    14. 

  14. #if (NMM_CORE==1)
    15. 

  15.      &          ,psim,p_phy,q10,q2,qfx,taux,tauy,qsfc,qshltr,qz0      &
    16. 

  16. #else
    17. 

  17.      &          ,psim,p_phy,q10,q2,qfx,qsfc,qshltr,qz0                &
    18. 

  18. #endif
    19. 

  19.      &          ,raincv,rho,sfcevp,sfcexc,sfcrunoff                   &
    20. 

  20.      &          ,smois,smstav,smstot,snoalb,snow,snowc,snowh,stepbl   &
    21. 

  21.      &          ,smcrel                                               &
    22. 

  22.      &          ,th10,th2,thz0,th_phy,tmn,tshltr,tsk,tslb             &
    23. 

  23.      &          ,tyr,tyra,tdly,tlag,lagday,nyear,nday,tmn_update,yr   &
    24. 

  24.      &          ,t_phy,u10,udrunoff,ust,uz0,u_frame,u_phy,v10,vegfra  &
    25. 

  25.      &          ,vz0,v_frame,v_phy,warm_rain,wspd,xice,xland,z,znt,zs &
    26. 

  26. #if (NMM_CORE==1)
    27. 

  27.      &          ,xicem,isice,iswater,ct,tke_pbl,sfenth                &
    28. 

  28. #else
    29. 

  29.      &          ,xicem,isice,iswater,ct,tke_pbl                       &
    30. 

  30. #endif
    31. 

  31.      &          ,albbck,embck,lh,sh2o,shdmax,shdmin,z0                &
    32. 

  32.      &          ,flqc,flhc,psfc,sst,sstsk,dtw,sst_update,sst_skin     &
    33. 

  33.      &          ,scm_force_skintemp,scm_force_flux,t2,emiss           &
    34. 

  34.      &          ,sf_sfclay_physics,sf_surface_physics,ra_lw_physics   & 
    35. 

  35.      &          ,mosaic_lu,mosaic_soil                                &
    36. 

  36.      &          ,landusef,soilctop,soilcbot,ra,rs,nlcat,nscat,vegf_px & ! PX-LSM
    37. 

  37.      &          ,snowncv, anal_interval, lai, pxlsm_smois_init        & ! PX-LSM
    38. 

  38.      &          ,pxlsm_soil_nudge                                     & ! PX-LSM
    39. 

  39.      &          ,idveg, iopt_crs, iopt_btr, iopt_run, iopt_sfc, iopt_frz, iopt_inf     &
    40. 

  40.      &          ,iopt_rad ,iopt_alb ,iopt_snf ,iopt_tbot, iopt_stc                     & 
    41. 

  41.      &          ,isnowxy  ,tvxy     ,tgxy     ,canicexy                                &
    42. 

  42.      &          ,canliqxy ,eahxy    ,tahxy    ,cmxy     ,chxy                          &
    43. 

  43.      &          ,fwetxy   ,sneqvoxy ,alboldxy ,qsnowxy  ,wslakexy ,zwtxy    ,waxy      &
    44. 

  44.      &          ,wtxy     ,tsnoxy   ,zsnsoxy  ,snicexy  ,snliqxy  ,lfmassxy ,rtmassxy  &
    45. 

  45.      &          ,stmassxy ,woodxy   ,stblcpxy ,fastcpxy ,xsaixy                        &
    46. 

  46.      &          ,tradxy   ,tsxy     ,neexy    ,gppxy    ,nppxy    ,fvegxy   ,qinxy     &
    47. 

  47.      &          ,runsfxy  ,runsbxy  ,ecanxy   ,edirxy   ,etranxy  ,fsaxy    ,firaxy    &
    48. 

  48.      &          ,aparxy   ,psnxy    ,savxy    ,sagxy                                   &
    49. 

  49.      &          ,fsnoxy                                                                &
    50. 

  50.      &          ,t2mvxy   ,t2mbxy   ,chstarxy ,rssunxy  ,rsshaxy  ,bgapxy   ,wgapxy    &    
    51. 

  51.      &          ,gapxy    ,tgvxy    ,tgbxy    ,q2mvxy   ,q2mbxy   ,chvxy    ,chbxy     &
    52. 

  52. #if ( EM_CORE==1)
    53. 

  53.      &          ,ch,tsq,qsq,cov                                       & ! MYNN
    54. 

  54. #endif
    55. 

  55.             !  Optional urban
    56. 

  56.      &          ,slope_rad,topo_shading,shadowmask                    & !I solar
    57. 

  57.      &          ,swnorm,slope,slp_azi                                 & !I solar
    58. 

  58.      &          ,declin,solcon,coszen,hrang,xlat_urb2d                & !I solar/urban
    59. 

  59.      &          ,num_roof_layers, num_wall_layers                     & !I urban
    60. 

  60.      &          ,num_road_layers, dzr, dzb, dzg                       & !I urban
    61. 

  61.      &          ,tr_urb2d,tb_urb2d,tg_urb2d,tc_urb2d,qc_urb2d         & !H urban
    62. 

  62.      &          ,uc_urb2d                                             & !H urban
    63. 

  63.      &          ,xxxr_urb2d,xxxb_urb2d,xxxg_urb2d,xxxc_urb2d          & !H urban
    64. 

  64.      &          ,trl_urb3d,tbl_urb3d,tgl_urb3d                        & !H urban
    65. 

  65.      &          ,sh_urb2d,lh_urb2d,g_urb2d,rn_urb2d,ts_urb2d          & !H urban
    66. 

  66.      &          ,frc_urb2d, utype_urb2d                               & !H urban
    67. 

  67.      &          ,cmr_sfcdif,chr_sfcdif,cmc_sfcdif,chc_sfcdif          &
    68. 

  68. !-----SSiB LSM (fds 06/2010)---------------------------------------------------
    69. 

  69.      &          ,alswvisdir, alswvisdif, alswnirdir, alswnirdif       & ! ssib
    70. 

  70.      &          ,swvisdir, swvisdif, swnirdir, swnirdif               & ! ssib
    71. 

  71.      &          ,ssib_br  ,ssib_fm  ,ssib_fh  ,ssib_cm  ,ssibxdd      & ! ssib
    72. 

  72.      &          ,ssib_lhf ,ssib_shf ,ssib_ghf ,ssib_egs ,ssib_eci     & ! ssib
    73. 

  73.      &          ,ssib_ect ,ssib_egi ,ssib_egt ,ssib_sdn ,ssib_sup     & ! ssib
    74. 

  74.      &          ,ssib_ldn ,ssib_lup ,ssib_wat ,ssib_shc ,ssib_shg     & ! ssib
    75. 

  75.      &          ,ssib_lai ,ssib_vcf ,ssib_z00 ,ssib_veg               & ! ssib
    76. 

  76.      &          ,ISNOW ,SWE ,SNOWDEN ,SNOWDEPTH ,TKAIR                & ! ssib-snow
    77. 

  77.      &          ,DZO1 ,WO1 ,TSSN1 ,TSSNO1 ,BWO1 ,BTO1                 & ! ssib-snow
    78. 

  78.      &          ,CTO1 ,FIO1 ,FLO1 ,BIO1 ,BLO1 ,HO1                    & ! ssib-snow
    79. 

  79.      &          ,DZO2 ,WO2 ,TSSN2 ,TSSNO2 ,BWO2 ,BTO2                 & ! ssib-snow
    80. 

  80.      &          ,CTO2 ,FIO2 ,FLO2 ,BIO2 ,BLO2 ,HO2                    & ! ssib-snow
    81. 

  81.      &          ,DZO3 ,WO3 ,TSSN3 ,TSSNO3 ,BWO3 ,BTO3                 & ! ssib-snow
    82. 

  82.      &          ,CTO3 ,FIO3 ,FLO3 ,BIO3 ,BLO3 ,HO3                    & ! ssib-snow
    83. 

  83.      &          ,DZO4 ,WO4 ,TSSN4 ,TSSNO4 ,BWO4 ,BTO4                 & ! ssib-snow
    84. 

  84.      &          ,CTO4 ,FIO4 ,FLO4 ,BIO4 ,BLO4 ,HO4                    & ! ssib-snow
    85. 

  85.      &          ,ra_sw_physics                                        & ! ssib 
    86. 

  86. !------------------------------------------------------------------------------
    87. 

  87.      &          , ids,ide,jds,jde,kds,kde                             &
    88. 

  88.      &          , ims,ime,jms,jme,kms,kme                             &
    89. 

  89.      &          , i_start,i_end,j_start,j_end,kts,kte,num_tiles       &
    90. 

  90.              !  Optional moisture tracers
    91. 

  91.      &           ,qv_curr, qc_curr, qr_curr                           &
    92. 

  92.      &           ,qi_curr, qs_curr, qg_curr                           &
    93. 

  93.              !  Optional moisture tracer flags
    94. 

  94.      &           ,f_qv,f_qc,f_qr                                      &
    95. 

  95.      &           ,f_qi,f_qs,f_qg                                      &
    96. 

  96.              !  Other optionals (more or less em specific)
    97. 

  97.      &          ,capg,hol,mol                                         &
    98. 

  98.      &          ,rainncv,rainshv,rainbl,regime,thc                    &
    99. 

  99.      &          ,qsg,qvg,qcg,soilt1,tsnav                             &
    100. 

  100.      &          ,smfr3d,keepfr3dflag,dew                              &
    101. 

  101.              !  Other optionals (more or less nmm specific)
    102. 

  102.      &          ,potevp,snopcx,soiltb,sr                              &
    103. 

  103.              !  Optional observation PX LSM surface nudging
    104. 

  104.      &          ,t2_ndg_old, q2_ndg_old, t2_ndg_new, q2_ndg_new       &
    105. 

  105.      &          ,sn_ndg_old, sn_ndg_new                               &
    106. 

  106.      &          ,t2obs, q2obs                                         &
    107. 

  107.              ! OPTIONAL, Required by TEMF surface layer 1/7/09 WA
    108. 

  108.      &          ,hd_temf,te_temf,fCor,exch_temf,wm_temf               &
    109. 

  109.              ! Required by ideal SCM surface layer 1/6/10 WA
    110. 

  110.      &          ,hfx_force,lh_force,tsk_force                       &
    111. 

  111.      &          ,hfx_force_tend,lh_force_tend,tsk_force_tend          &
    112. 

  112.              !  Optional observation nudging
    113. 

  113.      &          ,uratx,vratx,tratx                                    &
    114. 

  114.              !  Optional simple oml model
    115. 

  115.      &          ,omlcall,oml_hml0,oml_gamma                           &
    116. 

  116.      &          ,tml,t0ml,hml,h0ml,huml,hvml,f,tmoml                  &
    117. 

  117.      &          ,ustm,ck,cka,cd,cda,isftcflx,iz0tlnd                  &
    118. 

  118.      &         ,isurban, mminlu                                       &
    119. 

  119.      &          ,snotime                                              &
    120. 

  120.      &           ,rdlai2d                                             &
    121. 

  121.      &          ,usemonalb                                            &
    122. 

  122.      &          ,noahres                                              &
    123. 

  123.              !  Optional adaptive time step
    124. 

  124.      &          ,bldt,curr_secs,adapt_step_flag,bldtacttime           & 
    125. 

  125.          ! Optional urban with BEP
    126. 

  126.      &          ,sf_urban_physics,gmt,xlat,xlong,julday               &
    127. 

  127.      &          ,num_urban_layers                                     & !multi-layer urban
    128. 

  128.      &          ,trb_urb4d,tw1_urb4d,tw2_urb4d,tgb_urb4d              & !multi-layer urban
    129. 

  129.      &          ,tlev_urb3d,qlev_urb3d                                & !multi-layer urban
    130. 

  130.      &          ,tw1lev_urb3d,tw2lev_urb3d                            & !multi-layer urban
    131. 

  131.      &          ,tglev_urb3d,tflev_urb3d                              & !multi-layer urban
    132. 

  132.      &          ,sf_ac_urb3d,lf_ac_urb3d,cm_ac_urb3d                  & !multi-layer urban
    133. 

  133.      &          ,sfvent_urb3d,lfvent_urb3d                            & !multi-layer urban 
    134. 

  134.      &          ,sfwin1_urb3d,sfwin2_urb3d                            & !multi-layer urban       
    135. 

  135.      &          ,sfw1_urb3d,sfw2_urb3d,sfr_urb3d,sfg_urb3d            & !multi-layer urban
    136. 

  136.      &          ,a_u_bep,a_v_bep,a_t_bep,a_q_bep                      &
    137. 

  137.      &          ,b_u_bep,b_v_bep,b_t_bep,b_q_bep                      &
    138. 

  138.      &          ,sf_bep,vl_bep                                        &
    139. 

  139.      &          ,a_e_bep,b_e_bep,dlg_bep                              &
    140. 

  140.      &          ,dl_u_bep                                             &                          
    141. 

  141.      &          ,cldfra                                               & !ssib
    142. 

  142.          ! Optional urban Bep end
    143. 

  143.      &                                                             )
    144. 

  144.               
    145. 

  145. #if ( ! NMM_CORE == 1 )
    146. 

  146.    USE module_state_description, ONLY : SFCLAYSCHEME              &
    147. 

  147.                                        ,SFCLAYREVSCHEME           &
    148. 

  148.                                        ,MYJSFCSCHEME              &
    149. 

  149.                                        ,QNSESFCSCHEME             &
    150. 

  150.                                        ,GFSSFCSCHEME              &
    151. 

  151.                                        ,PXSFCSCHEME               &
    152. 

  152.                                        ,NOAHMPSCHEME              &
    153. 

  153.                                        ,TEMFSFCSCHEME             &
    154. 

  154.                                        ,IDEALSCMSFCSCHEME         &
    155. 

  155.                                        ,SLABSCHEME                &
    156. 

  156.                                        ,LSMSCHEME                 &
    157. 

  157.                                        ,RUCLSMSCHEME              &
    158. 

  158.                                        ,PXLSMSCHEME               &
    159. 

  159.                                        ,SSIBSCHEME                & !ssib
    160. 

  160.                                        ,MYNNSFCSCHEME             
    161. 

  161. #else
    162. 

  162.    USE module_state_description, ONLY : SFCLAYSCHEME              &
    163. 

  163.                                        ,SFCLAYREVSCHEME           &
    164. 

  164.                                        ,MYJSFCSCHEME              &
    165. 

  165.                                        ,QNSESFCSCHEME             &
    166. 

  166.                                        ,GFSSFCSCHEME              &
    167. 

  167.                                        ,PXSFCSCHEME               &
    168. 

  168.                                        ,NOAHMPSCHEME              &
    169. 

  169.                                        ,SLABSCHEME                &
    170. 

  170.                                        ,LSMSCHEME                 &
    171. 

  171.                                        ,RUCLSMSCHEME              &
    172. 

  172.                                        ,PXLSMSCHEME               &
    173. 

  173.                                        ,TEMFSFCSCHEME             &
    174. 

  174.                                        ,GFDLSFCSCHEME             &
    175. 

  175.                                        ,SSIBSCHEME                & ! ssib
    176. 

  176.                                        ,GFDLSLAB 
    177. 

  177. 

  178. 

  179. #endif
    180. 

  180.    USE module_model_constants
    181. 

  181. ! *** add new modules of schemes here
    182. 

  182. 

  183.    USE module_sf_sfclay
    184. 

  184.    USE module_sf_myjsfc
    185. 

  185.    USE module_sf_qnsesfc
    186. 

  186.    USE module_sf_gfs
    187. 

  187.    USE module_sf_noahdrv, only : lsm
    188. 

  188.    USE module_sf_noahmpdrv, only : noahmplsm
    189. 

  189.    USE module_sf_noah_seaice_drv
    190. 

  190.    USE module_sf_ssib  ! ssib
    191. 

  191.    USE module_sf_ruclsm
    192. 

  192.    USE module_sf_pxsfclay
    193. 

  193.    USE module_sf_pxlsm
    194. 

  194.    USE module_sf_temfsfclay
    195. 

  195.    USE module_sf_sfclayrev
    196. 

  196.    USE module_sf_noah_seaice_drv
    197. 

  197. #if ( EM_CORE==1)
    198. 

  198.    USE module_sf_mynn
    199. 

  199.    USE module_sf_oml
    200. 

  200.    USE module_sf_idealscmsfclay
    201. 

  201. #endif
    202. 

  202.    USE module_sf_scmflux       
    203. 

  203.    USE module_sf_scmskintemp
    204. 

  204. 

  205. #if ( NMM_CORE == 1 )
    206. 

  206.    USE module_sf_gfdl
    207. 

  207. #endif
    208. 

  208. 

  209.    USE module_sf_slab
    210. 

  210. !
    211. 

  211.    USE module_sf_sfcdiags
    212. 

  212.    USE module_sf_sfcdiags_ruclsm
    213. 

  213.    USE module_sf_sstskin
    214. 

  214.    USE module_sf_tmnupdate
    215. 

  215. !
    216. 

  216. !  This driver calls subroutines for the surface parameterizations.
    217. 

  217. !
    218. 

  218. !  surface layer: (between surface and pbl)
    219. 

  219. !      1. sfclay
    220. 

  220. !      2. myjsfc
    221. 

  221. !      7. Pleim surface layer
    222. 

  222. !      5. MYNN surface layer
    223. 

  223. !  surface: ground temp/lsm scheme:
    224. 

  224. !      1. slab
    225. 

  225. !      2. Noah LSM
    226. 

  226. !      7. Pleim-Xiu LSM
    227. 

  227. !     11. Revised sfclay (option 1) 
    228. 

  228. !
    229. 

  229. !  surface: ground temp/lsm scheme for urban:
    230. 

  230. !      2.  BEP
    231. 

  231. !
    232. 

  232. !  ocean mixed layer model
    233. 

  233. !      omlcall = 1
    234. 

  234. !------------------------------------------------------------------
    235. 

  235.    IMPLICIT NONE
    236. 

  236. !======================================================================
    237. 

  237. ! Grid structure in physics part of WRF
    238. 

  238. !----------------------------------------------------------------------
    239. 

  239. ! The horizontal velocities used in the physics are unstaggered
    240. 

  240. ! relative to temperature/moisture variables. All predicted
    241. 

  241. ! variables are carried at half levels except w, which is at full
    242. 

  242. ! levels. Some arrays with names (*8w) are at w (full) levels.
    243. 

  243. !
    244. 

  244. !----------------------------------------------------------------------
    245. 

  245. ! In WRF, kms (smallest number) is the bottom level and kme (largest
    246. 

  246. ! number) is the top level.  In your scheme, if 1 is at the top level,
    247. 

  247. ! then you have to reverse the order in the k direction.
    248. 

  248. !
    249. 

  249. !         kme      -   half level (no data at this level)
    250. 

  250. !         kme    ----- full level
    251. 

  251. !         kme-1    -   half level
    252. 

  252. !         kme-1  ----- full level
    253. 

  253. !         .
    254. 

  254. !         kms+2    -   half level
    255. 

  255. !         kms+2  ----- full level
    256. 

  256. !         kms+1    -   half level
    257. 

  257. !         kms+1  ----- full level
    258. 

  258. !         kms      -   half level
    259. 

  259. !         kms    ----- full level
    260. 

  260. !
    261. 

  261. !======================================================================
    262. 

  262. ! Definitions
    263. 

  263. !-----------
    264. 

  264. ! Theta      potential temperature (K)
    265. 

  265. ! Qv         water vapor mixing ratio (kg/kg)
    266. 

  266. ! Qc         cloud water mixing ratio (kg/kg)
    267. 

  267. ! Qr         rain water mixing ratio (kg/kg)
    268. 

  268. ! Qi         cloud ice mixing ratio (kg/kg)
    269. 

  269. ! Qs         snow mixing ratio (kg/kg)
    270. 

  270. !-----------------------------------------------------------------
    271. 

  271. !-- itimestep     number of time steps
    272. 

  272. !-- GLW           downward long wave flux at ground surface (W/m^2)
    273. 

  273. !-- GSW           net short wave flux at ground surface (W/m^2)
    274. 

  274. !-- SWDOWN        downward short wave flux at ground surface (W/m^2)
    275. 

  275. !-- EMISS         surface emissivity (between 0 and 1)
    276. 

  276. !-- TSK           surface temperature (K)
    277. 

  277. !-- TMN           soil temperature at lower boundary (K)
    278. 

  278. !-- TYR           annual mean surface temperature of previous year (K)
    279. 

  279. !-- TYRA          accumulated surface temperature in the current year (K)
    280. 

  280. !-- TLAG          mean surface temperature of previous 140 days (K)
    281. 

  281. !-- TDLY          accumulated daily mean surface temperature of the current day (K)
    282. 

  282. !-- XLAND         land mask (1 for land, 2 for water)
    283. 

  283. !-- ZNT           time-varying roughness length (m)
    284. 

  284. !-- Z0            background roughness length (m)
    285. 

  285. !-- MAVAIL        surface moisture availability (between 0 and 1)
    286. 

  286. !-- UST           u* in similarity theory (m/s)
    287. 

  287. !-- MOL           T* (similarity theory) (K)
    288. 

  288. !-- HOL           PBL height over Monin-Obukhov length
    289. 

  289. !-- PBLH          PBL height (m)
    290. 

  290. !-- CAPG          heat capacity for soil (J/K/m^3)
    291. 

  291. !-- THC           thermal inertia (Cal/cm/K/s^0.5)
    292. 

  292. !-- SNOWC         flag indicating snow coverage (1 for snow cover)
    293. 

  293. !-- HFX           net upward heat flux at the surface (W/m^2)
    294. 

  294. !-- QFX           net upward moisture flux at the surface (kg/m^2/s)
    295. 

  295. !-- TAUX          RHO*U**2 for ocean coupling
    296. 

  296. !-- TAUY          RHO*U**2 for ocean coupling
    297. 

  297. !-- LH            net upward latent heat flux at surface (W/m^2)
    298. 

  298. !-- REGIME        flag indicating PBL regime (stable, unstable, etc.)
    299. 

  299. !-- tke_pbl       turbulence kinetic energy from PBL schemes (m^2/s^2)
    300. 

  300. !-- akhs          sfc exchange coefficient of heat/moisture from MYJ
    301. 

  301. !-- akms          sfc exchange coefficient of momentum from MYJ
    302. 

  302. !-- thz0          potential temperature at roughness length (K)
    303. 

  303. !-- uz0           u wind component at roughness length (m/s)
    304. 

  304. !-- vz0           v wind component at roughness length (m/s)
    305. 

  305. !-- qsfc          specific humidity at lower boundary (kg/kg)
    306. 

  306. !-- uratx         ratio of u over u10 (Added for obs-nudging)
    307. 

  307. !-- vratx         ratio of v over v10 (Added for obs-nudging)
    308. 

  308. !-- tratx         ratio of t over th2 (Added for obs-nudging)
    309. 

  309. !-- u10           diagnostic 10-m u component from surface layer
    310. 

  310. !-- v10           diagnostic 10-m v component from surface layer
    311. 

  311. !-- th2           diagnostic 2-m theta from surface layer and lsm
    312. 

  312. !-- t2            diagnostic 2-m temperature from surface layer and lsm
    313. 

  313. !-- q2            diagnostic 2-m mixing ratio from surface layer and lsm
    314. 

  314. !-- tshltr        diagnostic 2-m theta from MYJ
    315. 

  315. !-- th10          diagnostic 10-m theta from MYJ
    316. 

  316. !-- qshltr        diagnostic 2-m specific humidity from MYJ
    317. 

  317. !-- q10           diagnostic 10-m specific humidity from MYJ
    318. 

  318. !-- lowlyr        index of lowest model layer above ground
    319. 

  319. !-- rr            dry air density (kg/m^3)
    320. 

  320. !-- u_phy         u-velocity interpolated to theta points (m/s)
    321. 

  321. !-- v_phy         v-velocity interpolated to theta points (m/s)
    322. 

  322. !-- th_phy        potential temperature (K)
    323. 

  323. !-- moist         moisture array (4D - last index is species) (kg/kg)
    324. 

  324. !-- p_phy         pressure (Pa)
    325. 

  325. !-- pi_phy        exner function (dimensionless)
    326. 

  326. !-- pshltr        diagnostic shelter (2m) pressure from MYJ (Pa)
    327. 

  327. !-- p8w           pressure at full levels (Pa)
    328. 

  328. !-- t_phy         temperature (K)
    329. 

  329. !-- dz8w          dz between full levels (m)
    330. 

  330. !-- z             height above sea level (m)
    331. 

  331. !-- DX            horizontal space interval (m)
    332. 

  332. !-- DT            time step (second)
    333. 

  333. !-- PSFC          pressure at the surface (Pa)
    334. 

  334. !-- SST           sea-surface temperature (K)
    335. 

  335. !-- SSTSK         skin sea-surface temperature (K)
    336. 

  336. !-- DTW           warm layer temp diff (K)
    337. 

  337. !-- TSLB
    338. 

  338. !-- ZS
    339. 

  339. !-- DZS
    340. 

  340. !-- num_soil_layers number of soil layer
    341. 

  341. !-- IFSNOW      ifsnow=1 for snow-cover effects
    342. 

  342. !-- omlcall       whether to call simple ocean mixed layer model from slab (1 = use oml)
    343. 

  343. !-- oml_hml0      initial mixed layer depth (if real-data not available, default 50 m)
    344. 

  344. !-- oml_gamma     lapse rate below mixed layer in ocean (default 0.14 K m-1)
    345. 

  345. !-- ck            enthalpy exchange coeff at 10 meters
    346. 

  346. !-- cd            momentum exchange coeff at 10 meters
    347. 

  347. !-- cka           enthalpy exchange coeff at the lowest model level
    348. 

  348. !-- cda           momentum exchange coeff at the lowest model level
    349. 

  349. !!!!!!!!!!!!!!
    350. 

  350. !
    351. 

  351. !
    352. 

  352. !-- LANDUSEF     Landuse fraction                      ! P-X LSM
    353. 

  353. !-- SOILCTOP     Top soil fraction                     ! P-X LSM
    354. 

  354. !-- SOILCBOT     Bottom soil fraction                  ! P-X LSM
    355. 

  355. !-- RA           Aerodynamic resistence                        ! P-X LSM
    356. 

  356. !-- RS           Stomatal resistence                   ! P-X LSM
    357. 

  357. !-- NLCAT        Number of landuse categories          ! P-X LSM
    358. 

  358. !-- NSCAT        Number of soil categories             ! P-X LSM
    359. 

  359. !-- ch - drag coefficient for heat/moisture            ! MYNN LSM
    360. 

  360. 

  361. !
    362. 

  362. !-- ids           start index for i in domain
    363. 

  363. !-- ide           end index for i in domain
    364. 

  364. !-- jds           start index for j in domain
    365. 

  365. !-- jde           end index for j in domain
    366. 

  366. !-- kds           start index for k in domain
    367. 

  367. !-- kde           end index for k in domain
    368. 

  368. !-- ims           start index for i in memory
    369. 

  369. !-- ime           end index for i in memory
    370. 

  370. !-- jms           start index for j in memory
    371. 

  371. !-- jme           end index for j in memory
    372. 

  372. !-- kms           start index for k in memory
    373. 

  373. !-- kme           end index for k in memory
    374. 

  374. !-- its           start index for i in tile
    375. 

  375. !-- ite           end index for i in tile
    376. 

  376. !-- jts           start index for j in tile
    377. 

  377. !-- jte           end index for j in tile
    378. 

  378. !-- kts           start index for k in tile
    379. 

  379. !-- kte           end index for k in tile
    380. 

  380. !
    381. 

  381. !******************************************************************
    382. 

  382. !------------------------------------------------------------------
    383. 

  383. 

  384.    INTEGER, INTENT(IN) ::                                             &
    385. 

  385.      &           ids,ide,jds,jde,kds,kde                              &
    386. 

  386.      &          ,ims,ime,jms,jme,kms,kme                              &
    387. 

  387.      &          ,kts,kte,num_tiles
    388. 

  388. 

  389.    INTEGER, INTENT(IN)::   FRACTIONAL_SEAICE
    390. 

  390.    INTEGER, INTENT(IN)::   SEAICE_ALBEDO_OPT
    391. 

  391. 

  392.    INTEGER, INTENT(IN)::   NLCAT, mosaic_lu, mosaic_soil
    393. 

  393.    INTEGER, INTENT(IN)::   NSCAT
    394. 

  394. 

  395.    INTEGER, INTENT(IN) :: sf_sfclay_physics, sf_surface_physics,      &
    396. 

  396.                           sf_urban_physics,ra_lw_physics,sst_update,  &
    397. 

  397.                           ra_sw_physics                    
    398. 

  398.    INTEGER, INTENT(IN),OPTIONAL :: sst_skin, tmn_update,              &
    399. 

  399.                                    scm_force_skintemp, scm_force_flux 
    400. 

  400. 

  401.    INTEGER, DIMENSION(num_tiles), INTENT(IN) ::                       &
    402. 

  402.      &           i_start,i_end,j_start,j_end
    403. 

  403. 

  404.    INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT )::  ISLTYP
    405. 

  405.    INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   IVGTYP
    406. 

  406.    INTEGER, DIMENSION( ims:ime , jms:jme ), INTENT(IN )::   LOWLYR
    407. 

  407.    INTEGER, INTENT(IN )::   IFSNOW
    408. 

  408.    INTEGER, INTENT(IN )::   ISFFLX
    409. 

  409.    INTEGER, INTENT(IN )::   ITIMESTEP
    410. 

  410.    INTEGER, INTENT(IN )::   NUM_SOIL_LAYERS
    411. 

  411.    REAL,    INTENT(IN ),OPTIONAL ::   JULIAN_in
    412. 

  412.    INTEGER, INTENT(IN )::   LAGDAY
    413. 

  413.    INTEGER, INTENT(IN )::   STEPBL
    414. 

  414.    INTEGER, INTENT(IN )::   ISICE
    415. 

  415.    INTEGER, INTENT(IN )::   ISWATER
    416. 

  416.    INTEGER, INTENT(IN ), OPTIONAL :: ISURBAN
    417. 

  417.    CHARACTER(LEN=*), INTENT(IN ), OPTIONAL :: MMINLU
    418. 

  418.    LOGICAL, INTENT(IN )::   WARM_RAIN
    419. 

  419.    LOGICAL, INTENT(IN)::   tice2tsk_if2cold
    420. 

  420.    INTEGER, INTENT(INOUT ),OPTIONAL ::   NYEAR
    421. 

  421.    REAL   , INTENT(INOUT ),OPTIONAL ::   NDAY
    422. 

  422.    INTEGER, INTENT(IN ),OPTIONAL ::   YR
    423. 

  423.    REAL , INTENT(IN )::   U_FRAME
    424. 

  424.    REAL , INTENT(IN )::   V_FRAME
    425. 

  425. #if (NMM_CORE==1)
    426. 

  426.    real , intent(IN )::   SFENTH
    427. 

  427. #endif
    428. 

  428.    REAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), INTENT(INOUT)::   SMOIS
    429. 

  429.    REAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), INTENT(INOUT)::   TSLB
    430. 

  430.    REAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), INTENT(OUT)  ::   SMCREL
    431. 

  431.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN )::   GLW
    432. 

  432.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   GSW,SWDOWN
    433. 

  433.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN )::   HT
    434. 

  434.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN )::   RAINCV
    435. 

  435.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   SST
    436. 

  436.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT ),OPTIONAL ::   SSTSK
    437. 

  437.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT ),OPTIONAL ::   DTW
    438. 

  438.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   TMN
    439. 

  439.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT ),OPTIONAL ::   TYR
    440. 

  440.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT ),OPTIONAL ::   TYRA
    441. 

  441.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT ),OPTIONAL ::   TDLY
    442. 

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

  443.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   VEGFRA
    444. 

  444. !------fds (06/2010)--------------------------
    445. 

  445.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   XICE
    446. 

  446. !---------------------------------------------
    447. 

  447.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   XLAND
    448. 

  448.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   XICEM
    449. 

  449.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   MAVAIL
    450. 

  450.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT)::   SNOALB
    451. 

  451.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   ACSNOW
    452. 

  452.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   SNOTIME
    453. 

  453.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   AKHS
    454. 

  454.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   AKMS
    455. 

  455.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   ALBEDO
    456. 

  456.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   CANWAT
    457. 

  457. 

  458.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   GRDFLX
    459. 

  459.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   HFX
    460. 

  460.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   RMOL
    461. 

  461.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   PBLH
    462. 

  462.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   Q2
    463. 

  463.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   QFX
    464. 

  464. #if (NMM_CORE==1)
    465. 

  465.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(OUT):: TAUX
    466. 

  466.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(OUT):: TAUY
    467. 

  467. #endif
    468. 

  468.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   QSFC
    469. 

  469.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   QZ0
    470. 

  470.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   SFCRUNOFF
    471. 

  471.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   SMSTAV
    472. 

  472.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   SMSTOT
    473. 

  473.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   SNOW
    474. 

  474.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   SNOWC
    475. 

  475.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   SNOWH
    476. 

  476.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   TH2
    477. 

  477.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   THZ0
    478. 

  478.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   TSK
    479. 

  479.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   UDRUNOFF
    480. 

  480.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   UST
    481. 

  481.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   UZ0
    482. 

  482.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   VZ0
    483. 

  483.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   WSPD
    484. 

  484.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT)::   ZNT
    485. 

  485. !-----fds (06/2010)---------------------------------------------
    486. 

  486.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_LHF ! SSiB output
    487. 

  487.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_SHF ! SSiB output
    488. 

  488.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_GHF ! SSiB output
    489. 

  489.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_EGS ! SSiB output
    490. 

  490.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_ECI ! SSiB output
    491. 

  491.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_ECT ! SSiB output
    492. 

  492.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_EGI ! SSiB output
    493. 

  493.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_EGT ! SSiB output
    494. 

  494.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_SDN ! SSiB output
    495. 

  495.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_SUP ! SSiB output
    496. 

  496.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_LDN ! SSiB output
    497. 

  497.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_LUP ! SSiB output
    498. 

  498.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_WAT ! SSiB output
    499. 

  499.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_SHC ! SSiB output
    500. 

  500.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_SHG ! SSiB output
    501. 

  501.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_LAI ! SSiB output
    502. 

  502.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_VCF ! SSiB output
    503. 

  503.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_Z00 ! SSiB output
    504. 

  504.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   SSIB_VEG ! SSiB output
    505. 

  505.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   ALSWVISDIR! SSiB
    506. 

  506.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   ALSWVISDIF! SSiB
    507. 

  507.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   ALSWNIRDIR! SSiB
    508. 

  508.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(OUT)::   ALSWNIRDIF! SSiB
    509. 

  509.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(IN)::      SWVISDIR! SSiB
    510. 

  510.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(IN)::      SWVISDIF! SSiB
    511. 

  511.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(IN)::      SWNIRDIR! SSiB
    512. 

  512.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(IN)::      SWNIRDIF! SSiB
    513. 

  513.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   SSiB_BR ! SSiB
    514. 

  514.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   SSiB_FM ! SSiB
    515. 

  515.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   SSiB_FH ! SSiB
    516. 

  516.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   SSiB_CM ! SSiB
    517. 

  517.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   SSiBXDD ! SSiB
    518. 

  518.    INTEGER, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT):: ISNOW  ! ssib-snow
    519. 

  519.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   SWE     ! ssib-snow
    520. 

  520.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   SNOWDEN ! ssib-snow
    521. 

  521.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT):: SNOWDEPTH ! ssib-snow
    522. 

  522.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   TKAIR   ! ssib-snow
    523. 

  523.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   DZO1    ! ssib-snow
    524. 

  524.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   WO1     ! ssib-snow
    525. 

  525.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   TSSN1   ! ssib-snow
    526. 

  526.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   TSSNO1  ! ssib-snow
    527. 

  527.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BWO1    ! ssib-snow
    528. 

  528.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BTO1    ! ssib-snow
    529. 

  529.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   CTO1    ! ssib-snow
    530. 

  530.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   FIO1    ! ssib-snow
    531. 

  531.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   FLO1    ! ssib-snow
    532. 

  532.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BIO1    ! ssib-snow
    533. 

  533.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BLO1    ! ssib-snow
    534. 

  534.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   HO1     ! ssib-snow
    535. 

  535.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   DZO2    ! ssib-snow
    536. 

  536.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   WO2     ! ssib-snow
    537. 

  537.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   TSSN2   ! ssib-snow
    538. 

  538.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   TSSNO2  ! ssib-snow
    539. 

  539.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BWO2    ! ssib-snow
    540. 

  540.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BTO2    ! ssib-snow
    541. 

  541.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   CTO2    ! ssib-snow
    542. 

  542.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   FIO2    ! ssib-snow
    543. 

  543.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   FLO2    ! ssib-snow
    544. 

  544.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BIO2    ! ssib-snow
    545. 

  545.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BLO2    ! ssib-snow
    546. 

  546.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   HO2     ! ssib-snow
    547. 

  547.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   DZO3    ! ssib-snow
    548. 

  548.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   WO3     ! ssib-snow
    549. 

  549.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   TSSN3   ! ssib-snow
    550. 

  550.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   TSSNO3  ! ssib-snow
    551. 

  551.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BWO3    ! ssib-snow
    552. 

  552.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BTO3    ! ssib-snow
    553. 

  553.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   CTO3    ! ssib-snow
    554. 

  554.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   FIO3    ! ssib-snow
    555. 

  555.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   FLO3    ! ssib-snow
    556. 

  556.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BIO3    ! ssib-snow
    557. 

  557.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BLO3    ! ssib-snow
    558. 

  558.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   HO3     ! ssib-snow
    559. 

  559.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   DZO4    ! ssib-snow
    560. 

  560.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   WO4     ! ssib-snow
    561. 

  561.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   TSSN4   ! ssib-snow
    562. 

  562.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   TSSNO4  ! ssib-snow
    563. 

  563.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BWO4    ! ssib-snow
    564. 

  564.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BTO4    ! ssib-snow
    565. 

  565.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   CTO4    ! ssib-snow
    566. 

  566.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   FIO4    ! ssib-snow
    567. 

  567.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   FLO4    ! ssib-snow
    568. 

  568.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BIO4    ! ssib-snow
    569. 

  569.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   BLO4    ! ssib-snow
    570. 

  570.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   HO4     ! ssib-snow
    571. 

  571. !----------------------------------------------------------
    572. 

  572.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   BR
    573. 

  573.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   CHKLOWQ
    574. 

  574.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   GZ1OZ0
    575. 

  575.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   PSHLTR
    576. 

  576.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   PSIH
    577. 

  577.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   PSIM
    578. 

  578.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   Q10
    579. 

  579.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   QSHLTR
    580. 

  580.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   TH10
    581. 

  581.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   TSHLTR
    582. 

  582.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   U10
    583. 

  583.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   V10
    584. 

  584.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(OUT)::   PSFC
    585. 

  585.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)::   ACSNOM
    586. 

  586.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)::   SFCEVP
    587. 

  587.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT),OPTIONAL ::   ACHFX
    588. 

  588.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT),OPTIONAL ::   ACLHF
    589. 

  589.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(INOUT),OPTIONAL ::   ACGRDFLX
    590. 

  590.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)::   SFCEXC
    591. 

  591.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)::   FLHC
    592. 

  592.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)::   FLQC
    593. 

  593.    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) ::   CT
    594. 

  594.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   DZ8W
    595. 

  595.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   P8W
    596. 

  596.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   PI_PHY
    597. 

  597.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   P_PHY
    598. 

  598.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   RHO
    599. 

  599.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   TH_PHY
    600. 

  600.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   T_PHY
    601. 

  601.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   U_PHY
    602. 

  602.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   V_PHY
    603. 

  603.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN )::   Z
    604. 

  604. 

  605.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::   TKE_PBL
    606. 

  606.    REAL, DIMENSION(1:num_soil_layers), INTENT(IN)::   DZS
    607. 

  607.    REAL, DIMENSION(1:num_soil_layers), INTENT(IN)::   ZS
    608. 

  608.    REAL, INTENT(IN )::   DT
    609. 

  609.    REAL, INTENT(IN )::   DX
    610. 

  610.    REAL,       INTENT(IN   ),OPTIONAL    ::     bldt
    611. 

  611.    REAL,       INTENT(IN   ),OPTIONAL    ::     curr_secs
    612. 

  612.    LOGICAL,    INTENT(IN   ),OPTIONAL    ::     adapt_step_flag
    613. 

  613.    REAL,       INTENT(INOUT),OPTIONAL    ::     bldtacttime  
    614. 

  614. 

  615. !  arguments for NCAR surface physics
    616. 

  616. 

  617.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT )::   ALBBCK  ! INOUT needed for NMM
    618. 

  618.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT )::   EMBCK
    619. 

  619.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT )::   LH
    620. 

  620.    REAL, DIMENSION( ims:ime , 1:num_soil_layers, jms:jme ), INTENT(INOUT)::   SH2O
    621. 

  621.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN )::   SHDMAX
    622. 

  622.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN )::   SHDMIN
    623. 

  623.    REAL, DIMENSION( ims:ime , jms:jme ), INTENT(INOUT )::   Z0
    624. 

  624. 

  625.    INTEGER, OPTIONAL, INTENT(IN) :: idveg, iopt_crs, iopt_btr, iopt_run, iopt_sfc, iopt_frz, iopt_inf, iopt_rad, iopt_alb, iopt_snf, iopt_tbot, iopt_stc
    626. 

  626.    INTEGER, OPTIONAL, DIMENSION(ims:ime , jms:jme), INTENT(INOUT) :: ISNOWXY
    627. 

  627.    REAL, OPTIONAL, DIMENSION(ims:ime ,-2:num_soil_layers, jms:jme), INTENT(INOUT) :: zsnsoxy
    628. 

  628.    REAL, OPTIONAL, DIMENSION(ims:ime ,-2:0,               jms:jme), INTENT(INOUT) :: tsnoxy, snicexy, snliqxy
    629. 

  629.    REAL, OPTIONAL, DIMENSION(ims:ime , jms:jme), INTENT(INOUT) :: tvxy, tgxy, canicexy, canliqxy, eahxy, tahxy, cmxy, chxy,               &
    630. 

  630.         fwetxy, sneqvoxy, alboldxy, qsnowxy, wslakexy, zwtxy, waxy, wtxy, lfmassxy, rtmassxy, stmassxy, woodxy, stblcpxy, fastcpxy,       &
    631. 

  631.         xsaixy, tradxy, tsxy, neexy, gppxy, nppxy, fvegxy, qinxy, runsfxy, runsbxy, ecanxy, edirxy, etranxy, fsaxy, firaxy,               &
    632. 

  632.         aparxy, psnxy, savxy, sagxy, fsnoxy, q2mvxy, q2mbxy
    633. 

  633.    REAL, OPTIONAL, DIMENSION(ims:ime , jms:jme), INTENT(INOUT) ::  t2mvxy   ,t2mbxy   ,chstarxy, rssunxy, rsshaxy, bgapxy,wgapxy,gapxy , &
    634. 

  634.         tgvxy ,tgbxy, chvxy, chbxy
    635. 

  635. 

  636. ! Variables for multi-layer UCM
    637. 

  637.    REAL, OPTIONAL, INTENT(IN  )   ::                                   GMT 
    638. 

  638.    INTEGER, OPTIONAL, INTENT(IN  ) ::                               JULDAY
    639. 

  639.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN   )        ::XLAT, XLONG
    640. 

  640.    INTEGER, INTENT(IN )::   NUM_URBAN_LAYERS
    641. 

  641.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: trb_urb4d
    642. 

  642.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tw1_urb4d
    643. 

  643.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tw2_urb4d
    644. 

  644.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tgb_urb4d
    645. 

  645.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tlev_urb3d
    646. 

  646.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: qlev_urb3d
    647. 

  647.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tw1lev_urb3d
    648. 

  648.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tw2lev_urb3d
    649. 

  649.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tglev_urb3d
    650. 

  650.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: tflev_urb3d
    651. 

  651.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: lf_ac_urb3d
    652. 

  652.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: sf_ac_urb3d
    653. 

  653.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: cm_ac_urb3d
    654. 

  654.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: sfvent_urb3d
    655. 

  655.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: lfvent_urb3d
    656. 

  656.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfwin1_urb3d
    657. 

  657.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfwin2_urb3d
    658. 

  658.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfw1_urb3d
    659. 

  659.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfw2_urb3d
    660. 

  660.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfr_urb3d
    661. 

  661.    REAL, OPTIONAL, DIMENSION( ims:ime, 1:num_urban_layers, jms:jme ), INTENT(INOUT) :: sfg_urb3d
    662. 

  662.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::a_u_bep   !Implicit momemtum component X-direction
    663. 

  663.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::a_v_bep   !Implicit momemtum component Y-direction
    664. 

  664.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::a_t_bep   !Implicit component pot. temperature
    665. 

  665.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::a_e_bep   !Implicit component TKE
    666. 

  666.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::a_q_bep   !Implicit component TKE
    667. 

  667.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::b_u_bep   !Explicit momentum component X-direction
    668. 

  668.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::b_v_bep   !Explicit momentum component Y-direction
    669. 

  669.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::b_t_bep   !Explicit component pot. temperature
    670. 

  670.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::b_e_bep   !Explicit component TKE
    671. 

  671.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::b_q_bep   !Explicit component TKE
    672. 

  672.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::vl_bep    !Fraction air volume in grid cell
    673. 

  673.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::dlg_bep   !Height above ground
    674. 

  674.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::sf_bep  !Fraction air at the face of grid cell
    675. 

  675.    REAL, OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) ::dl_u_bep  !Length scale
    676. 

  676. 

  677. ! Optional
    678. 

  678. !
    679. 

  679. !  arguments for Ocean Mixed Layer Model
    680. 

  680.    REAL, DIMENSION( ims:ime , jms:jme ), OPTIONAL, INTENT(INOUT )::   TML, T0ML, HML, H0ML, HUML, HVML
    681. 

  681.    REAL, DIMENSION( ims:ime , jms:jme ), OPTIONAL, INTENT(IN    )::   F, TMOML
    682. 

  682.    REAL, DIMENSION( ims:ime , jms:jme ), OPTIONAL, INTENT(OUT   )::   CK, CKA, CD, CDA, USTM
    683. 

  683. 

  684. #if ( EM_CORE==1)
    685. 

  685.    REAL, DIMENSION( ims:ime , jms:jme ), &
    686. 

  686.         &OPTIONAL, INTENT(INOUT   ):: ch
    687. 

  687.    
    688. 

  688.    REAL, DIMENSION( ims:ime , kms:kme, jms:jme ), &
    689. 

  689.         &OPTIONAL, INTENT(IN   ):: tsq,qsq,cov
    690. 

  690. #endif
    691. 

  691. 

  692. 

  693.    INTEGER, OPTIONAL, INTENT(IN )::   slope_rad, topo_shading
    694. 

  694.    INTEGER, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(IN):: shadowmask
    695. 

  695.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT):: swnorm
    696. 

  696.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(IN):: slope,slp_azi
    697. 

  697. 

  698.    INTEGER, OPTIONAL, INTENT(IN )::   ISFTCFLX,IZ0TLND
    699. 

  699.    INTEGER, OPTIONAL, INTENT(IN )::   OMLCALL
    700. 

  700.    REAL   , OPTIONAL, INTENT(IN )::   OML_HML0
    701. 

  701.    REAL   , OPTIONAL, INTENT(IN )::   OML_GAMMA
    702. 

  702. !
    703. 

  703. !  Observation nudging
    704. 

  704. !
    705. 

  705.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(OUT)::   uratx  !Added for obs-nudging
    706. 

  706.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(OUT)::   vratx  !Added for obs-nudging
    707. 

  707.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(OUT)::   tratx  !Added for obs-nudging
    708. 

  708. !
    709. 

  709. !  PX LSM Surface Grid Analysis nudging
    710. 

  710. !
    711. 

  711.    INTEGER, OPTIONAL, INTENT(IN)    :: pxlsm_smois_init, pxlsm_soil_nudge, ANAL_INTERVAL
    712. 

  712.    REAL, DIMENSION( ims:ime, NLCAT, jms:jme ) , OPTIONAL, INTENT(INOUT)::   LANDUSEF
    713. 

  713.    REAL, DIMENSION( ims:ime, NSCAT, jms:jme ) , OPTIONAL, INTENT(INOUT)::   SOILCTOP, SOILCBOT
    714. 

  714.    REAL, DIMENSION( ims:ime , jms:jme ), OPTIONAL, INTENT(INOUT)::   VEGF_PX
    715. 

  715.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   RA
    716. 

  716.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   RS
    717. 

  717.    REAL, DIMENSION( ims:ime, jms:jme ) , OPTIONAL, INTENT(INOUT)::   LAI
    718. 

  718.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(OUT)::   T2OBS
    719. 

  719.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(OUT)::   Q2OBS
    720. 

  720. 

  721.    REAL,       DIMENSION( ims:ime,  jms:jme ),                           &
    722. 

  722.                OPTIONAL, INTENT(INOUT)    ::      t2_ndg_old,            &
    723. 

  723.                                                   q2_ndg_old,            &
    724. 

  724.                                                   t2_ndg_new,            &
    725. 

  725.                                                   q2_ndg_new,            &
    726. 

  726.                                                   sn_ndg_old,            &
    727. 

  727.                                                   sn_ndg_new
    728. 

  728. !
    729. 

  729. !
    730. 

  730. ! Flags relating to the optional tendency arrays declared above
    731. 

  731. ! Models that carry the optional tendencies will provdide the
    732. 

  732. ! optional arguments at compile time; these flags all the model
    733. 

  733. ! to determine at run-time whether a particular tracer is in
    734. 

  734. ! use or not.
    735. 

  735. !
    736. 

  736.    LOGICAL, INTENT(IN), OPTIONAL ::                             &
    737. 

  737.                                                       f_qv      &
    738. 

  738.                                                      ,f_qc      &
    739. 

  739.                                                      ,f_qr      &
    740. 

  740.                                                      ,f_qi      &
    741. 

  741.                                                      ,f_qs      &
    742. 

  742.                                                      ,f_qg
    743. 

  743. 

  744.    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ),                 &
    745. 

  745.          OPTIONAL, INTENT(INOUT) ::                              &
    746. 

  746.                       ! optional moisture tracers
    747. 

  747.                       ! 2 time levels; if only one then use CURR
    748. 

  748.                       qv_curr, qc_curr, qr_curr                  &
    749. 

  749.                      ,qi_curr, qs_curr, qg_curr
    750. 

  750.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(IN)   ::   snowncv
    751. 

  751.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   capg
    752. 

  752.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   emiss
    753. 

  753.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   hol
    754. 

  754.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   mol
    755. 

  755.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   regime
    756. 

  756.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(IN )::     rainncv
    757. 

  757.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(IN )::     rainshv
    758. 

  758.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   RAINBL
    759. 

  759.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   t2
    760. 

  760.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(IN )::     thc
    761. 

  761.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   qsg
    762. 

  762.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   qvg
    763. 

  763.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   qcg
    764. 

  764.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   dew
    765. 

  765.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   soilt1
    766. 

  766.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   tsnav
    767. 

  767.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   potevp ! NMM LSM
    768. 

  768.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   snopcx ! NMM LSM
    769. 

  769.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   soiltb ! NMM LSM
    770. 

  770.    REAL, DIMENSION( ims:ime, jms:jme ), OPTIONAL, INTENT(INOUT)::   sr ! NMM and RUC LSM
    771. 

  771.    REAL, DIMENSION( ims:ime, 1:num_soil_layers, jms:jme ), OPTIONAL, INTENT(INOUT)::   smfr3d
    772. 

  772.    REAL, DIMENSION( ims:ime, 1:num_soil_layers, jms:jme ), OPTIONAL, INTENT(INOUT)::   keepfr3dflag
    773. 

  773. 

  774.    REAL, DIMENSION( ims:ime, jms:jme ) , INTENT(OUT), OPTIONAL  ::   NOAHRES
    775. 

  775. 

  776. ! Variables for TEMF surface layer
    777. 

  777.    REAL,OPTIONAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT) :: te_temf
    778. 

  778.    REAL,OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: hd_temf, exch_temf, wm_temf
    779. 

  779.    REAL,OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN   ) :: fCor
    780. 

  780. 

  781. ! Variables for ideal SCM surface layer
    782. 

  782.    REAL,OPTIONAL, INTENT(INOUT) :: hfx_force,lh_force,tsk_force
    783. 

  783.    REAL,OPTIONAL, INTENT(IN   ) :: hfx_force_tend,lh_force_tend,tsk_force_tend
    784. 

  784. 

  785. !  LOCAL  VAR
    786. 

  786. 

  787.    REAL,       DIMENSION( ims:ime, kms:kme, jms:jme ) ::v_phytmp
    788. 

  788.    REAL,       DIMENSION( ims:ime, kms:kme, jms:jme ) ::u_phytmp
    789. 

  789. 

  790.    REAL,       DIMENSION( ims:ime, jms:jme )          ::  ZOL
    791. 

  791. 

  792.    REAL,       DIMENSION( ims:ime, jms:jme )          ::          &
    793. 

  793.                                                              QGH, &
    794. 

  794.                                                              CHS, &
    795. 

  795.                                                              CPM, &
    796. 

  796.                                                             CHS2, &
    797. 

  797.                                                             CQS2
    798. 

  798. ! SSIB local variables
    799. 

  799.    REAL ZDIFF
    800. 

  800. !
    801. 

  801.    REAL    :: DTMIN,DTBL
    802. 

  802. !
    803. 

  803.    INTEGER :: i,J,K,NK,jj,ij
    804. 

  804.    INTEGER :: gfdl_ntsflg
    805. 

  805.    LOGICAL :: radiation, myj, frpcpn, isisfc
    806. 

  806.    LOGICAL, INTENT(in), OPTIONAL :: rdlai2d
    807. 

  807.    LOGICAL, INTENT(in), OPTIONAL :: usemonalb
    808. 

  808.    REAL    :: total_depth,mid_point_depth
    809. 

  809.    REAL    :: tconst,tprior,tnew,yrday,deltat
    810. 

  810.    REAL    :: SWSAVE
    811. 

  811.    REAL,       DIMENSION( ims:ime, jms:jme )          ::  GSWSAVE
    812. 

  812. !-------------------------------------------------
    813. 

  813. ! urban related variables are added to declaration
    814. 

  814. !-------------------------------------------------
    815. 

  815.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: CMR_SFCDIF
    816. 

  816.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: CHR_SFCDIF
    817. 

  817.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: CMC_SFCDIF
    818. 

  818.    REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: CHC_SFCDIF
    819. 

  819.      REAL, OPTIONAL, INTENT(IN) :: DECLIN, SOLCON
    820. 

  820.      REAL, OPTIONAL , DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: COSZEN
    821. 

  821.      REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: HRANG
    822. 

  822.      REAL, OPTIONAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN) :: XLAT_URB2D  !urban
    823. 

  823.      INTEGER,  INTENT(IN) :: num_roof_layers                         !urban
    824. 

  824.      INTEGER,  INTENT(IN) :: num_wall_layers                     …
    825.

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