/[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.42 by adcroft, Thu Nov 7 21:51:15 2002 UTC revision 1.89 by heimbach, Wed Apr 6 20:18:19 2005 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 65  C     Vertical gridding Line 52  C     Vertical gridding
52    
53  C     Horizontal gridding  C     Horizontal gridding
54  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
55        usingCartesianGrid = .TRUE.        usingCartesianGrid = .FALSE.
56        DO i=1,Nx        DO i=1,Nx
57         delX(i) = UNSET_RL         delX(i) = UNSET_RL
58        ENDDO        ENDDO
# 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    C     Coriolis map:
72          useConstantF  = .FALSE.
73          useBetaPlaneF = .FALSE.
74          useSphereF    = .FALSE.
75    
76    
77  C--   Set default "physical" parameters  C--   Set default "physical" parameters
78        DO K =1,Nr        DO K =1,Nr
79         tRef(K) = 30.D0 - FLOAT( K )         tRef(K) = 30.D0 - FLOAT( K )
80  Cml       sRef(K) = 35.D0  Cml       sRef(K) = 35.D0
81        ENDDO        ENDDO
       gravitySign         = 1.D0  
82        gravity             = 9.81D0        gravity             = 9.81D0
83        rhoNil              = 999.8D0        rhoNil              = 999.8D0
84        rhoConstFresh       = 999.8D0  C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)
85    C         (so that the default produces same results as before)
86    c     rhoConstFresh       = 999.8D0
87        f0                  = 1.D-4        f0                  = 1.D-4
88        beta                = 1.D-11        beta                = 1.D-11
89        omega               = 2.D0 * PI / ( 3600.D0 * 24.D0 )  C-    Always use 1 day in the past but should be 86164 (=86400*365.25/366.25)
90          rotationPeriod      = 86400. _d 0
91        viscAh              = 0.D3        viscAh              = 0.D3
92          viscAhGrid          = 0.D0
93          viscAhGridMin       = 0.D0
94          viscAhGridMax       = 1.D21
95          viscAhMax           = 1.D21
96          viscC2leith         = 0.D0
97          viscC2leithD        = 0.D0
98          viscC2smag          = 0.D0
99        viscAstrain         = 0.D3        viscAstrain         = 0.D3
100        viscAtension        = 0.D3        viscAtension        = 0.D3
101        diffKhT             = 0.D3        diffKhT             = 0.D3
# Line 102  Cml       sRef(K) = 35.D0 Line 105  Cml       sRef(K) = 35.D0
105        no_slip_bottom      = .TRUE.        no_slip_bottom      = .TRUE.
106        diffKrTDefault      = 0.D-3        diffKrTDefault      = 0.D-3
107        diffKrSDefault      = 0.D-3        diffKrSDefault      = 0.D-3
108          diffKrBL79surf      = 0.D0
109          diffKrBL79deep      = 0.D0
110          diffKrBL79scl       = 200.D0
111          diffKrBL79Ho        = -2000.D0
112        viscA4              = 0.D11        viscA4              = 0.D11
113          viscA4Grid          = 0.D0
114          viscA4GridMax       = 0.D0
115          viscA4GridMin       = 0.D0
116          viscA4Max           = 1.D21
117          viscC4leith         = 0.D0
118          viscC4leithD        = 0.D0
119        diffK4T             = 0.D11        diffK4T             = 0.D11
120        diffK4S             = 0.D11        diffK4S             = 0.D11
121        cosPower            = 0.        cosPower            = 0.
# Line 127  Cml      sBeta               = 7.4D-4 Line 140  Cml      sBeta               = 7.4D-4
140        saltStepping        = .TRUE.        saltStepping        = .TRUE.
141        saltAdvection       = .TRUE.        saltAdvection       = .TRUE.
142        saltForcing         = .TRUE.        saltForcing         = .TRUE.
       tr1Stepping         = .FALSE.  
143        metricTerms         = .TRUE.        metricTerms         = .TRUE.
144        useNHMTerms         = .TRUE.        useNHMTerms         = .FALSE.
145          useFullLeith        = .FALSE.
146        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
147        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
148          momImplVertAdv      = .FALSE.
149          tempImplVertAdv     = .FALSE.
150          saltImplVertAdv     = .FALSE.
151        nonHydrostatic      = .FALSE.        nonHydrostatic      = .FALSE.
152        quasiHydrostatic    = .FALSE.        quasiHydrostatic    = .FALSE.
153        globalFiles         = .FALSE.        globalFiles         = .FALSE.
154          useSingleCpuIO      = .FALSE.
155        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
156          useOldFreezing      = .FALSE.
157        ivdc_kappa          = 0.D0        ivdc_kappa          = 0.D0
158        groundAtK1          = .FALSE.        groundAtK1          = .FALSE.
159        bottomDragLinear    = 0.        bottomDragLinear    = 0.
160        bottomDragQuadratic = 0.        bottomDragQuadratic = 0.
161        usePickupBeforeC35    = .FALSE.        usePickupBeforeC35    = .FALSE.
162          usePickupBeforeC54    = .FALSE.
163        debugMode             = .FALSE.        debugMode             = .FALSE.
164        readPickupWithTracer  = .FALSE.        readPickupWithTracer  = .FALSE.
165        writePickupWithTracer = .FALSE.        writePickupWithTracer = .FALSE.
166        tempAdvScheme       = 2        tempAdvScheme       = 2
167        saltAdvScheme       = 2        saltAdvScheme       = 2
       tracerAdvScheme     = 2  
168        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
169          useCDscheme         = .FALSE.
170        useEnergyConservingCoriolis = .FALSE.        useEnergyConservingCoriolis = .FALSE.
171        useJamartWetPoints  = .FALSE.        useJamartWetPoints  = .FALSE.
172          useJamartMomAdv     = .FALSE.
173          SadournyCoriolis    = .FALSE.
174          upwindVorticity     = .FALSE.
175          highOrderVorticity  = .FALSE.
176          useAbsVorticity     = .FALSE.
177          debugLevel          = debLevA
178          inAdMode            = .FALSE.
179          inAdExact           = .TRUE.
180    
181  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
182        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
# Line 159  C--   Set (free)surface-related paramete Line 186  C--   Set (free)surface-related paramete
186        exactConserv        = .FALSE.        exactConserv        = .FALSE.
187        uniformLin_PhiSurf  = .TRUE.        uniformLin_PhiSurf  = .TRUE.
188        nonlinFreeSurf      = 0        nonlinFreeSurf      = 0
189        hFacInf             = 1.D0        hFacInf             = 0.2 _d 0
190        hFacSup             = 1.D0        hFacSup             = 2.0 _d 0
191          select_rStar        = 0
192        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
193        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
194        salt_EvPrRn = 0.        salt_EvPrRn = 0.
       trac_EvPrRn = UNSET_RL  
195    
196  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
197        atm_po =  1.D5        celsius2K = 273.16 _d 0
198        atm_cp =  1004.D0        atm_Po =  1. _d 5
199        atm_kappa = 2.D0 / 7.D0        atm_Cp = 1004. _d 0
200        Integr_GeoPot = 2        atm_Rd = UNSET_RL
201          atm_kappa = 2. _d 0 / 7. _d 0
202          atm_Rq = 0. _d 0
203          integr_GeoPot = 2
204          selectFindRoSurf = 0
205    
206  C--   Elliptic solver parameters  C--   Elliptic solver parameters
207        cg2dMaxIters       = 150        cg2dMaxIters       = 150
# Line 178  C--   Elliptic solver parameters Line 209  C--   Elliptic solver parameters
209        cg2dTargetResWunit = -1.        cg2dTargetResWunit = -1.
210        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
211        cg2dpcOffDFac      = 0.51D0        cg2dpcOffDFac      = 0.51D0
212          cg2dPreCondFreq    = 1
213        cg3dMaxIters       = 150        cg3dMaxIters       = 150
214        cg3dTargetResidual = 1.D-7        cg3dTargetResidual = 1.D-7
215        cg3dChkResFreq     = 1        cg3dChkResFreq     = 1
216    
217  C--   Time stepping parameters  C--   Time stepping parameters
218        deltaT            = 0.        deltaT            = 0. _d 0
219          deltaTmom         = 0. _d 0
220          deltaTfreesurf    = 0. _d 0
221          DO k=1,Nr
222            dTtracerLev(k)  = 0. _d 0
223          ENDDO
224          baseTime          = 0. _d 0
225        nIter0            = 0        nIter0            = 0
226        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
227          pickupSuff        = ' '
228        nTimeSteps        = 0        nTimeSteps        = 0
229        nEndIter          = nIter0+nTimeSteps        nEndIter          = nIter0+nTimeSteps
230        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
# Line 193  C--   Time stepping parameters Line 232  C--   Time stepping parameters
232        abEps             = 0.01        abEps             = 0.01
233        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
234        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
235          outputTypesInclusive = .FALSE.
236          pickup_read_mdsio = .TRUE.
237          pickup_write_mdsio= .TRUE.
238          pickup_write_immed= .FALSE.
239        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
240          adjDumpFreq       = deltaT*0
241        diagFreq          = deltaT*0        diagFreq          = deltaT*0
242          snapshot_mdsio    = .TRUE.
243        monitorFreq       = -1.        monitorFreq       = -1.
244          adjMonitorFreq    = 0.
245          monitor_stdio     = .TRUE.
246        taveFreq          = deltaT*0        taveFreq          = deltaT*0
247          timeave_mdsio     = .TRUE.
248        tave_lastIter     = 0.5 _d 0        tave_lastIter     = 0.5 _d 0
249        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
250        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
# Line 212  C--   Time stepping parameters Line 260  C--   Time stepping parameters
260        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
261        externForcingPeriod     = 0.        externForcingPeriod     = 0.
262        externForcingCycle      = 0.        externForcingCycle      = 0.
263          tCylIn             = 0.
264          tCylOut            = 20.
265          calendarDumps     = .FALSE.
266    
267  C--   Input files  C--   Input files
268        bathyFile       = ' '        bathyFile       = ' '
# Line 223  C--   Input files Line 274  C--   Input files
274        thetaClimFile   = ' '        thetaClimFile   = ' '
275        saltClimFile    = ' '        saltClimFile    = ' '
276        EmPmRfile       = ' '        EmPmRfile       = ' '
277          saltFluxFile    = ' '
278        surfQfile       = ' '        surfQfile       = ' '
279        surfQswfile     = ' '        surfQnetFile    = ' '
280          surfQswFile     = ' '
281        uVelInitFile    = ' '        uVelInitFile    = ' '
282        vVelInitFile    = ' '        vVelInitFile    = ' '
283        pSurfInitFile   = ' '        pSurfInitFile   = ' '
284        dQdTFile        = ' '        dQdTFile        = ' '
285        ploadFile       = ' '        ploadFile       = ' '
286          eddyTauxFile    = ' '
287          eddyTauyFile    = ' '
288          lambdaThetaFile = ' '
289          lambdaSaltFile  = ' '
290          mdsioLocalDir   = ' '
291    
292  C  C
293        RETURN        RETURN
Legend:
Removed from v.1.42  
changed lines
  Added in v.1.89
  ViewVC Help
Powered by ViewVC 1.1.22