84 |
& saltClimFile, buoyancyRelation, |
& saltClimFile, buoyancyRelation, |
85 |
& EmPmRfile, saltFluxFile, |
& EmPmRfile, saltFluxFile, |
86 |
& surfQfile, surfQnetFile, surfQswFile, |
& surfQfile, surfQnetFile, surfQswFile, |
87 |
|
& lambdaThetaFile, lambdaSaltFile, |
88 |
& uVelInitFile, vVelInitFile, pSurfInitFile, |
& uVelInitFile, vVelInitFile, pSurfInitFile, |
89 |
& dQdTfile, ploadFile, |
& dQdTfile, ploadFile, |
90 |
& eddyTauxFile, eddyTauyFile, |
& eddyTauxFile, eddyTauyFile, |
112 |
CHARACTER*(MAX_LEN_FNAM) ploadFile |
CHARACTER*(MAX_LEN_FNAM) ploadFile |
113 |
CHARACTER*(MAX_LEN_FNAM) eddyTauxFile |
CHARACTER*(MAX_LEN_FNAM) eddyTauxFile |
114 |
CHARACTER*(MAX_LEN_FNAM) eddyTauyFile |
CHARACTER*(MAX_LEN_FNAM) eddyTauyFile |
115 |
|
CHARACTER*(MAX_LEN_FNAM) lambdaThetaFile |
116 |
|
CHARACTER*(MAX_LEN_FNAM) lambdaSaltFile |
117 |
CHARACTER*(MAX_LEN_FNAM) mdsioLocalDir |
CHARACTER*(MAX_LEN_FNAM) mdsioLocalDir |
118 |
CHARACTER*(MAX_LEN_FNAM) the_run_name |
CHARACTER*(MAX_LEN_FNAM) the_run_name |
119 |
CHARACTER*(6) eosType |
CHARACTER*(6) eosType |
220 |
C if F, converts P+R-E to salt flux (no SL effect) |
C if F, converts P+R-E to salt flux (no SL effect) |
221 |
C useFullLeith :: Set to true to use full Leith viscosity (may be unstable |
C useFullLeith :: Set to true to use full Leith viscosity (may be unstable |
222 |
C on irregular grids) |
C on irregular grids) |
223 |
|
C useAnisotropicViscAGridMax :: Set to true to use Alistair's latest |
224 |
|
C anisotropic length scale. It is used only for maximum viscosity |
225 |
|
C calculations. Alistair recommends a value of viscA*GridMax=.25 |
226 |
C rigidLid :: Set to true to use rigid lid |
C rigidLid :: Set to true to use rigid lid |
227 |
C implicitFreeSurface :: Set to true to use implcit free surface |
C implicitFreeSurface :: Set to true to use implcit free surface |
228 |
C exactConserv :: Set to true to conserve exactly the total Volume |
C exactConserv :: Set to true to conserve exactly the total Volume |
250 |
C useMultiDimAdvec :: True if multi-dim advection is used at least once |
C useMultiDimAdvec :: True if multi-dim advection is used at least once |
251 |
C forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution |
C forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution |
252 |
C out off Adams-Bashforth time stepping. |
C out off Adams-Bashforth time stepping. |
253 |
|
C startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup |
254 |
C doThetaClimRelax :: Set true if relaxation to temperature |
C doThetaClimRelax :: Set true if relaxation to temperature |
255 |
C climatology is required. |
C climatology is required. |
256 |
C doSaltClimRelax :: Set true if relaxation to salinity |
C doSaltClimRelax :: Set true if relaxation to salinity |
296 |
& tempAdvection, tempForcing, |
& tempAdvection, tempForcing, |
297 |
& saltAdvection, saltForcing, |
& saltAdvection, saltForcing, |
298 |
& useRealFreshWaterFlux, |
& useRealFreshWaterFlux, |
299 |
& useFullLeith, |
& useFullLeith, useAnisotropicViscAGridMax, |
300 |
& rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf, |
& rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf, |
301 |
& momStepping, tempStepping, saltStepping, |
& momStepping, tempStepping, saltStepping, |
302 |
& metricTerms, usingSphericalPolarMTerms, useNHMTerms, |
& metricTerms, usingSphericalPolarMTerms, useNHMTerms, |
314 |
& usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr, |
& usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr, |
315 |
& nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO, |
& nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO, |
316 |
& allowFreezing, useOldFreezing, groundAtK1, |
& allowFreezing, useOldFreezing, groundAtK1, |
317 |
& usePickupBeforeC35, usePickupBeforeC54, debugMode, |
& usePickupBeforeC35, usePickupBeforeC54, startFromPickupAB2, |
|
& readPickupWithTracer, writePickupWithTracer, |
|
318 |
& pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed, |
& pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed, |
319 |
& timeave_mdsio, snapshot_mdsio, monitor_stdio, |
& timeave_mdsio, snapshot_mdsio, monitor_stdio, |
320 |
& outputTypesInclusive, |
& outputTypesInclusive, debugMode, |
321 |
& inAdMode, inAdTrue, inAdFalse, inAdExact, |
& inAdMode, inAdTrue, inAdFalse, inAdExact, |
322 |
& calendarDumps |
& calendarDumps |
323 |
|
|
341 |
LOGICAL saltAdvection |
LOGICAL saltAdvection |
342 |
LOGICAL saltForcing |
LOGICAL saltForcing |
343 |
LOGICAL useRealFreshWaterFlux |
LOGICAL useRealFreshWaterFlux |
344 |
LOGICAL useFullLeith |
LOGICAL useFullLeith, useAnisotropicViscAGridMax |
345 |
LOGICAL rigidLid |
LOGICAL rigidLid |
346 |
LOGICAL implicitFreeSurface |
LOGICAL implicitFreeSurface |
347 |
LOGICAL exactConserv |
LOGICAL exactConserv |
388 |
LOGICAL groundAtK1 |
LOGICAL groundAtK1 |
389 |
LOGICAL usePickupBeforeC35 |
LOGICAL usePickupBeforeC35 |
390 |
LOGICAL usePickupBeforeC54 |
LOGICAL usePickupBeforeC54 |
391 |
|
LOGICAL startFromPickupAB2 |
392 |
LOGICAL debugMode |
LOGICAL debugMode |
|
LOGICAL readPickupWithTracer |
|
|
LOGICAL writePickupWithTracer |
|
393 |
LOGICAL pickup_read_mdsio, pickup_write_mdsio |
LOGICAL pickup_read_mdsio, pickup_write_mdsio |
394 |
LOGICAL pickup_write_immed |
LOGICAL pickup_write_immed |
395 |
LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio |
LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio |
430 |
C rhoConstFresh :: Constant reference density for fresh water (rain) |
C rhoConstFresh :: Constant reference density for fresh water (rain) |
431 |
C tRef :: reference vertical profile for potential temperature |
C tRef :: reference vertical profile for potential temperature |
432 |
C sRef :: reference vertical profile for salinity/specific humidity |
C sRef :: reference vertical profile for salinity/specific humidity |
|
C startTime :: Start time for model ( s ) |
|
433 |
C phiMin :: Latitude of southern most cell face. |
C phiMin :: Latitude of southern most cell face. |
434 |
C thetaMin :: Longitude of western most cell face (this |
C thetaMin :: Longitude of western most cell face (this |
435 |
C is an "inert" parameter but it is included |
C is an "inert" parameter but it is included |
505 |
C freesurfFac :: Parameter to turn implicit free surface term on or off |
C freesurfFac :: Parameter to turn implicit free surface term on or off |
506 |
C freesurfac = 1. uses implicit free surface |
C freesurfac = 1. uses implicit free surface |
507 |
C freesurfac = 0. uses rigid lid |
C freesurfac = 0. uses rigid lid |
508 |
|
C abEps :: Adams-Bashforth-2 stabilizing weight |
509 |
|
C alph_AB :: Adams-Bashforth-3 primary factor |
510 |
|
C beta_AB :: Adams-Bashforth-3 secondary factor |
511 |
C implicSurfPress :: parameter of the Crank-Nickelson time stepping : |
C implicSurfPress :: parameter of the Crank-Nickelson time stepping : |
512 |
C Implicit part of Surface Pressure Gradient ( 0-1 ) |
C Implicit part of Surface Pressure Gradient ( 0-1 ) |
513 |
C implicDiv2Dflow :: parameter of the Crank-Nickelson time stepping : |
C implicDiv2Dflow :: parameter of the Crank-Nickelson time stepping : |
520 |
C hFacSup that control vanishing and creating levels |
C hFacSup that control vanishing and creating levels |
521 |
C tauCD :: CD scheme coupling timescale ( 1/s ) |
C tauCD :: CD scheme coupling timescale ( 1/s ) |
522 |
C rCD :: CD scheme normalised coupling parameter ( 0-1 ) |
C rCD :: CD scheme normalised coupling parameter ( 0-1 ) |
523 |
|
C baseTime :: model base time (time origin) = time @ iteration zero |
524 |
C startTime :: Starting time for this integration ( s ). |
C startTime :: Starting time for this integration ( s ). |
525 |
C endTime :: Ending time for this integration ( s ). |
C endTime :: Ending time for this integration ( s ). |
526 |
C chkPtFreq :: Frequency of rolling check pointing ( s ). |
C chkPtFreq :: Frequency of rolling check pointing ( s ). |
544 |
C step (of each taveFreq period) put in the time average. |
C step (of each taveFreq period) put in the time average. |
545 |
C (fraction for 1rst iter = 1 - tave_lastIter) |
C (fraction for 1rst iter = 1 - tave_lastIter) |
546 |
C tauThetaClimRelax :: Relaxation to climatology time scale ( s ). |
C tauThetaClimRelax :: Relaxation to climatology time scale ( s ). |
|
C lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ). |
|
547 |
C tauSaltClimRelax :: Relaxation to climatology time scale ( s ). |
C tauSaltClimRelax :: Relaxation to climatology time scale ( s ). |
|
C lambdaSaltClimRelax :: Inverse time scale for relaxation ( 1/s ). |
|
548 |
C latBandClimRelax :: latitude band where Relaxation to Clim. is applied, |
C latBandClimRelax :: latitude band where Relaxation to Clim. is applied, |
549 |
C i.e. where |yC| <= latBandClimRelax |
C i.e. where |yC| <= latBandClimRelax |
550 |
C externForcingPeriod :: Is the period of which forcing varies (eg. 1 month) |
C externForcingPeriod :: Is the period of which forcing varies (eg. 1 month) |
569 |
& cg2dpcOffDFac, cg3dTargetResidual, |
& cg2dpcOffDFac, cg3dTargetResidual, |
570 |
& delR, delRc, delX, delY, |
& delR, delRc, delX, delY, |
571 |
& deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock, |
& deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock, |
572 |
& abeps, startTime, |
& abEps, alph_AB, beta_AB, |
573 |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
& phiMin, thetaMin, rSphere, recip_RSphere, f0, beta, |
|
& fCori, fCoriG, fCoriCos, |
|
574 |
& viscAh, viscAhW, viscAhMax, |
& viscAh, viscAhW, viscAhMax, |
575 |
& viscAhGrid, viscAhGridMax, viscAhGridMin, |
& viscAhGrid, viscAhGridMax, viscAhGridMin, |
576 |
& viscC2leith, viscC2leithD, |
& viscC2leith, viscC2leithD, |
578 |
& viscAhD, viscAhZ, viscA4D, viscA4Z, |
& viscAhD, viscAhZ, viscA4D, viscA4Z, |
579 |
& viscA4, viscA4W, |
& viscA4, viscA4W, |
580 |
& viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin, |
& viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin, |
581 |
& viscC4leith, viscAr, viscAstrain, viscAtension, |
& viscC4leith, viscC4leithD, viscAr, viscAstrain, viscAtension, |
582 |
& diffKhT, diffK4T, diffKrNrT, |
& diffKhT, diffK4T, diffKrNrT, |
583 |
& diffKhS, diffK4S, diffKrNrS, |
& diffKhS, diffK4S, diffKrNrS, |
584 |
& diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho, |
& diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho, |
587 |
& gravity, recip_Gravity, gBaro, rhonil, recip_rhonil, |
& gravity, recip_Gravity, gBaro, rhonil, recip_rhonil, |
588 |
& recip_rhoConst, rhoConst, |
& recip_rhoConst, rhoConst, |
589 |
& rhoConstFresh, convertEmP2rUnit, tRef, sRef, |
& rhoConstFresh, convertEmP2rUnit, tRef, sRef, |
590 |
& endTime, chkPtFreq, pchkPtFreq, dumpFreq, adjDumpFreq, |
& baseTime, startTime, endTime, |
591 |
|
& chkPtFreq, pchkPtFreq, dumpFreq, adjDumpFreq, |
592 |
& diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq, |
& diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq, |
593 |
& afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom, |
& afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom, |
594 |
& cosPower, cAdjFreq, omega, rotationPeriod, |
& cosPower, cAdjFreq, omega, rotationPeriod, |
595 |
& tauThetaClimRelax, lambdaThetaClimRelax, |
& tauThetaClimRelax, |
596 |
& tauSaltClimRelax, lambdaSaltClimRelax, |
& tauSaltClimRelax, |
597 |
& tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax, |
& tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax, |
598 |
& externForcingCycle, externForcingPeriod, |
& externForcingCycle, externForcingPeriod, |
599 |
& convertFW2Salt, temp_EvPrRn, salt_EvPrRn, |
& convertFW2Salt, temp_EvPrRn, salt_EvPrRn, |
616 |
_RL deltaTmom |
_RL deltaTmom |
617 |
_RL dTtracerLev(Nr) |
_RL dTtracerLev(Nr) |
618 |
_RL deltaTfreesurf |
_RL deltaTfreesurf |
619 |
_RL abeps |
_RL abEps, alph_AB, beta_AB |
620 |
_RL phiMin |
_RL phiMin |
621 |
_RL thetaMin |
_RL thetaMin |
622 |
_RL rSphere |
_RL rSphere |
678 |
_RL convertEmP2rUnit |
_RL convertEmP2rUnit |
679 |
_RL tRef(Nr) |
_RL tRef(Nr) |
680 |
_RL sRef(Nr) |
_RL sRef(Nr) |
681 |
_RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL baseTime |
|
_RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS fCoriCos(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
682 |
_RL startTime |
_RL startTime |
683 |
_RL endTime |
_RL endTime |
684 |
_RL chkPtFreq |
_RL chkPtFreq |
701 |
_RL omega |
_RL omega |
702 |
_RL rotationPeriod |
_RL rotationPeriod |
703 |
_RL tauThetaClimRelax |
_RL tauThetaClimRelax |
|
_RL lambdaThetaClimRelax |
|
704 |
_RL tauSaltClimRelax |
_RL tauSaltClimRelax |
|
_RL lambdaSaltClimRelax |
|
705 |
_RL tauTr1ClimRelax |
_RL tauTr1ClimRelax |
706 |
_RL lambdaTr1ClimRelax |
_RL lambdaTr1ClimRelax |
707 |
_RL latBandClimRelax |
_RL latBandClimRelax |