chem/module_mixactivate_wrappers.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 ! MOSAIC module: see module_mosaic_driver.F for references and terms of use
   8 !**********************************************************************************
   9
  10 !----------------------------------------------------------------------
  11 ! This module contains interface wrapper routines to couple the aerosol
  12 ! modules with mixactivate in the physics directory. Due to compiling
  13 ! dependencies, these cannot be placed in module_mixactivate.
  14 !----------------------------------------------------------------------
  15
  16 MODULE module_mixactivate_wrappers
  17
  18 CONTAINS
  19
  20 !----------------------------------------------------------------------
  21 !----------------------------------------------------------------------
  22       subroutine mosaic_mixactivate (                        &
  23            id, ktau, dtstep, config_flags, idrydep_onoff,    &
  24            rho_phy, t_phy, w, cldfra, cldfra_old,            &
  25            ddvel, z, dz8w, p_at_w, t_at_w, exch_h,           &
  26            qv, qc, qi, qndrop3d, f_qc, f_qi, chem,           &
  27                ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource,      &
  28                qsrflx, &
  29            ids,ide, jds,jde, kds,kde,                        &
  30            ims,ime, jms,jme, kms,kme,                        &
  31            its,ite, jts,jte, kts,kte                         )
  32
  33     USE module_configure, only: grid_config_rec_type
  34         use module_state_description, only:  num_chem
  35         use module_data_mosaic_asect
  36         use module_mixactivate, only: mixactivate
  37
  38 ! wrapper to call mixactivate for mosaic description of aerosol
  39
  40         implicit none
  41
  42 !   subr arguments
  43         integer, intent(in) ::               &
  44          id, ktau,                       &
  45          ids, ide, jds, jde, kds, kde,   &
  46          ims, ime, jms, jme, kms, kme,   &
  47          its, ite, jts, jte, kts, kte,   &
  48          idrydep_onoff
  49
  50         real, intent(in) :: dtstep
  51
  52         real, intent(in),   &
  53                 dimension( ims:ime, kms:kme, jms:jme ) :: &
  54                 rho_phy, t_phy, w,   &
  55                 z, dz8w, p_at_w, t_at_w, exch_h
  56
  57         real, intent(inout),   &
  58                 dimension( ims:ime, kms:kme, jms:jme ) :: cldfra, cldfra_old
  59
  60         real, intent(in),   &
  61                 dimension( its:ite, jts:jte, num_chem ) :: ddvel
  62
  63         real, intent(in),   &
  64                 dimension( ims:ime, kms:kme, jms:jme ) :: &
  65                 qv, qc, qi
  66
  67     LOGICAL, intent(in) :: f_qc, f_qi
  68
  69         real, intent(inout),   &
  70                 dimension( ims:ime, kms:kme, jms:jme ) :: &
  71                 qndrop3d
  72
  73         real, intent(inout),   &
  74                 dimension( ims:ime, kms:kme, jms:jme, 1:num_chem ) :: &
  75                 chem
  76       real, intent(out), dimension(ims:ime,kms:kme,jms:jme) :: nsource,&
  77              ccn1,ccn2,ccn3,ccn4,ccn5,ccn6  ! number conc of aerosols activated at supersat
  78
  79         type(grid_config_rec_type), intent(in) :: config_flags
  80  real, intent(out) :: qsrflx(ims:ime, jms:jme, num_chem) ! wet deposition flux of aerosol
  81 ! local vars
  82 !       real qsrflx(ims:ime, jms:jme, num_chem) ! wet deposition flux of aerosol
  83         real sumhygro,sumvol
  84         integer i,j,k,l,m,n
  85         real hygro( its:ite, kts:kte, jts:jte, maxd_asize, maxd_atype ) ! bulk
  86
  87       qsrflx(:,:,:) = 0.0
  88
  89 ! calculate volume-weighted bulk hygroscopicity for each type and size
  90       do 100 j=jts,jte
  91       do 100 k=kts,kte
  92       do 100 i=its,ite
  93        do n=1,ntype_aer
  94        do m=1,nsize_aer(n)
  95                sumhygro=0.
  96                sumvol=0.
  97                do l=1,ncomp_aer(n)
  98                   sumhygro = sumhygro+hygro_aer(l,n)*   &
  99                    chem(i,k,j,massptr_aer(l,m,n,ai_phase))/dens_aer(l,n)
 100                   sumvol = sumvol+chem(i,k,j,massptr_aer(l,m,n,ai_phase))/dens_aer(l,n)
 101                end do ! comp
 102            hygro(i,k,j,m,n)=sumhygro/sumvol
 103         end do ! size
 104         end do ! type
 105   100 continue
 106
 107 ! check arguments of mixactivate for consistency between send, receive
 108 ! 06-nov-2005 rce - id & ktau added to arg list
 109       call mixactivate(  msectional, &
 110            chem, num_chem, qv, qc, qi, qndrop3d,   &
 111            t_phy, w, ddvel, idrydep_onoff,  &
 112            maxd_acomp, maxd_asize, maxd_atype, maxd_aphase,   &
 113            ncomp_aer, nsize_aer, ntype_aer, nphase_aer,  &
 114            numptr_aer, massptr_aer, dlo_sect, dhi_sect, sigmag_aer, dcen_sect,  &
 115            dens_aer, mw_aer,           &
 116            waterptr_aer, hygro,  ai_phase, cw_phase,                &
 117            ids,ide, jds,jde, kds,kde,                            &
 118            ims,ime, jms,jme, kms,kme,                            &
 119            its,ite, jts,jte, kts,kte,                            &
 120            rho_phy, z, dz8w, p_at_w, t_at_w, exch_h,      &
 121            cldfra, cldfra_old, qsrflx, &
 122                ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource,       &
 123                id, ktau, dtstep, &
 124            f_qc, f_qi               )
 125
 126       end subroutine mosaic_mixactivate
 127
 128
 129 !----------------------------------------------------------------------
 130 !----------------------------------------------------------------------
 131
 132       subroutine mosaic_mixactivate_init(                    &
 133            config_flags, chem, scalar,                       &
 134            chem_in_opt,                                      & !BSINGH(12/04/13):  Added for SAPRC 8 bin vbs
 135            ims,ime, jms,jme, kms,kme,                        &
 136            its,ite, jts,jte, kts,kte                         )
 137
 138       USE module_configure, only: grid_config_rec_type
 139       use module_state_description, only:  num_chem, num_scalar, p_qndrop
 140       use module_data_mosaic_asect
 141
 142         implicit none
 143
 144 ! subr arguments
 145       type(grid_config_rec_type), intent(in) :: config_flags
 146
 147       integer, intent(in) ::               &
 148            ims, ime, jms, jme, kms, kme,   &
 149            its, ite, jts, jte, kts, kte
 150       INTEGER,      INTENT(IN   ) :: chem_in_opt !BSINGH(12/04/13):  Added for SAPRC 8 bin vbs
 151       real, intent(inout),   &
 152            dimension( ims:ime, kms:kme, jms:jme, 1:num_chem ) :: &
 153            chem
 154
 155       real, intent(inout),   &
 156            dimension( ims:ime, kms:kme, jms:jme, 1:num_scalar ) :: &
 157            scalar
 158
 159       integer :: i, j, k, m, n, l
 160
 161       do j=jts,jte
 162          do k=kts,kte
 163             do i=its,ite
 164                scalar(i,k,j,p_qndrop) = 0.
 165             end do
 166          end do
 167       end do
 168
 169       if( cw_phase > 0 ) then   !sanity check in case using prognostic
 170                                 !drop number without aq. chemistry
 171          if (config_flags%chem_in_opt == 1) then !added by MS to transfer cloud borne species to interstitial species when using chem_in_opt=1
 172             do n=1,ntype_aer
 173                do m=1,nsize_aer(n)
 174                   chem(its:ite,kts:kte,jts:jte,numptr_aer(m,n,cw_phase)) = 0.
 175                   do l=1,ncomp_aer(n)
 176                      if( ai_phase > 0 ) then
 177                         ! add cloud borne aerosol mass to interstitial aerosol mass
 178                         chem(its:ite,kts:kte,jts:jte,massptr_aer(l,m,n,ai_phase))= &
 179                              chem(its:ite,kts:kte,jts:jte,massptr_aer(l,m,n,ai_phase)) + &
 180                              chem(its:ite,kts:kte,jts:jte,massptr_aer(l,m,n,cw_phase))
 181
 182                      endif ! for ai_phase>0
 183                      chem(its:ite,kts:kte,jts:jte,massptr_aer(l,m,n,cw_phase)) = 0.
 184                   end do              ! comp
 185                end do                 ! size
 186             end do
 187          else
 188             do n=1,ntype_aer
 189                do m=1,nsize_aer(n)
 190                   chem(its:ite,kts:kte,jts:jte,numptr_aer(m,n,cw_phase)) = 0.
 191                   do l=1,ncomp_aer(n)
 192                      chem(its:ite,kts:kte,jts:jte,massptr_aer(l,m,n,cw_phase)) = 0.
 193                   end do              ! comp
 194                end do                 ! size
 195             end do                 ! type
 196          endif !chem_in_opt
 197       end if
 198
 199       end subroutine mosaic_mixactivate_init
 200
 201
 202 !----------------------------------------------------------------------
 203 !----------------------------------------------------------------------
 204
 205
 206       subroutine sorgam_mixactivate (                        &
 207            id, ktau, dtstep, config_flags, idrydep_onoff,    &
 208            rho_phy, t_phy, w, cldfra, cldfra_old,            &
 209            ddvel, z, dz8w, p_at_w, t_at_w, exch_h,           &
 210            qv, qc, qi, qndrop3d, f_qc, f_qi, chem,           &
 211                ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource,      &
 212            ids,ide, jds,jde, kds,kde,                        &
 213            ims,ime, jms,jme, kms,kme,                        &
 214            its,ite, jts,jte, kts,kte                         )
 215
 216     USE module_configure, only: grid_config_rec_type
 217         use module_state_description, only:  num_chem
 218         use module_data_sorgam
 219         use module_mixactivate, only: mixactivate
 220
 221 ! wrapper to call mixactivate for sorgam description of aerosol
 222
 223         implicit none
 224
 225 !   subr arguments
 226         integer, intent(in) ::                  &
 227                 id, ktau,                       &
 228                 ids, ide, jds, jde, kds, kde,   &
 229                 ims, ime, jms, jme, kms, kme,   &
 230                 its, ite, jts, jte, kts, kte,   &
 231                 idrydep_onoff
 232
 233         real, intent(in) :: dtstep
 234
 235         real, intent(in),   &
 236                 dimension( ims:ime, kms:kme, jms:jme ) :: &
 237                 rho_phy, t_phy, w,   &
 238                 z, dz8w, p_at_w, t_at_w, exch_h
 239
 240         real, intent(inout),   &
 241                 dimension( ims:ime, kms:kme, jms:jme ) :: cldfra, cldfra_old
 242
 243         real, intent(in),   &
 244                 dimension( its:ite, jts:jte, num_chem ) :: ddvel
 245
 246         real, intent(in),   &
 247                 dimension( ims:ime, kms:kme, jms:jme ) :: &
 248                 qv, qc, qi
 249
 250     LOGICAL, intent(in) :: f_qc, f_qi
 251
 252         real, intent(inout),   &
 253                 dimension( ims:ime, kms:kme, jms:jme  ) :: &
 254                 qndrop3d
 255
 256         real, intent(inout),   &
 257                 dimension( ims:ime, kms:kme, jms:jme, 1:num_chem ) :: &
 258                 chem
 259       real, intent(out), dimension(ims:ime,kms:kme,jms:jme) :: nsource, &
 260              ccn1,ccn2,ccn3,ccn4,ccn5,ccn6  ! number conc of aerosols activated at supersat
 261
 262         type(grid_config_rec_type), intent(in) :: config_flags
 263
 264 ! local vars
 265         real qsrflx(ims:ime, jms:jme, num_chem) ! dry deposition flux of aerosol
 266         real sumhygro,sumvol
 267         integer i,j,k,l,m,n
 268         real hygro( its:ite, kts:kte, jts:jte,maxd_asize, maxd_atype )
 269
 270 ! calculate volume-weighted bulk hygroscopicity for each type and size
 271
 272       do 100 j=jts,jte
 273       do 100 k=kts,kte
 274       do 100 i=its,ite
 275        do n=1,ntype_aer
 276        do m=1,nsize_aer(n)
 277                sumhygro=0
 278                sumvol=0
 279                do l=1,ncomp_aer(n)
 280                   sumhygro = sumhygro+hygro_aer(l,n)*   &
 281                    chem(i,k,j,massptr_aer(l,m,n,ai_phase))/dens_aer(l,n)
 282                   sumvol = sumvol+chem(i,k,j,massptr_aer(l,m,n,ai_phase))/dens_aer(l,n)
 283                end do ! comp
 284                hygro(i,k,j,m,n)=sumhygro/sumvol
 285         end do ! size
 286         end do ! type
 287   100 continue
 288
 289
 290 ! check arguments of mixactivate for consistency between send, receive
 291 ! 06-nov-2005 rce - id & ktau added to arg list
 292       call mixactivate(  msectional, &
 293            chem, num_chem, qv, qc, qi, qndrop3d,   &
 294            t_phy, w, ddvel, idrydep_onoff,  &
 295            maxd_acomp, maxd_asize, maxd_atype, maxd_aphase,   &
 296            ncomp_aer, nsize_aer, ntype_aer, nphase_aer,  &
 297            numptr_aer, massptr_aer, dlo_sect, dhi_sect, sigmag_aer, dcen_sect,  &
 298            dens_aer, mw_aer,           &
 299            waterptr_aer, hygro,  ai_phase, cw_phase,                 &
 300            ids,ide, jds,jde, kds,kde,                            &
 301            ims,ime, jms,jme, kms,kme,                            &
 302            its,ite, jts,jte, kts,kte,                            &
 303            rho_phy, z, dz8w, p_at_w, t_at_w, exch_h,      &
 304            cldfra, cldfra_old, qsrflx,                      &
 305                ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource,       &
 306                id, ktau, dtstep, &
 307            f_qc, f_qi               )
 308
 309       end subroutine sorgam_mixactivate
 310
 311 !!! TUCCELLA
 312       subroutine soa_vbs_mixactivate (                       &
 313            id, ktau, dtstep, config_flags, idrydep_onoff,    &
 314            rho_phy, t_phy, w, cldfra, cldfra_old,            &
 315            ddvel, z, dz8w, p_at_w, t_at_w, exch_h,           &
 316            qv, qc, qi, qndrop3d, f_qc, f_qi, chem,           &
 317            ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource,      &
 318            ids,ide, jds,jde, kds,kde,                        &
 319            ims,ime, jms,jme, kms,kme,                        &
 320            its,ite, jts,jte, kts,kte                         )
 321
 322     USE module_configure, only: grid_config_rec_type
 323         use module_state_description, only:  num_chem
 324         use module_data_soa_vbs
 325         use module_mixactivate, only: mixactivate
 326
 327 ! wrapper to call mixactivate for sorgam description of aerosol
 328
 329         implicit none
 330
 331 !   subr arguments
 332         integer, intent(in) ::                  &
 333                 id, ktau,                       &
 334                 ids, ide, jds, jde, kds, kde,   &
 335                 ims, ime, jms, jme, kms, kme,   &
 336                 its, ite, jts, jte, kts, kte,   &
 337                 idrydep_onoff
 338
 339         real, intent(in) :: dtstep
 340
 341         real, intent(in),   &
 342                 dimension( ims:ime, kms:kme, jms:jme ) :: &
 343                 rho_phy, t_phy, w,   &
 344                 z, dz8w, p_at_w, t_at_w, exch_h
 345
 346         real, intent(inout),   &
 347                 dimension( ims:ime, kms:kme, jms:jme ) :: cldfra, cldfra_old
 348
 349         real, intent(in),   &
 350                 dimension( its:ite, jts:jte, num_chem ) :: ddvel
 351
 352         real, intent(in),   &
 353                 dimension( ims:ime, kms:kme, jms:jme ) :: &
 354                 qv, qc, qi
 355
 356     LOGICAL, intent(in) :: f_qc, f_qi
 357
 358         real, intent(inout),   &
 359                 dimension( ims:ime, kms:kme, jms:jme  ) :: &
 360                 qndrop3d
 361
 362         real, intent(inout),   &
 363                 dimension( ims:ime, kms:kme, jms:jme, 1:num_chem ) :: &
 364                 chem
 365       real, intent(out), dimension(ims:ime,kms:kme,jms:jme) :: nsource, &
 366              ccn1,ccn2,ccn3,ccn4,ccn5,ccn6  ! number conc of aerosols activated at supersat
 367
 368         type(grid_config_rec_type), intent(in) :: config_flags ! local vars
 369         real qsrflx(ims:ime, jms:jme, num_chem) ! dry deposition flux of aerosol
 370         real sumhygro,sumvol
 371         integer i,j,k,l,m,n
 372         real hygro( its:ite, kts:kte, jts:jte,maxd_asize, maxd_atype )
 373
 374 ! calculate volume-weighted bulk hygroscopicity for each type and size
 375
 376       do 100 j=jts,jte
 377       do 100 k=kts,kte
 378       do 100 i=its,ite
 379        do n=1,ntype_aer
 380        do m=1,nsize_aer(n)
 381                sumhygro=0
 382                sumvol=0
 383                do l=1,ncomp_aer(n)
 384                   sumhygro = sumhygro+hygro_aer(l,n)*   &
 385                              chem(i,k,j,massptr_aer(l,m,n,ai_phase))/dens_aer(l,n)
 386                   sumvol = sumvol+chem(i,k,j,massptr_aer(l,m,n,ai_phase))/dens_aer(l,n)
 387                end do ! comp
 388                hygro(i,k,j,m,n)=sumhygro/sumvol
 389         end do ! size
 390         end do ! type
 391   100 continue
 392
 393
 394 ! check arguments of mixactivate for consistency between send, receive
 395 ! 06-nov-2005 rce - id & ktau added to arg list
 396       call mixactivate(  msectional, &
 397            chem, num_chem, qv, qc, qi, qndrop3d,   &
 398            t_phy, w, ddvel, idrydep_onoff,  &
 399            maxd_acomp, maxd_asize, maxd_atype, maxd_aphase,   &
 400            ncomp_aer, nsize_aer, ntype_aer, nphase_aer,  &
 401            numptr_aer, massptr_aer, dlo_sect, dhi_sect, sigmag_aer, dcen_sect, &
 402            dens_aer, mw_aer,           &
 403            waterptr_aer, hygro,  ai_phase, cw_phase,                 &
 404            ids,ide, jds,jde, kds,kde,                            &
 405            ims,ime, jms,jme, kms,kme,                            &
 406            its,ite, jts,jte, kts,kte,                            &
 407            rho_phy, z, dz8w, p_at_w, t_at_w, exch_h,      &
 408            cldfra, cldfra_old, qsrflx,                      &
 409                ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource,       &
 410                id, ktau, dtstep, &
 411            f_qc, f_qi               )
 412
 413       end subroutine soa_vbs_mixactivate
 414
 415       subroutine sorgam_vbs_mixactivate (                        &
 416            id, ktau, dtstep, config_flags, idrydep_onoff,    &
 417            rho_phy, t_phy, w, cldfra, cldfra_old,            &
 418            ddvel, z, dz8w, p_at_w, t_at_w, exch_h,           &
 419            qv, qc, qi, qndrop3d, f_qc, f_qi, chem,           &
 420            ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource,      &
 421            ids,ide, jds,jde, kds,kde,                        &
 422            ims,ime, jms,jme, kms,kme,                        &
 423            its,ite, jts,jte, kts,kte                         )
 424
 425     USE module_configure, only: grid_config_rec_type
 426         use module_state_description, only:  num_chem
 427         use module_data_sorgam_vbs
 428         use module_mixactivate, only: mixactivate
 429
 430 ! wrapper to call mixactivate for sorgam description of aerosol
 431
 432         implicit none
 433
 434 !   subr arguments
 435         integer, intent(in) ::                  &
 436                 id, ktau,                       &
 437                 ids, ide, jds, jde, kds, kde,   &
 438                 ims, ime, jms, jme, kms, kme,   &
 439                 its, ite, jts, jte, kts, kte,   &
 440                 idrydep_onoff
 441
 442         real, intent(in) :: dtstep
 443
 444         real, intent(in),   &
 445                 dimension( ims:ime, kms:kme, jms:jme ) :: &
 446                 rho_phy, t_phy, w,   &
 447                 z, dz8w, p_at_w, t_at_w, exch_h
 448
 449         real, intent(inout),   &
 450                 dimension( ims:ime, kms:kme, jms:jme ) :: cldfra, cldfra_old
 451
 452         real, intent(in),   &
 453                 dimension( its:ite, jts:jte, num_chem ) :: ddvel
 454
 455         real, intent(in),   &
 456                 dimension( ims:ime, kms:kme, jms:jme ) :: &
 457                 qv, qc, qi
 458
 459     LOGICAL, intent(in) :: f_qc, f_qi
 460
 461         real, intent(inout),   &
 462                 dimension( ims:ime, kms:kme, jms:jme  ) :: &
 463                 qndrop3d
 464
 465         real, intent(inout),   &
 466                 dimension( ims:ime, kms:kme, jms:jme, 1:num_chem ) :: &
 467                 chem
 468       real, intent(out), dimension(ims:ime,kms:kme,jms:jme) :: nsource, &
 469              ccn1,ccn2,ccn3,ccn4,ccn5,ccn6  ! number conc of aerosols activated at supersat
 470
 471         type(grid_config_rec_type), intent(in) :: config_flags
 472
 473 ! local vars
 474         real qsrflx(ims:ime, jms:jme, num_chem) ! dry deposition flux of aerosol
 475         real sumhygro,sumvol
 476         integer i,j,k,l,m,n
 477         real hygro( its:ite, kts:kte, jts:jte,maxd_asize, maxd_atype )
 478
 479 ! calculate volume-weighted bulk hygroscopicity for each type and size
 480
 481       do 100 j=jts,jte
 482       do 100 k=kts,kte
 483       do 100 i=its,ite
 484        do n=1,ntype_aer
 485        do m=1,nsize_aer(n)
 486                sumhygro=0
 487                sumvol=0
 488                do l=1,ncomp_aer(n)
 489                   sumhygro = sumhygro+hygro_aer(l,n)*   &
 490                    chem(i,k,j,massptr_aer(l,m,n,ai_phase))/dens_aer(l,n)
 491                   sumvol = sumvol+chem(i,k,j,massptr_aer(l,m,n,ai_phase))/dens_aer(l,n)
 492                end do ! comp
 493                hygro(i,k,j,m,n)=sumhygro/sumvol
 494         end do ! size
 495         end do ! type
 496   100 continue
 497
 498
 499 ! check arguments of mixactivate for consistency between send, receive
 500 ! 06-nov-2005 rce - id & ktau added to arg list
 501       call mixactivate(  msectional, &
 502            chem, num_chem, qv, qc, qi, qndrop3d,   &
 503            t_phy, w, ddvel, idrydep_onoff,  &
 504            maxd_acomp, maxd_asize, maxd_atype, maxd_aphase,   &
 505            ncomp_aer, nsize_aer, ntype_aer, nphase_aer,  &
 506            numptr_aer, massptr_aer, dlo_sect, dhi_sect, sigmag_aer, dcen_sect,  &
 507            dens_aer, mw_aer,           &
 508            waterptr_aer, hygro,  ai_phase, cw_phase,                 &
 509            ids,ide, jds,jde, kds,kde,                            &
 510            ims,ime, jms,jme, kms,kme,                            &
 511            its,ite, jts,jte, kts,kte,                            &
 512            rho_phy, z, dz8w, p_at_w, t_at_w, exch_h,      &
 513            cldfra, cldfra_old, qsrflx,                      &
 514                ccn1, ccn2, ccn3, ccn4, ccn5, ccn6, nsource,       &
 515                id, ktau, dtstep, &
 516            f_qc, f_qi               )
 517
 518       end subroutine sorgam_vbs_mixactivate
 519
 520 END MODULE module_mixactivate_wrappers