33 |
C adTapeDir :: read-write checkpointing tape files from/to this |
C adTapeDir :: read-write checkpointing tape files from/to this |
34 |
C directory name instead of current dir. Conflicts |
C directory name instead of current dir. Conflicts |
35 |
C mdsioLocalDir, so only one of the two can be set. |
C mdsioLocalDir, so only one of the two can be set. |
|
C In contrast to mdsioLocalDir, if specified adTapeDir |
|
|
C must exist before the model starts. |
|
36 |
C tRefFile :: File containing reference Potential Temperat. tRef (1.D) |
C tRefFile :: File containing reference Potential Temperat. tRef (1.D) |
37 |
C sRefFile :: File containing reference salinity/spec.humid. sRef (1.D) |
C sRefFile :: File containing reference salinity/spec.humid. sRef (1.D) |
38 |
C rhoRefFile :: File containing reference density profile rhoRef (1.D) |
C rhoRefFile :: File containing reference density profile rhoRef (1.D) |
160 |
C nIter0 :: Start time-step number of for this run |
C nIter0 :: Start time-step number of for this run |
161 |
C nTimeSteps :: Number of timesteps to execute |
C nTimeSteps :: Number of timesteps to execute |
162 |
C nTimeSteps_l2 :: Number of inner timesteps to execute per timestep |
C nTimeSteps_l2 :: Number of inner timesteps to execute per timestep |
|
C writeStatePrec :: Precision used for writing model state. |
|
|
C writeBinaryPrec :: Precision used for writing binary files |
|
|
C readBinaryPrec :: Precision used for reading binary files |
|
163 |
C selectCoriMap :: select setting of Coriolis parameter map: |
C selectCoriMap :: select setting of Coriolis parameter map: |
164 |
C =0 f-Plane (Constant Coriolis, = f0) |
C =0 f-Plane (Constant Coriolis, = f0) |
165 |
C =1 Beta-Plane Coriolis (= f0 + beta.y) |
C =1 Beta-Plane Coriolis (= f0 + beta.y) |
199 |
C =0: average KE from grid center to U & V location |
C =0: average KE from grid center to U & V location |
200 |
C =1: use local velocity norm @ U & V location |
C =1: use local velocity norm @ U & V location |
201 |
C =2: same with wet-point averaging of other component |
C =2: same with wet-point averaging of other component |
202 |
|
C readBinaryPrec :: Precision used for reading binary files |
203 |
|
C writeStatePrec :: Precision used for writing model state. |
204 |
|
C writeBinaryPrec :: Precision used for writing binary files |
205 |
|
C rwSuffixType :: controls the format of the mds file suffix. |
206 |
|
C =0 (default): use iteration number (myIter, I10.10); |
207 |
|
C =1: 100*myTime (100th sec); =2: myTime (seconds); |
208 |
|
C =3: myTime/360 (10th of hr); =4: myTime/3600 (hours). |
209 |
C monitorSelect :: select group of variables to monitor |
C monitorSelect :: select group of variables to monitor |
210 |
C =1 : dynvars ; =2 : + vort ; =3 : + surface |
C =1 : dynvars ; =2 : + vort ; =3 : + surface |
211 |
C- debugLevel :: controls printing of algorithm intermediate results |
C- debugLevel :: controls printing of algorithm intermediate results |
212 |
C and statistics ; higher -> more writing |
C and statistics ; higher -> more writing |
213 |
|
C- plotLevel :: controls printing of field maps ; higher -> more flds |
214 |
|
|
215 |
COMMON /PARM_I/ |
COMMON /PARM_I/ |
216 |
& cg2dMaxIters, cg2dChkResFreq, |
& cg2dMaxIters, cg2dChkResFreq, |
218 |
& cg3dMaxIters, cg3dChkResFreq, |
& cg3dMaxIters, cg3dChkResFreq, |
219 |
& printResidualFreq, |
& printResidualFreq, |
220 |
& nIter0, nTimeSteps, nTimeSteps_l2, nEndIter, |
& nIter0, nTimeSteps, nTimeSteps_l2, nEndIter, |
|
& writeStatePrec, |
|
|
& writeBinaryPrec, readBinaryPrec, |
|
221 |
& selectCoriMap, |
& selectCoriMap, |
222 |
& selectSigmaCoord, |
& selectSigmaCoord, |
223 |
& nonlinFreeSurf, select_rStar, |
& nonlinFreeSurf, select_rStar, |
228 |
& saltAdvScheme, saltVertAdvScheme, |
& saltAdvScheme, saltVertAdvScheme, |
229 |
& selectKEscheme, selectVortScheme, |
& selectKEscheme, selectVortScheme, |
230 |
& selectBotDragQuadr, |
& selectBotDragQuadr, |
231 |
& monitorSelect, debugLevel |
& readBinaryPrec, writeBinaryPrec, writeStatePrec, |
232 |
|
& rwSuffixType, monitorSelect, debugLevel, plotLevel |
233 |
INTEGER cg2dMaxIters |
INTEGER cg2dMaxIters |
234 |
INTEGER cg2dChkResFreq |
INTEGER cg2dChkResFreq |
235 |
INTEGER cg2dPreCondFreq |
INTEGER cg2dPreCondFreq |
241 |
INTEGER nTimeSteps |
INTEGER nTimeSteps |
242 |
INTEGER nTimeSteps_l2 |
INTEGER nTimeSteps_l2 |
243 |
INTEGER nEndIter |
INTEGER nEndIter |
|
INTEGER writeStatePrec |
|
|
INTEGER writeBinaryPrec |
|
|
INTEGER readBinaryPrec |
|
244 |
INTEGER selectCoriMap |
INTEGER selectCoriMap |
245 |
INTEGER selectSigmaCoord |
INTEGER selectSigmaCoord |
246 |
INTEGER nonlinFreeSurf |
INTEGER nonlinFreeSurf |
255 |
INTEGER selectKEscheme |
INTEGER selectKEscheme |
256 |
INTEGER selectVortScheme |
INTEGER selectVortScheme |
257 |
INTEGER selectBotDragQuadr |
INTEGER selectBotDragQuadr |
258 |
|
INTEGER readBinaryPrec |
259 |
|
INTEGER writeStatePrec |
260 |
|
INTEGER writeBinaryPrec |
261 |
|
INTEGER rwSuffixType |
262 |
INTEGER monitorSelect |
INTEGER monitorSelect |
263 |
INTEGER debugLevel |
INTEGER debugLevel |
264 |
|
INTEGER plotLevel |
265 |
|
|
266 |
C-- COMMON /PARM_L/ Logical valued parameters used by the model. |
C-- COMMON /PARM_L/ Logical valued parameters used by the model. |
267 |
C- Coordinate + Grid params: |
C- Coordinate + Grid params: |
285 |
C deepAtmosphere :: deep model (drop the shallow-atmosphere approximation) |
C deepAtmosphere :: deep model (drop the shallow-atmosphere approximation) |
286 |
C setInterFDr :: set Interface depth (put cell-Center at the middle) |
C setInterFDr :: set Interface depth (put cell-Center at the middle) |
287 |
C setCenterDr :: set cell-Center depth (put Interface at the middle) |
C setCenterDr :: set cell-Center depth (put Interface at the middle) |
288 |
|
C useMin4hFacEdges :: set hFacW,hFacS as minimum of adjacent hFacC factor |
289 |
C- Momentum params: |
C- Momentum params: |
290 |
C no_slip_sides :: Impose "no-slip" at lateral boundaries. |
C no_slip_sides :: Impose "no-slip" at lateral boundaries. |
291 |
C no_slip_bottom :: Impose "no-slip" at bottom boundary. |
C no_slip_bottom :: Impose "no-slip" at bottom boundary. |
399 |
C monitor_stdio :: use stdio for monitor output |
C monitor_stdio :: use stdio for monitor output |
400 |
C dumpInitAndLast :: dumps model state to files at Initial (nIter0) |
C dumpInitAndLast :: dumps model state to files at Initial (nIter0) |
401 |
C & Last iteration, in addition multiple of dumpFreq iter. |
C & Last iteration, in addition multiple of dumpFreq iter. |
|
C printDomain :: controls printing of domain fields (bathy, hFac ...). |
|
402 |
|
|
403 |
COMMON /PARM_L/ |
COMMON /PARM_L/ |
404 |
& fluidIsAir, fluidIsWater, |
& fluidIsAir, fluidIsWater, |
405 |
& usingPCoords, usingZCoords, |
& usingPCoords, usingZCoords, |
406 |
& usingCartesianGrid, usingSphericalPolarGrid, rotateGrid, |
& usingCartesianGrid, usingSphericalPolarGrid, rotateGrid, |
407 |
& usingCylindricalGrid, usingCurvilinearGrid, hasWetCSCorners, |
& usingCylindricalGrid, usingCurvilinearGrid, hasWetCSCorners, |
408 |
& deepAtmosphere, setInterFDr, setCenterDr, |
& deepAtmosphere, setInterFDr, setCenterDr, useMin4hFacEdges, |
409 |
& no_slip_sides, no_slip_bottom, bottomVisc_pCell, useSmag3D, |
& no_slip_sides, no_slip_bottom, bottomVisc_pCell, useSmag3D, |
410 |
& useFullLeith, useStrainTensionVisc, useAreaViscLength, |
& useFullLeith, useStrainTensionVisc, useAreaViscLength, |
411 |
& momViscosity, momAdvection, momForcing, |
& momViscosity, momAdvection, momForcing, |
441 |
& pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed, |
& pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed, |
442 |
& writePickupAtEnd, |
& writePickupAtEnd, |
443 |
& timeave_mdsio, snapshot_mdsio, monitor_stdio, |
& timeave_mdsio, snapshot_mdsio, monitor_stdio, |
444 |
& outputTypesInclusive, dumpInitAndLast, |
& outputTypesInclusive, dumpInitAndLast |
|
& printDomain |
|
445 |
|
|
446 |
LOGICAL fluidIsAir |
LOGICAL fluidIsAir |
447 |
LOGICAL fluidIsWater |
LOGICAL fluidIsWater |
454 |
LOGICAL deepAtmosphere |
LOGICAL deepAtmosphere |
455 |
LOGICAL setInterFDr |
LOGICAL setInterFDr |
456 |
LOGICAL setCenterDr |
LOGICAL setCenterDr |
457 |
|
LOGICAL useMin4hFacEdges |
458 |
|
|
459 |
LOGICAL no_slip_sides |
LOGICAL no_slip_sides |
460 |
LOGICAL no_slip_bottom |
LOGICAL no_slip_bottom |
541 |
LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio |
LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio |
542 |
LOGICAL outputTypesInclusive |
LOGICAL outputTypesInclusive |
543 |
LOGICAL dumpInitAndLast |
LOGICAL dumpInitAndLast |
|
LOGICAL printDomain |
|
544 |
|
|
545 |
C-- COMMON /PARM_R/ "Real" valued parameters used by the model. |
C-- COMMON /PARM_R/ "Real" valued parameters used by the model. |
546 |
C cg2dTargetResidual |
C cg2dTargetResidual |
964 |
_RL cg2dMinColumnEps, pReleaseVisc, pReleaseDamp |
_RL cg2dMinColumnEps, pReleaseVisc, pReleaseDamp |
965 |
#endif |
#endif |
966 |
|
|
967 |
|
|
968 |
C-- COMMON /PARM_A/ Thermodynamics constants ? |
C-- COMMON /PARM_A/ Thermodynamics constants ? |
969 |
COMMON /PARM_A/ HeatCapacity_Cp |
COMMON /PARM_A/ HeatCapacity_Cp |
970 |
_RL HeatCapacity_Cp |
_RL HeatCapacity_Cp |