/[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.61 by adcroft, Wed Jun 2 13:23:55 2004 UTC revision 1.78 by dimitri, Mon Dec 13 17:24:55 2004 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7  CBOP  CBOP
8  C     !ROUTINE: SET_DEFAULTS  C     !ROUTINE: SET_DEFAULTS
9  C     !INTERFACE:  C     !INTERFACE:
# Line 11  C     !INTERFACE: Line 12  C     !INTERFACE:
12       O   hFacMinDrDefault, delRdefault, rkFacDefault,       O   hFacMinDrDefault, delRdefault, rkFacDefault,
13       I   myThid )       I   myThid )
14    
15  C     !DESCRIPTION: \bv  C     !DESCRIPTION:
16  C     *==========================================================*  C     Routine to set model "parameter defaults".
 C     | SUBROUTINE SET_DEFAULTS                                    
 C     | o Routine to set model "parameters"                        
 C     *==========================================================*  
 C     | Notes:                                                      
 C     | ======                                                      
 C     | The present version of this routine is a place-holder.      
 C     | A production version needs to handle parameters from an    
 C     | external file and possibly reading in some initial field    
 C     | values.                                                    
 C     *==========================================================*  
 C     \ev  
17    
18  C     !USES:  C     !USES:
19        IMPLICIT NONE        IMPLICIT NONE
 C     === Global variables ===  
20  #include "SIZE.h"  #include "SIZE.h"
21  #include "EEPARAMS.h"  #include "EEPARAMS.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
# Line 35  Cml#include "EOS.h" Line 24  Cml#include "EOS.h"
24  #include "GRID.h"  #include "GRID.h"
25    
26  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
 C     === Routine arguments ===  
27  C     myThid - Number of this instance of INI_PARMS  C     myThid - Number of this instance of INI_PARMS
28        INTEGER myThid        INTEGER myThid
29        _RL viscArDefault        _RL viscArDefault
# Line 46  C     myThid - Number of this instance o Line 34  C     myThid - Number of this instance o
34        _RS rkFacDefault        _RS rkFacDefault
35    
36  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
 C     === Local variables ===  
37  C     K, I, J - Loop counters  C     K, I, J - Loop counters
38        INTEGER K, I, J        INTEGER K, I, J
39  CEOP  CEOP
# Line 79  C     In spherical polar distances are i Line 66  C     In spherical polar distances are i
66        rSphere              = 6370. * 1.D3        rSphere              = 6370. * 1.D3
67  C     General curvilinear coordinate system  C     General curvilinear coordinate system
68        usingCurvilinearGrid = .FALSE.        usingCurvilinearGrid = .FALSE.
69    C     General cylindrical coordinate system
70          usingCylindricalGrid = .FALSE.
71    
72    
73  C--   Set default "physical" parameters  C--   Set default "physical" parameters
74        DO K =1,Nr        DO K =1,Nr
# Line 107  C-    Always use 1 day in the past but s Line 97  C-    Always use 1 day in the past but s
97        no_slip_bottom      = .TRUE.        no_slip_bottom      = .TRUE.
98        diffKrTDefault      = 0.D-3        diffKrTDefault      = 0.D-3
99        diffKrSDefault      = 0.D-3        diffKrSDefault      = 0.D-3
100          diffKrBL79surf      = 0.D0
101          diffKrBL79deep      = 0.D0
102          diffKrBL79scl       = 200.D0
103          diffKrBL79Ho        = -2000.D0
104        viscA4              = 0.D11        viscA4              = 0.D11
105        viscA4Grid          = 0.D0        viscA4Grid          = 0.D0
106          viscA4GridMax       = 0.D0
107          viscA4GridMin       = 0.D0
108        viscA4Max           = 1.D21        viscA4Max           = 1.D21
109        viscC4leith         = 0.D0        viscC4leith         = 0.D0
110        diffK4T             = 0.D11        diffK4T             = 0.D11
# Line 135  Cml      sBeta               = 7.4D-4 Line 131  Cml      sBeta               = 7.4D-4
131        saltStepping        = .TRUE.        saltStepping        = .TRUE.
132        saltAdvection       = .TRUE.        saltAdvection       = .TRUE.
133        saltForcing         = .TRUE.        saltForcing         = .TRUE.
       tr1Stepping         = .FALSE.  
134        metricTerms         = .TRUE.        metricTerms         = .TRUE.
135        useNHMTerms         = .FALSE.        useNHMTerms         = .FALSE.
136        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
# Line 154  Cml      sBeta               = 7.4D-4 Line 149  Cml      sBeta               = 7.4D-4
149        bottomDragLinear    = 0.        bottomDragLinear    = 0.
150        bottomDragQuadratic = 0.        bottomDragQuadratic = 0.
151        usePickupBeforeC35    = .FALSE.        usePickupBeforeC35    = .FALSE.
152          usePickupBeforeC54    = .FALSE.
153        debugMode             = .FALSE.        debugMode             = .FALSE.
154        readPickupWithTracer  = .FALSE.        readPickupWithTracer  = .FALSE.
155        writePickupWithTracer = .FALSE.        writePickupWithTracer = .FALSE.
156        tempAdvScheme       = 2        tempAdvScheme       = 2
157        saltAdvScheme       = 2        saltAdvScheme       = 2
       tracerAdvScheme     = 2  
158        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
159        useCDscheme         = .FALSE.        useCDscheme         = .FALSE.
160        useEnergyConservingCoriolis = .FALSE.        useEnergyConservingCoriolis = .FALSE.
161        useJamartWetPoints  = .FALSE.        useJamartWetPoints  = .FALSE.
162          useJamartMomAdv     = .FALSE.
163        SadournyCoriolis    = .FALSE.        SadournyCoriolis    = .FALSE.
164        upwindVorticity     = .FALSE.        upwindVorticity     = .FALSE.
165        highOrderVorticity  = .FALSE.        highOrderVorticity  = .FALSE.
166        useAbsVorticity     = .FALSE.        useAbsVorticity     = .FALSE.
167        debugLevel          = debLevA  C hack to allow deblevelzero, until debugmode flag is removed
168          debugLevel          = 999
169    C      debugLevel          = debLevA
170    
171  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
172        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
# Line 184  C--   Set (free)surface-related paramete Line 182  C--   Set (free)surface-related paramete
182        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
183        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
184        salt_EvPrRn = 0.        salt_EvPrRn = 0.
       trac_EvPrRn = UNSET_RL  
185    
186  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
187        celsius2K = 273.16 _d 0        celsius2K = 273.16 _d 0
# Line 202  C--   Elliptic solver parameters Line 199  C--   Elliptic solver parameters
199        cg2dTargetResWunit = -1.        cg2dTargetResWunit = -1.
200        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
201        cg2dpcOffDFac      = 0.51D0        cg2dpcOffDFac      = 0.51D0
202          cg2dPreCondFreq    = 1
203        cg3dMaxIters       = 150        cg3dMaxIters       = 150
204        cg3dTargetResidual = 1.D-7        cg3dTargetResidual = 1.D-7
205        cg3dChkResFreq     = 1        cg3dChkResFreq     = 1
206    
207  C--   Time stepping parameters  C--   Time stepping parameters
208        deltaT            = 0.        deltaT            = 0. _d 0
209          deltaTmom         = 0. _d 0
210          deltaTfreesurf    = 0. _d 0
211          DO k=1,Nr
212            dTtracerLev(k)  = 0. _d 0
213          ENDDO
214        nIter0            = 0        nIter0            = 0
215        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
216        pickupSuff        = ' '        pickupSuff        = ' '
# Line 218  C--   Time stepping parameters Line 221  C--   Time stepping parameters
221        abEps             = 0.01        abEps             = 0.01
222        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
223        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
224          outputTypesInclusive = .FALSE.
225          pickup_read_mdsio = .TRUE.
226          pickup_write_mdsio= .TRUE.
227          pickup_write_immed= .FALSE.
228        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
229        adjDumpFreq       = deltaT*0        adjDumpFreq       = deltaT*0
230        diagFreq          = deltaT*0        diagFreq          = deltaT*0
231          snapshot_mdsio    = .TRUE.
232        monitorFreq       = -1.        monitorFreq       = -1.
233          monitor_stdio     = .TRUE.
234        taveFreq          = deltaT*0        taveFreq          = deltaT*0
235          timeave_mdsio     = .TRUE.
236        tave_lastIter     = 0.5 _d 0        tave_lastIter     = 0.5 _d 0
237        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
238        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
# Line 238  C--   Time stepping parameters Line 248  C--   Time stepping parameters
248        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
249        externForcingPeriod     = 0.        externForcingPeriod     = 0.
250        externForcingCycle      = 0.        externForcingCycle      = 0.
251          tCyl              = 0.
252    
253  C--   Input files  C--   Input files
254        bathyFile       = ' '        bathyFile       = ' '
# Line 250  C--   Input files Line 261  C--   Input files
261        saltClimFile    = ' '        saltClimFile    = ' '
262        EmPmRfile       = ' '        EmPmRfile       = ' '
263        surfQfile       = ' '        surfQfile       = ' '
264        surfQswfile     = ' '        surfQnetFile    = ' '
265          surfQswFile     = ' '
266        uVelInitFile    = ' '        uVelInitFile    = ' '
267        vVelInitFile    = ' '        vVelInitFile    = ' '
268        pSurfInitFile   = ' '        pSurfInitFile   = ' '
Legend:
Removed from v.1.61  
changed lines
  Added in v.1.78
  ViewVC Help
Powered by ViewVC 1.1.22