chem/module_data_ISRPIA.F

   1       MODULE ISRPIA
   2 ! %W% %P% %G% %U%
   3 !***********************************************************************
   4 !***********************************************************************
   5 ! Development of this code was sponsored by EPRI, 3412 Hillview Ave.,  *
   6 ! Palo Alto, CA 94304 under Contract WO8221-01                         *
   7 !                                                                      *
   8 ! Developed by Yang Zhang, Betty Pun, and Christian Seigneur,          *
   9 ! Atmospheric and Environmental Research, Inc., 2682 Bishop Drive,     *
  10 ! Suite 120, San Ramon, CA 94583                                       *
  11 !
  12 ! Development of previously available modules are listed at the        *
  13 ! begining of the code of the corresponding module.  Some of these     *
  14 ! modules may be copyrighted                                           *
  15 !***********************************************************************
  16 ! RCS file, release, date & time of last delta, author, state, [and locker]
  17 ! $Header: /ncsu/volume2/Mac_Share3/yang/madrid_final/revision/code/july2002/models/CCTM/src/aero/aero_MADRID1_mebi/ISRPIA.EXT,v 1.1 2002/10/24 00:55:29 models3 Exp $
  18
  19 !=======================================================================
  20 !
  21 ! *** ISORROPIA PLUS CODE
  22 ! *** INCLUDE FILE 'ISRPIA.EXT'
  23 ! *** THIS FILE CONTAINS THE DECLARATIONS OF THE GLOBAL CONSTANTS
  24 !     AND VARIABLES.
  25 !
  26 ! *** COPYRIGHT 1996-98, UNIVERSITY OF MIAMI, CARNEGIE MELLON UNIVERSITY
  27 ! *** WRITTEN BY ATHANASIOS NENES
  28 !
  29 ! **************************************************************************************
  30 !  FUNCTION: Include file for the ISORROPIA thermodynamic module                       *
  31 !  PRECONDITION REQUIRED:  Aerosol Option: MADRID (Development and                     *
  32 !                                     Application of the Model of Aerosol              *
  33 !                               Dynamics, Reaction, Ionization and Dissolution)        *
  34 !  RETURN VALUES:  gaseous and PM species concentrations                               *
  35 !  KEY SUBROUTINES AND FUNCTIONS CALLED:                                               *
  36 !  REVISION HISTORY:                                                                   *
  37 !         Original coagsolv code provided by ATHANASIOS NENES, Caltech (now at         *
  38 !             Georgia Tech), 2000                                                      *
  39 !         Revised by Yang Zhang and Xiao-Ming Hu, NCSU, May-June 2005                  *
  40 !             to incorporate MADRID into WRF/Chem                                      *
  41 !         Revised by Ying Pan and Yang Zhang, NCSU, Oct. 2009                          *
  42 !             to couple CB05 with MADE/SORGAM-AQ-ISORROPIA                             *
  43 ! **************************************************************************************
  44
  45 !=======================================================================
  46 !
  47 !      USE OTHR
  48 !      USE CGEN
  49 !      USE SOLN
  50 !      USE IONSdata
  51 !      USE ZSRdata
  52 !      USE SALT
  53 !      USE GAS_MADRID
  54 !      USE EQUK
  55       USE module_data_isrpia_data
  56
  57       IMPLICIT DOUBLE PRECISION (A-H, O-Z)
  58        PARAMETER (NCOMP=5, NIONS=7, NGASAQ=3, NSLDS=9, NPAIR=13,   &
  59                   NZSR=100, NERRMX=25)
  60 !      PARAMETER (NIONS=7, NSLDS=9, NERRMX=25)
  61
  62       DOUBLE PRECISION  W(NCOMP), WAER(NCOMP), TEMP, RH, IPROB
  63       INTEGER METSTBL, NADJ
  64
  65       DOUBLE PRECISION  DRH2SO4,  DRNH42S4, DRNAHSO4, DRNACL,   DRNANO3,   &
  66                         DRNA2SO4, DRNH4HS4, DRLC,     DRNH4NO3, DRNH4CL
  67
  68       DOUBLE PRECISION :: DRMLCAB,DRMLCAS,DRMASAN,DRMG1,  DRMG2,     &
  69                     DRMG3,    DRMH1,    DRMH2,    DRMI1,    DRMI2,   &
  70                     DRMI3,    DRMQ1,    DRMR1,    DRMR2,    DRMR3,   &
  71                     DRMR4,    DRMR5,    DRMR6,    DRMR7,    DRMR8,   &
  72                     DRMR9,    DRMR10,   DRMR11,   DRMR12,   DRMR13
  73       INTEGER WFTYP
  74
  75 !     DOUBLE PRECISION CH2SO4,  CNH42S4, CNH4HS4, CNACL,   CNA2SO4, &
  76 !                      CNANO3,  CNH4NO3, CNH4CL,  CNAHSO4, CLC
  77
  78 !     DOUBLE PRECISION GNH3, GHNO3, GHCL
  79
  80 !     DOUBLE PRECISION XK1, XK2, XK3, XK4, XK5, XK6, XK7, XK8, XK9, XK10,   &
  81 !                      XK11, XK12, XK13, XK14, XKW, XK21, XK22, XK31, XK32, &
  82 !                      XK41, XK42
  83
  84       CHARACTER SCASE*15
  85       DOUBLE PRECISION SULRATW, SULRAT, SODRAT
  86
  87 !      DOUBLE PRECISION EPS, MAXIT, NSWEEP, NDIV, ICLACT
  88 !     DATA   EPS/1D-6/, MAXIT/100/,NSWEEP/4/,NDIV/5/
  89
  90       CHARACTER*40 ERRMSG(NERRMX)
  91       INTEGER   ERRSTK(NERRMX), NOFER
  92       LOGICAL   STKOFL
  93       DATA ERRSTK/NERRMX*0/,   ERRMSG/NERRMX*' '/,  NOFER/0/,   &
  94            STKOFL/.FALSE./
  95
  96 !     CHARACTER VERSION*15
  97 !     DOUBLE PRECISION GREAT, TINY, TINY2, ZERO, ONE
  98 !     DATA   &
  99 !          TINY/1D-20/, GREAT/1D10/, ZERO/0.0D0/, ONE/1.0D0/,   &
 100 !          TINY2/1D-11/
 101 !       DATA VERSION /'1.7 (03/26/07)'/
 102
 103 !
 104 ! *** INPUT VARIABLES **************************************************
 105 !
 106 !      INTEGER METSTBL
 107 !      COMMON /INPT/ W(NCOMP), WAER(NCOMP), TEMP, RH, IPROB, METSTBL
 108 !
 109 ! *** WATER ACTIVITIES OF PURE SALT SOLUTIONS **************************
 110 !
 111 !      COMMON /ZSR / AWAS(NZSR), AWSS(NZSR), AWAC(NZSR), AWSC(NZSR),
 112 !     &              AWAN(NZSR), AWSN(NZSR), AWSB(NZSR), AWAB(NZSR),
 113 !     &              AWSA(NZSR), AWLC(NZSR)
 114 !
 115 ! *** DELIQUESCENCE RELATIVE HUMIDITIES ********************************
 116 !
 117 !      INTEGER WFTYP
 118 !      COMMON /DRH / DRH2SO4,  DRNH42S4, DRNAHSO4, DRNACL,   DRNANO3,
 119 !     &              DRNA2SO4, DRNH4HS4, DRLC,     DRNH4NO3, DRNH4CL
 120 !      COMMON /MDRH/ DRMLCAB,  DRMLCAS,  DRMASAN,  DRMG1,    DRMG2,
 121 !     &              DRMG3,    DRMH1,    DRMH2,    DRMI1,    DRMI2,
 122 !     &              DRMI3,    DRMQ1,    DRMR1,    DRMR2,    DRMR3,
 123 !     &              DRMR4,    DRMR5,    DRMR6,    DRMR7,    DRMR8,
 124 !     &              DRMR9,    DRMR10,   DRMR11,   DRMR12,   DRMR13,
 125 !     &              WFTYP
 126 !
 127 ! *** VARIABLES FOR LIQUID AEROSOL PHASE *******************************
 128 !
 129 !      DOUBLE PRECISION MOLAL, MOLALR, M0
 130 !      REAL IONIC
 131 !      LOGICAL CALAOU, CALAIN, FRST, DRYF
 132 !      COMMON /IONS/ MOLAL(NIONS), MOLALR(NPAIR), GAMA(NPAIR), ZZ(NPAIR),
 133 !     &              Z(NIONS),     GAMOU(NPAIR),  GAMIN(NPAIR),M0(NPAIR),
 134 !     &              GASAQ(NGASAQ),
 135 !     &              EPSACT,       COH,           CHNO3,       CHCL,
 136 !     &              WATER,        IONIC,         IACALC,
 137 !     &              FRST,         CALAIN,        CALAOU,      DRYF
 138 !
 139 ! *** VARIABLES FOR SOLID AEROSOL PHASE ********************************
 140 !
 141 !      COMMON /SALT/ CH2SO4,  CNH42S4, CNH4HS4, CNACL,   CNA2SO4,
 142 !     &              CNANO3,  CNH4NO3, CNH4CL,  CNAHSO4, CLC
 143 !
 144 ! *** VARIABLES FOR GAS PHASE ******************************************
 145 !
 146 !      COMMON /GAS / GNH3, GHNO3, GHCL
 147 !
 148 ! *** EQUILIBRIUM CONSTANTS ********************************************
 149 !
 150 !      COMMON /EQUK/ XK1, XK2, XK3, XK4, XK5, XK6, XK7, XK8, XK9, XK10,
 151 !     &              XK11, XK12, XK13, XK14, XKW, XK21, XK22, XK31, XK32,
 152 !     &              XK41, XK42
 153 !
 154 ! *** MOLECULAR WEIGHTS ************************************************
 155 !
 156 !      DOUBLE PRECISION IMW
 157 !      COMMON /OTHR/ R, IMW(NIONS), WMW(NCOMP), SMW(NPAIR)
 158 !
 159 ! *** SOLUTION/INFO VARIABLES ******************************************
 160 !
 161 !      CHARACTER SCASE*15
 162 !      COMMON /CASE/ SULRATW, SULRAT, SODRAT, SCASE
 163 !
 164 !      COMMON /SOLN/ EPS, MAXIT, NSWEEP, NDIV, ICLACT
 165 !
 166 ! *** ERROR SYSTEM *****************************************************
 167 !
 168 !      CHARACTER ERRMSG*40
 169 !      INTEGER   ERRSTK, NOFER
 170 !      LOGICAL   STKOFL
 171 !      COMMON /EROR/ STKOFL, NOFER, ERRSTK(NERRMX), ERRMSG(NERRMX)
 172 !
 173 ! *** GENERIC VARIABLES ************************************************
 174 !
 175 !      CHARACTER VERSION*14
 176 !      COMMON /CGEN/ GREAT, TINY, TINY2, ZERO, ONE, VERSION
 177
 178 !
 179 ! *** END OF INCLUDE FILE **********************************************
 180       END MODULE ISRPIA
 181 !