C VERSION(8/6/96) C C >>>> PARAMETERS C cc PARAMETER (IM=147,JM=452,KB=11,KSL=10,KSMAX=10,LAYMAX=11) ! NYHOPS GRID PARAMETER (IM=200,JM=200,KB=3,KSL=3,KSMAX=3,LAYMAX=4) cc PARAMETER (IM=115,JM=76,KB=11,KSL=3,KSMAX=3,LAYMAX=4) ! JEM 2012 Grid cc PARAMETER (IM=115,JM=76,KB=11,KSL=10,KSMAX=10,LAYMAX=11) ! JEM 2016 Grid cc PARAMETER (IM=1069,JM=550,KB=11,KSL=2,KSMAX=2,LAYMAX=3) ! J-Bay 30m grid cc PARAMETER (IM=190,JM=5,KB=6,KSL=4,KSMAX=4,LAYMAX=5) cc PARAMETER (IM=210,JM=5,KB=6,KSL=4,KSMAX=4,LAYMAX=5) cc PARAMETER (IM=685,JM=219,KB=6,KSL=4,KSMAX=4,LAYMAX=5) C C >>>> OTHER PARAMETERS C PARAMETER (IMM1=IM-1,JMM1=JM-1,KBM1=KB-1) PARAMETER (IMM2=IM-2,JMM2=JM-2,KBM2=KB-2) PARAMETER (MAXWET=IM*JM) C C MAX NUMBER OF LAYERS IN CHEM BED MODEL C PARAMETER (NBEDMAX=100) C C REAL ISPI,ISP2I,KM,KH,KQ,L,NU INTEGER EPTS,VPTS,FPTS INTEGER IHISTM,ITRACM,EPTSM,VPTSM,FPTSM,EBCM,QBCM,DBCM,PSCM INTEGER EPCM PARAMETER (IHISTM=19600,ITRACM=19600) PARAMETER (EPTSM=500,VPTSM=500,FPTSM=500) PARAMETER (EBCM=2*(IM+JM),QBCM=200,DBCM=650,PSCM=100,NUMDAM=150) PARAMETER (EPCM=500) ! NG04232008 ELEVATION POTENTIAL CELLS C C RMarsooli_August2015, Mellor's wave model C wave parameter mmm = Number of angle segment + 2 parameter (mmm=24+2) C C********************************************************************* C C FOR PARTICLE TRACKING C C NSRCM = MAXIMUM NO. OF PARTICLE SOURCES C NPARTM = MAXIMUM NO. OF PARTICLES/RELEASE C NCONVM = MAXIMUM NO. OF PARTICLE RELEASES BEFORE CONVERSION TO C CONCENTRATION C C NOTE: MAX. NO. OF PARTICLES = NPARTM * NCONVM C PARAMETER (NSRCM=5,NPARTM=100,NCONVM=1000) C C********************************************************************* C LOGICAL CONSPLT,CONSTSR,CONSTRN,CDFOUT,CREADWBC,CREADPET LOGICAL CREADMET !RMarsooli_May2015 C CHARACTER*1 VSX,VSY,PTU,PTV,PTW,PTAM,PTS,PTT,PRHO, . PTQ2,PTL,PTKM,PTKH . ,PTHF,PTCD,PTUF,PTVF ! NG 03/04/2008 included . ,PTWAVE ! NG 01/05/2010 included CHARACTER*3 DEV CHARACTER*4 DIRFLX, SEDTYPE, OPTPSTR CHARACTER*7 TRACER C CHARACTER*5 PDEPFORM CHARACTER*7 SEDTRAN,CHEMTRAN,PARTICLE,SHERODE CHARACTER*7 VARIP0,VARIA0N,SEDLOAD,VARIBULK CHARACTER*8 WAVEDYN,HYDTYPE C CHARACTER*10 TOR,ADVECT,HORZMIX,VERTMIX,SCHEME CHARACTER OPTEXTC*4 CHARACTER*20 OPTEBC,OPTTSI,OPTMBC CHARACTER OPTEPC*8 ! NG05062008 NEW OPTION FOR ELEVATION POTENTIAL CELLS CHARACTER*7 WETLAND !RMarsooli_Jan2015 CHARACTER*10 MADSCHEME !RMarsooli_Feb2015 CHARACTER*8 VEGCDM,VEGCDWAVE !RMarsooli_Feb2015 C REAL*8 TSED(LAYMAX,IM,JM),EBTOT(LAYMAX,IM,JM),EBMAX(LAYMAX,IM,JM) REAL*8 TSED0IN(LAYMAX),EBCUR(LAYMAX,IM,JM),PSED1(LAYMAX,IM,JM), + P0(2),THICK(IM,JM),PSED2(LAYMAX,IM,JM), + PTOT(IM,JM),TSEDT(IM,JM),TSET0T,E(2,IM,JM),DD(2,IM,JM), + BEDTHI,BEDTH(2,IM,JM),TSED0(LAYMAX,IM,JM) REAL*8 FPBED(KSMAX,IM,JM),ACTLAY(2,IM,JM),CARMOR(IM,JM), + FALAY(KSMAX,IM,JM) C REAL*8 ARCTHIK(IM,JM),THSAVE(VPTSM) C REAL NIKURH,KHZSAVE,KMZSAVE,AMZSAVE C REAL*8 SEDMASS(NBEDMAX,IM,JM),CHEMMASS(NBEDMAX,IM,JM), + SEDDEP(IM,JM) REAL*8 QEPCNUDG(EPCM),QEPCNUDF(EPCM) ! NG04232008 C REAL*8 BVEG(IM,JM),HVEG(IM,JM),HVEGA(IM,JM), . CDVEG(IM,JM),CDVEGW(IM,JM),STIFF(IM,JM),NVEG(IM,JM) !RMarsooli_Jan2015 REAL*8 DVEG3D(IM,JM,KB),DVEG2D(IM,JM),IVEG3D(IM,JM,KB), . IVEG2D(IM,JM) !RMarsooli_Jan2015 REAL*8 DWALLU3D(IM,JM,KB),DWALLV3D(IM,JM,KB),SVF(IM,JM), . DWALLU2D(IM,JM),DWALLV2D(IM,JM) !RMarsooli_Jan2015 REAL*8 Z0B(IM,JM) REAL*8 CH_SEN(IM,JM),CE_LAT(IM,JM) !RMarsooli_May2015 REAL*8 CDwindN(IM,JM),CH_SENN(IM,JM),CE_LATN(IM,JM), . ustarwind(IM,JM),uNwind(IM,JM) REAL*8 SHEATFLUX(IM,JM),LHEATFLUX(IM,JM),RHEATFLUX(IM,JM), . TMOMEFLUX(IM,JM),RMOMEFLUX(IM,JM) !RMarsooli_July2015 INTEGER IDXVEG(MAXWET),IDYVEG(MAXWET) !RMarsooli_April2015 REAL*8 COSANG(IM,JM),SINANG(IM,JM) !RMarsooli REAL*8 DELX(IM,JM),DELY(IM,JM) !RMarsooli REAL*8 SMALLVAL !RMarsooli C----------------------------------------------------------------------- C RMarsooli_August2015 C Common variables associated with the Mellor's wave model integer & m,nwavemel,isplitw real & dthinv,dtw,dtw2,dtww,ent,enwnd,thtav,thtpws,Hs, & aren,areninv,x,y, & css,snn,thetaa,thetam,csm,snm, & enw,ens,ene,enn, & sigthw,sigths,sigthe,sigthn, & sigthf,sigth,sigthb,sigthav, & cth,cg,kth,kthD,kthav,kthavD, & kpD,cgx,cgy, & age,ageinv,tht_wnd, & sigp,sigpws,ustw,cdd, & cpp,kpp, cthth, & enb,en,enf, & qdiss,Hm,bdis1,bdis2, & bdisveg,vegQ, !RMarsooli & sdiss,swin,fspr, & swintot,sdistot,bdistot,ustp, & tps2d,ud,vd, C variables used by circulation model & tpx0,tpy0,ttx0,tty0, & tpx,tpy,tpzdist,ttx,tty, & ust,vst, & ucb,ucf,vcb,vcf, & Sxx,Sxy,Syy, & F1,F2,F3,F4 & ,WS2D !RMarsooli(Magnitude of wind speed) INTEGER NUMWMBC,IWM,JWM,IWMC,JWMC !RMarsooli c------------------------------------------------------------------------- DATA RHO0/1000./ CNG COMMON/TESTNG/WU(IM,JM),WV(IM,JM),WU1(IM,JM),WV1(IM,JM), CNG . WU2(IM,JM),WV2(IM,JM),T1WV,T2WV COMMON/BLKCON/ISPI,ISP2I,SMOTH,THETA,DAYI, . DTE,DTI,HPRNU,VPRNU,UMOL,GRAV,RAMP,NU,HORCON,ALPHA,TLAG,TIME, . THOUR,ISTART,IEND,INTX,IEXT,IPRINT,IPRTSTART,IAVGE,ISKILL, . IURUN,IUGRD,IUTAS,IUUAV,IURRS,IUPRT,IUPLT,IUTSR, . IUTRN,IUWRS,IUT90,IUT91,IUT92,IUT93,IUT94,IUT95,IUT96, . IUT103, . IUT191,IUT192,ISPLIT,IUT194, . CALCSWOBS,CALCCLD,WETMIN,WETEPS, . TOR,ADVECT,SCHEME,HORZMIX,VERTMIX,OPTEBC,OPTTSI,OPTMBC . ,OPTEXTC . ,HMINWAVE !RMarsooli_AUGUST_2015 C COMMON /BLKCON2/ TSOPTION,SEDTRAN, + PDEPFORM,CHEMTRAN,WAVEDYN,HYDTYPE,SEDTYPE, + VARIBULK, OPTPSTR C COMMON /BLKCON4/IUT401,IUT402,IUT403,IUT404,IUT405,IUT406 COMMON /BLKCON5/IUT501,IUT502,IUT503,IUT504,IUT505,IUT506 COMMON /BLKCON6/IUT601,IUT602,IUT603,IUT604,IUT605,IUT606 COMMON /BLKCON7/IUT701,IUT702,IUT703,IUT704,IUT705 COMMON /BLKCON8/IUT801 C C--------------- PRINT CONTROL ARRAYS ---------------------------------- C COMMON/BLKPRT/IROW,JROW COMMON/BLKCHR/ . DEV,VSX,VSY,PTU,PTV,PTW,PTAM,PTS,PTT, . PRHO,PTQ2,PTL,PTKM,PTKH . ,PTHF,PTCD,PTUF,PTVF ! NG 03/04/2008 included . ,PTWAVE ! NG 01/05/2010 included C C---------------- 1-D ARRAYS ------------------------------------------- C COMMON/BLK1D/ . DZR(KB),Z(KB),ZZ(KB),DZ(KB),DZZ(KB), . DPTHSL(KSL),TSI(KSL),SSI(KSL) C C---------------- 2-D ARRAYS ------------------------------------------- C COMMON/BLK2D/ . H(IM,JM),H1(IM,JM),H2(IM,JM),D(IM,JM),DT(IM,JM),ANG(IM,JM), . ART(IM,JM),ARU(IM,JM),ARV(IM,JM),CBC(IM,JM),DATUM(IM,JM), . DUM(IM,JM),DVM(IM,JM),FSM(IM,JM),COR(IM,JM),CURV42D(IM,JM), . TXSURF(IM,JM),TYSURF(IM,JM),TXBOT(IM,JM),TYBOT(IM,JM), . WTSURF(IM,JM),WSSURF(IM,JM),TPS(IM,JM),AAM2D(IM,JM), . UAF(IM,JM),UA(IM,JM),UAB(IM,JM),VAF(IM,JM),VA(IM,JM), . VAB(IM,JM),ELF(IM,JM),EL(IM,JM),ELB(IM,JM),ATELB(IM,JM), . ETF(IM,JM),ET(IM,JM),ETB(IM,JM),FLUXUA(IM,JM),FLUXVA(IM,JM), . EGF(IM,JM),EGB(IM,JM),UTF(IM,JM),UTB(IM,JM), . VTF(IM,JM),VTB(IM,JM), . C11(IM,JM),C12(IM,JM),C21(IM,JM),C22(IM,JM), . XGRID(IM,JM),YGRID(IM,JM),XCORNER(IM+1,JM+1),YCORNER(IM+1,JM+1), !csvv080326 for 2d heat . VARBF(IM,JM),VARHF(IM,JM),VARUF(IM,JM),VARMN(IM,JM), . EXTC1(IM,JM) . ,DUMWAD(IM,JM),DVMWAD(IM,JM) ! NG05192010 . ,DHMWAD(IM,JM),FSMRAD(IM,JM) !RMarsooli_June2015 C COMMON/ICE/ ! NG 2011 . VARWIF(IM,JM),TXICE(IM,JM),TYICE(IM,JM), . IUT196, ICEFRICTION, . T1I, T2I, . THICKICE(IM,JM),DTHICKICE(IM,JM,2), . AREAICE(IM,JM),DAREAICE(IM,JM,2) COMMON/VEG/ !RMarsooli_Jan2015 . BVEG,HVEG,HVEGA,CDVEG,CDVEGW,STIFF,CDINER,CDQ2,CDQ2L, . DVEG3D,DVEG2D,VEGCDM,NUMVEG,IDXVEG,IDYVEG,VEGCDWAVE, . IUVEG,NVEG,IVEG3D,IVEG2D,WETLAND COMMON/WALL/ . DWALLU3D,DWALLV3D,DWALLU2D,DWALLV2D,CDWALL !RMarsooli_Jan2015 COMMON/URBAN/ . SVF !RMarsooli_April2015 COMMON/HORVISBACK/ . HORBACK !RMarsooli_Jan2015 COMMON/COAREDATA/ . ZU_COARE,ZT_COARE,ZQ_COARE,ZS_COARE, . ZI_COARE,COAREVERSION,CH_SEN,CE_LAT, . CDwindN,CH_SENN,CE_LATN,ustarwind,uNwind, . JWARM_COARE,JCOOL_COARE,JWAVE_COARE,NCOARE,CREADMET !RMarsooli_April2015 COMMON/SURFLUXES/ . SHEATFLUX,LHEATFLUX,RHEATFLUX,TMOMEFLUX,RMOMEFLUX !RMarsooli_July2015 COMMON/MESHVAR/COSANG,SINANG,SMALLVAL,DELX,DELY !RMarsooli C C---------------- 3-D ARRAYS ------------------------------------------- C COMMON/BLK3D/ . A(IM,JM,KB),C(IM,JM,KB),VH(IM,JM,KB),VHP(IM,JM,KB), . PROD(IM,JM,KB),DTEF(IM,JM,KB),XMFL3D(IM,JM,KB), . YMFL3D(IM,JM,KB),UF(IM,JM,KB),VF(IM,JM,KB),KM(IM,JM,KB), . KH(IM,JM,KB),KQ(IM,JM,KB),L(IM,JM,KB),Q2(IM,JM,KB), . Q2B(IM,JM,KB),AAM(IM,JM,KB),Q2L(IM,JM,KB),Q2LB(IM,JM,KB), . U(IM,JM,KB),UB(IM,JM,KB),W(IM,JM,KB),WR(IM,JM,KB), . V(IM,JM,KB),VB(IM,JM,KB),T(IM,JM,KB),TB(IM,JM,KB), . S(IM,JM,KB),SB(IM,JM,KB),RHO(IM,JM,KB),RMEAN(IM,JM,KB), . TMEAN(IM,JM,KB),SMEAN(IM,JM,KB), . AAMAX(IM,JM,KB),AAMAY(IM,JM,KB) C C----------- ARCHIVING AND SKILL ASSESSMENT ARRAYS --------------------- C COMMON/CONS/ . CONSPLT,CONSTSR,CONSTRN,CDFOUT COMMON/ARCH/DEI, . ARCU(IM,JM,KB),ARCV(IM,JM,KB),ARCUX(IM,JM,KB), . ARCVX(IM,JM,KB),ARCS(IM,JM,KB),ARCQ2(IM,JM,KB), . ARCT(IM,JM,KB),ARCW(IM,JM,KB),ARCKH(IM,JM,KB), . ARCET(IM,JM),ARCAM(IM,JM,KB),ARCKM(IM,JM,KB), . IHIST(IHISTM,2),JHM, . ARCPATM(IM,JM),ARCQPREC(IM,JM),ARCQEVAP(IM,JM), . ARCCLD(IM,JM),ARCAIRT(IM,JM),ARCRHUM(IM,JM), ! NG04092008 . ARCSWOBS(IM,JM),ARCWSX2D(IM,JM),ARCWSY2D(IM,JM), ! NG04092008 . ARCWVPD(IM,JM),ARCWVHT(IM,JM),ARCWVDR(IM,JM), ! NG04142008 . ARCWVSP(IM,JM),ARCWVCP(IM,JM),ARCWVEN(IM,JM), ! RM05312016 . ARCWVST(IM,JM,mmm),ARCWVCT(IM,JM,mmm),ARCWVET(IM,JM,mmm) ! RM05312016 COMMON/TRANS/FLTWT, . ITRAC(ITRACM,2),NPLPF,IWET COMMON/SKILL/ . ESAVE(EPTSM),DZSAVE(VPTSM),CCFLUX(FPTSM,KB), . UZSAVE(VPTSM,KB),VZSAVE(VPTSM,KB),SZSAVE(VPTSM,KB), . TZSAVE(VPTSM,KB),KMZSAVE(VPTSM,KB),KHZSAVE(VPTSM,KB), . WZSAVE(VPTSM,KB),AMZSAVE(VPTSM,KB), . SKILLI,AREA,EMI,APEI, . EPTS,VPTS,FPTS,TXSAVE(EPTSM),TYSAVE(EPTSM), . QPRECSAVE(EPTSM),QEVAPSAVE(EPTSM),PATMSAVE(EPTSM), . CLDSAVE(EPTSM),AIRTSAVE(EPTSM),RHUMSAVE(EPTSM), ! NG04092008 . SWOBSSAVE(EPTSM),WSX2DSAVE(EPTSM),WSY2DSAVE(EPTSM), ! NG04092008 . WVPDSAVE(EPTSM),WVHTSAVE(EPTSM),WVDRSAVE(EPTSM), ! NG04142008 . WVSPSAVE(EPTSM),WVCPSAVE(EPTSM),WVENSAVE(EPTSM), ! RM05312016 . WVSTSAVE(EPTSM,mmm),WVCTSAVE(EPTSM,mmm),WVETSAVE(EPTSM,mmm), ! RM05312016 . INXIE(EPTSM),INXJE(EPTSM),INXIV(VPTSM),INXJV(VPTSM), . ISFLX(FPTSM),JSFLX(FPTSM),NFLXE(FPTSM), . CS1FLUX(FPTSM,KB),CS2FLUX(FPTSM,KB) COMMON/SKICHR/ DIRFLX(FPTSM) C C---------------------- BOUNDARY VALUE ARRAYS -------------------------- C COMMON/EBC/ . AMP(EBCM,6),PHASE(EBCM,6),PERIOD(6),EMEAN(EBCM),VTRAN(EBCM), . TBDRYSL(EBCM,KSL),SBDRYSL(EBCM,KSL), . NUMEBC,IETA(EBCM),JETA(EBCM),ICON(EBCM),JCON(EBCM), . NUMECR,ICHARD(EBCM),JCHARD(EBCM),NCHARD(EBCM), ! CNG 2012 . NUMEBCTR,ITRED(EBCM),JTRED(EBCM),ITREC(EBCM),JTREC(EBCM), . NUMEBCSE,ISEED(EBCM),JSEED(EBCM),ISEEC(EBCM),JSEEC(EBCM), . NUMEBCCH,ICHED(EBCM),JCHED(EBCM),ICHEC(EBCM),JCHEC(EBCM), . ITAS(EBCM),JTAS(EBCM) COMMON/QBC/ VQDIST(QBCM,KBM1), . NUMQBC,IQD(QBCM),JQD(QBCM),IQC(QBCM),JQC(QBCM), . NUMQBCTR,ITRQD(QBCM),JTRQD(QBCM),ITRQC(QBCM),JTRQC(QBCM), . NUMQBCSE,ISEQD(QBCM),JSEQD(QBCM),ISEQC(QBCM),JSEQC(QBCM), . NUMQBCCH,ICHQD(QBCM),JCHQD(QBCM),ICHQC(QBCM),JCHQC(QBCM) COMMON/DBC/ VDDIST(DBCM,KBM1),VDDISTSE(DBCM,KBM1), . VDDISTCH(DBCM,KBM1), . NUMDBC,NUMDBC1,NUMDBC2,IDD(DBCM),JDD(DBCM), . NUMDBCSE,NUMDBCCH, . IDDSE(DBCM),JDDSE(DBCM),IDDCH(DBCM),JDDCH(DBCM) COMMON /PSC/NUMPSTR,IPSTR(PSCM),JPSTR(PSCM),KPSTR(PSCM) C C---------------------- FORCING ARRAYS --------------------------------- C COMMON/FORCING1/ . T1E,T2E,EBDRY(EBCM),DEBDRY(EBCM,2), . T1Q,T2Q,QDIS (QBCM),DQDIS (QBCM,2),TDIS (QBCM),DTDIS (QBCM,2), . SDIS (QBCM),DSDIS (QBCM,2), . T1D,T2D,QDIFF(DBCM),DQDIFF(DBCM,2),TDIFF(DBCM),DTDIFF(DBCM,2), . SDIFF(DBCM),DSDIFF(DBCM,2),T1D2,T2D2,T1D2CON,T2D2CON, . T1M,T2M,QPREC,DQPREC(2),QEVAP,DQEVAP(2),HFLUX,DHFLUX(2), . TX,DTX(2),TY,DTY(2),WDSMAX,PATMOPT,DHFLX2D(IM,JM,2), . T1W,T2W,TX2D(IM,JM),DTX2D(IM,JM,2),TY2D(IM,JM),DTY2D(IM,JM,2), . T1TS,T2TS,TBDRY(EBCM,KBM1),DTBDRY(EBCM,KBM1,2), . SBDRY(EBCM,KBM1),DSBDRY(EBCM,KBM1,2), . TDIF(DBCM,KBM1),SDIF(DBCM,KBM1), . PATM(IM,JM),DPATM(IM,JM,2), . T1EX,T2EX,DEXTC(IM,JM,2) COMMON/FORCING2/ + T1WAVE,T2WAVE,DWVPD(IM,JM,2),DWVHT(IM,JM,2),DWVDR(IM,JM,2), + T1SEDM,T2SEDM,T1CON,T2CON,T1CHMM,T2CHMM, + T1SEDS,T2SEDS,T1CHMS,T2CHMS, + T1SEDOBM,T2SEDOBM,T1CONOB,T2CONOB,T1CHMOBM,T2CHMOBM, + T1SEDOBS,T2SEDOBS,T1CHMOBS,T2CHMOBS, + T1DCON,T2DCON,T1SEDDM,T2SEDDM,T1SEDDS,T2SEDDS, + T1CHMDM,T2CHMDM,T1CHMDS,T2CHMDS CNG04232008 INCLUDE ELEVATION POTENTIAL CELLS COMMON/FORCING3/ + IUT65,NUMEPC,IEPC(EPCM),JEPC(EPCM),EPCMEAN(EPCM), + T1EPC,T2EPC,TNUDGEPC(EPCM),EPCOR(EPCM),DEPCOR(EPCM,2), + QEPCNUDG,TQEPCIN(EPCM,KBM1),SQEPCIN(EPCM,KBM1), + TQEPC(EPCM,KBM1),SQEPC(EPCM,KBM1),OPTEPC,QEPCNUDF C C----------------------------------------------------------------------- C COMMON /BAR/ DULEV,ICLOSE(100), IOPEN(100), ICL, IOP C C-----------------MODEL CALCULATED HEAT FLUX, jeff JI, 8/31/94 C COMMON /MHFLUX/ J1JDAY,J2JDAY,J1YR,CLDX !(6/21/02) * ,IYR,IMO,IDA,IHR,IHOUR,ALAT,ALON,CLDFRC,SDAY,EDAY,TR,EXTC * ,DAIRTM(2),DRELHU(2),DBAROP(2),DSWOBS(2),DWSX(2),DWSY(2) * ,CLOUD(2), SWRAD(IM,JM), RADP(IM,JM,KB) * ,EXTCOEF(2), ZEROS(IM,JM), WNDSH, REFL * ,WNDSH2D(IM,JM),WNDSH_VEG(IM,JM) !RMarsooli2015 * ,UTCSHIFT ! NG04232008 * ,CREADPET ! NG09022014 For 2D Precipitation and Evaporation synop_pet COMMON /MHFLUX2D/ T1P,T2P,QPREC2D(IM,JM),DQPREC2D(IM,JM,2) * ,QEVAP2D(IM,JM),DQEVAP2D(IM,JM,2) * ,AIRTM2D(IM,JM),DAIRTM2D(IM,JM,2) * ,RELHU2D(IM,JM),DRELHU2D(IM,JM,2) * ,SWOBS2D(IM,JM),DSWOBS2D(IM,JM,2) * ,CLDFRC2D(IM,JM),DCLD2D(IM,JM,2) * ,EXTC2D(IM,JM), DEXTC2D(IM,JM,2) * ,WSX2D(IM,JM),DWSX2D(IM,JM,2) * ,WSY2D(IM,JM),DWSY2D(IM,JM,2) CNG04092008 This should use PATM -> * ,BAROP2D(IM,JM),DBAROP2D(IM,JM,2) C COMMON /TRACER0/TRACER COMMON /TRACER1/CONC1(IM,JM,KB),CONC1B(IM,JM,KB), + WCSURF(IM,JM),CDIFF1(DBCM),DCDIFF1(DBCM,2),ARCC1(IM,JM,KB), + CBDRY1(EBCM,KBM1),CDIS1(QBCM),DCDIS1(QBCM,2),CMEAN1(IM,JM,KB), + IUT98,IUT99,CBDRYSL1(EBCM,KSL),DCBDRY1(EBCM,KBM1,2), + CSI1(KSL),C1ZSAVE(VPTSM,KB),CDIF1(DBCM,KBM1), & NUMDBCTR,NUMDBCTR1,NUMDBCTR2,ITRDD(DBCM),JTRDD(DBCM) !hli COMMON /TRACER2/T1PSTR,T2PSTR,DCPSTR(PSCM,2),CPSTR(PSCM) Cqa Added decay rate function of T,S,Io and Vs by Quamrul QA 11/12/01 COMMON /TRACER3/CONDRAT,THETAT,VSRATE,ASOL,ASAL, +DKRATE(IM,JM,KB),DKD(IM,JM,KB),DKI(IM,JM,KB),DKS(IM,JM,KB) C C----------------------------------------------------------------------- C RMarsooli_August2015 C Common variables associated with the Mellor's wave model common/blkwave/ & dthinv, & ent(im,jm),enwnd(im,jm),thtav(im,jm),thtpws(im,jm),Hs(im,jm), & aren(im,jm),areninv(im,jm),x(im),y(jm), & css(mmm),snn(mmm),thetaa(mmm),thetam(mmm),csm(mmm),snm(mmm), & enw(jm,mmm),ens(im,mmm),ene(jm,mmm),enn(im,mmm), & sigthw(jm,mmm),sigths(im,mmm),sigthe(jm,mmm),sigthn(im,mmm), & sigthf(im,jm,mmm),sigth(im,jm,mmm),sigthb(im,jm,mmm), & sigthav(im,jm), & cth(im,jm,mmm),cg(im,jm,mmm),kth(im,jm,mmm),kthD(im,jm,mmm), & kpD(im,jm),kthav(im,jm),kthavD(im,jm), & cgx(im,jm,mmm),cgy(im,jm,mmm), & age(im,jm),ageinv(im,jm),tht_wnd(im,jm), & sigp(im,jm),sigpws(im,jm),ustw(im,jm),cdd(im,jm), & cpp(im,jm),kpp(im,jm), cthth(im,jm,mmm), & enb(im,jm,mmm),en(im,jm,mmm),enf(im,jm,mmm), & qdiss(im,jm),Hm(im,jm),bdis1(im,jm,mmm),bdis2(im,jm,mmm), & bdisveg(im,jm,mmm),vegQ(im,jm), !RMarsooli & sdiss(im,jm,mmm),swin(im,jm,mmm),fspr(im,jm,mmm), & swintot(im,jm),sdistot(im,jm),bdistot(im,jm), & ustp(im,jm), & tps2d(im,jm),ud(im,jm),vd(im,jm), C variables used by circulation model & tpx0(im,jm),tpy0(im,jm),ttx0(im,jm),tty0(im,jm), & tpx(im,jm,kb),tpy(im,jm,kb),tpzdist(im,jm,kb), CRM & ttx(im,jm,kb),tty(im,jm,kb), & ust(im,jm,kb),vst(im,jm,kb), & ucb(im,jm,kb),ucf(im,jm,kb),vcb(im,jm,kb),vcf(im,jm,kb), & Sxx(im,jm,kb),Sxy(im,jm,kb),Syy(im,jm,kb), & F1(kb),F2(kb),F3(kb),F4(kb) c RMarsooli, other new variables & ,WS2D(im,jm) !RMarsooli & ,NUMWMBC & ,IWM(EBCM),JWM(EBCM),IWMC(EBCM),JWMC(EBCM) & ,WHINITIAL,WPINITIAL !RMarsooli C C---------------- SEDIMENT TRANSPORT ----------------------------------- C COMMON /SED1/A0IN,ALPHAS,WS2,TCRDEP,TCRDP2,RESEXP,N24CNT, + EXPM,IP0,IA0,XNSBED,IBED,UCR2,DENNON,DENCOH,BPOR,SUSARM, + KSED,WS1MIN C COMMON /SED2/IBMSK(IM,JM),LAYER(LAYMAX,IM,JM) C COMMON /SED3/COUT(IM,JM,KB),REXP(IM,JM),CSED(20),CBED(IM,JM), . CSED1(IM,JM,KB),CSED2(IM,JM,KB),WCT1BOT(IM,JM),WCT2BOT(IM,JM), . CSED1B(IM,JM,KB),CSED2B(IM,JM,KB), . TAU(IM,JM,KB),TAUCR(LAYMAX),TCRIT(LAYMAX,IM,JM), . FTIME2(LAYMAX),TAUMAX(LAYMAX,IM,JM), . ASED(LAYMAX,IM,JM),A0(LAYMAX,IM,JM),TAUCUR(LAYMAX,IM,JM), . WS1(6250000),TAUD2(LAYMAX,2000),FTIME(LAYMAX), . WSET1(IM,JM,KB),WSET2(IM,JM,KB),UCRVAR(IM,JM), . CSI(KSMAX),CBDRYSL(KSMAX,EBCM,KSL),CBDRY(KSMAX,EBCM,KBM1), . DCDIS(KSMAX,QBCM,2),DCBDRY(KSMAX,EBCM,KBM1,2), . CDIS(KSMAX,QBCM),CDIFF(KSMAX,DBCM) C C FOR SGI COMPILING C COMMON /SED5/TSED COMMON /SED6/EBTOT COMMON /SED7/EBMAX COMMON /SED8/TSED0IN COMMON /SED8A/TSED0 COMMON /SED9/EBCUR COMMON /SED10/PSED1 COMMON /SED11/PSED2 COMMON /SED12/P0 COMMON /SED13/THICK COMMON /SED14/PTOT COMMON /SED15/TSEDT COMMON /SED16/TSET0T COMMON /SED17/E COMMON /SED18/DD COMMON /SED19/BEDTHI COMMON /SED20/BEDTH COMMON /SED21/FPBED COMMON /SED22/ACTLAY COMMON /SED23/CARMOR COMMON /SED24/FALAY COMMON /SED25/WINDSP,WINDIR,DWDS(2),DWDD(2),PI COMMON /SED26/FETCH(36,IM,JM),HMEAN(36,IM,JM),Z0B,Z0BCOH COMMON /SED27/XWCOR(IM,JM),YWCOR(IM,JM) COMMON /SED28/ARCWS1(IM,JM,KB),ARCWS2(IM,JM,KB), + WS1SAVE(VPTSM,KB),WS2SAVE(VPTSM,KB) COMMON /SED30/ARCSED1(IM,JM,KB),ARCSED2(IM,JM,KB), + ARCTAU(IM,JM,KB), + C1SAVE(VPTSM,KB),C2SAVE(VPTSM,KB), + TAUSAVE(VPTSM,KB) COMMON /SED31/ARCTHIK COMMON /SED32/THSAVE COMMON /SED33/SXADV(KSMAX,IM,JM,KB),SYADV(KSMAX,IM,JM,KB), + SZADV(KSMAX,IM,JM,KB),SXDIF(KSMAX,IM,JM,KB), + SYDIF(KSMAX,IM,JM,KB),SZDIF(KSMAX,IM,JM,KB) COMMON /SED34/ADEP,DEPEXP,PDEP(20000) COMMON /SED35/DETA(IM,JM) COMMON /SED36/USHEAR(IM,JM),QBAR(IM,JM),UBAR(IM,JM,KB), + VBAR(IM,JM,KB) COMMON /SED37/WVHT(IM,JM),WVPD(IM,JM),WVDR(IM,JM),CURANG(IM,JM) COMMON /SED38/CSDIFF(KSMAX,DBCM),DCSDIFF(KSMAX,DBCM,2) COMMON /SED39/REXP35(IM,JM),REXP59(IM,JM), + ASED35(LAYMAX,IM,JM),A035(LAYMAX,IM,JM), + ASED59(LAYMAX,IM,JM),A059(LAYMAX,IM,JM) COMMON /SED40/SEDTHK(IM,JM) COMMON /ZWAVE/Z0WAVE,BFRIC,BFCOH,NWAVE,NWAVECNT,WIFRIC ! NG FOR ICE C C CHEMICAL TRANSPORT C COMMON /CHEM1/CHEM1(IM,JM,KB),CHEM2(IM,JM,KB), + CHEM1B(IM,JM,KB),CHEM2B(IM,JM,KB), + CHEMBOT1(IM,JM),CHEMBOT2(IM,JM),CHEMDIS(2,QBCM), + CHEMI(2),CHEMBDRY(2,EBCM,KBM1), + PBDRY(2,EBCM,KBM1),DPDIS(2,QBCM,2), + DPBDRY(2,EBCM,KBM1,2),PDIS(2,QBCM), + DPDIFF(2,DBCM,2),PDIFF(2,DBCM) COMMON /CHEM2/NCHEMLAY,NACTLAY,CHEMTHIK,CHEMACT, + CBEDCHEM(NBEDMAX,IM,JM) COMMON /CHEM3/SEDMASS COMMON /CHEM4/CHEMMASS COMMON /CHEM5/SEDDEP COMMON /CHEM6/ARCHEM1(IM,JM,KB),ARCHEM2(IM,JM,KB), + ARCPBED(NBEDMAX,IM,JM),ARCDEP(IM,JM), + P1SAVE(VPTSM,KB),P2SAVE(VPTSM,KB), + PBEDSAVE(VPTSM,NBEDMAX) C C MODIFIED BY C.K.Z. ON 12/10/97 C C FOR SHORELINE EROSION OPTION C COMMON /SERODE1/SHERODE COMMON /SERODE2/NSERODE,KSE(20),ISERODE(20,50),JSERODE(20,50) COMMON /SERODE3/TSERODE(10000),SEMASS(10000,20) COMMON /SERODE4/SHERFLX(IM,JM) C C FOR PARTICLE TRACKING C COMMON /PTRACK0/PARTICLE C COMMON /PTRACK1/NFREQ,NPART,NCONV,IRELST,NPCLASS,NSOURCE, + NGRADELOOP C COMMON /PTRACK2/CP(IM,JM,KB) C COMMON /PTRACK3/XP(NSRCM,NPARTM,NCONVM),YP(NSRCM,NPARTM,NCONVM), + ZP(NSRCM,NPARTM,NCONVM),IP(NSRCM,NPARTM,NCONVM), + JP(NSRCM,NPARTM,NCONVM),KP(NSRCM,NPARTM,NCONVM), + INOUT(NSRCM,NPARTM,NCONVM) C COMMON /PTRACK4/PZERO1(DBCM,KBM1),PZERO2(QBCM),PZERO3(EBCM,KBM1) C COMMON /PTRACK5/XCOR(IM+1,JM+1),YCOR(IM+1,JM+1),H1P(IM,JM), + H2P(IM,JM) C COMMON /PTRACK6/IPART(NPARTM),JPART(NPARTM),KPART(NPARTM) C COMMON /PTRACK7/ISOURCE(NSRCM),JSOURCE(NSRCM),KSOURCE(NSRCM) C C FOR WAVE MODEL BASED ON SCHWAB MODEL BY QA C CNG COMMON /WAVEBLK1/CO(IM,JM),XMOM(IM,JM),YMOM(IM,JM), COMMON /WAVEBLK1/FSMW(IM,JM),CO(IM,JM),XMOM(IM,JM),YMOM(IM,JM), 2 SIG(IM,JM),CS(IM,JM),SN(IM,JM),NR,WN(IM,JM),CREADWBC C C SPECIFY WAVE PARAMTERS ON OPEN BOUNDARIES (hli,03/19/04) C COMMON /WAVEBLK3/WHBDRY(EBCM),WPBDRY(EBCM),WDBDRY(EBCM), & WHBFLX(EBCM),WPBFLX(EBCM),WDBFLX(EBCM), & DWHBFLX(EBCM,2),DWPBFLX(EBCM,2),DWDBFLX(EBCM,2), & IUT195,IUT802 C C C FOR DIFFERENT TYPE OF OPEN BOUNDARY CONDITIONS BY Quamrul QA 12/28/98 C CHARACTER*7 BCTYPE COMMON /OBCBLK/BCTYPE C END OF COMDECK