************************************************************************ ** ** ** CE-QUAL-W2-INC ** ** ** ** The INCLUDE file for CE-QUAL-W2 ** ** ** ** Version 3.0 ** ** March 1998 ** ** ** ** Developed by: Thomas M. Cole ** ** Water Quality Modeling Group ** ** U.S. Army Corps of Engineers ** ** Waterways Experiment Station ** ** Vicksburg, Mississippi 39180 ** ** Phone number: (601) 634-3283 ** ** ** ** Modified to work with version 3 by Jim Bowen UNC Charlotte ** ** ** ************************************************************************ ************************************************************************ ** ** ** This file contains the PARAMETER specifications for CE-QUAL-W2 ** ** and the EQUIVALENCEs between the constituent array and the var- ** ** iables used in the kinetics subroutines. ** ** ** ** The PARAMETERS are defined as follows: ** ** ** ** CONFN - control filename ** ** IMP - total number of segments ** ** IMC - total number of segments ** ** KMP - total number of layers ** ** KMC - total number of layers ** ** NAP - total number of algal groups ** ** NBP - total number of branches ** ** NCP - total number of constituents ** ** NDP - total number of derived constituents ** ** NFP - total number of kinetic flux variables ** ** NHP - total number of hydrodynamic variables ** ** NPP - Maximum number of dates for I/O parameters ** ** (see manual, page __) ** ** NRP - ? ** NSP - total number of outlet structures ** ** NTP - total number of tributaries ** ** NWP - total number of withdrawals ** ** ** ** Notes: ** ** ** ** 1. The variables IMC and KMC may be set to 1 in order to save ** ** memory if constituents are not being modeled. If consti- ** ** tuents are modeled, then IMC and KMC should be set equal ** ** to IMP and KMP. ** ** 2. Do not change the value of NCP unless more water quality ** ** state variables have been added to the model. ** ** 3. No parameters may be set to zero. They must be set to at ** ** least 1. ** ** ** ************************************************************************ ***** Type declarations REAL LDOM, LDOMSS REAL LPOM, LPOMSS REAL LDOMD, LPOMD REAL LRDOMD, LRPOMD REAL LPOMAP, LPOMNS REAL LDOMAP REAL NH4, NH4D, NH4SS REAL NO3, NO3D, NO3SS, NO3SED REAL NH4AR, NH4AG, NH4AP, NH4DOM, NH4POM, NH4OM, NH4SD, . NH4SR REAL NO3AG REAL KF, KFS REAL NVIOL INTEGER SNP, TSR, PRF, VPL, CPL, FLX, SPR, . WSF, XTSR, OPT !jb03/98 CHARACTER CONFN*72 ***** PARAMETER declarations ********************* ** This include file was created by Jim Bowen (jdbowen@uncc.edu) 10/98 * * Application: Neuse Estuary, System N9 * * 66 segments, 20 layers, 2 branches, 12 tributaries * ********************** PARAMETER (CONFN='w2_con.npt') PARAMETER (IMP=66,KMP=20,IMC=66,KMC=20,NAP=3,NBP=2,NCP=22,NDP=14, . NFP=61,NHP=17,NPP=100,NRP=1,NSP=1,NTP=12,NWP=1) ***** DIMENSION declarations DIMENSION CSSK(KMC,IMC,NCP), C1(KMC,IMC,NCP), C2(KMC,IMC,NCP), . CD(KMC,IMC,NDP), KF(KMC,IMC,NFP), KFS(KMC,IMC,NFP), . AF(KMC,IMC,NAP,5), HYD(KMP,IMP,NHP), CPH(KMC,IMC,4) ***** I/O units DIMENSION OPT(NRP,9) DIMENSION SNP(NRP), TSR(NRP), PRF(NRP), . VPL(NRP), CPL(NRP), SPR(NRP), . FLX(NRP), WSF(NRP), XTSR(NRP) !jb03/98 ***** Hydrodynamic variables DIMENSION U(KMP,IMP), W(KMP,IMP), T2(KMP,IMP), . AZ(KMP,IMP), RHO(KMP,IMP), ST(KMP,IMP), . SB(KMP,IMP), NVIOL(KMP,IMP), VSH(KMP,IMP), . ADMX(KMP,IMP), DM(KMP,IMP), ADMZ(KMP,IMP), . HDG(KMP,IMP), HPG(KMP,IMP), RI(KMP,IMP), !JB06/98 . WSH2(KMP,IMP), BVF2(kmp,imp) !JB06/98 ***** State variables DIMENSION TDS(KMC,IMC), COL(KMC,IMC), SS(KMC,IMC), . NH4(KMC,IMC), NO3(KMC,IMC), PO4(KMC,IMC), . FE(KMC,IMC), DSI(KMC,IMC), PSI(KMC,IMC), . LDOM(KMC,IMC), RDOM(KMC,IMC), LPOM(KMC,IMC), . RPOM(KMC,IMC), CBOD(KMC,IMC), DO(KMC,IMC), . TIC(KMC,IMC), ALK(KMC,IMC), ALG(KMC,IMC,NAP) ***** State variable sources/sinks DIMENSION COLSS(KMC,IMC), AGESS(KMC,IMC), SSSS(KMC,IMC), . NH4SS(KMC,IMC), NO3SS(KMC,IMC), PO4SS(KMC,IMC), . FESS(KMC,IMC), DSISS(KMC,IMC), PSISS(KMC,IMC), . LDOMSS(KMC,IMC), RDOMSS(KMC,IMC), LPOMSS(KMC,IMC), . RPOMSS(KMC,IMC), CBODSS(KMC,IMC), DOSS(KMC,IMC), . TICSS(KMC,IMC), ASS(KMC,IMC,NAP) ***** Derived constituents DIMENSION PH(KMC,IMC), CO2(KMC,IMC), HCO3(KMC,IMC), . CO3(KMC,IMC) DIMENSION TN(KMC,IMC), TP(KMC,IMC) DIMENSION DON(KMC,IMC), DOP(KMC,IMC), DOC(KMC,IMC) DIMENSION PON(KMC,IMC), POP(KMC,IMC), POC(KMC,IMC) DIMENSION TON(KMC,IMC), TOP(KMC,IMC), TOC(KMC,IMC) DIMENSION APR(KMC,IMC), CHLA(KMC,IMC), ATOT(KMC,IMC) ***** Inorganic suspended solids fluxes DIMENSION SSSI(KMC,IMC), SSSO(KMC,IMC) ***** Phosphorus kinetic fluxes DIMENSION PO4AR(KMC,IMC), PO4AG(KMC,IMC), PO4AP(KMC,IMC), . PO4POM(KMC,IMC), PO4DOM(KMC,IMC), PO4OM(KMC,IMC), . PO4SD(KMC,IMC), PO4SR(KMC,IMC), PO4NS(KMC,IMC) ***** Ammonia kinetic fluxes DIMENSION NH4AR(KMC,IMC), NH4AG(KMC,IMC), NH4AP(KMC,IMC), . NH4POM(KMC,IMC), NH4DOM(KMC,IMC), NH4OM(KMC,IMC), . NH4SD(KMC,IMC), NH4SR(KMC,IMC), NH4D(KMC,IMC) ***** Nitrate/nitrite fluxes DIMENSION NO3AG(KMC,IMC), NO3D(KMC,IMC), NO3SED(KMC,IMC) ***** Silica fluxes DIMENSION DSIAG(KMC,IMC), DSID(KMC,IMC), DSISD(KMC,IMC), . DSISR(KMC,IMC), DSIS(KMC,IMC) DIMENSION PSIAM(KMC,IMC), PSID(KMC,IMC), PSINS(KMC,IMC) ***** Iron fluxes DIMENSION FENS(KMC,IMC), FESR(KMC,IMC) ***** Organic matter fluxes DIMENSION LDOMAP(KMC,IMC), LDOMD(KMC,IMC), LRDOMD(KMC,IMC), . RDOMD(KMC,IMC) DIMENSION LPOMAP(KMC,IMC), LPOMD(KMC,IMC), LRPOMD(KMC,IMC), . RPOMD(KMC,IMC), LPOMNS(KMC,IMC), RPOMNS(KMC,IMC) DIMENSION CBODD(KMC,IMC) ***** Algal fluxes DIMENSION AGR(KMC,IMC,NAP), ARR(KMC,IMC,NAP), AER(KMC,IMC,NAP), . AMR(KMC,IMC,NAP), ASR(KMC,IMC,NAP) ***** Dissolved oxygen fluxes DIMENSION DOAP(KMC,IMC), DOAR(KMC,IMC), DODOM(KMC,IMC), . DOPOM(KMC,IMC), DOOM(KMC,IMC), DONIT(KMC,IMC), . DOSED(KMC,IMC), DOSOD(KMC,IMC), DOBOD(KMC,IMC), . DOAE(KMC,IMC) ***** Total inorganic carbon DIMENSION TICAP(KMC,IMC) ***** Sediment fluxes DIMENSION SEDD(KMC,IMC), SODD(KMC,IMC), SEDAS(KMC,IMC), . SEDOMS(KMC,IMC), SEDNS(KMC,IMC), SEDB(KMC,IMC) !JB10/98 ***** State variable EQUIVALENCEs EQUIVALENCE (TDS(1,1), C2(1,1,3)), . (COL(1,1), C2(1,1,4)), . (SS(1,1), C2(1,1,5)), . (PO4(1,1), C2(1,1,6)), . (NH4(1,1), C2(1,1,7)), . (NO3(1,1), C2(1,1,8)), . (DSI(1,1), C2(1,1,9)), . (PSI(1,1), C2(1,1,10)), . (FE(1,1), C2(1,1,11)), . (LDOM(1,1), C2(1,1,12)), . (RDOM(1,1), C2(1,1,13)), . (LPOM(1,1), C2(1,1,14)), . (RPOM(1,1), C2(1,1,15)), . (CBOD(1,1), C2(1,1,16)), . (ALG(1,1,1), C2(1,1,17)), . (ALG(1,1,2), C2(1,1,18)), . (ALG(1,1,3), C2(1,1,19)), . (DO(1,1), C2(1,1,20)), . (TIC(1,1), C2(1,1,21)), . (ALK(1,1), C2(1,1,22)) ***** I/O unit EQUIVALENCEs EQUIVALENCE (SNP(1), OPT(1,1)), . (TSR(1), OPT(1,2)), . (PRF(1), OPT(1,3)), . (VPL(1), OPT(1,4)), . (CPL(1), OPT(1,5)), . (SPR(1), OPT(1,6)), . (FLX(1), OPT(1,7)), . (WSF(1), OPT(1,8)), . (XTSR(1), OPT(1,9)) ***** State variable sources/sinks EQUIVALENCEs EQUIVALENCE (AGESS(1,1), CSSK(1,1,2)), . (COLSS(1,1), CSSK(1,1,4)), . (SSSS(1,1), CSSK(1,1,5)), . (PO4SS(1,1), CSSK(1,1,6)), . (NH4SS(1,1), CSSK(1,1,7)), . (NO3SS(1,1), CSSK(1,1,8)), . (DSISS(1,1), CSSK(1,1,9)), . (PSISS(1,1), CSSK(1,1,10)), . (FESS(1,1), CSSK(1,1,11)), . (LDOMSS(1,1), CSSK(1,1,12)), . (RDOMSS(1,1), CSSK(1,1,13)), . (LPOMSS(1,1), CSSK(1,1,14)), . (RPOMSS(1,1), CSSK(1,1,15)), . (CBODSS(1,1), CSSK(1,1,16)), . (ASS(1,1,1), CSSK(1,1,17)), . (ASS(1,1,2), CSSK(1,1,18)), . (ASS(1,1,3), CSSK(1,1,19)), . (DOSS(1,1), CSSK(1,1,20)), . (TICSS(1,1), CSSK(1,1,21)) ***** pH related EQUIVALENCEs EQUIVALENCE (PH(1,1), CPH(1,1,1)), . (CO2(1,1), CPH(1,1,2)), . (HCO3(1,1), CPH(1,1,3)), . (CO3(1,1), CPH(1,1,4)) ***** Derived variable EQUIVALENCEs EQUIVALENCE (DOC(1,1), CD(1,1,1)), . (POC(1,1), CD(1,1,2)), . (TOC(1,1), CD(1,1,3)), . (DON(1,1), CD(1,1,4)), . (PON(1,1), CD(1,1,5)), . (TON(1,1), CD(1,1,6)), . (TN(1,1), CD(1,1,7)), . (DOP(1,1), CD(1,1,8)), . (POP(1,1), CD(1,1,9)), . (TOP(1,1), CD(1,1,10)), . (TP(1,1), CD(1,1,11)), . (APR(1,1), CD(1,1,12)), . (CHLA(1,1), CD(1,1,13)), . (ATOT(1,1), CD(1,1,14)) ***** Kinetic flux EQUIVALENCEs EQUIVALENCE (SSSI(1,1), KF(1,1,1)), . (SSSO(1,1), KF(1,1,2)) EQUIVALENCE (PO4AR(1,1), KF(1,1,3)), . (PO4AG(1,1), KF(1,1,4)), . (PO4AP(1,1), KF(1,1,5)), . (PO4POM(1,1), KF(1,1,6)), . (PO4DOM(1,1), KF(1,1,7)), . (PO4OM(1,1), KF(1,1,8)), . (PO4SD(1,1), KF(1,1,9)), . (PO4SR(1,1), KF(1,1,10)), . (PO4NS(1,1), KF(1,1,11)) EQUIVALENCE (NH4D(1,1), KF(1,1,12)), . (NH4AR(1,1), KF(1,1,13)), . (NH4AG(1,1), KF(1,1,14)), . (NH4AP(1,1), KF(1,1,15)), . (NH4POM(1,1), KF(1,1,16)), . (NH4DOM(1,1), KF(1,1,17)), . (NH4OM(1,1), KF(1,1,18)), . (NH4SD(1,1), KF(1,1,19)), . (NH4SR(1,1), KF(1,1,20)) EQUIVALENCE (NO3D(1,1), KF(1,1,21)), . (NO3AG(1,1), KF(1,1,22)), . (NO3SED(1,1), KF(1,1,23)) EQUIVALENCE (DSIAG(1,1), KF(1,1,24)), . (DSID(1,1), KF(1,1,25)), . (DSISD(1,1), KF(1,1,26)), . (DSISR(1,1), KF(1,1,27)), . (DSIS(1,1), KF(1,1,28)), . (PSIAM(1,1), KF(1,1,29)), . (PSINS(1,1), KF(1,1,30)), . (PSID(1,1), KF(1,1,31)) EQUIVALENCE (FENS(1,1), KF(1,1,32)), . (FESR(1,1), KF(1,1,33)) EQUIVALENCE (LDOMD(1,1), KF(1,1,34)), . (LRDOMD(1,1), KF(1,1,35)), . (RDOMD(1,1), KF(1,1,36)), . (LDOMAP(1,1), KF(1,1,37)), . (LPOMD(1,1), KF(1,1,38)), . (LRPOMD(1,1), KF(1,1,39)), . (RPOMD(1,1), KF(1,1,40)), . (LPOMAP(1,1), KF(1,1,41)), . (LPOMNS(1,1), KF(1,1,42)), . (RPOMNS(1,1), KF(1,1,43)) EQUIVALENCE (CBODD(1,1), KF(1,1,44)) EQUIVALENCE (DOAP(1,1), KF(1,1,45)), . (DOAR(1,1), KF(1,1,46)), . (DOPOM(1,1), KF(1,1,47)), . (DODOM(1,1), KF(1,1,48)), . (DOOM(1,1), KF(1,1,49)), . (DONIT(1,1), KF(1,1,50)), . (DOBOD(1,1), KF(1,1,51)), . (DOAE(1,1), KF(1,1,52)), . (DOSED(1,1), KF(1,1,53)), . (DOSOD(1,1), KF(1,1,54)) EQUIVALENCE (TICAP(1,1), KF(1,1,55)) EQUIVALENCE (SEDD(1,1), KF(1,1,56)), . (SEDAS(1,1), KF(1,1,57)), . (SEDOMS(1,1), KF(1,1,58)), . (SEDNS(1,1), KF(1,1,59)), . (SODD(1,1), KF(1,1,60)), . (SEDB(1,1), KF(1,1,61)) !JB10/98 ***** Algal rate EQUIVALENCEs EQUIVALENCE (AGR(1,1,1), AF(1,1,1,1)), . (ARR(1,1,1), AF(1,1,1,2)), . (AER(1,1,1), AF(1,1,1,3)), . (AMR(1,1,1), AF(1,1,1,4)), . (ASR(1,1,1), AF(1,1,1,5)) ***** Hydrodynamic EQUIVALENCEs EQUIVALENCE (NVIOL(1,1), HYD(1,1,1)), . (U(1,1), HYD(1,1,2)), . (W(1,1), HYD(1,1,3)), . (T2(1,1), HYD(1,1,4)), . (RHO(1,1), HYD(1,1,5)), . (AZ(1,1), HYD(1,1,6)), . (VSH(1,1), HYD(1,1,7)), . (ST(1,1), HYD(1,1,8)), . (SB(1,1), HYD(1,1,9)), . (ADMX(1,1), HYD(1,1,10)), . (DM(1,1), HYD(1,1,11)), . (HDG(1,1), HYD(1,1,12)), . (ADMZ(1,1), HYD(1,1,13)), . (HPG(1,1), HYD(1,1,14)), . (RI(1,1), HYD(1,1,15)), . (WSH2(1,1), HYD(1,1,16)), . (BVF2(1,1), HYD(1,1,17)) ***** COMMON declarations COMMON /GCNSTC/ HYD, C1, C2, CD, CSSK, KF, KFS, AF, OPT