PageRenderTime 52ms CodeModel.GetById 14ms RepoModel.GetById 1ms app.codeStats 0ms

/wrfv2_fire/phys/module_physics_addtendc.F

http://github.com/jbeezley/wrf-fire
FORTRAN Legacy | 2022 lines | 1180 code | 394 blank | 448 comment | 100 complexity | 624d41128f439e094353cf29568747e3 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: PHYSICS
    2. 

  2. !
    3. 

  3. ! note: this module really belongs in the dyn_em directory since it is 
    4. 

  4. !       specific only to the EM core. Leaving here for now, with an 
    5. 

  5. !       #if ( EM_CORE == 1 ) directive. JM 20031201
    6. 

  6. !
    7. 

  7. 

  8. !  This MODULE holds the routines which are used to perform updates of the
    9. 

  9. !  model C-grid tendencies with physics A-grid tendencies
    10. 

  10. !  The module consolidates code that was (up to v1.2) duplicated in 
    11. 

  11. !  module_em and module_rk and in
    12. 

  12. !  module_big_step_utilities.F and module_big_step_utilities_em.F
    13. 

  13. 

  14. !  This MODULE CONTAINS the following routines:
    15. 

  15. !  update_phy_ten, phy_ra_ten, phy_bl_ten, phy_cu_ten, advance_ppt,
    16. 

  16. !  add_a2a, add_a2c_u, and add_a2c_v
    17. 

  17. 

  18. 

  19. MODULE module_physics_addtendc
    20. 

  20. 

  21. #if ( EM_CORE == 1 )
    22. 

  22. 

  23.    USE module_state_description
    24. 

  24.    USE module_configure
    25. 

  25. 

  26. CONTAINS
    27. 

  27. 

  28. SUBROUTINE update_phy_ten(rph_tendf,rt_tendf,ru_tendf,rv_tendf,moist_tendf, &
    29. 

  29.                       scalar_tendf,mu_tendf,                                &
    30. 

  30.                       RTHRATEN,RTHBLTEN,RTHCUTEN,RTHSHTEN,                  &
    31. 

  31.                       RUBLTEN,RUCUTEN,RUSHTEN,                              &
    32. 

  32.                       RVBLTEN,RVCUTEN,RVSHTEN,                              &
    33. 

  33.                       RQVBLTEN,RQCBLTEN,RQIBLTEN,                           &
    34. 

  34.                       RQVCUTEN,RQCCUTEN,RQRCUTEN,RQICUTEN,RQSCUTEN,         &
    35. 

  35.                       RQVSHTEN,RQCSHTEN,RQRSHTEN,RQISHTEN,RQSSHTEN,RQGSHTEN,&
    36. 

  36.                       RUNDGDTEN,RVNDGDTEN,RTHNDGDTEN,RPHNDGDTEN,            &
    37. 

  37.                       RQVNDGDTEN,RMUNDGDTEN,                                &
    38. 

  38.                       rthfrten,rqvfrten,                                    & !fire
    39. 

  39.                       n_moist,n_scalar,config_flags,rk_step,adv_moist_cond, &
    40. 

  40.                       ids, ide, jds, jde, kds, kde,                         &
    41. 

  41.                       ims, ime, jms, jme, kms, kme,                         &
    42. 

  42.                       its, ite, jts, jte, kts, kte                          )
    43. 

  43. !-------------------------------------------------------------------
    44. 

  44.    IMPLICIT NONE
    45. 

  45. !-------------------------------------------------------------------
    46. 

  46. 

  47.    TYPE(grid_config_rec_type   ) ,   INTENT(IN   ) :: config_flags
    48. 

  48. 

  49.    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde,   &
    50. 

  50.                                    ims, ime, jms, jme, kms, kme,   &
    51. 

  51.                                    its, ite, jts, jte, kts, kte,   &
    52. 

  52.                                    n_moist,n_scalar,rk_step
    53. 

  53. 

  54.    LOGICAL , INTENT(IN)        :: adv_moist_cond
    55. 

  55. 

  56.    REAL , DIMENSION(ims:ime , kms:kme, jms:jme),INTENT(INOUT) ::   &
    57. 

  57.                                                          ru_tendf, &
    58. 

  58.                                                          rv_tendf, &
    59. 

  59.                                                          rt_tendf, &
    60. 

  60.                                                          rph_tendf
    61. 

  61. 

  62.    REAL , DIMENSION(ims:ime , jms:jme),INTENT(INOUT) ::  mu_tendf
    63. 

  63. 

  64.    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_moist),           &
    65. 

  65.           INTENT(INOUT)     ::                        moist_tendf
    66. 

  66. 

  67.    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_scalar),           &
    68. 

  68.           INTENT(INOUT)     ::                        scalar_tendf
    69. 

  69. 

  70.    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
    71. 

  71.                                                        RTHRATEN, &
    72. 

  72.                                                        RTHBLTEN, &
    73. 

  73.                                                        RTHCUTEN, &
    74. 

  74.                                                        RTHSHTEN, &
    75. 

  75.                                                         RUBLTEN, &
    76. 

  76.                                                         RUCUTEN, &
    77. 

  77.                                                         RUSHTEN, &
    78. 

  78.                                                         RVBLTEN, &
    79. 

  79.                                                         RVCUTEN, &
    80. 

  80.                                                         RVSHTEN, &
    81. 

  81.                                                        RQVBLTEN, &
    82. 

  82.                                                        RQCBLTEN, &
    83. 

  83.                                                        RQIBLTEN, &
    84. 

  84.                                                        RQVCUTEN, &
    85. 

  85.                                                        RQCCUTEN, &
    86. 

  86.                                                        RQRCUTEN, &
    87. 

  87.                                                        RQICUTEN, &
    88. 

  88.                                                        RQSCUTEN, &
    89. 

  89.                                                        RQVSHTEN, &
    90. 

  90.                                                        RQCSHTEN, &
    91. 

  91.                                                        RQRSHTEN, &
    92. 

  92.                                                        RQISHTEN, &
    93. 

  93.                                                        RQSSHTEN, &
    94. 

  94.                                                        RQGSHTEN, &
    95. 

  95.                                                      RTHNDGDTEN, &
    96. 

  96.                                                      RPHNDGDTEN, &
    97. 

  97.                                                      RQVNDGDTEN, &
    98. 

  98.                                                       RUNDGDTEN, &
    99. 

  99.                                                       RVNDGDTEN
    100. 

  100. 

  101.    REAL, DIMENSION(ims:ime, jms:jme), INTENT(IN  ) :: RMUNDGDTEN
    102. 

  102. 

  103.    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   & ! fire
    104. 

  104.                                                        rthfrten, &
    105. 

  105.                                                        rqvfrten   
    106. 

  106. !------------------------------------------------------------------
    107. 

  107. !  set up loop bounds for this grid's boundary conditions
    108. 

    109. 

  109.    if (config_flags%ra_lw_physics .gt. 0 .or.                  &
    110. 

  110.        config_flags%ra_sw_physics .gt. 0)                      &
    111. 

  111.       CALL phy_ra_ten(config_flags,rt_tendf,RTHRATEN,          &
    112. 

  112.                       ids, ide, jds, jde, kds, kde,            &
    113. 

  113.                       ims, ime, jms, jme, kms, kme,            &
    114. 

  114.                       its, ite, jts, jte, kts, kte             )
    115. 

  115. 

  116.    if (config_flags%bl_pbl_physics .gt. 0)                     &
    117. 

  117.       CALL phy_bl_ten(config_flags,rk_step,n_moist,n_scalar,   &
    118. 

  118.                       rt_tendf,ru_tendf,rv_tendf,moist_tendf,  &
    119. 

  119.                       scalar_tendf,adv_moist_cond,             &
    120. 

  120.                       RTHBLTEN,RUBLTEN,RVBLTEN,                &
    121. 

  121.                       RQVBLTEN,RQCBLTEN,RQIBLTEN,              &
    122. 

  122.                       ids, ide, jds, jde, kds, kde,            &
    123. 

  123.                       ims, ime, jms, jme, kms, kme,            &
    124. 

  124.                       its, ite, jts, jte, kts, kte             )
    125. 

  125. 

  126.    if (config_flags%cu_physics .gt. 0)                         &
    127. 

  127.       CALL phy_cu_ten(config_flags,rk_step,n_moist,n_scalar,   &
    128. 

  128.                       rt_tendf,ru_tendf,rv_tendf,              &
    129. 

  129.                       RUCUTEN,RVCUTEN,RTHCUTEN,                &
    130. 

  130.                       RQVCUTEN,RQCCUTEN,RQRCUTEN,              &
    131. 

  131.                       RQICUTEN,RQSCUTEN,moist_tendf,           &
    132. 

  132.                       scalar_tendf,adv_moist_cond,             &
    133. 

  133.                       ids, ide, jds, jde, kds, kde,            &
    134. 

  134.                       ims, ime, jms, jme, kms, kme,            &
    135. 

  135.                       its, ite, jts, jte, kts, kte             )
    136. 

  136. 

  137.    if (config_flags%shcu_physics .gt. 0)                       &
    138. 

  138.       CALL phy_shcu_ten(config_flags,rk_step,n_moist,          &
    139. 

  139.                       rt_tendf,ru_tendf,rv_tendf,              &
    140. 

  140.                       RUSHTEN,RVSHTEN,RTHSHTEN,                &
    141. 

  141.                       RQVSHTEN,RQCSHTEN,RQRSHTEN,              &
    142. 

  142.                       RQISHTEN,RQSSHTEN,RQGSHTEN,moist_tendf,  &
    143. 

  143.                       ids, ide, jds, jde, kds, kde,            &
    144. 

  144.                       ims, ime, jms, jme, kms, kme,            &
    145. 

  145.                       its, ite, jts, jte, kts, kte             )
    146. 

  146. 

  147.    if (config_flags%grid_fdda .gt. 0)                          &
    148. 

  148.       CALL phy_fg_ten(config_flags,rk_step,n_moist,            &
    149. 

  149.                       rph_tendf,rt_tendf,ru_tendf,rv_tendf,    &
    150. 

  150.                       mu_tendf, moist_tendf,                   &
    151. 

  151.                       RUNDGDTEN,RVNDGDTEN,RTHNDGDTEN,          &
    152. 

  152.                       RPHNDGDTEN,RQVNDGDTEN,RMUNDGDTEN,        &
    153. 

  153.                       ids, ide, jds, jde, kds, kde,            &
    154. 

  154.                       ims, ime, jms, jme, kms, kme,            &
    155. 

  155.                       its, ite, jts, jte, kts, kte             )
    156. 

  156. 

  157.    if (config_flags%ifire .gt. 0)                              & ! fire
    158. 

  158.       CALL phy_fr_ten(config_flags,rk_step,n_moist,            &
    159. 

  159.                       rt_tendf,ru_tendf,rv_tendf,              &
    160. 

  160.                       mu_tendf, moist_tendf,                   &
    161. 

  161.                       rthfrten,rqvfrten,                       &
    162. 

  162.                       ids, ide, jds, jde, kds, kde,            &
    163. 

  163.                       ims, ime, jms, jme, kms, kme,            &
    164. 

  164.                       its, ite, jts, jte, kts, kte             )
    165. 

  165. 

  166. END SUBROUTINE update_phy_ten
    167. 

  167. 

  168. !=================================================================
    169. 

  169. SUBROUTINE phy_ra_ten(config_flags,rt_tendf,RTHRATEN,            &
    170. 

  170.                       ids, ide, jds, jde, kds, kde,              &
    171. 

  171.                       ims, ime, jms, jme, kms, kme,              &
    172. 

  172.                       its, ite, jts, jte, kts, kte               )
    173. 

  173. !-----------------------------------------------------------------
    174. 

  174.    IMPLICIT NONE
    175. 

  175. !-----------------------------------------------------------------
    176. 

  176.    TYPE(grid_config_rec_type  ) , INTENT(IN   ) :: config_flags
    177. 

  177. 

  178.    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde, &
    179. 

  179.                                    ims, ime, jms, jme, kms, kme, &
    180. 

  180.                                    its, ite, jts, jte, kts, kte
    181. 

  181. 

  182.    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
    183. 

  183.                                                        RTHRATEN
    184. 

  184. 

  185.    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) ::  &
    186. 

  186.                                                        rt_tendf
    187. 

  187. 

  188. ! LOCAL VARS
    189. 

  189. 

  190.    INTEGER :: i,j,k
    191. 

  191. 

  192.    CALL add_a2a(rt_tendf,RTHRATEN,config_flags,                  &
    193. 

  193.                 ids,ide, jds, jde, kds, kde,                     &
    194. 

  194.                 ims, ime, jms, jme, kms, kme,                    &
    195. 

  195.                 its, ite, jts, jte, kts, kte                     )
    196. 

  196. 

  197. END SUBROUTINE phy_ra_ten
    198. 

  198. 

  199. !=================================================================
    200. 

  200. SUBROUTINE phy_bl_ten(config_flags,rk_step,n_moist,n_scalar,     &
    201. 

  201.                       rt_tendf,ru_tendf,rv_tendf,moist_tendf,    &
    202. 

  202.                       scalar_tendf,adv_moist_cond,               &
    203. 

  203.                       RTHBLTEN,RUBLTEN,RVBLTEN,                  & 
    204. 

  204.                       RQVBLTEN,RQCBLTEN,RQIBLTEN,                &
    205. 

  205.                       ids, ide, jds, jde, kds, kde,              &
    206. 

  206.                       ims, ime, jms, jme, kms, kme,              &
    207. 

  207.                       its, ite, jts, jte, kts, kte               )
    208. 

  208. !-----------------------------------------------------------------
    209. 

  209.    IMPLICIT NONE
    210. 

  210. !-----------------------------------------------------------------
    211. 

  211.    TYPE(grid_config_rec_type) ,     INTENT(IN   ) :: config_flags
    212. 

  212. 

  213.    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde, &
    214. 

  214.                                    ims, ime, jms, jme, kms, kme, &
    215. 

  215.                                    its, ite, jts, jte, kts, kte, &
    216. 

  216.                                    n_moist, n_scalar, rk_step
    217. 

  217. 

  218.    LOGICAL , INTENT(IN)     :: adv_moist_cond
    219. 

  219. 

  220.    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_moist),         &
    221. 

  221.           INTENT(INOUT)     ::                      moist_tendf
    222. 

  222. 

  223.    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_scalar),         &
    224. 

  224.           INTENT(INOUT)     ::                      scalar_tendf
    225. 

  225. 

  226.    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   & 
    227. 

  227.                                                        RTHBLTEN, &
    228. 

  228.                                                         RUBLTEN, &
    229. 

  229.                                                         RVBLTEN, &
    230. 

  230.                                                        RQVBLTEN, &
    231. 

  231.                                                        RQCBLTEN, &
    232. 

  232.                                                        RQIBLTEN
    233. 

  233. 

  234.    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT)::   &
    235. 

  235.                                                        rt_tendf, &
    236. 

  236.                                                        ru_tendf, &
    237. 

  237.                                                        rv_tendf
    238. 

  238. ! LOCAL VARS
    239. 

  239. 

  240.    INTEGER :: i,j,k,IBGN,IEND,JBGN,JEND
    241. 

  241. 

  242. !-----------------------------------------------------------------
    243. 

  243. 

  244.    SELECT CASE(config_flags%bl_pbl_physics)
    245. 

  245. 

  246.       CASE (YSUSCHEME)
    247. 

  247. 

  248.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    249. 

  249.                 ids,ide, jds, jde, kds, kde,                     &
    250. 

  250.                 ims, ime, jms, jme, kms, kme,                    &
    251. 

  251.                 its, ite, jts, jte, kts, kte                     )
    252. 

  252. 

  253.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    254. 

  254.                 ids,ide, jds, jde, kds, kde,                     &
    255. 

  255.                 ims, ime, jms, jme, kms, kme,                    &
    256. 

  256.                 its, ite, jts, jte, kts, kte                     )
    257. 

  257. 

  258.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    259. 

  259.                 ids,ide, jds, jde, kds, kde,                     &
    260. 

  260.                 ims, ime, jms, jme, kms, kme,                    &
    261. 

  261.                 its, ite, jts, jte, kts, kte                     )
    262. 

  262. 

  263.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    264. 

  264.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    265. 

  265.                 config_flags,                                    &
    266. 

  266.                 ids,ide, jds, jde, kds, kde,                     &
    267. 

  267.                 ims, ime, jms, jme, kms, kme,                    &
    268. 

  268.                 its, ite, jts, jte, kts, kte                     )
    269. 

  269. 

  270.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    271. 

  271.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    272. 

  272.                 config_flags,                                    &
    273. 

  273.                 ids,ide, jds, jde, kds, kde,                     &
    274. 

  274.                 ims, ime, jms, jme, kms, kme,                    &
    275. 

  275.                 its, ite, jts, jte, kts, kte                     )
    276. 

  276.      
    277. 

  277.         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
    278. 

  278.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
    279. 

  279.                 config_flags,                                    &
    280. 

  280.                 ids,ide, jds, jde, kds, kde,                     &
    281. 

  281.                 ims, ime, jms, jme, kms, kme,                    &
    282. 

  282.                 its, ite, jts, jte, kts, kte                     )
    283. 

  283. 

  284.        IF(.not. adv_moist_cond)THEN
    285. 

  285. 

  286.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    287. 

  287.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
    288. 

  288.                 config_flags,                                    &
    289. 

  289.                 ids,ide, jds, jde, kds, kde,                     &
    290. 

  290.                 ims, ime, jms, jme, kms, kme,                    &
    291. 

  291.                 its, ite, jts, jte, kts, kte                     )
    292. 

  292.      
    293. 

  293.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    294. 

  294.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
    295. 

  295.                 config_flags,                                    &
    296. 

  296.                 ids,ide, jds, jde, kds, kde,                     &
    297. 

  297.                 ims, ime, jms, jme, kms, kme,                    &
    298. 

  298.                 its, ite, jts, jte, kts, kte                     )
    299. 

  299.        ENDIF
    300. 

  300. 

  301.       CASE (MRFSCHEME)
    302. 

  302. 

  303.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    304. 

  304.                 ids,ide, jds, jde, kds, kde,                     &
    305. 

  305.                 ims, ime, jms, jme, kms, kme,                    &
    306. 

  306.                 its, ite, jts, jte, kts, kte                     )
    307. 

  307. 

  308.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    309. 

  309.                 ids,ide, jds, jde, kds, kde,                     &
    310. 

  310.                 ims, ime, jms, jme, kms, kme,                    &
    311. 

  311.                 its, ite, jts, jte, kts, kte                     )
    312. 

  312. 

  313.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    314. 

  314.                 ids,ide, jds, jde, kds, kde,                     &
    315. 

  315.                 ims, ime, jms, jme, kms, kme,                    &
    316. 

  316.                 its, ite, jts, jte, kts, kte                     )
    317. 

  317. 

  318.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    319. 

  319.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    320. 

  320.                 config_flags,                                    &
    321. 

  321.                 ids,ide, jds, jde, kds, kde,                     &
    322. 

  322.                 ims, ime, jms, jme, kms, kme,                    &
    323. 

  323.                 its, ite, jts, jte, kts, kte                     )
    324. 

  324. 

  325.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    326. 

  326.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    327. 

  327.                 config_flags,                                    &
    328. 

  328.                 ids,ide, jds, jde, kds, kde,                     &
    329. 

  329.                 ims, ime, jms, jme, kms, kme,                    &
    330. 

  330.                 its, ite, jts, jte, kts, kte                     )
    331. 

  331.      
    332. 

  332.         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
    333. 

  333.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
    334. 

  334.                 config_flags,                                    &
    335. 

  335.                 ids,ide, jds, jde, kds, kde,                     &
    336. 

  336.                 ims, ime, jms, jme, kms, kme,                    &
    337. 

  337.                 its, ite, jts, jte, kts, kte                     )
    338. 

  338. 

  339.        IF(.not. adv_moist_cond)THEN
    340. 

  340. 

  341.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    342. 

  342.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
    343. 

  343.                 config_flags,                                    &
    344. 

  344.                 ids,ide, jds, jde, kds, kde,                     &
    345. 

  345.                 ims, ime, jms, jme, kms, kme,                    &
    346. 

  346.                 its, ite, jts, jte, kts, kte                     )
    347. 

  347.      
    348. 

  348.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    349. 

  349.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
    350. 

  350.                 config_flags,                                    &
    351. 

  351.                 ids,ide, jds, jde, kds, kde,                     &
    352. 

  352.                 ims, ime, jms, jme, kms, kme,                    &
    353. 

  353.                 its, ite, jts, jte, kts, kte                     )
    354. 

  354.        ENDIF
    355. 

  355. 

  356.       CASE (ACMPBLSCHEME)
    357. 

  357. 

  358.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    359. 

  359.                 ids,ide, jds, jde, kds, kde,                     &
    360. 

  360.                 ims, ime, jms, jme, kms, kme,                    &
    361. 

  361.                 its, ite, jts, jte, kts, kte                     )
    362. 

  362. 

  363.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    364. 

  364.                 ids,ide, jds, jde, kds, kde,                     &
    365. 

  365.                 ims, ime, jms, jme, kms, kme,                    &
    366. 

  366.                 its, ite, jts, jte, kts, kte                     )
    367. 

  367. 

  368.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    369. 

  369.                 ids,ide, jds, jde, kds, kde,                     &
    370. 

  370.                 ims, ime, jms, jme, kms, kme,                    &
    371. 

  371.                 its, ite, jts, jte, kts, kte                     )
    372. 

  372. 

  373.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    374. 

  374.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    375. 

  375.                 config_flags,                                    &
    376. 

  376.                 ids,ide, jds, jde, kds, kde,                     &
    377. 

  377.                 ims, ime, jms, jme, kms, kme,                    &
    378. 

  378.                 its, ite, jts, jte, kts, kte                     )
    379. 

  379. 

  380.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    381. 

  381.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    382. 

  382.                 config_flags,                                    &
    383. 

  383.                 ids,ide, jds, jde, kds, kde,                     &
    384. 

  384.                 ims, ime, jms, jme, kms, kme,                    &
    385. 

  385.                 its, ite, jts, jte, kts, kte                     )
    386. 

  386.      
    387. 

  387.         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
    388. 

  388.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
    389. 

  389.                 config_flags,                                    &
    390. 

  390.                 ids,ide, jds, jde, kds, kde,                     &
    391. 

  391.                 ims, ime, jms, jme, kms, kme,                    &
    392. 

  392.                 its, ite, jts, jte, kts, kte                     )
    393. 

  393. 

  394.        IF(.not. adv_moist_cond)THEN
    395. 

  395. 

  396.         if (P_QT .ge. PARAM_FIRST_SCALAR)THEN
    397. 

  397.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
    398. 

  398.                 config_flags,                                    &
    399. 

  399.                 ids,ide, jds, jde, kds, kde,                     &
    400. 

  400.                 ims, ime, jms, jme, kms, kme,                    &
    401. 

  401.                 its, ite, jts, jte, kts, kte                     )
    402. 

  402. 

  403.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
    404. 

  404.                 config_flags,                                    &
    405. 

  405.                 ids,ide, jds, jde, kds, kde,                     &
    406. 

  406.                 ims, ime, jms, jme, kms, kme,                    &
    407. 

  407.                 its, ite, jts, jte, kts, kte                     )
    408. 

  408.         ENDIF
    409. 

  409.      
    410. 

  410.        ENDIF
    411. 

  411. 

  412.       CASE (MYJPBLSCHEME)
    413. 

  413. 

  414.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    415. 

  415.                 ids,ide, jds, jde, kds, kde,                     &
    416. 

  416.                 ims, ime, jms, jme, kms, kme,                    &
    417. 

  417.                 its, ite, jts, jte, kts, kte                     )
    418. 

  418. 

  419.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    420. 

  420.                 ids,ide, jds, jde, kds, kde,                     &
    421. 

  421.                 ims, ime, jms, jme, kms, kme,                    &
    422. 

  422.                 its, ite, jts, jte, kts, kte                     )
    423. 

  423. 

  424.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    425. 

  425.                 ids,ide, jds, jde, kds, kde,                     &
    426. 

  426.                 ims, ime, jms, jme, kms, kme,                    &
    427. 

  427.                 its, ite, jts, jte, kts, kte                     )
    428. 

  428. 

  429.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    430. 

  430.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    431. 

  431.                 config_flags,                                    &
    432. 

  432.                 ids,ide, jds, jde, kds, kde,                     &
    433. 

  433.                 ims, ime, jms, jme, kms, kme,                    &
    434. 

  434.                 its, ite, jts, jte, kts, kte                     )
    435. 

  435. 

  436.        IF(.not. adv_moist_cond)THEN
    437. 

  437. 

  438.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    439. 

  439.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
    440. 

  440.                 config_flags,                                    &
    441. 

  441.                 ids,ide, jds, jde, kds, kde,                     &
    442. 

  442.                 ims, ime, jms, jme, kms, kme,                    &
    443. 

  443.                 its, ite, jts, jte, kts, kte                     )
    444. 

  444.      
    445. 

  445.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    446. 

  446.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
    447. 

  447.                 config_flags,                                    &
    448. 

  448.                 ids,ide, jds, jde, kds, kde,                     &
    449. 

  449.                 ims, ime, jms, jme, kms, kme,                    &
    450. 

  450.                 its, ite, jts, jte, kts, kte                     )
    451. 

  451.      
    452. 

  452. !       if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    453. 

  453. !          CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQSBLTEN,  &
    454. 

  454. !               config_flags,                                    &
    455. 

  455. !               ids,ide, jds, jde, kds, kde,                     &
    456. 

  456. !               ims, ime, jms, jme, kms, kme,                    &
    457. 

  457. !               its, ite, jts, jte, kts, kte                     )
    458. 

  458. !    
    459. 

  459. !       if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    460. 

  460. !          CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQRBLTEN,  &
    461. 

  461. !               config_flags,                                    &
    462. 

  462. !               ids,ide, jds, jde, kds, kde,                     &
    463. 

  463. !               ims, ime, jms, jme, kms, kme,                    &
    464. 

  464. !               its, ite, jts, jte, kts, kte                     )
    465. 

  465. !    
    466. 

  466. !       if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    467. 

  467. !          CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQGBLTEN,  &
    468. 

  468. !               config_flags,                                    &
    469. 

  469. !               ids,ide, jds, jde, kds, kde,                     &
    470. 

  470. !               ims, ime, jms, jme, kms, kme,                    &
    471. 

  471. !               its, ite, jts, jte, kts, kte                     )
    472. 

  472.      
    473. 

  473.        ELSE
    474. 

  474. 

  475.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    476. 

  476.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    477. 

  477.                 config_flags,                                    &
    478. 

  478.                 ids,ide, jds, jde, kds, kde,                     &
    479. 

  479.                 ims, ime, jms, jme, kms, kme,                    &
    480. 

  480.                 its, ite, jts, jte, kts, kte                     )
    481. 

  481. 

  482.         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
    483. 

  483.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
    484. 

  484.                 config_flags,                                    &
    485. 

  485.                 ids,ide, jds, jde, kds, kde,                     &
    486. 

  486.                 ims, ime, jms, jme, kms, kme,                    &
    487. 

  487.                 its, ite, jts, jte, kts, kte                     )
    488. 

  488.      
    489. 

  489. !       if (P_QS .ge. PARAM_FIRST_SCALAR)                                         &
    490. 

  490. !          CALL add_a2a(moist_tendf(ims,kms,jms,P_QS),RQSBLTEN,  &
    491. 

  491. !               config_flags,                                    &
    492. 

  492. !               ids,ide, jds, jde, kds, kde,                     &
    493. 

  493. !               ims, ime, jms, jme, kms, kme,                    &
    494. 

  494. !               its, ite, jts, jte, kts, kte                     )
    495. 

  495. !    
    496. 

  496. !       if (P_QR .ge. PARAM_FIRST_SCALAR)                                         &
    497. 

  497. !          CALL add_a2a(moist_tendf(ims,kms,jms,P_QR),RQRBLTEN,  &
    498. 

  498. !               config_flags,                                    &
    499. 

  499. !               ids,ide, jds, jde, kds, kde,                     &
    500. 

  500. !               ims, ime, jms, jme, kms, kme,                    &
    501. 

  501. !               its, ite, jts, jte, kts, kte                     )
    502. 

  502. !    
    503. 

  503. !       if (P_QG .ge. PARAM_FIRST_SCALAR)                                         &
    504. 

  504. !          CALL add_a2a(moist_tendf(ims,kms,jms,P_QG),RQGBLTEN,  &
    505. 

  505. !               config_flags,                                    &
    506. 

  506. !               ids,ide, jds, jde, kds, kde,                     &
    507. 

  507. !               ims, ime, jms, jme, kms, kme,                    &
    508. 

  508. !               its, ite, jts, jte, kts, kte                     )
    509. 

  509.      
    510. 

  510.        ENDIF
    511. 

  511. 

  512.       CASE (QNSEPBLSCHEME,QNSEPBL09SCHEME)
    513. 

  513. 

  514.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    515. 

  515.                 ids,ide, jds, jde, kds, kde,                     &
    516. 

  516.                 ims, ime, jms, jme, kms, kme,                    &
    517. 

  517.                 its, ite, jts, jte, kts, kte                     )
    518. 

  518. 

  519.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    520. 

  520.                 ids,ide, jds, jde, kds, kde,                     &
    521. 

  521.                 ims, ime, jms, jme, kms, kme,                    &
    522. 

  522.                 its, ite, jts, jte, kts, kte                     )
    523. 

  523. 

  524.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    525. 

  525.                 ids,ide, jds, jde, kds, kde,                     &
    526. 

  526.                 ims, ime, jms, jme, kms, kme,                    &
    527. 

  527.                 its, ite, jts, jte, kts, kte                     )
    528. 

  528. 

  529.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    530. 

  530.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    531. 

  531.                 config_flags,                                    &
    532. 

  532.                 ids,ide, jds, jde, kds, kde,                     &
    533. 

  533.                 ims, ime, jms, jme, kms, kme,                    &
    534. 

  534.                 its, ite, jts, jte, kts, kte                     )
    535. 

  535. 

  536.        IF(.not. adv_moist_cond)THEN
    537. 

  537. 

  538.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    539. 

  539.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
    540. 

  540.                 config_flags,                                    &
    541. 

  541.                 ids,ide, jds, jde, kds, kde,                     &
    542. 

  542.                 ims, ime, jms, jme, kms, kme,                    &
    543. 

  543.                 its, ite, jts, jte, kts, kte                     )
    544. 

  544.      
    545. 

  545.        ELSE
    546. 

  546. 

  547.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    548. 

  548.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    549. 

  549.                 config_flags,                                    &
    550. 

  550.                 ids,ide, jds, jde, kds, kde,                     &
    551. 

  551.                 ims, ime, jms, jme, kms, kme,                    &
    552. 

  552.                 its, ite, jts, jte, kts, kte                     )
    553. 

  553. 

  554.        ENDIF
    555. 

  555. 

  556.       CASE (GFSSCHEME)
    557. 

  557.                                                                                                                                         
    558. 

  558.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    559. 

  559.                 ids,ide, jds, jde, kds, kde,                     &
    560. 

  560.                 ims, ime, jms, jme, kms, kme,                    &
    561. 

  561.                 its, ite, jts, jte, kts, kte                     )
    562. 

  562.                                                                                                                                         
    563. 

  563.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    564. 

  564.                 ids,ide, jds, jde, kds, kde,                     &
    565. 

  565.                 ims, ime, jms, jme, kms, kme,                    &
    566. 

  566.                 its, ite, jts, jte, kts, kte                     )
    567. 

  567.                                                                                                                                         
    568. 

  568.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    569. 

  569.                 ids,ide, jds, jde, kds, kde,                     &
    570. 

  570.                 ims, ime, jms, jme, kms, kme,                    &
    571. 

  571.                 its, ite, jts, jte, kts, kte                     )
    572. 

  572.                                                                                                                                         
    573. 

  573.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    574. 

  574.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    575. 

  575.                 config_flags,                                    &
    576. 

  576.                 ids,ide, jds, jde, kds, kde,                     &
    577. 

  577.                 ims, ime, jms, jme, kms, kme,                    &
    578. 

  578.                 its, ite, jts, jte, kts, kte                     )
    579. 

  579.                                                                                                                                         
    580. 

  580.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    581. 

  581.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    582. 

  582.                 config_flags,                                    &
    583. 

  583.                 ids,ide, jds, jde, kds, kde,                     &
    584. 

  584.                 ims, ime, jms, jme, kms, kme,                    &
    585. 

  585.                 its, ite, jts, jte, kts, kte                     )
    586. 

  586.                                                                                                                                         
    587. 

  587.         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
    588. 

  588.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
    589. 

  589.                 config_flags,                                    &
    590. 

  590.                 ids,ide, jds, jde, kds, kde,                     &
    591. 

  591.                 ims, ime, jms, jme, kms, kme,                    &
    592. 

  592.                 its, ite, jts, jte, kts, kte                     )
    593. 

  593. 

  594.        IF(.not. adv_moist_cond)THEN
    595. 

  595. 

  596.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    597. 

  597.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
    598. 

  598.                 config_flags,                                    &
    599. 

  599.                 ids,ide, jds, jde, kds, kde,                     &
    600. 

  600.                 ims, ime, jms, jme, kms, kme,                    &
    601. 

  601.                 its, ite, jts, jte, kts, kte                     )
    602. 

  602.      
    603. 

  603.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    604. 

  604.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
    605. 

  605.                 config_flags,                                    &
    606. 

  606.                 ids,ide, jds, jde, kds, kde,                     &
    607. 

  607.                 ims, ime, jms, jme, kms, kme,                    &
    608. 

  608.                 its, ite, jts, jte, kts, kte                     )
    609. 

  609.        ENDIF
    610. 

  610. 

  611.       CASE (MYNNPBLSCHEME2,MYNNPBLSCHEME3)
    612. 

  612. 

  613.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    614. 

  614.                 ids,ide, jds, jde, kds, kde,                     &
    615. 

  615.                 ims, ime, jms, jme, kms, kme,                    &
    616. 

  616.                 its, ite, jts, jte, kts, kte                     )
    617. 

  617. 

  618.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    619. 

  619.                 ids,ide, jds, jde, kds, kde,                     &
    620. 

  620.                 ims, ime, jms, jme, kms, kme,                    &
    621. 

  621.                 its, ite, jts, jte, kts, kte                     )
    622. 

  622. 

  623.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    624. 

  624.                 ids,ide, jds, jde, kds, kde,                     &
    625. 

  625.                 ims, ime, jms, jme, kms, kme,                    &
    626. 

  626.                 its, ite, jts, jte, kts, kte                     )
    627. 

  627. 

  628.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    629. 

  629.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    630. 

  630.                 config_flags,                                    &
    631. 

  631.                 ids,ide, jds, jde, kds, kde,                     &
    632. 

  632.                 ims, ime, jms, jme, kms, kme,                    &
    633. 

  633.                 its, ite, jts, jte, kts, kte                     )
    634. 

  634. 

  635.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    636. 

  636.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    637. 

  637.                 config_flags,                                    &
    638. 

  638.                 ids,ide, jds, jde, kds, kde,                     &
    639. 

  639.                 ims, ime, jms, jme, kms, kme,                    &
    640. 

  640.                 its, ite, jts, jte, kts, kte                     )
    641. 

  641.      
    642. 

  642.        IF(.not. adv_moist_cond)THEN
    643. 

  643. 

  644.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    645. 

  645.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
    646. 

  646.                 config_flags,                                    &
    647. 

  647.                 ids,ide, jds, jde, kds, kde,                     &
    648. 

  648.                 ims, ime, jms, jme, kms, kme,                    &
    649. 

  649.                 its, ite, jts, jte, kts, kte                     )
    650. 

  650.      
    651. 

  651.        ENDIF
    652. 

  652. 

  653.        CASE (BOULACSCHEME)
    654. 

  654. 

  655.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    656. 

  656.                 ids,ide, jds, jde, kds, kde,                     &
    657. 

  657.                 ims, ime, jms, jme, kms, kme,                    &
    658. 

  658.                 its, ite, jts, jte, kts, kte                     )
    659. 

  659. 

  660.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    661. 

  661.                 ids,ide, jds, jde, kds, kde,                     &
    662. 

  662.                 ims, ime, jms, jme, kms, kme,                    &
    663. 

  663.                 its, ite, jts, jte, kts, kte                     )
    664. 

  664. 

  665.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    666. 

  666.                 ids,ide, jds, jde, kds, kde,                     &
    667. 

  667.                 ims, ime, jms, jme, kms, kme,                    &
    668. 

  668.                 its, ite, jts, jte, kts, kte                     )
    669. 

  669. 

  670.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    671. 

  671.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    672. 

  672.                 config_flags,                                    &
    673. 

  673.                 ids,ide, jds, jde, kds, kde,                     &
    674. 

  674.                 ims, ime, jms, jme, kms, kme,                    &
    675. 

  675.                 its, ite, jts, jte, kts, kte                     )
    676. 

  676. 

  677.        IF(.not. adv_moist_cond)THEN
    678. 

  678. 

  679.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    680. 

  680.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
    681. 

  681.                 config_flags,                                    &
    682. 

  682.                 ids,ide, jds, jde, kds, kde,                     &
    683. 

  683.                 ims, ime, jms, jme, kms, kme,                    &
    684. 

  684.                 its, ite, jts, jte, kts, kte                     )
    685. 

  685. 

  686.        ELSE
    687. 

  687. 

  688.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    689. 

  689.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    690. 

  690.                 config_flags,                                    &
    691. 

  691.                 ids,ide, jds, jde, kds, kde,                     &
    692. 

  692.                 ims, ime, jms, jme, kms, kme,                    &
    693. 

  693.                 its, ite, jts, jte, kts, kte                     )
    694. 

  694.        ENDIF
    695. 

  695. 

  696.        CASE (CAMUWPBLSCHEME)
    697. 

  697.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    698. 

  698.                 ids,ide, jds, jde, kds, kde,                     &
    699. 

  699.                 ims, ime, jms, jme, kms, kme,                    &
    700. 

  700.                 its, ite, jts, jte, kts, kte                     )
    701. 

  701. 

  702.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    703. 

  703.                 ids,ide, jds, jde, kds, kde,                     &
    704. 

  704.                 ims, ime, jms, jme, kms, kme,                    &
    705. 

  705.                 its, ite, jts, jte, kts, kte                     )
    706. 

  706. 

  707.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    708. 

  708.                 ids,ide, jds, jde, kds, kde,                     &
    709. 

  709.                 ims, ime, jms, jme, kms, kme,                    &
    710. 

  710.                 its, ite, jts, jte, kts, kte                     )
    711. 

  711. 

  712.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    713. 

  713.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    714. 

  714.                 config_flags,                                    &
    715. 

  715.                 ids,ide, jds, jde, kds, kde,                     &
    716. 

  716.                 ims, ime, jms, jme, kms, kme,                    &
    717. 

  717.                 its, ite, jts, jte, kts, kte                     )
    718. 

  718. 

  719.        IF(.not. adv_moist_cond)THEN
    720. 

  720. 

  721.         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
    722. 

  722.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
    723. 

  723.                 config_flags,                                    &
    724. 

  724.                 ids,ide, jds, jde, kds, kde,                     &
    725. 

  725.                 ims, ime, jms, jme, kms, kme,                    &
    726. 

  726.                 its, ite, jts, jte, kts, kte                     )
    727. 

  727. 

  728.        ELSE
    729. 

  729. 

  730.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    731. 

  731.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    732. 

  732.                 config_flags,                                    &
    733. 

  733.                 ids,ide, jds, jde, kds, kde,                     &
    734. 

  734.                 ims, ime, jms, jme, kms, kme,                    &
    735. 

  735.                 its, ite, jts, jte, kts, kte                     )
    736. 

  736.        ENDIF
    737. 

  737. 

  738. 

  739.       CASE (TEMFPBLSCHEME)
    740. 

  740.            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
    741. 

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

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

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

  744. 

  745.            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
    746. 

  746.                 ids,ide, jds, jde, kds, kde,                     &
    747. 

  747.                 ims, ime, jms, jme, kms, kme,                    &
    748. 

  748.                 its, ite, jts, jte, kts, kte                     )
    749. 

  749. 

  750.            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
    751. 

  751.                 ids,ide, jds, jde, kds, kde,                     &
    752. 

  752.                 ims, ime, jms, jme, kms, kme,                    &
    753. 

  753.                 its, ite, jts, jte, kts, kte                     )
    754. 

  754. 

  755.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    756. 

  756.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
    757. 

  757.                 config_flags,                                    &
    758. 

  758.                 ids,ide, jds, jde, kds, kde,                     &
    759. 

  759.                 ims, ime, jms, jme, kms, kme,                    &
    760. 

  760.                 its, ite, jts, jte, kts, kte                     )
    761. 

  761. 

  762.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    763. 

  763.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
    764. 

  764.                 config_flags,                                    &
    765. 

  765.                 ids,ide, jds, jde, kds, kde,                     &
    766. 

  766.                 ims, ime, jms, jme, kms, kme,                    &
    767. 

  767.                 its, ite, jts, jte, kts, kte                     )
    768. 

  768.      
    769. 

  769.         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
    770. 

  770.            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
    771. 

  771.                 config_flags,                                    &
    772. 

  772.                 ids,ide, jds, jde, kds, kde,                     &
    773. 

  773.                 ims, ime, jms, jme, kms, kme,                    &
    774. 

  774.                 its, ite, jts, jte, kts, kte                     )
    775. 

  775. 

  776.       CASE DEFAULT
    777. 

  777. 

  778.        print*,'phy_bl_ten: The pbl scheme does not exist'
    779. 

  779. 

  780.    END SELECT
    781. 

  781. 

  782. END SUBROUTINE phy_bl_ten
    783. 

  783. 

  784. !=================================================================
    785. 

  785. SUBROUTINE phy_cu_ten(config_flags,rk_step,n_moist,n_scalar,     &
    786. 

  786.                       rt_tendf,ru_tendf,rv_tendf,                &
    787. 

  787.                       RUCUTEN,RVCUTEN,RTHCUTEN,                  &
    788. 

  788.                       RQVCUTEN,RQCCUTEN,RQRCUTEN,                &
    789. 

  789.                       RQICUTEN,RQSCUTEN,moist_tendf,             &
    790. 

  790.                       scalar_tendf,adv_moist_cond,               &
    791. 

  791.                       ids, ide, jds, jde, kds, kde,              &
    792. 

  792.                       ims, ime, jms, jme, kms, kme,              &
    793. 

  793.                       its, ite, jts, jte, kts, kte               )
    794. 

  794. !-----------------------------------------------------------------
    795. 

  795.    IMPLICIT NONE
    796. 

  796. !-----------------------------------------------------------------
    797. 

  797.    TYPE(grid_config_rec_type  ) , INTENT(IN   ) :: config_flags
    798. 

  798. 

  799.    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde, &
    800. 

  800.                                    ims, ime, jms, jme, kms, kme, &
    801. 

  801.                                    its, ite, jts, jte, kts, kte, &
    802. 

  802.                                    n_moist, n_scalar, rk_step
    803. 

  803. 

  804.    LOGICAL , INTENT(IN)     :: adv_moist_cond
    805. 

  805. 

  806.    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_moist),         &
    807. 

  807.           INTENT(INOUT)     ::                      moist_tendf
    808. 

  808. 

  809.    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_scalar),         &
    810. 

  810.           INTENT(INOUT)     ::                      scalar_tendf
    811. 

  811. 

  812.    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
    813. 

  813.                                                         RUCUTEN, &
    814. 

  814.                                                         RVCUTEN, &
    815. 

  815.                                                        RTHCUTEN, &
    816. 

  816.                                                        RQVCUTEN, &
    817. 

  817.                                                        RQCCUTEN, &
    818. 

  818.                                                        RQRCUTEN, &
    819. 

  819.                                                        RQICUTEN, &
    820. 

  820.                                                        RQSCUTEN
    821. 

  821. 

  822.    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) ::  &
    823. 

  823.                                                        rt_tendf, &
    824. 

  824.                                                        ru_tendf, &
    825. 

  825.                                                        rv_tendf
    826. 

  826. 

  827. 

  828. ! LOCAL VARS
    829. 

  829. 

  830.    INTEGER :: i,j,k
    831. 

  831. 

  832.    SELECT CASE (config_flags%cu_physics)   
    833. 

  833. 

  834.    CASE (KFSCHEME)
    835. 

  835.         CALL add_a2a(rt_tendf,RTHCUTEN,config_flags,             &
    836. 

  836.                 ids,ide, jds, jde, kds, kde,                     &
    837. 

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

  838.                 its, ite, jts, jte, kts, kte                     )
    839. 

  839. 

  840.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    841. 

  841.         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
    842. 

  842.                 config_flags,                                    &
    843. 

  843.                 ids,ide, jds, jde, kds, kde,                     &
    844. 

  844.                 ims, ime, jms, jme, kms, kme,                    &
    845. 

  845.                 its, ite, jts, jte, kts, kte                     )
    846. 

  846. 

  847.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    848. 

  848.         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCCUTEN,     &
    849. 

  849.                 config_flags,                                    &
    850. 

  850.                 ids,ide, jds, jde, kds, kde,                     &
    851. 

  851.                 ims, ime, jms, jme, kms, kme,                    &
    852. 

  852.                 its, ite, jts, jte, kts, kte                     )
    853. 

  853. 

  854.         if (P_QR .ge. PARAM_FIRST_SCALAR)                                         &
    855. 

  855.         CALL add_a2a(moist_tendf(ims,kms,jms,P_QR),RQRCUTEN,     &
    856. 

  856.                 config_flags,                                    &
    857. 

  857.                 ids,ide, jds, jde, kds, kde,                     &
    858. 

  858.                 ims, ime, jms, jme, kms, kme,                    &
    859. 

  859.                 its, ite, jts, jte, kts, kte                     )
    860. 

  860. 

  861.         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
    862. 

  862.         CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQICUTEN,     &
    863. 

  863.                 config_flags,                                    &
    864. 

  864.                 ids,ide, jds, jde, kds, kde,                     &
    865. 

  865.                 ims, ime, jms, jme, kms, kme,                    &
    866. 

  866.                 its, ite, jts, jte, kts, kte                     )
    867. 

  867. 

  868.         if (P_QS .ge. PARAM_FIRST_SCALAR)                                         &
    869. 

  869.         CALL add_a2a(moist_tendf(ims,kms,jms,P_QS),RQSCUTEN,     &
    870. 

  870.                 config_flags,                                    &
    871. 

  871.                 ids,ide, jds, jde, kds, kde,                     &
    872. 

  872.                 ims, ime, jms, jme, kms, kme,                    &
    873. 

  873.                 its, ite, jts, jte, kts, kte                     )
    874. 

  874. 

  875.        IF(.not. adv_moist_cond)THEN
    876. 

  876. 

  877.         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
    878. 

  878.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCCUTEN, &
    879. 

  879.                 config_flags,                                    &
    880. 

  880.                 ids,ide, jds, jde, kds, kde,                     &
    881. 

  881.                 ims, ime, jms, jme, kms, kme,                    &
    882. 

  882.                 its, ite, jts, jte, kts, kte                     )
    883. 

  883. 

  884.         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
    885. 

  885.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQRCUTEN, &
    886. 

  886.                 config_flags,                                    &
    887. 

  887.                 ids,ide, jds, jde, kds, kde,                     &
    888. 

  888.                 ims, ime, jms, jme, kms, kme,                    &
    889. 

  889.                 its, ite, jts, jte, kts, kte                     )
    890. 

  890. 

  891.         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
    892. 

  892.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQICUTEN, &
    893. 

  893.                 config_flags,                                    &
    894. 

  894.                 ids,ide, jds, jde, kds, kde,                     &
    895. 

  895.                 ims, ime, jms, jme, kms, kme,                    &
    896. 

  896.                 its, ite, jts, jte, kts, kte                     )
    897. 

  897. 

  898.         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
    899. 

  899.            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQSCUTEN, &
    900. 

  900.                 config_flags,                                    &
    901. 

  901.                 ids,ide, jds, jde, kds, kde,                     &
    902. 

  902.                 ims, ime, jms, jme, kms, kme,                    &
    903. 

  903.                 its, ite, jts, jte, kts, kte                     )
    904. 

  904. 

  905.        ENDIF
    906. 

  906. 

  907.    CASE (BMJSCHEME)
    908. 

  908.         CALL add_a2a(rt_tendf,RTHCUTEN,                          &
    909. 

  909.                 config_flags,                                    &
    910. 

  910.                 ids,ide, jds, jde, kds, kde,                     &
    911. 

  911.                 ims, ime, jms, jme, kms, kme,                    &
    912. 

  912.                 its, ite, jts, jte, kts, kte                     )
    913. 

  913. 

  914.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    915. 

  915.         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
    916. 

  916.                 config_flags,                                    &
    917. 

  917.                 ids,ide, jds, jde, kds, kde,                     &
    918. 

  918.                 ims, ime, jms, jme, kms, kme,                    &
    919. 

  919.                 its, ite, jts, jte, kts, kte                     )
    920. 

  920. 

  921.    CASE (KFETASCHEME)
    922. 

  922.         CALL add_a2a(rt_tendf,RTHCUTEN,config_flags,             &
    923. 

  923.                 ids,ide, jds, jde, kds, kde,                     &
    924. 

  924.                 ims, ime, jms, jme, kms, kme,                    &
    925. 

  925.                 its, ite, jts, jte, kts, kte                     )
    926. 

  926. 

  927.         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
    928. 

  928.         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
    929. 

  929.                 config_flags,                                    &
    930. 

  930.                 ids,ide, jds, jde, kds, kde,                     &
    931. 

  931.                 ims, ime, jms, jme, kms, kme,                    &
    932. 

  932.                 its, ite, jts, jte, kts, kte                     )
    933. 

  933. 

  934.         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
    935. 

  935.         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCCUTEN,     &
    936. 

  936.                 config_flags,                                    &
    937. 

  937.                 ids,ide, jds, jde, kds, kde,                     &
    938. 

  938.       …
    939.

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