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 routine zeroes the fast microphysics sinks and sources, 
       6             : !! at one spatial point per call.
       7             : !!
       8             : !! @author Andy Ackerman
       9             : !! @version Oct-1997
      10   213271518 : subroutine zeromicro(carma, cstate, iz, rc)
      11             : 
      12             :   ! types
      13             :   use carma_precision_mod
      14             :   use carma_enums_mod
      15             :   use carma_constants_mod
      16             :   use carma_types_mod
      17             :   use carmastate_mod
      18             :   use carma_mod
      19             : 
      20             :   implicit none
      21             : 
      22             :   type(carma_type), intent(in)         :: carma   !! the carma object
      23             :   type(carmastate_type), intent(inout) :: cstate  !! the carma state object
      24             :   integer, intent(in)                  :: iz      !! vertical index
      25             :   integer, intent(inout)               :: rc      !! return code, negative indicates failure
      26             :   
      27             : 
      28             :   ! Set production terms and loss rates due to nucleation, growth,
      29             :   ! and evaporation to zero.  Also set index of smallest bin nuceleated
      30             :   ! during time step equal to <NBIN> first time through spatial loop.
      31             :   
      32   213271518 :   if (do_grow) then
      33             : 
      34   213271518 :     phprod         = 0._f
      35   213271518 :     rlprod         = 0._f
      36  9170675274 :     dtpart(iz,:,:) = 0._f
      37             : 
      38   639814554 :     if (NGAS > 0) gasprod(:) = 0._f
      39             : 
      40 49478992176 :     rhompe(:, :)  = 0._f
      41 49478992176 :     rnucpe(:,:)   = 0._f
      42 49478992176 :     growpe(:,:)   = 0._f
      43 49478992176 :     evappe(:,:)   = 0._f
      44 18554622066 :     rnuclg(:,:,:) = 0._f
      45  9170675274 :     growlg(:,:)   = 0._f
      46  9170675274 :     evaplg(:,:)   = 0._f
      47             : 
      48             :   end if
      49             : 
      50             :   ! Return to caller with fast microphysics sinks and sources zeroed.
      51   213271518 :   return
      52   213271518 : end