/[MITgcm]/MITgcm/model/src/set_defaults.F
ViewVC logotype

Diff of /MITgcm/model/src/set_defaults.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.116 by heimbach, Wed Jun 7 01:55:13 2006 UTC revision 1.156 by jmc, Wed Dec 7 23:31:15 2011 UTC
# Line 9  C     !ROUTINE: SET_DEFAULTS Line 9  C     !ROUTINE: SET_DEFAULTS
9  C     !INTERFACE:  C     !INTERFACE:
10        SUBROUTINE SET_DEFAULTS(        SUBROUTINE SET_DEFAULTS(
11       O   viscArDefault, diffKrTDefault, diffKrSDefault,       O   viscArDefault, diffKrTDefault, diffKrSDefault,
12       O   hFacMinDrDefault, delRdefault, rkFacDefault,       O   hFacMinDrDefault, delRdefault,
13       I   myThid )       I   myThid )
14    
15  C     !DESCRIPTION:  C     !DESCRIPTION:
# Line 21  C     !USES: Line 21  C     !USES:
21  #include "EEPARAMS.h"  #include "EEPARAMS.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
23  Cml#include "EOS.h"  Cml#include "EOS.h"
24  #include "GRID.h"  c#include "GRID.h"
25    
26  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
27  C     myThid - Number of this instance of INI_PARMS  C     myThid :: my Thread Id. Number
28        INTEGER myThid        INTEGER myThid
29        _RL viscArDefault        _RL viscArDefault
30        _RL diffKrTDefault        _RL diffKrTDefault
31        _RL diffKrSDefault        _RL diffKrSDefault
32        _RL hFacMinDrDefault        _RL hFacMinDrDefault
33        _RL delRDefault(Nr)        _RL delRDefault(Nr)
       _RS rkFacDefault  
34    
35  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
36  C     K, I, J - Loop counters  C     i, j, k :: Loop counters
37        INTEGER K, I, J        INTEGER i, j, k
38  CEOP  CEOP
39    
40  C--   Grid parameters  C--   Grid parameters
41  C     Vertical gridding  C-    Vertical gridding
42        rkFacDefault        = 1. _d 0        delRFile            = ' '
43        horiVertRatio       = 1. _d 0        delRcFile           = ' '
44          hybSigmFile         = ' '
45        Ro_SeaLevel         = 0.        Ro_SeaLevel         = 0.
46          rSigmaBnd           = UNSET_RL
47          selectSigmaCoord    = 0
48        DO k=1,Nr        DO k=1,Nr
49         delRdefault(k)     = 0.         delRdefault(k)     = 0.
50        ENDDO        ENDDO
51        DO k=1,Nr+1        DO k=1,Nr+1
52         delRc(k)           = UNSET_RL         delRc(k)           = UNSET_RL
53        ENDDO        ENDDO
54    C-    vertical profile
55          tRefFile            = ' '
56          sRefFile            = ' '
57          rhoRefFile          = ' '
58    
59  C     Horizontal gridding  C-    Horizontal gridding
60        delXFile            = ' '        delXFile            = ' '
61        delYFile            = ' '        delYFile            = ' '
62        horizGridFile       = ' '        horizGridFile       = ' '
63  C     In cartesian coords distances are in metres        deepAtmosphere      = .FALSE.
64        usingCartesianGrid  = .FALSE.        xgOrigin            = UNSET_RL
65          ygOrigin            = UNSET_RL
66        DO i=1,Nx        DO i=1,Nx
67         delX(i)            = UNSET_RL         delX(i)            = UNSET_RL
68        ENDDO        ENDDO
69        DO j=1,Ny        DO j=1,Ny
70         delY(j)            = UNSET_RL         delY(j)            = UNSET_RL
71        ENDDO        ENDDO
72    C     In cartesian coords distances are in metres
73          usingCartesianGrid  = .FALSE.
74  C     In spherical polar distances are in degrees  C     In spherical polar distances are in degrees
75        usingSphericalPolarGrid = .FALSE.        usingSphericalPolarGrid = .FALSE.
76        phiMin              = 0.0        rSphere             = UNSET_RL
       thetaMin            = 0.  
       rSphere             = 6370. _d 3  
77  C     General curvilinear coordinate system  C     General curvilinear coordinate system
78        usingCurvilinearGrid= .FALSE.        usingCurvilinearGrid= .FALSE.
79          radius_fromHorizGrid= UNSET_RL
80  C     General cylindrical coordinate system  C     General cylindrical coordinate system
81        usingCylindricalGrid= .FALSE.        usingCylindricalGrid= .FALSE.
82  C     Coriolis map:  C     Coriolis map:
83        useConstantF        = .FALSE.        selectCoriMap       = -1
84        useBetaPlaneF       = .FALSE.        use3dCoriolis       = .TRUE.
85        useSphereF          = .FALSE.  C     grid rotation
86          rotateGrid          = .FALSE.
87          phiEuler            = 0. _d 0
88          thetaEuler          = 0. _d 0
89          psiEuler            = 0. _d 0
90    
91  C--   Set default "physical" parameters  C--   Set default "physical" parameters
92        nh_Am2              = 1. _d 0        nh_Am2              = 1. _d 0
# Line 86  C         (so that the default produces Line 97  C         (so that the default produces
97  c     rhoConstFresh       = 999.8 _d 0  c     rhoConstFresh       = 999.8 _d 0
98        f0                  = 1. _d -4        f0                  = 1. _d -4
99        beta                = 1. _d -11        beta                = 1. _d -11
100          fPrime              = 0. _d 0
101  C-    Earth rotation period is 86400*365.25/366.25 (use to be 1.day)  C-    Earth rotation period is 86400*365.25/366.25 (use to be 1.day)
102        rotationPeriod      = 86164. _d 0        rotationPeriod      = 86164. _d 0
103        viscAh              = 0. _d 3        viscAh              = 0. _d 3
# Line 112  C-    Earth rotation period is 86400*365 Line 124  C-    Earth rotation period is 86400*365
124        diffKrBL79deep      = 0. _d 0        diffKrBL79deep      = 0. _d 0
125        diffKrBL79scl       = 200. _d 0        diffKrBL79scl       = 200. _d 0
126        diffKrBL79Ho        = -2000. _d 0        diffKrBL79Ho        = -2000. _d 0
127          BL79LatVary         = 30
128          diffKrBLEQsurf      = UNSET_RL
129          diffKrBLEQdeep      = UNSET_RL
130          diffKrBLEQscl       = UNSET_RL
131          diffKrBLEQHo        = UNSET_RL
132        viscA4              = 0. _d 11        viscA4              = 0. _d 11
133        viscA4Grid          = 0. _d 0        viscA4Grid          = 0. _d 0
134        viscA4GridMax       = 1. _d 21        viscA4GridMax       = 1. _d 21
# Line 133  Cml      sBeta               = 7.4 _d -4 Line 150  Cml      sBeta               = 7.4 _d -4
150        hFacMinDrDefault    = 0. _d 0        hFacMinDrDefault    = 0. _d 0
151        implicitIntGravWave = .FALSE.        implicitIntGravWave = .FALSE.
152        staggerTimeStep     = .FALSE.        staggerTimeStep     = .FALSE.
153          doResetHFactors     = .FALSE.
154        momViscosity        = .TRUE.        momViscosity        = .TRUE.
155        momAdvection        = .TRUE.        momAdvection        = .TRUE.
156        momForcing          = .TRUE.        momForcing          = .TRUE.
# Line 161  Cml      sBeta               = 7.4 _d -4 Line 179  Cml      sBeta               = 7.4 _d -4
179        globalFiles         = .FALSE.        globalFiles         = .FALSE.
180        useSingleCpuIO      = .FALSE.        useSingleCpuIO      = .FALSE.
181        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
182          allowInteriorFreezing = .FALSE.
183        useOldFreezing      = .FALSE.        useOldFreezing      = .FALSE.
184        ivdc_kappa          = 0. _d 0        ivdc_kappa          = 0. _d 0
185        usePickupBeforeC35    = .FALSE.        hMixCriteria        = -.8 _d 0
186          dRhoSmall           = 1. _d -6
187          hMixSmooth          = 0. _d 0
188        usePickupBeforeC54    = .FALSE.        usePickupBeforeC54    = .FALSE.
       debugMode             = .FALSE.  
189        tempAdvScheme       = 2        tempAdvScheme       = 2
190        saltAdvScheme       = 2        saltAdvScheme       = 2
191        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
192          useMultiDimAdvec    = .FALSE.
193        useCDscheme         = .FALSE.        useCDscheme         = .FALSE.
194        useEnergyConservingCoriolis = .FALSE.        useEnergyConservingCoriolis = .FALSE.
195        useJamartWetPoints  = .FALSE.        useJamartWetPoints  = .FALSE.
196        useJamartMomAdv     = .FALSE.        useJamartMomAdv     = .FALSE.
197        SadournyCoriolis    = .FALSE.        selectVortScheme    = UNSET_I
198        upwindVorticity     = .FALSE.        upwindVorticity     = .FALSE.
199        highOrderVorticity  = .FALSE.        highOrderVorticity  = .FALSE.
200        useAbsVorticity     = .FALSE.        useAbsVorticity     = .FALSE.
201        upwindShear         = .FALSE.        upwindShear         = .FALSE.
202        selectKEscheme      = 0        selectKEscheme      = 0
203        debugLevel          = debLevA        IF ( debugMode ) THEN
204            debugLevel        = debLevD
205          ELSE
206            debugLevel        = debLevB
207    #ifdef ALLOW_AUTODIFF
208            debugLevel        = debLevA
209    #endif
210          ENDIF
211        inAdMode            = .FALSE.        inAdMode            = .FALSE.
212        inAdExact           = .TRUE.        inAdExact           = .TRUE.
213    
# Line 187  C--   Set (free)surface-related paramete Line 215  C--   Set (free)surface-related paramete
215        implicitFreeSurface = .FALSE.        implicitFreeSurface = .FALSE.
216        rigidLid            = .FALSE.        rigidLid            = .FALSE.
217        implicSurfPress     = 1. _d 0        implicSurfPress     = 1. _d 0
218        implicDiv2DFlow     = 1. _d 0        implicDiv2Dflow     = 1. _d 0
219        exactConserv        = .FALSE.        exactConserv        = .FALSE.
220          linFSConserveTr     = .FALSE.
221        uniformLin_PhiSurf  = .TRUE.        uniformLin_PhiSurf  = .TRUE.
222        nonlinFreeSurf      = 0        nonlinFreeSurf      = 0
223        hFacInf             = 0.2 _d 0        hFacInf             = 0.2 _d 0
224        hFacSup             = 2.0 _d 0        hFacSup             = 2.0 _d 0
225        select_rStar        = 0        select_rStar        = 0
226          selectNHfreeSurf    = 0
227          selectAddFluid      = 0
228        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
229        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
230        salt_EvPrRn = 0.        salt_EvPrRn = 0.
231          temp_addMass = UNSET_RL
232          salt_addMass = UNSET_RL
233        balanceEmPmR        = .FALSE.        balanceEmPmR        = .FALSE.
234        balanceQnet         = .FALSE.        balanceQnet         = .FALSE.
235        balancePrintMean    = .FALSE.        balancePrintMean    = .FALSE.
# Line 221  C--   Elliptic solver parameters Line 254  C--   Elliptic solver parameters
254        cg3dMaxIters       = 150        cg3dMaxIters       = 150
255        cg3dTargetResidual = 1. _d -7        cg3dTargetResidual = 1. _d -7
256        cg3dChkResFreq     = 1        cg3dChkResFreq     = 1
257          useSRCGSolver      = .FALSE.
258    
259  C--   Time stepping parameters  C--   Time stepping parameters
260        deltaT            = 0. _d 0        deltaT            = 0. _d 0
# Line 230  C--   Time stepping parameters Line 264  C--   Time stepping parameters
264          dTtracerLev(k)  = 0. _d 0          dTtracerLev(k)  = 0. _d 0
265        ENDDO        ENDDO
266        baseTime          = 0. _d 0        baseTime          = 0. _d 0
267        nIter0            = 0        nIter0            = -1
268        startTime         = deltaT*float(nIter0)        startTime         = UNSET_RL
269        pickupSuff        = ' '        pickupSuff        = ' '
270          pickupStrictlyMatch = .TRUE.
271        nTimeSteps        = 0        nTimeSteps        = 0
272        nEndIter          = nIter0+nTimeSteps        nEndIter          = 0
273        endTime           = deltaT*float(nEndIter)        endTime           = 0. _d 0
274        momForcingOutAB   = UNSET_I        momForcingOutAB   = UNSET_I
275        tracForcingOutAB  = UNSET_I        tracForcingOutAB  = UNSET_I
276        momDissip_In_AB   = .TRUE.        momDissip_In_AB   = .TRUE.
# Line 256  C--   Time stepping parameters Line 291  C--   Time stepping parameters
291        pickup_read_mdsio = .TRUE.        pickup_read_mdsio = .TRUE.
292        pickup_write_mdsio= .TRUE.        pickup_write_mdsio= .TRUE.
293        pickup_write_immed= .FALSE.        pickup_write_immed= .FALSE.
294          writePickupAtEnd  = .TRUE.
295        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
296        adjDumpFreq       = deltaT*0        adjDumpFreq       = deltaT*0
297        diagFreq          = deltaT*0        diagFreq          = deltaT*0
# Line 263  C--   Time stepping parameters Line 299  C--   Time stepping parameters
299        snapshot_mdsio    = .TRUE.        snapshot_mdsio    = .TRUE.
300        monitorFreq       = -1.        monitorFreq       = -1.
301        adjMonitorFreq    = 0.        adjMonitorFreq    = 0.
302          monitorSelect     = UNSET_I
303        monitor_stdio     = .TRUE.        monitor_stdio     = .TRUE.
304        taveFreq          = deltaT*0        taveFreq          = deltaT*0
305        timeave_mdsio     = .TRUE.        timeave_mdsio     = .TRUE.
# Line 270  C--   Time stepping parameters Line 307  C--   Time stepping parameters
307        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
308        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
309        readBinaryPrec    = precFloat32        readBinaryPrec    = precFloat32
       nCheckLev         = 1  
       checkPtSuff(1)    = 'ckptA'  
       checkPtSuff(2)    = 'ckptB'  
310        cAdjFreq          =  0. _d 0        cAdjFreq          =  0. _d 0
311        tauCD             =  0. _d 0        tauCD             =  0. _d 0
312        tauThetaClimRelax =  0. _d 0        tauThetaClimRelax =  0. _d 0
313        tauSaltClimRelax  =  0. _d 0        tauSaltClimRelax  =  0. _d 0
       tauTr1ClimRelax   =  0. _d 0  
314        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
315        externForcingPeriod     = 0.        externForcingPeriod     = 0.
316        externForcingCycle      = 0.        externForcingCycle      = 0.
# Line 287  C--   Time stepping parameters Line 320  C--   Time stepping parameters
320  C--   Input files  C--   Input files
321        bathyFile       = ' '        bathyFile       = ' '
322        topoFile        = ' '        topoFile        = ' '
       shelfIceFile    = ' '  
323        hydrogSaltFile  = ' '        hydrogSaltFile  = ' '
324        hydrogThetaFile = ' '        hydrogThetaFile = ' '
325          maskIniTemp     = .TRUE.
326          maskIniSalt     = .TRUE.
327          checkIniTemp    = .TRUE.
328          checkIniSalt    = .TRUE.
329          diffKrFile      = ' '
330          viscAhDfile     = ' '
331          viscAhZfile     = ' '
332          viscA4Dfile     = ' '
333          viscA4Zfile     = ' '
334        zonalWindFile   = ' '        zonalWindFile   = ' '
335        meridWindFile   = ' '        meridWindFile   = ' '
336        thetaClimFile   = ' '        thetaClimFile   = ' '
# Line 304  C--   Input files Line 345  C--   Input files
345        pSurfInitFile   = ' '        pSurfInitFile   = ' '
346        dQdTFile        = ' '        dQdTFile        = ' '
347        ploadFile       = ' '        ploadFile       = ' '
348        eddyTauxFile    = ' '        eddyPsiXFile    = ' '
349        eddyTauyFile    = ' '        eddyPsiYFile    = ' '
350        lambdaThetaFile = ' '        lambdaThetaFile = ' '
351        lambdaSaltFile  = ' '        lambdaSaltFile  = ' '
352        mdsioLocalDir   = ' '        mdsioLocalDir   = ' '
353          adTapeDir       = ' '
354          the_run_name    = ' '
355    
356        RETURN        RETURN
357        END        END
Legend:
Removed from v.1.116  
changed lines
  Added in v.1.156
  ViewVC Help
Powered by ViewVC 1.1.22