wrfv2_fire/phys/module_physics_addtendc.F

   1 !WRF:MODEL_LAYER: PHYSICS
   2 !
   3 ! note: this module really belongs in the dyn_em directory since it is
   4 !       specific only to the EM core. Leaving here for now, with an
   5 !       #if ( EM_CORE == 1 ) directive. JM 20031201
   6 !
   7
   8 !  This MODULE holds the routines which are used to perform updates of the
   9 !  model C-grid tendencies with physics A-grid tendencies
  10 !  The module consolidates code that was (up to v1.2) duplicated in
  11 !  module_em and module_rk and in
  12 !  module_big_step_utilities.F and module_big_step_utilities_em.F
  13
  14 !  This MODULE CONTAINS the following routines:
  15 !  update_phy_ten, phy_ra_ten, phy_bl_ten, phy_cu_ten, advance_ppt,
  16 !  add_a2a, add_a2c_u, and add_a2c_v
  17
  18
  19 MODULE module_physics_addtendc
  20
  21 #if ( EM_CORE == 1 )
  22
  23    USE module_state_description
  24    USE module_configure
  25
  26 CONTAINS
  27
  28 SUBROUTINE update_phy_ten(rph_tendf,rt_tendf,ru_tendf,rv_tendf,moist_tendf, &
  29                       scalar_tendf,mu_tendf,                                &
  30                       RTHRATEN,RTHBLTEN,RTHCUTEN,RTHSHTEN,                  &
  31                       RUBLTEN,RUCUTEN,RUSHTEN,                              &
  32                       RVBLTEN,RVCUTEN,RVSHTEN,                              &
  33                       RQVBLTEN,RQCBLTEN,RQIBLTEN,                           &
  34                       RQVCUTEN,RQCCUTEN,RQRCUTEN,RQICUTEN,RQSCUTEN,         &
  35                       RQVSHTEN,RQCSHTEN,RQRSHTEN,RQISHTEN,RQSSHTEN,RQGSHTEN,&
  36                       RUNDGDTEN,RVNDGDTEN,RTHNDGDTEN,RPHNDGDTEN,            &
  37                       RQVNDGDTEN,RMUNDGDTEN,                                &
  38                       rthfrten,rqvfrten,                                    & !fire
  39                       n_moist,n_scalar,config_flags,rk_step,adv_moist_cond, &
  40                       ids, ide, jds, jde, kds, kde,                         &
  41                       ims, ime, jms, jme, kms, kme,                         &
  42                       its, ite, jts, jte, kts, kte                          )
  43 !-------------------------------------------------------------------
  44    IMPLICIT NONE
  45 !-------------------------------------------------------------------
  46
  47    TYPE(grid_config_rec_type   ) ,   INTENT(IN   ) :: config_flags
  48
  49    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde,   &
  50                                    ims, ime, jms, jme, kms, kme,   &
  51                                    its, ite, jts, jte, kts, kte,   &
  52                                    n_moist,n_scalar,rk_step
  53
  54    LOGICAL , INTENT(IN)        :: adv_moist_cond
  55
  56    REAL , DIMENSION(ims:ime , kms:kme, jms:jme),INTENT(INOUT) ::   &
  57                                                          ru_tendf, &
  58                                                          rv_tendf, &
  59                                                          rt_tendf, &
  60                                                          rph_tendf
  61
  62    REAL , DIMENSION(ims:ime , jms:jme),INTENT(INOUT) ::  mu_tendf
  63
  64    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_moist),           &
  65           INTENT(INOUT)     ::                        moist_tendf
  66
  67    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_scalar),           &
  68           INTENT(INOUT)     ::                        scalar_tendf
  69
  70    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
  71                                                        RTHRATEN, &
  72                                                        RTHBLTEN, &
  73                                                        RTHCUTEN, &
  74                                                        RTHSHTEN, &
  75                                                         RUBLTEN, &
  76                                                         RUCUTEN, &
  77                                                         RUSHTEN, &
  78                                                         RVBLTEN, &
  79                                                         RVCUTEN, &
  80                                                         RVSHTEN, &
  81                                                        RQVBLTEN, &
  82                                                        RQCBLTEN, &
  83                                                        RQIBLTEN, &
  84                                                        RQVCUTEN, &
  85                                                        RQCCUTEN, &
  86                                                        RQRCUTEN, &
  87                                                        RQICUTEN, &
  88                                                        RQSCUTEN, &
  89                                                        RQVSHTEN, &
  90                                                        RQCSHTEN, &
  91                                                        RQRSHTEN, &
  92                                                        RQISHTEN, &
  93                                                        RQSSHTEN, &
  94                                                        RQGSHTEN, &
  95                                                      RTHNDGDTEN, &
  96                                                      RPHNDGDTEN, &
  97                                                      RQVNDGDTEN, &
  98                                                       RUNDGDTEN, &
  99                                                       RVNDGDTEN
 100
 101    REAL, DIMENSION(ims:ime, jms:jme), INTENT(IN  ) :: RMUNDGDTEN
 102
 103    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   & ! fire
 104                                                        rthfrten, &
 105                                                        rqvfrten
 106 !------------------------------------------------------------------
 107 !  set up loop bounds for this grid's boundary conditions
 108
 109    if (config_flags%ra_lw_physics .gt. 0 .or.                  &
 110        config_flags%ra_sw_physics .gt. 0)                      &
 111       CALL phy_ra_ten(config_flags,rt_tendf,RTHRATEN,          &
 112                       ids, ide, jds, jde, kds, kde,            &
 113                       ims, ime, jms, jme, kms, kme,            &
 114                       its, ite, jts, jte, kts, kte             )
 115
 116    if (config_flags%bl_pbl_physics .gt. 0)                     &
 117       CALL phy_bl_ten(config_flags,rk_step,n_moist,n_scalar,   &
 118                       rt_tendf,ru_tendf,rv_tendf,moist_tendf,  &
 119                       scalar_tendf,adv_moist_cond,             &
 120                       RTHBLTEN,RUBLTEN,RVBLTEN,                &
 121                       RQVBLTEN,RQCBLTEN,RQIBLTEN,              &
 122                       ids, ide, jds, jde, kds, kde,            &
 123                       ims, ime, jms, jme, kms, kme,            &
 124                       its, ite, jts, jte, kts, kte             )
 125
 126    if (config_flags%cu_physics .gt. 0)                         &
 127       CALL phy_cu_ten(config_flags,rk_step,n_moist,n_scalar,   &
 128                       rt_tendf,ru_tendf,rv_tendf,              &
 129                       RUCUTEN,RVCUTEN,RTHCUTEN,                &
 130                       RQVCUTEN,RQCCUTEN,RQRCUTEN,              &
 131                       RQICUTEN,RQSCUTEN,moist_tendf,           &
 132                       scalar_tendf,adv_moist_cond,             &
 133                       ids, ide, jds, jde, kds, kde,            &
 134                       ims, ime, jms, jme, kms, kme,            &
 135                       its, ite, jts, jte, kts, kte             )
 136
 137    if (config_flags%shcu_physics .gt. 0)                       &
 138       CALL phy_shcu_ten(config_flags,rk_step,n_moist,          &
 139                       rt_tendf,ru_tendf,rv_tendf,              &
 140                       RUSHTEN,RVSHTEN,RTHSHTEN,                &
 141                       RQVSHTEN,RQCSHTEN,RQRSHTEN,              &
 142                       RQISHTEN,RQSSHTEN,RQGSHTEN,moist_tendf,  &
 143                       ids, ide, jds, jde, kds, kde,            &
 144                       ims, ime, jms, jme, kms, kme,            &
 145                       its, ite, jts, jte, kts, kte             )
 146
 147    if (config_flags%grid_fdda .gt. 0)                          &
 148       CALL phy_fg_ten(config_flags,rk_step,n_moist,            &
 149                       rph_tendf,rt_tendf,ru_tendf,rv_tendf,    &
 150                       mu_tendf, moist_tendf,                   &
 151                       RUNDGDTEN,RVNDGDTEN,RTHNDGDTEN,          &
 152                       RPHNDGDTEN,RQVNDGDTEN,RMUNDGDTEN,        &
 153                       ids, ide, jds, jde, kds, kde,            &
 154                       ims, ime, jms, jme, kms, kme,            &
 155                       its, ite, jts, jte, kts, kte             )
 156
 157    if (config_flags%ifire .gt. 0)                              & ! fire
 158       CALL phy_fr_ten(config_flags,rk_step,n_moist,            &
 159                       rt_tendf,ru_tendf,rv_tendf,              &
 160                       mu_tendf, moist_tendf,                   &
 161                       rthfrten,rqvfrten,                       &
 162                       ids, ide, jds, jde, kds, kde,            &
 163                       ims, ime, jms, jme, kms, kme,            &
 164                       its, ite, jts, jte, kts, kte             )
 165
 166 END SUBROUTINE update_phy_ten
 167
 168 !=================================================================
 169 SUBROUTINE phy_ra_ten(config_flags,rt_tendf,RTHRATEN,            &
 170                       ids, ide, jds, jde, kds, kde,              &
 171                       ims, ime, jms, jme, kms, kme,              &
 172                       its, ite, jts, jte, kts, kte               )
 173 !-----------------------------------------------------------------
 174    IMPLICIT NONE
 175 !-----------------------------------------------------------------
 176    TYPE(grid_config_rec_type  ) , INTENT(IN   ) :: config_flags
 177
 178    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde, &
 179                                    ims, ime, jms, jme, kms, kme, &
 180                                    its, ite, jts, jte, kts, kte
 181
 182    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
 183                                                        RTHRATEN
 184
 185    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) ::  &
 186                                                        rt_tendf
 187
 188 ! LOCAL VARS
 189
 190    INTEGER :: i,j,k
 191
 192    CALL add_a2a(rt_tendf,RTHRATEN,config_flags,                  &
 193                 ids,ide, jds, jde, kds, kde,                     &
 194                 ims, ime, jms, jme, kms, kme,                    &
 195                 its, ite, jts, jte, kts, kte                     )
 196
 197 END SUBROUTINE phy_ra_ten
 198
 199 !=================================================================
 200 SUBROUTINE phy_bl_ten(config_flags,rk_step,n_moist,n_scalar,     &
 201                       rt_tendf,ru_tendf,rv_tendf,moist_tendf,    &
 202                       scalar_tendf,adv_moist_cond,               &
 203                       RTHBLTEN,RUBLTEN,RVBLTEN,                  &
 204                       RQVBLTEN,RQCBLTEN,RQIBLTEN,                &
 205                       ids, ide, jds, jde, kds, kde,              &
 206                       ims, ime, jms, jme, kms, kme,              &
 207                       its, ite, jts, jte, kts, kte               )
 208 !-----------------------------------------------------------------
 209    IMPLICIT NONE
 210 !-----------------------------------------------------------------
 211    TYPE(grid_config_rec_type) ,     INTENT(IN   ) :: config_flags
 212
 213    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde, &
 214                                    ims, ime, jms, jme, kms, kme, &
 215                                    its, ite, jts, jte, kts, kte, &
 216                                    n_moist, n_scalar, rk_step
 217
 218    LOGICAL , INTENT(IN)     :: adv_moist_cond
 219
 220    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_moist),         &
 221           INTENT(INOUT)     ::                      moist_tendf
 222
 223    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_scalar),         &
 224           INTENT(INOUT)     ::                      scalar_tendf
 225
 226    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
 227                                                        RTHBLTEN, &
 228                                                         RUBLTEN, &
 229                                                         RVBLTEN, &
 230                                                        RQVBLTEN, &
 231                                                        RQCBLTEN, &
 232                                                        RQIBLTEN
 233
 234    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT)::   &
 235                                                        rt_tendf, &
 236                                                        ru_tendf, &
 237                                                        rv_tendf
 238 ! LOCAL VARS
 239
 240    INTEGER :: i,j,k,IBGN,IEND,JBGN,JEND
 241
 242 !-----------------------------------------------------------------
 243
 244    SELECT CASE(config_flags%bl_pbl_physics)
 245
 246       CASE (YSUSCHEME)
 247
 248            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 249                 ids,ide, jds, jde, kds, kde,                     &
 250                 ims, ime, jms, jme, kms, kme,                    &
 251                 its, ite, jts, jte, kts, kte                     )
 252
 253            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 254                 ids,ide, jds, jde, kds, kde,                     &
 255                 ims, ime, jms, jme, kms, kme,                    &
 256                 its, ite, jts, jte, kts, kte                     )
 257
 258            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 259                 ids,ide, jds, jde, kds, kde,                     &
 260                 ims, ime, jms, jme, kms, kme,                    &
 261                 its, ite, jts, jte, kts, kte                     )
 262
 263         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 264            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 265                 config_flags,                                    &
 266                 ids,ide, jds, jde, kds, kde,                     &
 267                 ims, ime, jms, jme, kms, kme,                    &
 268                 its, ite, jts, jte, kts, kte                     )
 269
 270         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 271            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 272                 config_flags,                                    &
 273                 ids,ide, jds, jde, kds, kde,                     &
 274                 ims, ime, jms, jme, kms, kme,                    &
 275                 its, ite, jts, jte, kts, kte                     )
 276
 277         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
 278            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
 279                 config_flags,                                    &
 280                 ids,ide, jds, jde, kds, kde,                     &
 281                 ims, ime, jms, jme, kms, kme,                    &
 282                 its, ite, jts, jte, kts, kte                     )
 283
 284        IF(.not. adv_moist_cond)THEN
 285
 286         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 287            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
 288                 config_flags,                                    &
 289                 ids,ide, jds, jde, kds, kde,                     &
 290                 ims, ime, jms, jme, kms, kme,                    &
 291                 its, ite, jts, jte, kts, kte                     )
 292
 293         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 294            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
 295                 config_flags,                                    &
 296                 ids,ide, jds, jde, kds, kde,                     &
 297                 ims, ime, jms, jme, kms, kme,                    &
 298                 its, ite, jts, jte, kts, kte                     )
 299        ENDIF
 300
 301       CASE (MRFSCHEME)
 302
 303            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 304                 ids,ide, jds, jde, kds, kde,                     &
 305                 ims, ime, jms, jme, kms, kme,                    &
 306                 its, ite, jts, jte, kts, kte                     )
 307
 308            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 309                 ids,ide, jds, jde, kds, kde,                     &
 310                 ims, ime, jms, jme, kms, kme,                    &
 311                 its, ite, jts, jte, kts, kte                     )
 312
 313            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 314                 ids,ide, jds, jde, kds, kde,                     &
 315                 ims, ime, jms, jme, kms, kme,                    &
 316                 its, ite, jts, jte, kts, kte                     )
 317
 318         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 319            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 320                 config_flags,                                    &
 321                 ids,ide, jds, jde, kds, kde,                     &
 322                 ims, ime, jms, jme, kms, kme,                    &
 323                 its, ite, jts, jte, kts, kte                     )
 324
 325         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 326            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 327                 config_flags,                                    &
 328                 ids,ide, jds, jde, kds, kde,                     &
 329                 ims, ime, jms, jme, kms, kme,                    &
 330                 its, ite, jts, jte, kts, kte                     )
 331
 332         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
 333            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
 334                 config_flags,                                    &
 335                 ids,ide, jds, jde, kds, kde,                     &
 336                 ims, ime, jms, jme, kms, kme,                    &
 337                 its, ite, jts, jte, kts, kte                     )
 338
 339        IF(.not. adv_moist_cond)THEN
 340
 341         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 342            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
 343                 config_flags,                                    &
 344                 ids,ide, jds, jde, kds, kde,                     &
 345                 ims, ime, jms, jme, kms, kme,                    &
 346                 its, ite, jts, jte, kts, kte                     )
 347
 348         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 349            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
 350                 config_flags,                                    &
 351                 ids,ide, jds, jde, kds, kde,                     &
 352                 ims, ime, jms, jme, kms, kme,                    &
 353                 its, ite, jts, jte, kts, kte                     )
 354        ENDIF
 355
 356       CASE (ACMPBLSCHEME)
 357
 358            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 359                 ids,ide, jds, jde, kds, kde,                     &
 360                 ims, ime, jms, jme, kms, kme,                    &
 361                 its, ite, jts, jte, kts, kte                     )
 362
 363            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 364                 ids,ide, jds, jde, kds, kde,                     &
 365                 ims, ime, jms, jme, kms, kme,                    &
 366                 its, ite, jts, jte, kts, kte                     )
 367
 368            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 369                 ids,ide, jds, jde, kds, kde,                     &
 370                 ims, ime, jms, jme, kms, kme,                    &
 371                 its, ite, jts, jte, kts, kte                     )
 372
 373         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 374            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 375                 config_flags,                                    &
 376                 ids,ide, jds, jde, kds, kde,                     &
 377                 ims, ime, jms, jme, kms, kme,                    &
 378                 its, ite, jts, jte, kts, kte                     )
 379
 380         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 381            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 382                 config_flags,                                    &
 383                 ids,ide, jds, jde, kds, kde,                     &
 384                 ims, ime, jms, jme, kms, kme,                    &
 385                 its, ite, jts, jte, kts, kte                     )
 386
 387         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
 388            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
 389                 config_flags,                                    &
 390                 ids,ide, jds, jde, kds, kde,                     &
 391                 ims, ime, jms, jme, kms, kme,                    &
 392                 its, ite, jts, jte, kts, kte                     )
 393
 394        IF(.not. adv_moist_cond)THEN
 395
 396         if (P_QT .ge. PARAM_FIRST_SCALAR)THEN
 397            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
 398                 config_flags,                                    &
 399                 ids,ide, jds, jde, kds, kde,                     &
 400                 ims, ime, jms, jme, kms, kme,                    &
 401                 its, ite, jts, jte, kts, kte                     )
 402
 403            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
 404                 config_flags,                                    &
 405                 ids,ide, jds, jde, kds, kde,                     &
 406                 ims, ime, jms, jme, kms, kme,                    &
 407                 its, ite, jts, jte, kts, kte                     )
 408         ENDIF
 409
 410        ENDIF
 411
 412       CASE (MYJPBLSCHEME)
 413
 414            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 415                 ids,ide, jds, jde, kds, kde,                     &
 416                 ims, ime, jms, jme, kms, kme,                    &
 417                 its, ite, jts, jte, kts, kte                     )
 418
 419            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 420                 ids,ide, jds, jde, kds, kde,                     &
 421                 ims, ime, jms, jme, kms, kme,                    &
 422                 its, ite, jts, jte, kts, kte                     )
 423
 424            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 425                 ids,ide, jds, jde, kds, kde,                     &
 426                 ims, ime, jms, jme, kms, kme,                    &
 427                 its, ite, jts, jte, kts, kte                     )
 428
 429         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 430            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 431                 config_flags,                                    &
 432                 ids,ide, jds, jde, kds, kde,                     &
 433                 ims, ime, jms, jme, kms, kme,                    &
 434                 its, ite, jts, jte, kts, kte                     )
 435
 436        IF(.not. adv_moist_cond)THEN
 437
 438         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 439            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
 440                 config_flags,                                    &
 441                 ids,ide, jds, jde, kds, kde,                     &
 442                 ims, ime, jms, jme, kms, kme,                    &
 443                 its, ite, jts, jte, kts, kte                     )
 444
 445         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 446            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
 447                 config_flags,                                    &
 448                 ids,ide, jds, jde, kds, kde,                     &
 449                 ims, ime, jms, jme, kms, kme,                    &
 450                 its, ite, jts, jte, kts, kte                     )
 451
 452 !       if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 453 !          CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQSBLTEN,  &
 454 !               config_flags,                                    &
 455 !               ids,ide, jds, jde, kds, kde,                     &
 456 !               ims, ime, jms, jme, kms, kme,                    &
 457 !               its, ite, jts, jte, kts, kte                     )
 458 !
 459 !       if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 460 !          CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQRBLTEN,  &
 461 !               config_flags,                                    &
 462 !               ids,ide, jds, jde, kds, kde,                     &
 463 !               ims, ime, jms, jme, kms, kme,                    &
 464 !               its, ite, jts, jte, kts, kte                     )
 465 !
 466 !       if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 467 !          CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQGBLTEN,  &
 468 !               config_flags,                                    &
 469 !               ids,ide, jds, jde, kds, kde,                     &
 470 !               ims, ime, jms, jme, kms, kme,                    &
 471 !               its, ite, jts, jte, kts, kte                     )
 472
 473        ELSE
 474
 475         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 476            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 477                 config_flags,                                    &
 478                 ids,ide, jds, jde, kds, kde,                     &
 479                 ims, ime, jms, jme, kms, kme,                    &
 480                 its, ite, jts, jte, kts, kte                     )
 481
 482         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
 483            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
 484                 config_flags,                                    &
 485                 ids,ide, jds, jde, kds, kde,                     &
 486                 ims, ime, jms, jme, kms, kme,                    &
 487                 its, ite, jts, jte, kts, kte                     )
 488
 489 !       if (P_QS .ge. PARAM_FIRST_SCALAR)                                         &
 490 !          CALL add_a2a(moist_tendf(ims,kms,jms,P_QS),RQSBLTEN,  &
 491 !               config_flags,                                    &
 492 !               ids,ide, jds, jde, kds, kde,                     &
 493 !               ims, ime, jms, jme, kms, kme,                    &
 494 !               its, ite, jts, jte, kts, kte                     )
 495 !
 496 !       if (P_QR .ge. PARAM_FIRST_SCALAR)                                         &
 497 !          CALL add_a2a(moist_tendf(ims,kms,jms,P_QR),RQRBLTEN,  &
 498 !               config_flags,                                    &
 499 !               ids,ide, jds, jde, kds, kde,                     &
 500 !               ims, ime, jms, jme, kms, kme,                    &
 501 !               its, ite, jts, jte, kts, kte                     )
 502 !
 503 !       if (P_QG .ge. PARAM_FIRST_SCALAR)                                         &
 504 !          CALL add_a2a(moist_tendf(ims,kms,jms,P_QG),RQGBLTEN,  &
 505 !               config_flags,                                    &
 506 !               ids,ide, jds, jde, kds, kde,                     &
 507 !               ims, ime, jms, jme, kms, kme,                    &
 508 !               its, ite, jts, jte, kts, kte                     )
 509
 510        ENDIF
 511
 512       CASE (QNSEPBLSCHEME)
 513
 514            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 515                 ids,ide, jds, jde, kds, kde,                     &
 516                 ims, ime, jms, jme, kms, kme,                    &
 517                 its, ite, jts, jte, kts, kte                     )
 518
 519            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 520                 ids,ide, jds, jde, kds, kde,                     &
 521                 ims, ime, jms, jme, kms, kme,                    &
 522                 its, ite, jts, jte, kts, kte                     )
 523
 524            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 525                 ids,ide, jds, jde, kds, kde,                     &
 526                 ims, ime, jms, jme, kms, kme,                    &
 527                 its, ite, jts, jte, kts, kte                     )
 528
 529         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 530            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 531                 config_flags,                                    &
 532                 ids,ide, jds, jde, kds, kde,                     &
 533                 ims, ime, jms, jme, kms, kme,                    &
 534                 its, ite, jts, jte, kts, kte                     )
 535
 536        IF(.not. adv_moist_cond)THEN
 537
 538         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 539            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
 540                 config_flags,                                    &
 541                 ids,ide, jds, jde, kds, kde,                     &
 542                 ims, ime, jms, jme, kms, kme,                    &
 543                 its, ite, jts, jte, kts, kte                     )
 544
 545        ELSE
 546
 547         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 548            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 549                 config_flags,                                    &
 550                 ids,ide, jds, jde, kds, kde,                     &
 551                 ims, ime, jms, jme, kms, kme,                    &
 552                 its, ite, jts, jte, kts, kte                     )
 553
 554        ENDIF
 555
 556       CASE (GFSSCHEME)
 557
 558            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 559                 ids,ide, jds, jde, kds, kde,                     &
 560                 ims, ime, jms, jme, kms, kme,                    &
 561                 its, ite, jts, jte, kts, kte                     )
 562
 563            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 564                 ids,ide, jds, jde, kds, kde,                     &
 565                 ims, ime, jms, jme, kms, kme,                    &
 566                 its, ite, jts, jte, kts, kte                     )
 567
 568            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 569                 ids,ide, jds, jde, kds, kde,                     &
 570                 ims, ime, jms, jme, kms, kme,                    &
 571                 its, ite, jts, jte, kts, kte                     )
 572
 573         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 574            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 575                 config_flags,                                    &
 576                 ids,ide, jds, jde, kds, kde,                     &
 577                 ims, ime, jms, jme, kms, kme,                    &
 578                 its, ite, jts, jte, kts, kte                     )
 579
 580         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 581            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 582                 config_flags,                                    &
 583                 ids,ide, jds, jde, kds, kde,                     &
 584                 ims, ime, jms, jme, kms, kme,                    &
 585                 its, ite, jts, jte, kts, kte                     )
 586
 587         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
 588            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
 589                 config_flags,                                    &
 590                 ids,ide, jds, jde, kds, kde,                     &
 591                 ims, ime, jms, jme, kms, kme,                    &
 592                 its, ite, jts, jte, kts, kte                     )
 593
 594        IF(.not. adv_moist_cond)THEN
 595
 596         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 597            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
 598                 config_flags,                                    &
 599                 ids,ide, jds, jde, kds, kde,                     &
 600                 ims, ime, jms, jme, kms, kme,                    &
 601                 its, ite, jts, jte, kts, kte                     )
 602
 603         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 604            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQIBLTEN,  &
 605                 config_flags,                                    &
 606                 ids,ide, jds, jde, kds, kde,                     &
 607                 ims, ime, jms, jme, kms, kme,                    &
 608                 its, ite, jts, jte, kts, kte                     )
 609        ENDIF
 610
 611       CASE (MYNNPBLSCHEME2,MYNNPBLSCHEME3)
 612
 613            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 614                 ids,ide, jds, jde, kds, kde,                     &
 615                 ims, ime, jms, jme, kms, kme,                    &
 616                 its, ite, jts, jte, kts, kte                     )
 617
 618            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 619                 ids,ide, jds, jde, kds, kde,                     &
 620                 ims, ime, jms, jme, kms, kme,                    &
 621                 its, ite, jts, jte, kts, kte                     )
 622
 623            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 624                 ids,ide, jds, jde, kds, kde,                     &
 625                 ims, ime, jms, jme, kms, kme,                    &
 626                 its, ite, jts, jte, kts, kte                     )
 627
 628         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 629            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 630                 config_flags,                                    &
 631                 ids,ide, jds, jde, kds, kde,                     &
 632                 ims, ime, jms, jme, kms, kme,                    &
 633                 its, ite, jts, jte, kts, kte                     )
 634
 635         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 636            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 637                 config_flags,                                    &
 638                 ids,ide, jds, jde, kds, kde,                     &
 639                 ims, ime, jms, jme, kms, kme,                    &
 640                 its, ite, jts, jte, kts, kte                     )
 641
 642        IF(.not. adv_moist_cond)THEN
 643
 644         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 645            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
 646                 config_flags,                                    &
 647                 ids,ide, jds, jde, kds, kde,                     &
 648                 ims, ime, jms, jme, kms, kme,                    &
 649                 its, ite, jts, jte, kts, kte                     )
 650
 651        ENDIF
 652
 653        CASE (BOULACSCHEME)
 654
 655            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 656                 ids,ide, jds, jde, kds, kde,                     &
 657                 ims, ime, jms, jme, kms, kme,                    &
 658                 its, ite, jts, jte, kts, kte                     )
 659
 660            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 661                 ids,ide, jds, jde, kds, kde,                     &
 662                 ims, ime, jms, jme, kms, kme,                    &
 663                 its, ite, jts, jte, kts, kte                     )
 664
 665            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 666                 ids,ide, jds, jde, kds, kde,                     &
 667                 ims, ime, jms, jme, kms, kme,                    &
 668                 its, ite, jts, jte, kts, kte                     )
 669
 670         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 671            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 672                 config_flags,                                    &
 673                 ids,ide, jds, jde, kds, kde,                     &
 674                 ims, ime, jms, jme, kms, kme,                    &
 675                 its, ite, jts, jte, kts, kte                     )
 676
 677        IF(.not. adv_moist_cond)THEN
 678
 679         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 680            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
 681                 config_flags,                                    &
 682                 ids,ide, jds, jde, kds, kde,                     &
 683                 ims, ime, jms, jme, kms, kme,                    &
 684                 its, ite, jts, jte, kts, kte                     )
 685
 686        ELSE
 687
 688         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 689            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 690                 config_flags,                                    &
 691                 ids,ide, jds, jde, kds, kde,                     &
 692                 ims, ime, jms, jme, kms, kme,                    &
 693                 its, ite, jts, jte, kts, kte                     )
 694        ENDIF
 695
 696        CASE (CAMUWPBLSCHEME)
 697            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 698                 ids,ide, jds, jde, kds, kde,                     &
 699                 ims, ime, jms, jme, kms, kme,                    &
 700                 its, ite, jts, jte, kts, kte                     )
 701
 702            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 703                 ids,ide, jds, jde, kds, kde,                     &
 704                 ims, ime, jms, jme, kms, kme,                    &
 705                 its, ite, jts, jte, kts, kte                     )
 706
 707            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 708                 ids,ide, jds, jde, kds, kde,                     &
 709                 ims, ime, jms, jme, kms, kme,                    &
 710                 its, ite, jts, jte, kts, kte                     )
 711
 712         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 713            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 714                 config_flags,                                    &
 715                 ids,ide, jds, jde, kds, kde,                     &
 716                 ims, ime, jms, jme, kms, kme,                    &
 717                 its, ite, jts, jte, kts, kte                     )
 718
 719        IF(.not. adv_moist_cond)THEN
 720
 721         if (P_QT .ge. PARAM_FIRST_SCALAR)                                         &
 722            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCBLTEN,  &
 723                 config_flags,                                    &
 724                 ids,ide, jds, jde, kds, kde,                     &
 725                 ims, ime, jms, jme, kms, kme,                    &
 726                 its, ite, jts, jte, kts, kte                     )
 727
 728        ELSE
 729
 730         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 731            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 732                 config_flags,                                    &
 733                 ids,ide, jds, jde, kds, kde,                     &
 734                 ims, ime, jms, jme, kms, kme,                    &
 735                 its, ite, jts, jte, kts, kte                     )
 736        ENDIF
 737
 738
 739       CASE (TEMFPBLSCHEME)
 740            CALL add_a2a(rt_tendf,RTHBLTEN,config_flags,          &
 741                 ids,ide, jds, jde, kds, kde,                     &
 742                 ims, ime, jms, jme, kms, kme,                    &
 743                 its, ite, jts, jte, kts, kte                     )
 744
 745            CALL add_a2c_u(ru_tendf,RUBLTEN,config_flags,         &
 746                 ids,ide, jds, jde, kds, kde,                     &
 747                 ims, ime, jms, jme, kms, kme,                    &
 748                 its, ite, jts, jte, kts, kte                     )
 749
 750            CALL add_a2c_v(rv_tendf,RVBLTEN,config_flags,         &
 751                 ids,ide, jds, jde, kds, kde,                     &
 752                 ims, ime, jms, jme, kms, kme,                    &
 753                 its, ite, jts, jte, kts, kte                     )
 754
 755         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 756            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVBLTEN,  &
 757                 config_flags,                                    &
 758                 ids,ide, jds, jde, kds, kde,                     &
 759                 ims, ime, jms, jme, kms, kme,                    &
 760                 its, ite, jts, jte, kts, kte                     )
 761
 762         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 763            CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCBLTEN,  &
 764                 config_flags,                                    &
 765                 ids,ide, jds, jde, kds, kde,                     &
 766                 ims, ime, jms, jme, kms, kme,                    &
 767                 its, ite, jts, jte, kts, kte                     )
 768
 769         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
 770            CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQIBLTEN,  &
 771                 config_flags,                                    &
 772                 ids,ide, jds, jde, kds, kde,                     &
 773                 ims, ime, jms, jme, kms, kme,                    &
 774                 its, ite, jts, jte, kts, kte                     )
 775
 776       CASE DEFAULT
 777
 778        print*,'phy_bl_ten: The pbl scheme does not exist'
 779
 780    END SELECT
 781
 782 END SUBROUTINE phy_bl_ten
 783
 784 !=================================================================
 785 SUBROUTINE phy_cu_ten(config_flags,rk_step,n_moist,n_scalar,     &
 786                       rt_tendf,ru_tendf,rv_tendf,                &
 787                       RUCUTEN,RVCUTEN,RTHCUTEN,                  &
 788                       RQVCUTEN,RQCCUTEN,RQRCUTEN,                &
 789                       RQICUTEN,RQSCUTEN,moist_tendf,             &
 790                       scalar_tendf,adv_moist_cond,               &
 791                       ids, ide, jds, jde, kds, kde,              &
 792                       ims, ime, jms, jme, kms, kme,              &
 793                       its, ite, jts, jte, kts, kte               )
 794 !-----------------------------------------------------------------
 795    IMPLICIT NONE
 796 !-----------------------------------------------------------------
 797    TYPE(grid_config_rec_type  ) , INTENT(IN   ) :: config_flags
 798
 799    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde, &
 800                                    ims, ime, jms, jme, kms, kme, &
 801                                    its, ite, jts, jte, kts, kte, &
 802                                    n_moist, n_scalar, rk_step
 803
 804    LOGICAL , INTENT(IN)     :: adv_moist_cond
 805
 806    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_moist),         &
 807           INTENT(INOUT)     ::                      moist_tendf
 808
 809    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_scalar),         &
 810           INTENT(INOUT)     ::                      scalar_tendf
 811
 812    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
 813                                                         RUCUTEN, &
 814                                                         RVCUTEN, &
 815                                                        RTHCUTEN, &
 816                                                        RQVCUTEN, &
 817                                                        RQCCUTEN, &
 818                                                        RQRCUTEN, &
 819                                                        RQICUTEN, &
 820                                                        RQSCUTEN
 821
 822    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) ::  &
 823                                                        rt_tendf, &
 824                                                        ru_tendf, &
 825                                                        rv_tendf
 826
 827
 828 ! LOCAL VARS
 829
 830    INTEGER :: i,j,k
 831
 832    SELECT CASE (config_flags%cu_physics)
 833
 834    CASE (KFSCHEME)
 835         CALL add_a2a(rt_tendf,RTHCUTEN,config_flags,             &
 836                 ids,ide, jds, jde, kds, kde,                     &
 837                 ims, ime, jms, jme, kms, kme,                    &
 838                 its, ite, jts, jte, kts, kte                     )
 839
 840         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 841         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
 842                 config_flags,                                    &
 843                 ids,ide, jds, jde, kds, kde,                     &
 844                 ims, ime, jms, jme, kms, kme,                    &
 845                 its, ite, jts, jte, kts, kte                     )
 846
 847         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 848         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCCUTEN,     &
 849                 config_flags,                                    &
 850                 ids,ide, jds, jde, kds, kde,                     &
 851                 ims, ime, jms, jme, kms, kme,                    &
 852                 its, ite, jts, jte, kts, kte                     )
 853
 854         if (P_QR .ge. PARAM_FIRST_SCALAR)                                         &
 855         CALL add_a2a(moist_tendf(ims,kms,jms,P_QR),RQRCUTEN,     &
 856                 config_flags,                                    &
 857                 ids,ide, jds, jde, kds, kde,                     &
 858                 ims, ime, jms, jme, kms, kme,                    &
 859                 its, ite, jts, jte, kts, kte                     )
 860
 861         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
 862         CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQICUTEN,     &
 863                 config_flags,                                    &
 864                 ids,ide, jds, jde, kds, kde,                     &
 865                 ims, ime, jms, jme, kms, kme,                    &
 866                 its, ite, jts, jte, kts, kte                     )
 867
 868         if (P_QS .ge. PARAM_FIRST_SCALAR)                                         &
 869         CALL add_a2a(moist_tendf(ims,kms,jms,P_QS),RQSCUTEN,     &
 870                 config_flags,                                    &
 871                 ids,ide, jds, jde, kds, kde,                     &
 872                 ims, ime, jms, jme, kms, kme,                    &
 873                 its, ite, jts, jte, kts, kte                     )
 874
 875        IF(.not. adv_moist_cond)THEN
 876
 877         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
 878            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCCUTEN, &
 879                 config_flags,                                    &
 880                 ids,ide, jds, jde, kds, kde,                     &
 881                 ims, ime, jms, jme, kms, kme,                    &
 882                 its, ite, jts, jte, kts, kte                     )
 883
 884         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
 885            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQRCUTEN, &
 886                 config_flags,                                    &
 887                 ids,ide, jds, jde, kds, kde,                     &
 888                 ims, ime, jms, jme, kms, kme,                    &
 889                 its, ite, jts, jte, kts, kte                     )
 890
 891         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
 892            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQICUTEN, &
 893                 config_flags,                                    &
 894                 ids,ide, jds, jde, kds, kde,                     &
 895                 ims, ime, jms, jme, kms, kme,                    &
 896                 its, ite, jts, jte, kts, kte                     )
 897
 898         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
 899            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQSCUTEN, &
 900                 config_flags,                                    &
 901                 ids,ide, jds, jde, kds, kde,                     &
 902                 ims, ime, jms, jme, kms, kme,                    &
 903                 its, ite, jts, jte, kts, kte                     )
 904
 905        ENDIF
 906
 907    CASE (BMJSCHEME)
 908         CALL add_a2a(rt_tendf,RTHCUTEN,                          &
 909                 config_flags,                                    &
 910                 ids,ide, jds, jde, kds, kde,                     &
 911                 ims, ime, jms, jme, kms, kme,                    &
 912                 its, ite, jts, jte, kts, kte                     )
 913
 914         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 915         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
 916                 config_flags,                                    &
 917                 ids,ide, jds, jde, kds, kde,                     &
 918                 ims, ime, jms, jme, kms, kme,                    &
 919                 its, ite, jts, jte, kts, kte                     )
 920
 921    CASE (KFETASCHEME)
 922         CALL add_a2a(rt_tendf,RTHCUTEN,config_flags,             &
 923                 ids,ide, jds, jde, kds, kde,                     &
 924                 ims, ime, jms, jme, kms, kme,                    &
 925                 its, ite, jts, jte, kts, kte                     )
 926
 927         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
 928         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
 929                 config_flags,                                    &
 930                 ids,ide, jds, jde, kds, kde,                     &
 931                 ims, ime, jms, jme, kms, kme,                    &
 932                 its, ite, jts, jte, kts, kte                     )
 933
 934         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
 935         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCCUTEN,     &
 936                 config_flags,                                    &
 937                 ids,ide, jds, jde, kds, kde,                     &
 938                 ims, ime, jms, jme, kms, kme,                    &
 939                 its, ite, jts, jte, kts, kte                     )
 940
 941         if (P_QR .ge. PARAM_FIRST_SCALAR)                                         &
 942         CALL add_a2a(moist_tendf(ims,kms,jms,P_QR),RQRCUTEN,     &
 943                 config_flags,                                    &
 944                 ids,ide, jds, jde, kds, kde,                     &
 945                 ims, ime, jms, jme, kms, kme,                    &
 946                 its, ite, jts, jte, kts, kte                     )
 947
 948         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
 949         CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQICUTEN,     &
 950                 config_flags,                                    &
 951                 ids,ide, jds, jde, kds, kde,                     &
 952                 ims, ime, jms, jme, kms, kme,                    &
 953                 its, ite, jts, jte, kts, kte                     )
 954
 955         if (P_QS .ge. PARAM_FIRST_SCALAR)                                         &
 956         CALL add_a2a(moist_tendf(ims,kms,jms,P_QS),RQSCUTEN,     &
 957                 config_flags,                                    &
 958                 ids,ide, jds, jde, kds, kde,                     &
 959                 ims, ime, jms, jme, kms, kme,                    &
 960                 its, ite, jts, jte, kts, kte                     )
 961
 962        IF(.not. adv_moist_cond)THEN
 963
 964         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
 965            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCCUTEN, &
 966                 config_flags,                                    &
 967                 ids,ide, jds, jde, kds, kde,                     &
 968                 ims, ime, jms, jme, kms, kme,                    &
 969                 its, ite, jts, jte, kts, kte                     )
 970
 971         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
 972            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQRCUTEN, &
 973                 config_flags,                                    &
 974                 ids,ide, jds, jde, kds, kde,                     &
 975                 ims, ime, jms, jme, kms, kme,                    &
 976                 its, ite, jts, jte, kts, kte                     )
 977
 978         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
 979            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQICUTEN, &
 980                 config_flags,                                    &
 981                 ids,ide, jds, jde, kds, kde,                     &
 982                 ims, ime, jms, jme, kms, kme,                    &
 983                 its, ite, jts, jte, kts, kte                     )
 984
 985         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
 986            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQSCUTEN, &
 987                 config_flags,                                    &
 988                 ids,ide, jds, jde, kds, kde,                     &
 989                 ims, ime, jms, jme, kms, kme,                    &
 990                 its, ite, jts, jte, kts, kte                     )
 991
 992        ENDIF
 993    CASE (GDSCHEME, G3SCHEME)
 994         CALL add_a2a(rt_tendf,RTHCUTEN,config_flags,             &
 995                 ids,ide, jds, jde, kds, kde,                     &
 996                 ims, ime, jms, jme, kms, kme,                    &
 997                 its, ite, jts, jte, kts, kte                     )
 998
 999         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
1000         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
1001                 config_flags,                                    &
1002                 ids,ide, jds, jde, kds, kde,                     &
1003                 ims, ime, jms, jme, kms, kme,                    &
1004                 its, ite, jts, jte, kts, kte                     )
1005
1006         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
1007         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCCUTEN,     &
1008                 config_flags,                                    &
1009                 ids,ide, jds, jde, kds, kde,                     &
1010                 ims, ime, jms, jme, kms, kme,                    &
1011                 its, ite, jts, jte, kts, kte                     )
1012
1013         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
1014         CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQICUTEN,     &
1015                 config_flags,                                    &
1016                 ids,ide, jds, jde, kds, kde,                     &
1017                 ims, ime, jms, jme, kms, kme,                    &
1018                 its, ite, jts, jte, kts, kte                     )
1019
1020        IF(.not. adv_moist_cond)THEN
1021
1022         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1023            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCCUTEN, &
1024                 config_flags,                                    &
1025                 ids,ide, jds, jde, kds, kde,                     &
1026                 ims, ime, jms, jme, kms, kme,                    &
1027                 its, ite, jts, jte, kts, kte                     )
1028
1029         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1030            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQICUTEN, &
1031                 config_flags,                                    &
1032                 ids,ide, jds, jde, kds, kde,                     &
1033                 ims, ime, jms, jme, kms, kme,                    &
1034                 its, ite, jts, jte, kts, kte                     )
1035
1036        ENDIF
1037
1038    CASE (NSASSCHEME)
1039         CALL add_a2a(rt_tendf,RTHCUTEN,config_flags,             &
1040                 ids,ide, jds, jde, kds, kde,                     &
1041                 ims, ime, jms, jme, kms, kme,                    &
1042                 its, ite, jts, jte, kts, kte                     )
1043
1044         CALL add_a2c_u(ru_tendf,RUCUTEN,config_flags,            &
1045                 ids,ide, jds, jde, kds, kde,                     &
1046                 ims, ime, jms, jme, kms, kme,                    &
1047                 its, ite, jts, jte, kts, kte                     )
1048
1049         CALL add_a2c_v(rv_tendf,RVCUTEN,config_flags,            &
1050                 ids,ide, jds, jde, kds, kde,                     &
1051                 ims, ime, jms, jme, kms, kme,                    &
1052                 its, ite, jts, jte, kts, kte                     )
1053
1054         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
1055         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
1056                 config_flags,                                    &
1057                 ids,ide, jds, jde, kds, kde,                     &
1058                 ims, ime, jms, jme, kms, kme,                    &
1059                 its, ite, jts, jte, kts, kte                     )
1060
1061         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
1062         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCCUTEN,     &
1063                 config_flags,                                    &
1064                 ids,ide, jds, jde, kds, kde,                     &
1065                 ims, ime, jms, jme, kms, kme,                    &
1066                 its, ite, jts, jte, kts, kte                     )
1067
1068         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
1069         CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQICUTEN,     &
1070                 config_flags,                                    &
1071                 ids,ide, jds, jde, kds, kde,                     &
1072                 ims, ime, jms, jme, kms, kme,                    &
1073                 its, ite, jts, jte, kts, kte                     )
1074
1075        IF(.not. adv_moist_cond)THEN
1076
1077         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1078            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCCUTEN, &
1079                 config_flags,                                    &
1080                 ids,ide, jds, jde, kds, kde,                     &
1081                 ims, ime, jms, jme, kms, kme,                    &
1082                 its, ite, jts, jte, kts, kte                     )
1083
1084         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1085            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQICUTEN, &
1086                 config_flags,                                    &
1087                 ids,ide, jds, jde, kds, kde,                     &
1088                 ims, ime, jms, jme, kms, kme,                    &
1089                 its, ite, jts, jte, kts, kte                     )
1090
1091        ENDIF
1092
1093    CASE (SASSCHEME,OSASSCHEME)
1094         CALL add_a2a(rt_tendf,RTHCUTEN,config_flags,             &
1095                 ids,ide, jds, jde, kds, kde,                     &
1096                 ims, ime, jms, jme, kms, kme,                    &
1097                 its, ite, jts, jte, kts, kte                     )
1098
1099         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
1100         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
1101                 config_flags,                                    &
1102                 ids,ide, jds, jde, kds, kde,                     &
1103                 ims, ime, jms, jme, kms, kme,                    &
1104                 its, ite, jts, jte, kts, kte                     )
1105
1106         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
1107         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCCUTEN,     &
1108                 config_flags,                                    &
1109                 ids,ide, jds, jde, kds, kde,                     &
1110                 ims, ime, jms, jme, kms, kme,                    &
1111                 its, ite, jts, jte, kts, kte                     )
1112
1113         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
1114         CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQICUTEN,     &
1115                 config_flags,                                    &
1116                 ids,ide, jds, jde, kds, kde,                     &
1117                 ims, ime, jms, jme, kms, kme,                    &
1118                 its, ite, jts, jte, kts, kte                     )
1119
1120        IF(.not. adv_moist_cond)THEN
1121
1122         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1123            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCCUTEN, &
1124                 config_flags,                                    &
1125                 ids,ide, jds, jde, kds, kde,                     &
1126                 ims, ime, jms, jme, kms, kme,                    &
1127                 its, ite, jts, jte, kts, kte                     )
1128
1129         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1130            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQICUTEN, &
1131                 config_flags,                                    &
1132                 ids,ide, jds, jde, kds, kde,                     &
1133                 ims, ime, jms, jme, kms, kme,                    &
1134                 its, ite, jts, jte, kts, kte                     )
1135
1136        ENDIF
1137
1138    CASE (CAMZMSCHEME)
1139         CALL add_a2c_u(ru_tendf,RUCUTEN,config_flags,            &
1140                 ids,ide, jds, jde, kds, kde,                     &
1141                 ims, ime, jms, jme, kms, kme,                    &
1142                 its, ite, jts, jte, kts, kte                     )
1143
1144         CALL add_a2c_v(rv_tendf,RVCUTEN,config_flags,            &
1145                 ids,ide, jds, jde, kds, kde,                     &
1146                 ims, ime, jms, jme, kms, kme,                    &
1147                 its, ite, jts, jte, kts, kte                     )
1148
1149         CALL add_a2a(rt_tendf,RTHCUTEN,config_flags,             &
1150                 ids,ide, jds, jde, kds, kde,                     &
1151                 ims, ime, jms, jme, kms, kme,                    &
1152                 its, ite, jts, jte, kts, kte                     )
1153
1154         if (P_QV .ge. PARAM_FIRST_SCALAR)                        &
1155         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
1156                 config_flags,                                    &
1157                 ids,ide, jds, jde, kds, kde,                     &
1158                 ims, ime, jms, jme, kms, kme,                    &
1159                 its, ite, jts, jte, kts, kte                     )
1160
1161         if (P_QC .ge. PARAM_FIRST_SCALAR)                        &
1162         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCCUTEN,     &
1163                 config_flags,                                    &
1164                 ids,ide, jds, jde, kds, kde,                     &
1165                 ims, ime, jms, jme, kms, kme,                    &
1166                 its, ite, jts, jte, kts, kte                     )
1167
1168         if (P_QI .ge. PARAM_FIRST_SCALAR)                        &
1169         CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQICUTEN,     &
1170                 config_flags,                                    &
1171                 ids,ide, jds, jde, kds, kde,                     &
1172                 ims, ime, jms, jme, kms, kme,                    &
1173                 its, ite, jts, jte, kts, kte                     )
1174        IF(.not. adv_moist_cond)THEN
1175
1176         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1177            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCCUTEN, &
1178                 config_flags,                                    &
1179                 ids,ide, jds, jde, kds, kde,                     &
1180                 ims, ime, jms, jme, kms, kme,                    &
1181                 its, ite, jts, jte, kts, kte                     )
1182
1183         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1184            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQICUTEN, &
1185                 config_flags,                                    &
1186                 ids,ide, jds, jde, kds, kde,                     &
1187                 ims, ime, jms, jme, kms, kme,                    &
1188                 its, ite, jts, jte, kts, kte                     )
1189
1190        ENDIF
1191
1192    CASE (TIEDTKESCHEME)
1193         CALL add_a2a(rt_tendf,RTHCUTEN,config_flags,             &
1194                 ids,ide, jds, jde, kds, kde,                     &
1195                 ims, ime, jms, jme, kms, kme,                    &
1196                 its, ite, jts, jte, kts, kte                     )
1197
1198         CALL add_a2c_u(ru_tendf,RUCUTEN,config_flags,            &
1199                 ids,ide, jds, jde, kds, kde,                     &
1200                 ims, ime, jms, jme, kms, kme,                    &
1201                 its, ite, jts, jte, kts, kte                     )
1202
1203         CALL add_a2c_v(rv_tendf,RVCUTEN,config_flags,            &
1204                 ids,ide, jds, jde, kds, kde,                     &
1205                 ims, ime, jms, jme, kms, kme,                    &
1206                 its, ite, jts, jte, kts, kte                     )
1207
1208         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
1209         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVCUTEN,     &
1210                 config_flags,                                    &
1211                 ids,ide, jds, jde, kds, kde,                     &
1212                 ims, ime, jms, jme, kms, kme,                    &
1213                 its, ite, jts, jte, kts, kte                     )
1214
1215         if (P_QC .ge. PARAM_FIRST_SCALAR)                                         &
1216         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCCUTEN,     &
1217                 config_flags,                                    &
1218                 ids,ide, jds, jde, kds, kde,                     &
1219                 ims, ime, jms, jme, kms, kme,                    &
1220                 its, ite, jts, jte, kts, kte                     )
1221
1222         if (P_QI .ge. PARAM_FIRST_SCALAR)                                         &
1223         CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQICUTEN,     &
1224                 config_flags,                                    &
1225                 ids,ide, jds, jde, kds, kde,                     &
1226                 ims, ime, jms, jme, kms, kme,                    &
1227                 its, ite, jts, jte, kts, kte                     )
1228
1229        IF(.not. adv_moist_cond)THEN
1230
1231         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1232            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQCCUTEN, &
1233                 config_flags,                                    &
1234                 ids,ide, jds, jde, kds, kde,                     &
1235                 ims, ime, jms, jme, kms, kme,                    &
1236                 its, ite, jts, jte, kts, kte                     )
1237
1238         if (P_QT .ge. PARAM_FIRST_SCALAR)                        &
1239            CALL add_a2a(scalar_tendf(ims,kms,jms,P_QT),RQICUTEN, &
1240                 config_flags,                                    &
1241                 ids,ide, jds, jde, kds, kde,                     &
1242                 ims, ime, jms, jme, kms, kme,                    &
1243                 its, ite, jts, jte, kts, kte                     )
1244
1245        ENDIF
1246
1247    CASE DEFAULT
1248
1249    END SELECT
1250
1251 END SUBROUTINE phy_cu_ten
1252
1253 !=================================================================
1254 SUBROUTINE phy_shcu_ten(config_flags,rk_step,n_moist,            &
1255                       rt_tendf,ru_tendf,rv_tendf,                &
1256                       RUSHTEN,RVSHTEN,RTHSHTEN,                  &
1257                       RQVSHTEN,RQCSHTEN,RQRSHTEN,                &
1258                       RQISHTEN,RQSSHTEN,RQGSHTEN,moist_tendf,    &
1259                       ids, ide, jds, jde, kds, kde,              &
1260                       ims, ime, jms, jme, kms, kme,              &
1261                       its, ite, jts, jte, kts, kte               )
1262 !-----------------------------------------------------------------
1263    IMPLICIT NONE
1264 !-----------------------------------------------------------------
1265    TYPE(grid_config_rec_type  ) , INTENT(IN   ) :: config_flags
1266
1267    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde, &
1268                                    ims, ime, jms, jme, kms, kme, &
1269                                    its, ite, jts, jte, kts, kte, &
1270                                    n_moist, rk_step
1271
1272    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_moist),         &
1273           INTENT(INOUT)     ::                      moist_tendf
1274
1275    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
1276                                                         RUSHTEN, &
1277                                                         RVSHTEN, &
1278                                                        RTHSHTEN, &
1279                                                        RQVSHTEN, &
1280                                                        RQCSHTEN, &
1281                                                        RQRSHTEN, &
1282                                                        RQISHTEN, &
1283                                                        RQSSHTEN, &
1284                                                        RQGSHTEN
1285
1286    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT) ::  &
1287                                                        rt_tendf, &
1288                                                        ru_tendf, &
1289                                                        rv_tendf
1290
1291 ! LOCAL VARS
1292
1293    INTEGER :: i,j,k
1294
1295    SELECT CASE (config_flags%shcu_physics)
1296
1297    CASE (CAMUWSHCUSCHEME)
1298         CALL add_a2c_u(ru_tendf,RUSHTEN,config_flags,            &
1299                 ids,ide, jds, jde, kds, kde,                     &
1300                 ims, ime, jms, jme, kms, kme,                    &
1301                 its, ite, jts, jte, kts, kte                     )
1302
1303         CALL add_a2c_v(rv_tendf,RVSHTEN,config_flags,            &
1304                 ids,ide, jds, jde, kds, kde,                     &
1305                 ims, ime, jms, jme, kms, kme,                    &
1306                 its, ite, jts, jte, kts, kte                     )
1307
1308         CALL add_a2a(rt_tendf,RTHSHTEN,config_flags,             &
1309                 ids,ide, jds, jde, kds, kde,                     &
1310                 ims, ime, jms, jme, kms, kme,                    &
1311                 its, ite, jts, jte, kts, kte                     )
1312
1313         if (P_QV .ge. PARAM_FIRST_SCALAR)                        &
1314         CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVSHTEN,     &
1315                 config_flags,                                    &
1316                 ids,ide, jds, jde, kds, kde,                     &
1317                 ims, ime, jms, jme, kms, kme,                    &
1318                 its, ite, jts, jte, kts, kte                     )
1319
1320         if (P_QC .ge. PARAM_FIRST_SCALAR)                        &
1321         CALL add_a2a(moist_tendf(ims,kms,jms,P_QC),RQCSHTEN,     &
1322                 config_flags,                                    &
1323                 ids,ide, jds, jde, kds, kde,                     &
1324                 ims, ime, jms, jme, kms, kme,                    &
1325                 its, ite, jts, jte, kts, kte                     )
1326
1327         if (P_QR .ge. PARAM_FIRST_SCALAR)                        &
1328         CALL add_a2a(moist_tendf(ims,kms,jms,P_QR),RQRSHTEN,     &
1329                 config_flags,                                    &
1330                 ids,ide, jds, jde, kds, kde,                     &
1331                 ims, ime, jms, jme, kms, kme,                    &
1332                 its, ite, jts, jte, kts, kte                     )
1333
1334         if (P_QI .ge. PARAM_FIRST_SCALAR)                        &
1335         CALL add_a2a(moist_tendf(ims,kms,jms,P_QI),RQISHTEN,     &
1336                 config_flags,                                    &
1337                 ids,ide, jds, jde, kds, kde,                     &
1338                 ims, ime, jms, jme, kms, kme,                    &
1339                 its, ite, jts, jte, kts, kte                     )
1340
1341         if (P_QS .ge. PARAM_FIRST_SCALAR)                        &
1342         CALL add_a2a(moist_tendf(ims,kms,jms,P_QS),RQSSHTEN,     &
1343                 config_flags,                                    &
1344                 ids,ide, jds, jde, kds, kde,                     &
1345                 ims, ime, jms, jme, kms, kme,                    &
1346                 its, ite, jts, jte, kts, kte                     )
1347
1348         if (P_QG .ge. PARAM_FIRST_SCALAR)                        &
1349         CALL add_a2a(moist_tendf(ims,kms,jms,P_QG),RQGSHTEN,     &
1350                 config_flags,                                    &
1351                 ids,ide, jds, jde, kds, kde,                     &
1352                 ims, ime, jms, jme, kms, kme,                    &
1353                 its, ite, jts, jte, kts, kte                     )
1354
1355    CASE DEFAULT
1356
1357    END SELECT
1358
1359
1360 END SUBROUTINE phy_shcu_ten
1361
1362 !=================================================================
1363 SUBROUTINE phy_fg_ten(config_flags,rk_step,n_moist,            &
1364                       rph_tendf,rt_tendf,ru_tendf,rv_tendf,    &
1365                       mu_tendf, moist_tendf,                   &
1366                       RUNDGDTEN,RVNDGDTEN,RTHNDGDTEN,          &
1367                       RPHNDGDTEN,RQVNDGDTEN,RMUNDGDTEN,        &
1368                       ids, ide, jds, jde, kds, kde,              &
1369                       ims, ime, jms, jme, kms, kme,              &
1370                       its, ite, jts, jte, kts, kte               )
1371 !-----------------------------------------------------------------
1372    IMPLICIT NONE
1373 !-----------------------------------------------------------------
1374    TYPE(grid_config_rec_type) ,     INTENT(IN   ) :: config_flags
1375
1376    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde, &
1377                                    ims, ime, jms, jme, kms, kme, &
1378                                    its, ite, jts, jte, kts, kte, &
1379                                    n_moist, rk_step
1380
1381    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_moist),         &
1382           INTENT(INOUT)     ::                      moist_tendf
1383
1384    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
1385                                                        RTHNDGDTEN, &
1386                                                        RPHNDGDTEN, &
1387                                                         RUNDGDTEN, &
1388                                                         RVNDGDTEN, &
1389                                                        RQVNDGDTEN
1390
1391    REAL, DIMENSION(ims:ime, jms:jme), INTENT(IN  ) ::  RMUNDGDTEN
1392
1393    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT)::   &
1394                                                        rph_tendf,&
1395                                                        rt_tendf, &
1396                                                        ru_tendf, &
1397                                                        rv_tendf
1398
1399    REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT)::  mu_tendf
1400
1401 ! LOCAL VARS
1402
1403    INTEGER :: i,j,k,IBGN,IEND,JBGN,JEND
1404
1405 !-----------------------------------------------------------------
1406
1407    SELECT CASE(config_flags%grid_fdda)
1408
1409       CASE (PSUFDDAGD)
1410
1411            CALL add_a2a(rt_tendf,RTHNDGDTEN,config_flags,        &
1412                 ids,ide, jds, jde, kds, kde,                     &
1413                 ims, ime, jms, jme, kms, kme,                    &
1414                 its, ite, jts, jte, kts, kte                     )
1415
1416 ! note fdda u and v tendencies are staggered
1417            CALL add_c2c_u(ru_tendf,RUNDGDTEN,config_flags,       &
1418                 ids,ide, jds, jde, kds, kde,                     &
1419                 ims, ime, jms, jme, kms, kme,                    &
1420                 its, ite, jts, jte, kts, kte                     )
1421
1422            CALL add_c2c_v(rv_tendf,RVNDGDTEN,config_flags,       &
1423                 ids,ide, jds, jde, kds, kde,                     &
1424                 ims, ime, jms, jme, kms, kme,                    &
1425                 its, ite, jts, jte, kts, kte                     )
1426
1427            CALL add_a2a(mu_tendf,RMUNDGDTEN,config_flags,      &
1428                 ids,ide, jds, jde, kds, kds,                     &
1429                 ims, ime, jms, jme, kms, kms,                    &
1430                 its, ite, jts, jte, kts, kts                     )
1431
1432         if (P_QV .ge. PARAM_FIRST_SCALAR)                                         &
1433            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),RQVNDGDTEN,  &
1434                 config_flags,                                    &
1435                 ids,ide, jds, jde, kds, kde,                     &
1436                 ims, ime, jms, jme, kms, kme,                    &
1437                 its, ite, jts, jte, kts, kte                     )
1438
1439       CASE (SPNUDGING)
1440
1441 ! note fdda u and v tendencies are staggered
1442            CALL add_c2c_u(ru_tendf,RUNDGDTEN,config_flags,       &
1443                 ids,ide, jds, jde, kds, kde,                     &
1444                 ims, ime, jms, jme, kms, kme,                    &
1445                 its, ite, jts, jte, kts, kte                     )
1446
1447            CALL add_c2c_v(rv_tendf,RVNDGDTEN,config_flags,       &
1448                 ids,ide, jds, jde, kds, kde,                     &
1449                 ims, ime, jms, jme, kms, kme,                    &
1450                 its, ite, jts, jte, kts, kte                     )
1451
1452            CALL add_a2a(rt_tendf,RTHNDGDTEN,config_flags,        &
1453                 ids,ide, jds, jde, kds, kde,                     &
1454                 ims, ime, jms, jme, kms, kme,                    &
1455                 its, ite, jts, jte, kts, kte                     )
1456
1457            CALL add_a2a_ph(rph_tendf,RPHNDGDTEN,config_flags,        &
1458                 ids,ide, jds, jde, kds, kde,                     &
1459                 ims, ime, jms, jme, kms, kme,                    &
1460                 its, ite, jts, jte, kts, kte                     )
1461
1462       CASE DEFAULT
1463
1464    END SELECT
1465
1466 END SUBROUTINE phy_fg_ten
1467
1468 !=================================================================
1469 SUBROUTINE phy_fr_ten(config_flags,rk_step,n_moist,            &
1470                       rt_tendf,ru_tendf,rv_tendf,              &
1471                       mu_tendf, moist_tendf,                   &
1472                       rthfrten,rqvfrten,                       &
1473                       ids, ide, jds, jde, kds, kde,              &
1474                       ims, ime, jms, jme, kms, kme,              &
1475                       its, ite, jts, jte, kts, kte               )
1476 !-----------------------------------------------------------------
1477    USE module_state_description, ONLY :                         &
1478                    FIRE_SFIRE
1479 !-----------------------------------------------------------------
1480    IMPLICIT NONE
1481 !-----------------------------------------------------------------
1482    TYPE(grid_config_rec_type) ,     INTENT(IN   ) :: config_flags
1483
1484    INTEGER , INTENT(IN)        ::  ids, ide, jds, jde, kds, kde, &
1485                                    ims, ime, jms, jme, kms, kme, &
1486                                    its, ite, jts, jte, kts, kte, &
1487                                    n_moist, rk_step
1488
1489    REAL , DIMENSION(ims:ime, kms:kme, jms:jme, n_moist),         &
1490           INTENT(INOUT)     ::                      moist_tendf
1491
1492    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(IN  ) ::   &
1493                                                        rthfrten, &
1494                                                        rqvfrten
1495
1496    REAL, DIMENSION(ims:ime, kms:kme, jms:jme), INTENT(INOUT)::   &
1497                                                        rt_tendf, &
1498                                                        ru_tendf, &
1499                                                        rv_tendf
1500
1501    REAL, DIMENSION(ims:ime, jms:jme), INTENT(INOUT)::  mu_tendf
1502
1503 ! LOCAL VARS
1504
1505    INTEGER :: i,j,k,IBGN,IEND,JBGN,JEND
1506
1507 !-----------------------------------------------------------------
1508
1509    SELECT CASE(config_flags%ifire)
1510
1511       CASE (FIRE_SFIRE)
1512
1513            CALL add_a2a(rt_tendf,rthfrten,                       &
1514                 config_flags,                                    &
1515                 ids,ide, jds, jde, kds, kde,                     &
1516                 ims, ime, jms, jme, kms, kme,                    &
1517                 its, ite, jts, jte, kts, kte                     )
1518
1519            CALL add_a2a(moist_tendf(ims,kms,jms,P_QV),rqvfrten,  &
1520                 config_flags,                                    &
1521                 ids,ide, jds, jde, kds, kde,                     &
1522                 ims, ime, jms, jme, kms, kme,                    &
1523                 its, ite, jts, jte, kts, kte                     )
1524
1525       CASE DEFAULT
1526
1527    END SELECT
1528
1529 END SUBROUTINE phy_fr_ten
1530
1531 !----------------------------------------------------------------------
1532 SUBROUTINE advance_ppt(RTHCUTEN,RQVCUTEN,RQCCUTEN,RQRCUTEN,           &
1533                      RQICUTEN,RQSCUTEN,                               &
1534                      RAINC,RAINCV,RAINSH,PRATEC,PRATESH,  &
1535                      NCA, HTOP,HBOT,CUTOP,CUBOT,                      &
1536                      CUPPT, DT, config_flags,                         &
1537                      ids,ide, jds,jde, kds,kde,                       &
1538                      ims,ime, jms,jme, kms,kme,                       &
1539                      its,ite, jts,jte, kts,kte                        )
1540 !----------------------------------------------------------------------
1541    USE module_state_description
1542    USE module_cu_kf
1543    USE module_cu_kfeta
1544 !----------------------------------------------------------------------
1545    IMPLICIT NONE
1546 !----------------------------------------------------------------------
1547    TYPE(grid_config_rec_type),  INTENT(IN   )    :: config_flags
1548
1549    INTEGER,      INTENT(IN   )    ::                             &
1550                                       ids,ide, jds,jde, kds,kde, &
1551                                       ims,ime, jms,jme, kms,kme, &
1552                                       its,ite, jts,jte, kts,kte
1553
1554
1555    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ),                 &
1556          INTENT(INOUT)  ::                             RTHCUTEN, &
1557                                                        RQVCUTEN, &
1558                                                        RQCCUTEN, &
1559                                                        RQRCUTEN, &
1560                                                        RQICUTEN, &
1561                                                        RQSCUTEN
1562
1563    REAL, DIMENSION( ims:ime , jms:jme ),                         &
1564           INTENT(INOUT) ::                                RAINC, &
1565                                                          RAINSH, &
1566                                                          RAINCV, &
1567                                                          PRATEC, &
1568                                                         PRATESH, &
1569                                                             NCA, &
1570                                                            HTOP, &
1571                                                            HBOT, &
1572                                                           CUTOP, &
1573                                                           CUBOT, &
1574                                                           CUPPT
1575    REAL, INTENT(IN) ::                                       DT
1576
1577 ! LOCAL  VAR
1578
1579    INTEGER :: i,j,k,i_start,i_end,j_start,j_end,k_start,k_end
1580    INTEGER :: NCUTOP, NCUBOT
1581
1582 !-----------------------------------------------------------------
1583
1584    IF (config_flags%cu_physics .eq. 0) return
1585
1586 ! SET START AND END POINTS FOR TILES
1587
1588    i_start = its
1589    i_end   = min( ite,ide-1 )
1590    j_start = jts
1591    j_end   = min( jte,jde-1 )
1592 !
1593 !  IF( config_flags%nested .or. config_flags%specified ) THEN
1594 !    i_start = max( its,ids+1 )
1595 !    i_end   = min( ite,ide-2 )
1596 !    j_start = max( jts,jds+1 )
1597 !    j_end   = min( jte,jde-2 )
1598 !  ENDIF
1599 !
1600    k_start = kts
1601    k_end = min( kte, kde-1 )
1602
1603 ! Update total cumulus scheme precipitation
1604
1605 ! in mm
1606
1607    DO J = j_start,j_end
1608    DO i = i_start,i_end
1609       RAINC(I,J)  = RAINC(I,J) + PRATEC(I,J)*DT
1610       RAINSH(I,J) = RAINSH(I,J) + PRATESH(I,J)*DT
1611       CUPPT(I,J)  = CUPPT(I,J) + (PRATEC(I,J)+PRATESH(I,J))*DT/1000.
1612    ENDDO
1613    ENDDO
1614
1615    SELECT CASE (config_flags%cu_physics)
1616
1617    CASE (KFSCHEME)
1618
1619         DO J = j_start,j_end
1620         DO i = i_start,i_end
1621
1622            IF ( NCA(I,J) .GT. 0 ) THEN
1623
1624               IF ( NINT(NCA(I,J) / DT) .le. 0 ) THEN
1625
1626               ! set tendency to zero
1627 !                PRATEC(I,J)=0.
1628 !                RAINCV(I,J)=0.
1629                  DO k = k_start,k_end
1630                     RTHCUTEN(i,k,j)=0.
1631                     RQVCUTEN(i,k,j)=0.
1632                     RQCCUTEN(i,k,j)=0.
1633                     RQRCUTEN(i,k,j)=0.
1634                     if (P_QI .ge. PARAM_FIRST_SCALAR) RQICUTEN(i,k,j)=0.
1635                     if (P_QS .ge. PARAM_FIRST_SCALAR) RQSCUTEN(i,k,j)=0.
1636                  ENDDO
1637               ENDIF
1638
1639               NCA(I,J)=NCA(I,J)-DT ! Decrease NCA
1640
1641            ENDIF
1642 !
1643         ENDDO
1644         ENDDO
1645
1646    CASE (BMJSCHEME, CAMZMSCHEME)
1647
1648         DO J = j_start,j_end
1649         DO i = i_start,i_end
1650
1651 ! HTOP, HBOT FOR GFDL RADIATION
1652            NCUTOP=NINT(CUTOP(I,J))
1653            NCUBOT=NINT(CUBOT(I,J))
1654            IF(NCUTOP>1.AND.NCUTOP<KDE)THEN
1655              HTOP(I,J)=MAX(CUTOP(I,J),HTOP(I,J))
1656            ENDIF
1657            IF(NCUBOT>0.AND.NCUBOT<KDE)THEN
1658              HBOT(I,J)=MIN(CUBOT(I,J),HBOT(I,J))
1659            ENDIF
1660
1661         ENDDO
1662         ENDDO
1663
1664    CASE (KFETASCHEME)
1665
1666         DO J = j_start,j_end
1667         DO i = i_start,i_end
1668
1669 ! HTOP, HBOT FOR GFDL RADIATION
1670            NCUTOP=NINT(CUTOP(I,J))
1671            NCUBOT=NINT(CUBOT(I,J))
1672            IF(NCUTOP>1.AND.NCUTOP<KDE)THEN
1673              HTOP(I,J)=MAX(CUTOP(I,J),HTOP(I,J))
1674            ENDIF
1675            IF(NCUBOT>0.AND.NCUBOT<KDE)THEN
1676              HBOT(I,J)=MIN(CUBOT(I,J),HBOT(I,J))
1677            ENDIF
1678
1679            IF ( NCA(I,J) .GT. 0 ) THEN
1680
1681
1682               IF ( NINT(NCA(I,J) / DT) .LE. 1 ) THEN
1683
1684               ! set tendency to zero
1685 !                PRATEC(I,J)=0.
1686 !                RAINCV(I,J)=0.
1687                  DO k = k_start,k_end
1688                     RTHCUTEN(i,k,j)=0.
1689                     RQVCUTEN(i,k,j)=0.
1690                     RQCCUTEN(i,k,j)=0.
1691                     RQRCUTEN(i,k,j)=0.
1692                     if (P_QI .ge. PARAM_FIRST_SCALAR) RQICUTEN(i,k,j)=0.
1693                     if (P_QS .ge. PARAM_FIRST_SCALAR) RQSCUTEN(i,k,j)=0.
1694                  ENDDO
1695               ENDIF
1696
1697               NCA(I,J)=NCA(I,J)-DT ! Decrease NCA
1698 !              NCA(I,J)=NCA(I,J)-1. ! Decrease NCA
1699
1700            ENDIF
1701 !
1702         ENDDO
1703         ENDDO
1704
1705    CASE DEFAULT
1706
1707    END SELECT
1708
1709 END SUBROUTINE advance_ppt
1710
1711 SUBROUTINE add_a2a(lvar,rvar,config_flags,                  &
1712                    ids,ide, jds, jde, kds, kde,             &
1713                    ims, ime, jms, jme, kms, kme,            &
1714                    its, ite, jts, jte, kts, kte             )
1715 !------------------------------------------------------------
1716    IMPLICIT NONE
1717 !------------------------------------------------------------
1718    TYPE(grid_config_rec_type),  INTENT(IN) :: config_flags
1719
1720    INTEGER , INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
1721                               ims, ime, jms, jme, kms, kme, &
1722                               its, ite, jts, jte, kts, kte
1723
1724    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(IN   ) ::&
1725                                                       rvar
1726    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(INOUT) ::&
1727                                                       lvar
1728
1729 ! LOCAL VARS
1730    INTEGER :: i,j,k,i_start,i_end,j_start,j_end,ktf
1731
1732    i_start = its
1733    i_end   = MIN(ite,ide-1)
1734    j_start = jts
1735    j_end   = MIN(jte,jde-1)
1736    ktf = min(kte,kde-1)
1737
1738    IF ( config_flags%specified .or. &
1739         config_flags%nested) i_start = MAX(ids+1,its)
1740    IF ( config_flags%specified .or. &
1741         config_flags%nested) i_end   = MIN(ide-2,ite)
1742    IF ( config_flags%specified .or. &
1743         config_flags%nested) j_start = MAX(jds+1,jts)
1744    IF ( config_flags%specified .or. &
1745         config_flags%nested) j_end   = MIN(jde-2,jte)
1746       IF ( config_flags%periodic_x ) i_start = its
1747       IF ( config_flags%periodic_x ) i_end = MIN( ite, ide-1 )
1748
1749    DO j = j_start,j_end
1750    DO k = kts,ktf
1751    DO i = i_start,i_end
1752       lvar(i,k,j) = lvar(i,k,j) + rvar(i,k,j)
1753    ENDDO
1754    ENDDO
1755    ENDDO
1756
1757 END SUBROUTINE add_a2a
1758
1759 SUBROUTINE add_a2a_ph(lvar,rvar,config_flags,                  &
1760                    ids,ide, jds, jde, kds, kde,             &
1761                    ims, ime, jms, jme, kms, kme,            &
1762                    its, ite, jts, jte, kts, kte             )
1763 !------------------------------------------------------------
1764    IMPLICIT NONE
1765 !------------------------------------------------------------
1766    TYPE(grid_config_rec_type),  INTENT(IN) :: config_flags
1767
1768    INTEGER , INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
1769                               ims, ime, jms, jme, kms, kme, &
1770                               its, ite, jts, jte, kts, kte
1771
1772    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(IN   ) ::&
1773                                                       rvar
1774    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(INOUT) ::&
1775                                                       lvar
1776
1777 ! LOCAL VARS
1778    INTEGER :: i,j,k,i_start,i_end,j_start,j_end
1779
1780    i_start = its
1781    i_end   = MIN(ite,ide-1)
1782    j_start = jts
1783    j_end   = MIN(jte,jde-1)
1784
1785    IF ( config_flags%specified .or. &
1786         config_flags%nested) i_start = MAX(ids+1,its)
1787    IF ( config_flags%specified .or. &
1788         config_flags%nested) i_end   = MIN(ide-2,ite)
1789    IF ( config_flags%specified .or. &
1790         config_flags%nested) j_start = MAX(jds+1,jts)
1791    IF ( config_flags%specified .or. &
1792         config_flags%nested) j_end   = MIN(jde-2,jte)
1793       IF ( config_flags%periodic_x ) i_start = its
1794       IF ( config_flags%periodic_x ) i_end = MIN( ite, ide-1 )
1795
1796    DO j = j_start,j_end
1797    DO k = kts,kte
1798    DO i = i_start,i_end
1799       lvar(i,k,j) = lvar(i,k,j) + rvar(i,k,j)
1800    ENDDO
1801    ENDDO
1802    ENDDO
1803
1804 END SUBROUTINE add_a2a_ph
1805
1806 !------------------------------------------------------------
1807 SUBROUTINE add_a2c_u(lvar,rvar,config_flags,                &
1808                    ids,ide, jds, jde, kds, kde,             &
1809                    ims, ime, jms, jme, kms, kme,            &
1810                    its, ite, jts, jte, kts, kte             )
1811 !------------------------------------------------------------
1812 !------------------------------------------------------------
1813    IMPLICIT NONE
1814 !------------------------------------------------------------
1815
1816    TYPE(grid_config_rec_type),  INTENT(IN   )    :: config_flags
1817
1818    INTEGER , INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
1819                               ims, ime, jms, jme, kms, kme, &
1820                               its, ite, jts, jte, kts, kte
1821
1822    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(IN   ) ::&
1823                                                       rvar
1824    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(INOUT) ::&
1825                                                       lvar
1826
1827 ! LOCAL VARS
1828
1829    INTEGER :: i,j,k,i_start,i_end,j_start,j_end,ktf
1830
1831    ktf=min(kte,kde-1)
1832
1833    i_start = its
1834    i_end   = ite
1835    j_start = jts
1836    j_end   = MIN(jte,jde-1)
1837
1838    IF ( config_flags%specified .or. &
1839         config_flags%nested) i_start = MAX(ids+1,its)
1840    IF ( config_flags%specified .or. &
1841         config_flags%nested) i_end   = MIN(ide-1,ite)
1842    IF ( config_flags%specified .or. &
1843         config_flags%nested) j_start = MAX(jds+1,jts)
1844    IF ( config_flags%specified .or. &
1845         config_flags%nested) j_end   = MIN(jde-2,jte)
1846       IF ( config_flags%periodic_x ) i_start = its
1847       IF ( config_flags%periodic_x ) i_end = ite
1848
1849    DO j = j_start,j_end
1850    DO k = kts,ktf
1851    DO i = i_start,i_end
1852       lvar(i,k,j) = lvar(i,k,j) + &
1853                        0.5*(rvar(i,k,j)+rvar(i-1,k,j))
1854    ENDDO
1855    ENDDO
1856    ENDDO
1857
1858 END SUBROUTINE add_a2c_u
1859
1860 !------------------------------------------------------------
1861 SUBROUTINE add_a2c_v(lvar,rvar,config_flags,                &
1862                    ids,ide, jds, jde, kds, kde,             &
1863                    ims, ime, jms, jme, kms, kme,            &
1864                    its, ite, jts, jte, kts, kte             )
1865 !------------------------------------------------------------
1866 !------------------------------------------------------------
1867    IMPLICIT NONE
1868 !------------------------------------------------------------
1869
1870    TYPE(grid_config_rec_type),  INTENT(IN   )    :: config_flags
1871
1872    INTEGER , INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
1873                               ims, ime, jms, jme, kms, kme, &
1874                               its, ite, jts, jte, kts, kte
1875
1876    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(IN   ) ::&
1877                                                       rvar
1878    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(INOUT) ::&
1879                                                       lvar
1880
1881 ! LOCAL VARS
1882
1883    INTEGER :: i,j,k,i_start,i_end,j_start,j_end,ktf
1884
1885    ktf=min(kte,kde-1)
1886
1887    i_start = its
1888    i_end   = MIN(ite,ide-1)
1889    j_start = jts
1890    j_end   = jte
1891
1892    IF ( config_flags%specified .or. &
1893         config_flags%nested) i_start = MAX(ids+1,its)
1894    IF ( config_flags%specified .or. &
1895         config_flags%nested) i_end   = MIN(ide-2,ite)
1896    IF ( config_flags%specified .or. &
1897         config_flags%nested) j_start = MAX(jds+1,jts)
1898    IF ( config_flags%specified .or. &
1899         config_flags%nested) j_end   = MIN(jde-1,jte)
1900       IF ( config_flags%periodic_x ) i_start = its
1901       IF ( config_flags%periodic_x ) i_end = MIN( ite, ide-1 )
1902
1903    DO j = j_start,j_end
1904    DO k = kts,kte
1905    DO i = i_start,i_end
1906       lvar(i,k,j) = lvar(i,k,j) + &
1907                      0.5*(rvar(i,k,j)+rvar(i,k,j-1))
1908    ENDDO
1909    ENDDO
1910    ENDDO
1911
1912 END SUBROUTINE add_a2c_v
1913
1914 !------------------------------------------------------------
1915 SUBROUTINE add_c2c_u(lvar,rvar,config_flags,                &
1916                    ids,ide, jds, jde, kds, kde,             &
1917                    ims, ime, jms, jme, kms, kme,            &
1918                    its, ite, jts, jte, kts, kte             )
1919 !------------------------------------------------------------
1920 !------------------------------------------------------------
1921    IMPLICIT NONE
1922 !------------------------------------------------------------
1923
1924    TYPE(grid_config_rec_type),  INTENT(IN   )    :: config_flags
1925
1926    INTEGER , INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
1927                               ims, ime, jms, jme, kms, kme, &
1928                               its, ite, jts, jte, kts, kte
1929
1930    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(IN   ) ::&
1931                                                       rvar
1932    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(INOUT) ::&
1933                                                       lvar
1934
1935 ! LOCAL VARS
1936
1937    INTEGER :: i,j,k,i_start,i_end,j_start,j_end,ktf
1938
1939    ktf=min(kte,kde-1)
1940
1941    i_start = its
1942    i_end   = ite
1943    j_start = jts
1944    j_end   = MIN(jte,jde-1)
1945
1946
1947    IF ( config_flags%specified .or. &
1948         config_flags%nested) i_start = MAX(ids+1,its)
1949    IF ( config_flags%specified .or. &
1950         config_flags%nested) i_end   = MIN(ide-1,ite)
1951    IF ( config_flags%specified .or. &
1952         config_flags%nested) j_start = MAX(jds+1,jts)
1953    IF ( config_flags%specified .or. &
1954         config_flags%nested) j_end   = MIN(jde-2,jte)
1955
1956 !  write(*,'(a,6i4)') 'call c2cu, i_start, i_end, j_start, j_end=', i_start, i_end, j_start, j_end
1957
1958    DO j = j_start,j_end
1959    DO k = kts,ktf
1960    DO i = i_start,i_end
1961       lvar(i,k,j) = lvar(i,k,j) + rvar(i,k,j)
1962    ENDDO
1963    ENDDO
1964    ENDDO
1965
1966 END SUBROUTINE add_c2c_u
1967
1968 SUBROUTINE add_c2c_v(lvar,rvar,config_flags,                &
1969                    ids,ide, jds, jde, kds, kde,             &
1970                    ims, ime, jms, jme, kms, kme,            &
1971                    its, ite, jts, jte, kts, kte             )
1972 !------------------------------------------------------------
1973 !------------------------------------------------------------
1974    IMPLICIT NONE
1975 !------------------------------------------------------------
1976
1977    TYPE(grid_config_rec_type),  INTENT(IN   )    :: config_flags
1978
1979    INTEGER , INTENT(IN   ) :: ids, ide, jds, jde, kds, kde, &
1980                               ims, ime, jms, jme, kms, kme, &
1981                               its, ite, jts, jte, kts, kte
1982
1983    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(IN   ) ::&
1984                                                       rvar
1985    REAL, DIMENSION(ims:ime,kms:kme,jms:jme),INTENT(INOUT) ::&
1986                                                       lvar
1987
1988 ! LOCAL VARS
1989
1990    INTEGER :: i,j,k,i_start,i_end,j_start,j_end,ktf
1991
1992    ktf=min(kte,kde-1)
1993
1994    i_start = its
1995    i_end   = MIN(ite,ide-1)
1996    j_start = jts
1997    j_end   = jte
1998
1999    IF ( config_flags%specified .or. &
2000         config_flags%nested) i_start = MAX(ids+1,its)
2001    IF ( config_flags%specified .or. &
2002         config_flags%nested) i_end   = MIN(ide-2,ite)
2003    IF ( config_flags%specified .or. &
2004         config_flags%nested) j_start = MAX(jds+1,jts)
2005    IF ( config_flags%specified .or. &
2006         config_flags%nested) j_end   = MIN(jde-1,jte)
2007
2008 !  write(*,'(a,6i4)') 'call c2cv, i_start, i_end, j_start, j_end=', i_start, i_end, j_start, j_end
2009
2010    DO j = j_start,j_end
2011    DO k = kts,kte
2012    DO i = i_start,i_end
2013       lvar(i,k,j) = lvar(i,k,j) + rvar(i,k,j)
2014    ENDDO
2015    ENDDO
2016    ENDDO
2017
2018 END SUBROUTINE add_c2c_v
2019
2020 #endif
2021
2022 END MODULE module_physics_addtendc