wrfv2_fire/chem/module_cbmz_addemiss.F

   1 !**************************************************************************
   2 ! This computer software was prepared by Battelle Memorial Institute, hereinafter
   3 ! the Contractor, under Contract No. DE-AC05-76RL0 1830 with the Department of
   4 ! Energy (DOE). NEITHER THE GOVERNMENT NOR THE CONTRACTOR MAKES ANY WARRANTY,
   5 ! EXPRESS OR IMPLIED, OR ASSUMES ANY LIABILITY FOR THE USE OF THIS SOFTWARE.
   6 !
   7 ! CBMZ module: see module_cbmz.F for references and terms of use
   8 !**********************************************************************************
   9
  10 MODULE module_cbmz_addemiss
  11 !WRF:MODEL_LAYER:CHEMICS
  12
  13
  14
  15    integer, parameter :: cbmz_addemiss_masscheck = -1
  16                        ! only do emissions masscheck calcs when this is positive
  17
  18
  19
  20 CONTAINS
  21
  22
  23
  24 !----------------------------------------------------------------------
  25    subroutine cbmz_addemiss_anthro( id, dtstep, dz8w, config_flags,       &
  26                rho_phy, chem,emis_ant,alt,                                &
  27                ids,ide, jds,jde, kds,kde,                                 &
  28                ims,ime, jms,jme, kms,kme,                                 &
  29                its,ite, jts,jte, kts,kte                                  )
  30 !
  31 ! adds emissions for cbmz trace gas species
  32 ! (i.e., emissions tendencies over time dtstep are applied
  33 ! to the trace gas concentrations)
  34 !
  35
  36   USE module_configure
  37   USE module_state_description
  38   USE module_data_radm2
  39
  40   IMPLICIT NONE
  41
  42    TYPE(grid_config_rec_type),  INTENT(IN   )    :: config_flags
  43
  44    INTEGER,      INTENT(IN   ) :: id,                                      &
  45                                   ids,ide, jds,jde, kds,kde,               &
  46                                   ims,ime, jms,jme, kms,kme,               &
  47                                   its,ite, jts,jte, kts,kte
  48
  49    REAL, INTENT(IN   ) ::    dtstep
  50
  51 ! trace species mixing ratios (gases=ppm)
  52    REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ),                 &
  53          INTENT(INOUT ) ::   chem
  54 !
  55 ! emissions arrays (v.1: ppm m/min; v.2: mole km^-2 hr^-1)
  56 !
  57 !   REAL, DIMENSION( ims:ime, kms:kme, jms:jme ),                          &
  58    REAL, DIMENSION( ims:ime, kms:config_flags%kemit, jms:jme,num_emis_ant),&
  59          INTENT(IN ) ::                                                    &
  60                          emis_ant
  61 ! layer thickness (m)
  62    REAL,  DIMENSION( ims:ime , kms:kme , jms:jme )         ,               &
  63           INTENT(IN   ) ::   dz8w, rho_phy,alt
  64
  65 ! local variables
  66    integer :: i,j,k
  67    real, parameter :: efact1 = 1.0/60.0
  68    real :: conv,conv3
  69    double precision :: chem_sum(num_chem)
  70
  71
  72 ! do mass check initial calc
  73       if (cbmz_addemiss_masscheck > 0) call addemiss_masscheck(               &
  74                id, config_flags, 1, 'cbmz_addemiss',                          &
  75                dtstep, efact1, dz8w, chem, chem_sum,                          &
  76                ids,ide, jds,jde, kds,kde,                                     &
  77                ims,ime, jms,jme, kms,kme,                                     &
  78                its,ite, jts,jte, kts,kte,                                     &
  79                21,                                                            &
  80                emis_ant(ims,kms,jms,p_e_so2),emis_ant(ims,kms,jms,p_e_no),    &
  81                emis_ant(ims,kms,jms,p_e_co),emis_ant(ims,kms,jms,p_e_eth),    &
  82                emis_ant(ims,kms,jms,p_e_hc3),emis_ant(ims,kms,jms,p_e_hc5),   &
  83                emis_ant(ims,kms,jms,p_e_hc8),emis_ant(ims,kms,jms,p_e_xyl),   &
  84                emis_ant(ims,kms,jms,p_e_ol2),emis_ant(ims,kms,jms,p_e_olt),   &
  85                emis_ant(ims,kms,jms,p_e_oli),emis_ant(ims,kms,jms,p_e_tol),   &
  86                emis_ant(ims,kms,jms,p_e_csl),emis_ant(ims,kms,jms,p_e_hcho),  &
  87                emis_ant(ims,kms,jms,p_e_ald),emis_ant(ims,kms,jms,p_e_ket),   &
  88                emis_ant(ims,kms,jms,p_e_ora2),emis_ant(ims,kms,jms,p_e_nh3),  &
  89                emis_ant(ims,kms,jms,p_e_no2),emis_ant(ims,kms,jms,p_e_ch3oh), &
  90                emis_ant(ims,kms,jms,p_e_c2h5oh))
  91
  92
  93
  94 !
  95 ! add emissions
  96 !
  97       do 100 j=jts,jte
  98       do 100 i=its,ite
  99
 100       DO k=kts,min(config_flags%kemit,kte)
 101 !v1 units:        conv = dtstep/(dz8w(i,k,j)*60.)
 102 !v2 units:
 103         conv = 4.828e-4/rho_phy(i,k,j)*dtstep/(dz8w(i,k,j)*60.)
 104         conv3 = (dtstep/dz8w(i,k,j))*alt(i,k,j)*28/250*1e-3
 105
 106 #if (defined(CHEM_DBG_I) && defined(CHEM_DBG_J) && defined(CHEM_DBG_K))
 107        if( (i <= CHEM_DBG_I .and. i >= CHEM_DBG_I) .and. &
 108            (j <= CHEM_DBG_J .and. j >= CHEM_DBG_J) .and. &
 109            (k <= CHEM_DBG_K .and. k >= CHEM_DBG_K)  ) then
 110           print*
 111           print*,"Converted emissions for CBMZ:"
 112           print*,"e_csl=",emis_ant(i,k,j,p_e_csl)*conv
 113           print*,"e_so2=",emis_ant(i,k,j,p_e_so2)*conv
 114           print*,"e_no=",emis_ant(i,k,j,p_e_no)*conv
 115           print*,"e_ald=",emis_ant(i,k,j,p_e_ald)*conv
 116           print*,"e_hcho=",emis_ant(i,k,j,p_e_hcho)*conv
 117           print*,"e_ora2=",emis_ant(i,k,j,p_e_ora2)*conv
 118           print*,"e_nh3=",emis_ant(i,k,j,p_e_nh3)*conv
 119           print*,"e_hc3=",emis_ant(i,k,j,p_e_hc3)*conv
 120           print*,"e_hc5=",emis_ant(i,k,j,p_e_hc5)*conv
 121           print*,"e_hc8=",emis_ant(i,k,j,p_e_hc8)*conv
 122           print*,"e_eth=",emis_ant(i,k,j,p_e_eth)*conv
 123           print*,"e_co=",emis_ant(i,k,j,p_e_co)*conv
 124           print*,"e_ol2=",emis_ant(i,k,j,p_e_ol2)*conv
 125           print*,"e_olt=",emis_ant(i,k,j,p_e_olt)*conv
 126           print*,"e_oli=",emis_ant(i,k,j,p_e_oli)*conv
 127           print*,"e_tol=",emis_ant(i,k,j,p_e_tol)*conv
 128           print*,"e_xyl=",emis_ant(i,k,j,p_e_xyl)*conv
 129           print*,"e_ket=",emis_ant(i,k,j,p_e_ket)*conv
 130        end if
 131 #endif
 132
 133         chem(i,k,j,p_csl)  =  chem(i,k,j,p_csl)                        &
 134                          +emis_ant(i,k,j,p_e_csl)*conv
 135         chem(i,k,j,p_so2)  = chem(i,k,j,p_so2)                         &
 136                          +emis_ant(i,k,j,p_e_so2)*conv
 137         chem(i,k,j,p_no)   = chem(i,k,j,p_no)                          &
 138                          +emis_ant(i,k,j,p_e_no)*conv
 139         chem(i,k,j,p_ald)  = chem(i,k,j,p_ald)                         &
 140                          +emis_ant(i,k,j,p_e_ald)*conv
 141         chem(i,k,j,p_hcho) = chem(i,k,j,p_hcho)                        &
 142                          +emis_ant(i,k,j,p_e_hcho)*conv
 143         chem(i,k,j,p_ora2)  = chem(i,k,j,p_ora2)                       &
 144                          +emis_ant(i,k,j,p_e_ora2)*conv
 145         chem(i,k,j,p_nh3)  = chem(i,k,j,p_nh3)                         &
 146                          +emis_ant(i,k,j,p_e_nh3)*conv
 147         chem(i,k,j,p_eth)  = chem(i,k,j,p_eth)                         &
 148                          +emis_ant(i,k,j,p_e_eth)*conv
 149         chem(i,k,j,p_co)  = chem(i,k,j,p_co)                           &
 150                          +emis_ant(i,k,j,p_e_co)*conv
 151         chem(i,k,j,p_ol2)  = chem(i,k,j,p_ol2)                         &
 152                          +emis_ant(i,k,j,p_e_ol2)*conv
 153         chem(i,k,j,p_olt)  = chem(i,k,j,p_olt)                         &
 154                          +emis_ant(i,k,j,p_e_olt)*conv
 155         chem(i,k,j,p_oli)  = chem(i,k,j,p_oli)                         &
 156                          +emis_ant(i,k,j,p_e_oli)*conv
 157         chem(i,k,j,p_tol)  = chem(i,k,j,p_tol)                         &
 158                          +emis_ant(i,k,j,p_e_tol)*conv
 159         chem(i,k,j,p_xyl)  = chem(i,k,j,p_xyl)                         &
 160                          +emis_ant(i,k,j,p_e_xyl)*conv
 161         chem(i,k,j,p_ket)  =  chem(i,k,j,p_ket)                        &
 162                          +emis_ant(i,k,j,p_e_ket)*conv
 163
 164          chem_select_2 : SELECT CASE( config_flags%chem_opt )
 165
 166          END SELECT  chem_select_2
 167
 168
 169
 170
 171 ! when emissions input file is "radm2sorg" variety, calc par emissions as a
 172 !   combination of the anthropogenic emissions for radm2 primary voc species
 173         if ( (config_flags%emiss_inpt_opt == EMISS_INPT_DEFAULT) .or.   &
 174              (config_flags%emiss_inpt_opt == EMISS_INPT_PNNL_RS) ) then
 175             chem(i,k,j,p_par) = chem(i,k,j,p_par)             &
 176                 + conv*                                       &
 177                   ( 0.4*emis_ant(i,k,j,p_e_ald) + 2.9*emis_ant(i,k,j,p_e_hc3)       &
 178                   + 4.8*emis_ant(i,k,j,p_e_hc5) + 7.9*emis_ant(i,k,j,p_e_hc8)       &
 179                   + 0.9*emis_ant(i,k,j,p_e_ket) + 2.8*emis_ant(i,k,j,p_e_oli)       &
 180                   + 1.8*emis_ant(i,k,j,p_e_olt) + 1.0*emis_ant(i,k,j,p_e_ora2) )
 181
 182 ! when emissions input file is "cbmzmosaic" variety,
 183 !   the par emissions are read into e_hc5
 184 !   and there are emissions for other species
 185         else
 186             chem(i,k,j,p_par) = chem(i,k,j,p_par)             &
 187                 + conv*emis_ant(i,k,j,p_e_hc5)
 188
 189             chem(i,k,j,p_no2) = chem(i,k,j,p_no2)             &
 190                 + conv*emis_ant(i,k,j,p_e_no2)
 191             chem(i,k,j,p_ch3oh)  = chem(i,k,j,p_ch3oh)        &
 192                 + conv*emis_ant(i,k,j,p_e_ch3oh)
 193             chem(i,k,j,p_c2h5oh) = chem(i,k,j,p_c2h5oh)       &
 194                 + conv*emis_ant(i,k,j,p_e_c2h5oh)
 195         end if
 196
 197       END DO
 198  100  continue
 199
 200
 201 ! do mass check final calc
 202       if (cbmz_addemiss_masscheck > 0) call addemiss_masscheck(               &
 203                id, config_flags, 2, 'cbmz_addemiss',                          &
 204                dtstep, efact1, dz8w, chem, chem_sum,                          &
 205                ids,ide, jds,jde, kds,kde,                                     &
 206                ims,ime, jms,jme, kms,kme,                                     &
 207                its,ite, jts,jte, kts,kte,                                     &
 208                21,                                                            &
 209                emis_ant(ims,kms,jms,p_e_so2),emis_ant(ims,kms,jms,p_e_no),    &
 210                emis_ant(ims,kms,jms,p_e_co),emis_ant(ims,kms,jms,p_e_eth),    &
 211                emis_ant(ims,kms,jms,p_e_hc3),emis_ant(ims,kms,jms,p_e_hc5),   &
 212                emis_ant(ims,kms,jms,p_e_hc8),emis_ant(ims,kms,jms,p_e_xyl),   &
 213                emis_ant(ims,kms,jms,p_e_ol2),emis_ant(ims,kms,jms,p_e_olt),   &
 214                emis_ant(ims,kms,jms,p_e_oli),emis_ant(ims,kms,jms,p_e_tol),   &
 215                emis_ant(ims,kms,jms,p_e_csl),emis_ant(ims,kms,jms,p_e_hcho),  &
 216                emis_ant(ims,kms,jms,p_e_ald),emis_ant(ims,kms,jms,p_e_ket),   &
 217                emis_ant(ims,kms,jms,p_e_ora2),emis_ant(ims,kms,jms,p_e_nh3),  &
 218                emis_ant(ims,kms,jms,p_e_no2),emis_ant(ims,kms,jms,p_e_ch3oh), &
 219                emis_ant(ims,kms,jms,p_e_c2h5oh))
 220
 221
 222    END subroutine cbmz_addemiss_anthro
 223
 224
 225
 226 !----------------------------------------------------------------------
 227   subroutine cbmz_addemiss_bio( id, dtstep, dz8w, config_flags,       &
 228         rho_phy, chem, e_bio, ne_area, e_iso,                         &
 229         ids,ide, jds,jde, kds,kde,                                    &
 230         ims,ime, jms,jme, kms,kme,                                    &
 231         its,ite, jts,jte, kts,kte                                     )
 232
 233   USE module_configure
 234   USE module_state_description
 235   USE module_data_radm2
 236   USE module_aerosols_sorgam
 237
 238   IMPLICIT NONE
 239
 240 ! subr arguments
 241    TYPE(grid_config_rec_type),  INTENT(IN   )    :: config_flags
 242
 243    INTEGER,      INTENT(IN   ) :: id, ne_area,                             &
 244                                   ids,ide, jds,jde, kds,kde,               &
 245                                   ims,ime, jms,jme, kms,kme,               &
 246                                   its,ite, jts,jte, kts,kte
 247
 248    REAL,      INTENT(IN   ) ::    dtstep
 249
 250    REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ),                 &
 251          INTENT(INOUT ) ::        chem
 252
 253    REAL, DIMENSION( ims:ime, jms:jme,ne_area ),                            &
 254          INTENT(IN ) ::           e_bio
 255
 256    REAL, DIMENSION( ims:ime, kms:config_flags%kemit, jms:jme ),            &
 257          INTENT(IN ) ::           e_iso
 258
 259    REAL,  DIMENSION( ims:ime , kms:kme , jms:jme )         ,               &
 260           INTENT(IN   ) ::        dz8w, rho_phy
 261
 262
 263 ! local variables
 264    integer i,j,k,n
 265    real, parameter :: efact1 = 1.0/60.0
 266    double precision :: chem_sum(num_chem)
 267
 268
 269 !
 270 ! apply gunther online biogenic gas emissions when bio_emiss_opt == GUNTHER1
 271 !
 272    if (config_flags%bio_emiss_opt == GUNTHER1) then
 273
 274       if (cbmz_addemiss_masscheck > 0) call addemiss_masscheck(         &
 275                id, config_flags, 1, 'cbmz_addemiss_bioaa',              &
 276                dtstep, efact1, dz8w, chem, chem_sum,                    &
 277                ids,ide, jds,jde, kds,kde,                               &
 278                ims,ime, jms,jme, kms,kme,                               &
 279                its,ite, jts,jte, kts,kte,                               &
 280                13,                                                      &
 281                e_bio(ims,jms,lald),  e_bio(ims,jms,lhc3),               &
 282                e_bio(ims,jms,lhc5),  e_bio(ims,jms,lhc8),               &
 283                e_bio(ims,jms,lhcho), e_bio(ims,jms,liso),               &
 284                e_bio(ims,jms,lket),  e_bio(ims,jms,lno),                &
 285                e_bio(ims,jms,loli),  e_bio(ims,jms,lolt),               &
 286                e_bio(ims,jms,lora1), e_bio(ims,jms,lora2),              &
 287                e_bio(ims,jms,lxyl),                                     &
 288                e_bio(ims,jms,lxyl),  e_bio(ims,jms,lxyl),               &
 289                e_bio(ims,jms,lxyl),  e_bio(ims,jms,lxyl),               &
 290                e_bio(ims,jms,lxyl),  e_bio(ims,jms,lxyl),               &
 291                e_bio(ims,jms,lxyl),  e_bio(ims,jms,lxyl)                )
 292
 293       do j=jts,jte
 294       do i=its,ite
 295         chem(i,kts,j,p_so2) = chem(i,kts,j,p_so2)    &
 296                           + e_bio(i,j,lso2)/(dz8w(i,kts,j)*60.)*dtstep
 297         chem(i,kts,j,p_sulf) = chem(i,kts,j,p_sulf)    &
 298                           + e_bio(i,j,lsulf)/(dz8w(i,kts,j)*60.)*dtstep
 299         chem(i,kts,j,p_no2) = chem(i,kts,j,p_no2)    &
 300                           + e_bio(i,j,lno2)/(dz8w(i,kts,j)*60.)*dtstep
 301         chem(i,kts,j,p_no) = chem(i,kts,j,p_no)    &
 302                           + e_bio(i,j,lno)/(dz8w(i,kts,j)*60.)*dtstep
 303         chem(i,kts,j,p_o3) = chem(i,kts,j,p_o3)    &
 304                           + e_bio(i,j,lo3)/(dz8w(i,kts,j)*60.)*dtstep
 305         chem(i,kts,j,p_hno3) = chem(i,kts,j,p_hno3)    &
 306                           + e_bio(i,j,lhno3)/(dz8w(i,kts,j)*60.)*dtstep
 307         chem(i,kts,j,p_h2o2) = chem(i,kts,j,p_h2o2)    &
 308                           + e_bio(i,j,lh2o2)/(dz8w(i,kts,j)*60.)*dtstep
 309         chem(i,kts,j,p_ald) = chem(i,kts,j,p_ald)    &
 310                           + e_bio(i,j,lald)/(dz8w(i,kts,j)*60.)*dtstep
 311         chem(i,kts,j,p_hcho) = chem(i,kts,j,p_hcho)    &
 312                           + e_bio(i,j,lhcho)/(dz8w(i,kts,j)*60.)*dtstep
 313         chem(i,kts,j,p_op1) = chem(i,kts,j,p_op1)    &
 314                           + e_bio(i,j,lop1)/(dz8w(i,kts,j)*60.)*dtstep
 315         chem(i,kts,j,p_op2) = chem(i,kts,j,p_op2)    &
 316                           + e_bio(i,j,lop2)/(dz8w(i,kts,j)*60.)*dtstep
 317         chem(i,kts,j,p_ora1) = chem(i,kts,j,p_ora1)    &
 318                           + e_bio(i,j,lora1)/(dz8w(i,kts,j)*60.)*dtstep
 319         chem(i,kts,j,p_ora2) = chem(i,kts,j,p_ora2)    &
 320                           + e_bio(i,j,lora2)/(dz8w(i,kts,j)*60.)*dtstep
 321         chem(i,kts,j,p_nh3) = chem(i,kts,j,p_nh3)    &
 322                           + e_bio(i,j,lnh3)/(dz8w(i,kts,j)*60.)*dtstep
 323         chem(i,kts,j,p_n2o5) = chem(i,kts,j,p_n2o5)    &
 324                           + e_bio(i,j,ln2o5)/(dz8w(i,kts,j)*60.)*dtstep
 325         chem(i,kts,j,p_no2) = chem(i,kts,j,p_no2)    &
 326                           + e_bio(i,j,lno2)/(dz8w(i,kts,j)*60.)*dtstep
 327         chem(i,kts,j,p_pan) = chem(i,kts,j,p_pan)    &
 328                           + e_bio(i,j,lpan)/(dz8w(i,kts,j)*60.)*dtstep
 329         chem(i,kts,j,p_eth) = chem(i,kts,j,p_eth)    &
 330                           + e_bio(i,j,leth)/(dz8w(i,kts,j)*60.)*dtstep
 331         chem(i,kts,j,p_co) = chem(i,kts,j,p_co)    &
 332                           + e_bio(i,j,lco)/(dz8w(i,kts,j)*60.)*dtstep
 333         chem(i,kts,j,p_ol2) = chem(i,kts,j,p_ol2)    &
 334                           + e_bio(i,j,lol2)/(dz8w(i,kts,j)*60.)*dtstep
 335         chem(i,kts,j,p_olt) = chem(i,kts,j,p_olt)    &
 336                           + e_bio(i,j,lolt)/(dz8w(i,kts,j)*60.)*dtstep
 337         chem(i,kts,j,p_oli) = chem(i,kts,j,p_oli)    &
 338                           + e_bio(i,j,loli)/(dz8w(i,kts,j)*60.)*dtstep
 339         chem(i,kts,j,p_tol) = chem(i,kts,j,p_tol)    &
 340                           + e_bio(i,j,ltol)/(dz8w(i,kts,j)*60.)*dtstep
 341         chem(i,kts,j,p_xyl) = chem(i,kts,j,p_xyl)    &
 342                           + e_bio(i,j,lxyl)/(dz8w(i,kts,j)*60.)*dtstep
 343         chem(i,kts,j,p_hono) = chem(i,kts,j,p_hono)    &
 344                           + e_bio(i,j,lhono)/(dz8w(i,kts,j)*60.)*dtstep
 345         chem(i,kts,j,p_hno4) = chem(i,kts,j,p_hno4)    &
 346                           + e_bio(i,j,lhno4)/(dz8w(i,kts,j)*60.)*dtstep
 347         chem(i,kts,j,p_ket) = chem(i,kts,j,p_ket)    &
 348                           + e_bio(i,j,lket)/(dz8w(i,kts,j)*60.)*dtstep
 349         chem(i,kts,j,p_mgly) = chem(i,kts,j,p_mgly)    &
 350                           + e_bio(i,j,lmgly)/(dz8w(i,kts,j)*60.)*dtstep
 351         chem(i,kts,j,p_onit) = chem(i,kts,j,p_onit)    &
 352                           + e_bio(i,j,lonit)/(dz8w(i,kts,j)*60.)*dtstep
 353         chem(i,kts,j,p_csl) = chem(i,kts,j,p_csl)    &
 354                           + e_bio(i,j,lcsl)/(dz8w(i,kts,j)*60.)*dtstep
 355         chem(i,kts,j,p_iso) = chem(i,kts,j,p_iso)    &
 356                           + e_bio(i,j,liso)/(dz8w(i,kts,j)*60.)*dtstep
 357       end do
 358       end do
 359
 360 ! calc par emissions as a combination of the biogenic emissions
 361 ! for radm2 primary voc species
 362       do j = jts, jte
 363       do i = its, ite
 364          chem(i,kts,j,p_par)  =  chem(i,kts,j,p_par)               &
 365              + (dtstep/(dz8w(i,kts,j)*60.))*                       &
 366                ( 0.4*e_bio(i,j,lald) + 2.9*e_bio(i,j,lhc3)         &
 367                + 4.8*e_bio(i,j,lhc5) + 7.9*e_bio(i,j,lhc8)         &
 368                + 0.9*e_bio(i,j,lket) + 2.8*e_bio(i,j,loli)         &
 369                + 1.8*e_bio(i,j,lolt) + 1.0*e_bio(i,j,lora2)        )
 370       end do
 371       end do
 372
 373       if (cbmz_addemiss_masscheck > 0) call addemiss_masscheck(         &
 374                id, config_flags, 2, 'cbmz_addemiss_bioaa',              &
 375                dtstep, efact1, dz8w, chem, chem_sum,                    &
 376                ids,ide, jds,jde, kds,kde,                               &
 377                ims,ime, jms,jme, kms,kme,                               &
 378                its,ite, jts,jte, kts,kte,                               &
 379                13,                                                      &
 380                e_bio(ims,jms,lald),  e_bio(ims,jms,lhc3),               &
 381                e_bio(ims,jms,lhc5),  e_bio(ims,jms,lhc8),               &
 382                e_bio(ims,jms,lhcho), e_bio(ims,jms,liso),               &
 383                e_bio(ims,jms,lket),  e_bio(ims,jms,lno),                &
 384                e_bio(ims,jms,loli),  e_bio(ims,jms,lolt),               &
 385                e_bio(ims,jms,lora1), e_bio(ims,jms,lora2),              &
 386                e_bio(ims,jms,lxyl),                                     &
 387                e_bio(ims,jms,lxyl),  e_bio(ims,jms,lxyl),               &
 388                e_bio(ims,jms,lxyl),  e_bio(ims,jms,lxyl),               &
 389                e_bio(ims,jms,lxyl),  e_bio(ims,jms,lxyl),               &
 390                e_bio(ims,jms,lxyl),  e_bio(ims,jms,lxyl)                )
 391
 392    end if
 393
 394
 395 !
 396 ! apply offline isoprene emissions when bio_emiss_opt /= GUNTHER1
 397 !
 398    if (config_flags%bio_emiss_opt /= GUNTHER1) then
 399
 400       if (cbmz_addemiss_masscheck > 0) call addemiss_masscheck(            &
 401                id, config_flags, 1, 'cbmz_addemiss_biobb',                 &
 402                dtstep, efact1, dz8w, chem, chem_sum,                       &
 403                ids,ide, jds,jde, kds,kde,                                  &
 404                ims,ime, jms,jme, kms,kme,                                  &
 405                its,ite, jts,jte, kts,kte,                                  &
 406                1,                                                          &
 407                e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,                  &
 408                e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,                  &
 409                e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,e_iso                   )
 410
 411       do j = jts, jte
 412       do k = kts, min(config_flags%kemit,kte)
 413       do i = its, ite
 414          chem(i,k,j,p_iso) = chem(i,k,j,p_iso) + e_iso(i,k,j)              &
 415               *4.828e-4/rho_phy(i,k,j)*(dtstep/(dz8w(i,k,j)*60.))
 416       end do
 417       end do
 418       end do
 419
 420       if (cbmz_addemiss_masscheck > 0) call addemiss_masscheck(            &
 421                id, config_flags, 2, 'cbmz_addemiss_biobb',                 &
 422                dtstep, efact1, dz8w, chem, chem_sum,                       &
 423                ids,ide, jds,jde, kds,kde,                                  &
 424                ims,ime, jms,jme, kms,kme,                                  &
 425                its,ite, jts,jte, kts,kte,                                  &
 426                1,                                                          &
 427                e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,                  &
 428                e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,                  &
 429                e_iso,e_iso,e_iso,e_iso,e_iso,e_iso,e_iso                   )
 430
 431    end if
 432
 433
 434    END subroutine cbmz_addemiss_bio
 435
 436
 437 END MODULE module_cbmz_addemiss
 438
 439
 440