/[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.79 by jmc, Mon Jan 3 23:16:32 2005 UTC revision 1.93 by jmc, Thu Jun 9 15:53:55 2005 UTC
# Line 68  C     General curvilinear coordinate sys Line 68  C     General curvilinear coordinate sys
68        usingCurvilinearGrid = .FALSE.        usingCurvilinearGrid = .FALSE.
69  C     General cylindrical coordinate system  C     General cylindrical coordinate system
70        usingCylindricalGrid = .FALSE.        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
# Line 75  C--   Set default "physical" parameters Line 79  C--   Set default "physical" parameters
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
82          nh_Am2              = 1.D0
83        gravity             = 9.81D0        gravity             = 9.81D0
84        rhoNil              = 999.8D0        rhoNil              = 999.8D0
85  C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)  C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)
# Line 86  C-    Always use 1 day in the past but s Line 91  C-    Always use 1 day in the past but s
91        rotationPeriod      = 86400. _d 0        rotationPeriod      = 86400. _d 0
92        viscAh              = 0.D3        viscAh              = 0.D3
93        viscAhGrid          = 0.D0        viscAhGrid          = 0.D0
94          viscAhGridMin       = 0.D0
95          viscAhGridMax       = 1.D21
96        viscAhMax           = 1.D21        viscAhMax           = 1.D21
97        viscC2leith         = 0.D0        viscC2leith         = 0.D0
98          viscC2leithD        = 0.D0
99          viscC2smag          = 0.D0
100        viscAstrain         = 0.D3        viscAstrain         = 0.D3
101        viscAtension        = 0.D3        viscAtension        = 0.D3
102        diffKhT             = 0.D3        diffKhT             = 0.D3
# Line 103  C-    Always use 1 day in the past but s Line 112  C-    Always use 1 day in the past but s
112        diffKrBL79Ho        = -2000.D0        diffKrBL79Ho        = -2000.D0
113        viscA4              = 0.D11        viscA4              = 0.D11
114        viscA4Grid          = 0.D0        viscA4Grid          = 0.D0
115        viscA4GridMax       = 0.D0        viscA4GridMax       = 1.D21
116        viscA4GridMin       = 0.D0        viscA4GridMin       = 0.D0
117        viscA4Max           = 1.D21        viscA4Max           = 1.D21
118        viscC4leith         = 0.D0        viscC4leith         = 0.D0
119          viscC4leithD        = 0.D0
120        diffK4T             = 0.D11        diffK4T             = 0.D11
121        diffK4S             = 0.D11        diffK4S             = 0.D11
122        cosPower            = 0.        cosPower            = 0.
# Line 133  Cml      sBeta               = 7.4D-4 Line 143  Cml      sBeta               = 7.4D-4
143        saltForcing         = .TRUE.        saltForcing         = .TRUE.
144        metricTerms         = .TRUE.        metricTerms         = .TRUE.
145        useNHMTerms         = .FALSE.        useNHMTerms         = .FALSE.
146          useFullLeith        = .FALSE.
147          useAnisotropicViscAGridMax = .FALSE.
148        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
149        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
150        momImplVertAdv      = .FALSE.        momImplVertAdv      = .FALSE.
# Line 151  Cml      sBeta               = 7.4D-4 Line 163  Cml      sBeta               = 7.4D-4
163        usePickupBeforeC35    = .FALSE.        usePickupBeforeC35    = .FALSE.
164        usePickupBeforeC54    = .FALSE.        usePickupBeforeC54    = .FALSE.
165        debugMode             = .FALSE.        debugMode             = .FALSE.
       readPickupWithTracer  = .FALSE.  
       writePickupWithTracer = .FALSE.  
166        tempAdvScheme       = 2        tempAdvScheme       = 2
167        saltAdvScheme       = 2        saltAdvScheme       = 2
168        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
# Line 164  Cml      sBeta               = 7.4D-4 Line 174  Cml      sBeta               = 7.4D-4
174        upwindVorticity     = .FALSE.        upwindVorticity     = .FALSE.
175        highOrderVorticity  = .FALSE.        highOrderVorticity  = .FALSE.
176        useAbsVorticity     = .FALSE.        useAbsVorticity     = .FALSE.
177          upwindShear         = .FALSE.
178        debugLevel          = debLevA        debugLevel          = debLevA
179          inAdMode            = .FALSE.
180          inAdExact           = .TRUE.
181    
182  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
183        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
# Line 209  C--   Time stepping parameters Line 222  C--   Time stepping parameters
222        DO k=1,Nr        DO k=1,Nr
223          dTtracerLev(k)  = 0. _d 0          dTtracerLev(k)  = 0. _d 0
224        ENDDO        ENDDO
225          baseTime          = 0. _d 0
226        nIter0            = 0        nIter0            = 0
227        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
228        pickupSuff        = ' '        pickupSuff        = ' '
# Line 217  C--   Time stepping parameters Line 231  C--   Time stepping parameters
231        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
232        forcing_In_AB     = .TRUE.        forcing_In_AB     = .TRUE.
233        abEps             = 0.01        abEps             = 0.01
234    #ifdef ALLOW_ADAMSBASHFORTH_3
235          alph_AB           = 0.5 _d 0
236          beta_AB           = 5. _d 0 / 12. _d 0
237          startFromPickupAB2= .FALSE.
238    #else
239          alph_AB           = UNSET_RL
240          beta_AB           = UNSET_RL
241          startFromPickupAB2= .TRUE.
242    #endif
243        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
244        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
245        outputTypesInclusive = .FALSE.        outputTypesInclusive = .FALSE.
# Line 228  C--   Time stepping parameters Line 251  C--   Time stepping parameters
251        diagFreq          = deltaT*0        diagFreq          = deltaT*0
252        snapshot_mdsio    = .TRUE.        snapshot_mdsio    = .TRUE.
253        monitorFreq       = -1.        monitorFreq       = -1.
254          adjMonitorFreq    = 0.
255        monitor_stdio     = .TRUE.        monitor_stdio     = .TRUE.
256        taveFreq          = deltaT*0        taveFreq          = deltaT*0
257        timeave_mdsio     = .TRUE.        timeave_mdsio     = .TRUE.
# Line 246  C--   Time stepping parameters Line 270  C--   Time stepping parameters
270        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
271        externForcingPeriod     = 0.        externForcingPeriod     = 0.
272        externForcingCycle      = 0.        externForcingCycle      = 0.
273        tCyl              = 0.        tCylIn             = 0.
274          tCylOut            = 20.
275          calendarDumps     = .FALSE.
276    
277  C--   Input files  C--   Input files
278        bathyFile       = ' '        bathyFile       = ' '
# Line 258  C--   Input files Line 284  C--   Input files
284        thetaClimFile   = ' '        thetaClimFile   = ' '
285        saltClimFile    = ' '        saltClimFile    = ' '
286        EmPmRfile       = ' '        EmPmRfile       = ' '
287          saltFluxFile    = ' '
288        surfQfile       = ' '        surfQfile       = ' '
289        surfQnetFile    = ' '        surfQnetFile    = ' '
290        surfQswFile     = ' '        surfQswFile     = ' '
# Line 266  C--   Input files Line 293  C--   Input files
293        pSurfInitFile   = ' '        pSurfInitFile   = ' '
294        dQdTFile        = ' '        dQdTFile        = ' '
295        ploadFile       = ' '        ploadFile       = ' '
296          eddyTauxFile    = ' '
297          eddyTauyFile    = ' '
298          lambdaThetaFile = ' '
299          lambdaSaltFile  = ' '
300        mdsioLocalDir   = ' '        mdsioLocalDir   = ' '
301    
302  C  C
Legend:
Removed from v.1.79  
changed lines
  Added in v.1.93
  ViewVC Help
Powered by ViewVC 1.1.22