Line data    Source code

       1             : ! Include shortname defintions, so that the F77 code does not have to be modified to
       2             : ! reference the CARMA structure.
       3             : #include "carma_globaer.h"
       4             : 
       5             : !! This determines the maximum particle concentration for each group in each
       6             : !! gridbox. This can be used to make calculations more efficient by skipping
       7             : !! calculations when concentrations are low
       8             : !!
       9             : !! @author Chuck Bardeen
      10             : !! @version Nov 2009
      11  1565790798 : subroutine maxconc(carma, cstate, iz, rc)
      12             : 
      13             :         ! types
      14             :   use carma_precision_mod
      15             :   use carma_enums_mod
      16             :   use carma_constants_mod
      17             :   use carma_types_mod
      18             :   use carmastate_mod
      19             :   use carma_mod
      20             : 
      21             :   implicit none
      22             : 
      23             :   type(carma_type), intent(in)         :: carma   !! the carma object
      24             :   type(carmastate_type), intent(inout) :: cstate  !! the carma state object
      25             :   integer, intent(in)                  :: iz      !! z index
      26             :   integer, intent(inout)               :: rc      !! return code, negative indicates failure
      27             : 
      28             :         ! Locals
      29             :         integer       :: igrp
      30             :         integer       :: iep
      31             : 
      32             : 
      33             :   ! Find maximum particle concentration for each spatial grid box
      34             :   !  (in units of cm^-3)
      35  4697372394 :   do igrp = 1,NGROUP
      36  3131581596 :     iep = ienconc(igrp)
      37             : 
      38 68894795112 :     pconmax(iz,igrp) = maxval(pc(iz,:,iep))
      39             : 
      40           0 :     pconmax(iz,igrp) = pconmax(iz,igrp) &
      41  4697372394 :                             / zmet(iz)
      42             :   enddo  ! igrp
      43             : 
      44  1565790798 :   return
      45  1565790798 : end