/[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.70 by edhill, Fri Sep 10 12:19:30 2004 UTC revision 1.103 by baylor, Thu Sep 22 14:00:59 2005 UTC
# Line 51  C     Vertical gridding Line 51  C     Vertical gridding
51        ENDDO        ENDDO
52    
53  C     Horizontal gridding  C     Horizontal gridding
54          delXFile        = ' '
55          delYFile        = ' '
56          horizGridFile   = ' '
57  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
58        usingCartesianGrid = .FALSE.        usingCartesianGrid = .FALSE.
59        DO i=1,Nx        DO i=1,Nx
# Line 68  C     General curvilinear coordinate sys Line 71  C     General curvilinear coordinate sys
71        usingCurvilinearGrid = .FALSE.        usingCurvilinearGrid = .FALSE.
72  C     General cylindrical coordinate system  C     General cylindrical coordinate system
73        usingCylindricalGrid = .FALSE.        usingCylindricalGrid = .FALSE.
74    C     Coriolis map:
75          useConstantF  = .FALSE.
76          useBetaPlaneF = .FALSE.
77          useSphereF    = .FALSE.
78    
79    
80  C--   Set default "physical" parameters  C--   Set default "physical" parameters
# Line 75  C--   Set default "physical" parameters Line 82  C--   Set default "physical" parameters
82         tRef(K) = 30.D0 - FLOAT( K )         tRef(K) = 30.D0 - FLOAT( K )
83  Cml       sRef(K) = 35.D0  Cml       sRef(K) = 35.D0
84        ENDDO        ENDDO
85          nh_Am2              = 1.D0
86        gravity             = 9.81D0        gravity             = 9.81D0
87        rhoNil              = 999.8D0        rhoNil              = 999.8D0
88  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 94  C-    Always use 1 day in the past but s
94        rotationPeriod      = 86400. _d 0        rotationPeriod      = 86400. _d 0
95        viscAh              = 0.D3        viscAh              = 0.D3
96        viscAhGrid          = 0.D0        viscAhGrid          = 0.D0
97          viscAhGridMin       = 0.D0
98          viscAhGridMax       = 1.D21
99        viscAhMax           = 1.D21        viscAhMax           = 1.D21
100          viscAhReMax         = 0.D0
101        viscC2leith         = 0.D0        viscC2leith         = 0.D0
102        viscAstrain         = 0.D3        viscC2leithD        = 0.D0
103        viscAtension        = 0.D3        viscC2smag          = 0.D0
104        diffKhT             = 0.D3        diffKhT             = 0.D3
105        diffKhS             = 0.D3        diffKhS             = 0.D3
106        viscArDefault       = 0.D-3        viscArDefault       = 0.D-3
# Line 103  C-    Always use 1 day in the past but s Line 114  C-    Always use 1 day in the past but s
114        diffKrBL79Ho        = -2000.D0        diffKrBL79Ho        = -2000.D0
115        viscA4              = 0.D11        viscA4              = 0.D11
116        viscA4Grid          = 0.D0        viscA4Grid          = 0.D0
117          viscA4GridMax       = 1.D21
118          viscA4GridMin       = 0.D0
119        viscA4Max           = 1.D21        viscA4Max           = 1.D21
120          viscA4ReMax         = 0.D0
121        viscC4leith         = 0.D0        viscC4leith         = 0.D0
122          viscC4leithD        = 0.D0
123          viscC4smag          = 0.D0
124        diffK4T             = 0.D11        diffK4T             = 0.D11
125        diffK4S             = 0.D11        diffK4S             = 0.D11
126        cosPower            = 0.        cosPower            = 0.
# Line 129  Cml      sBeta               = 7.4D-4 Line 145  Cml      sBeta               = 7.4D-4
145        saltStepping        = .TRUE.        saltStepping        = .TRUE.
146        saltAdvection       = .TRUE.        saltAdvection       = .TRUE.
147        saltForcing         = .TRUE.        saltForcing         = .TRUE.
       tr1Stepping         = .FALSE.  
148        metricTerms         = .TRUE.        metricTerms         = .TRUE.
149        useNHMTerms         = .FALSE.        useNHMTerms         = .FALSE.
150          useFullLeith        = .FALSE.
151          useAreaViscLength   = .FALSE.
152          useStrainTensionVisc = .FALSE.
153        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
154        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
155        momImplVertAdv      = .FALSE.        momImplVertAdv      = .FALSE.
# Line 144  Cml      sBeta               = 7.4D-4 Line 162  Cml      sBeta               = 7.4D-4
162        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
163        useOldFreezing      = .FALSE.        useOldFreezing      = .FALSE.
164        ivdc_kappa          = 0.D0        ivdc_kappa          = 0.D0
       groundAtK1          = .FALSE.  
165        bottomDragLinear    = 0.        bottomDragLinear    = 0.
166        bottomDragQuadratic = 0.        bottomDragQuadratic = 0.
167        usePickupBeforeC35    = .FALSE.        usePickupBeforeC35    = .FALSE.
168        usePickupBeforeC54    = .FALSE.        usePickupBeforeC54    = .FALSE.
169        debugMode             = .FALSE.        debugMode             = .FALSE.
       readPickupWithTracer  = .FALSE.  
       writePickupWithTracer = .FALSE.  
170        tempAdvScheme       = 2        tempAdvScheme       = 2
171        saltAdvScheme       = 2        saltAdvScheme       = 2
       tracerAdvScheme     = 2  
172        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
173        useCDscheme         = .FALSE.        useCDscheme         = .FALSE.
174        useEnergyConservingCoriolis = .FALSE.        useEnergyConservingCoriolis = .FALSE.
# Line 164  Cml      sBeta               = 7.4D-4 Line 178  Cml      sBeta               = 7.4D-4
178        upwindVorticity     = .FALSE.        upwindVorticity     = .FALSE.
179        highOrderVorticity  = .FALSE.        highOrderVorticity  = .FALSE.
180        useAbsVorticity     = .FALSE.        useAbsVorticity     = .FALSE.
181          upwindShear         = .FALSE.
182        debugLevel          = debLevA        debugLevel          = debLevA
183          inAdMode            = .FALSE.
184          inAdExact           = .TRUE.
185    
186  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
187        implicitFreeSurface = .TRUE.        implicitFreeSurface = .FALSE.
188        rigidLid            = .FALSE.        rigidLid            = .FALSE.
189        implicSurfPress     = 1.D0        implicSurfPress     = 1.D0
190        implicDiv2DFlow     = 1.D0        implicDiv2DFlow     = 1.D0
# Line 180  C--   Set (free)surface-related paramete Line 197  C--   Set (free)surface-related paramete
197        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
198        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
199        salt_EvPrRn = 0.        salt_EvPrRn = 0.
       trac_EvPrRn = UNSET_RL  
200    
201  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
202        celsius2K = 273.16 _d 0        celsius2K = 273.16 _d 0
# Line 204  C--   Elliptic solver parameters Line 220  C--   Elliptic solver parameters
220        cg3dChkResFreq     = 1        cg3dChkResFreq     = 1
221    
222  C--   Time stepping parameters  C--   Time stepping parameters
223        deltaT            = 0.        deltaT            = 0. _d 0
224          deltaTmom         = 0. _d 0
225          deltaTfreesurf    = 0. _d 0
226          DO k=1,Nr
227            dTtracerLev(k)  = 0. _d 0
228          ENDDO
229          baseTime          = 0. _d 0
230        nIter0            = 0        nIter0            = 0
231        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
232        pickupSuff        = ' '        pickupSuff        = ' '
# Line 213  C--   Time stepping parameters Line 235  C--   Time stepping parameters
235        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
236        forcing_In_AB     = .TRUE.        forcing_In_AB     = .TRUE.
237        abEps             = 0.01        abEps             = 0.01
238    #ifdef ALLOW_ADAMSBASHFORTH_3
239          alph_AB           = 0.5 _d 0
240          beta_AB           = 5. _d 0 / 12. _d 0
241          startFromPickupAB2= .FALSE.
242    #else
243          alph_AB           = UNSET_RL
244          beta_AB           = UNSET_RL
245          startFromPickupAB2= .TRUE.
246    #endif
247        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
248        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
249          outputTypesInclusive = .FALSE.
250        pickup_read_mdsio = .TRUE.        pickup_read_mdsio = .TRUE.
251        pickup_write_mdsio= .TRUE.        pickup_write_mdsio= .TRUE.
252        pickup_write_immed= .FALSE.        pickup_write_immed= .FALSE.
253        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
254        adjDumpFreq       = deltaT*0        adjDumpFreq       = deltaT*0
255        diagFreq          = deltaT*0        diagFreq          = deltaT*0
256          dumpInitAndLast   = .TRUE.
257        snapshot_mdsio    = .TRUE.        snapshot_mdsio    = .TRUE.
258        monitorFreq       = -1.        monitorFreq       = -1.
259        monitor_mdsio     = .TRUE.        adjMonitorFreq    = 0.
260          monitor_stdio     = .TRUE.
261        taveFreq          = deltaT*0        taveFreq          = deltaT*0
262        timeave_mdsio     = .TRUE.        timeave_mdsio     = .TRUE.
263        tave_lastIter     = 0.5 _d 0        tave_lastIter     = 0.5 _d 0
# Line 241  C--   Time stepping parameters Line 275  C--   Time stepping parameters
275        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
276        externForcingPeriod     = 0.        externForcingPeriod     = 0.
277        externForcingCycle      = 0.        externForcingCycle      = 0.
278        tCyl              = 0.        tCylIn             = 0.
279          tCylOut            = 20.
280          calendarDumps     = .FALSE.
281    
282  C--   Input files  C--   Input files
283        bathyFile       = ' '        bathyFile       = ' '
# Line 253  C--   Input files Line 289  C--   Input files
289        thetaClimFile   = ' '        thetaClimFile   = ' '
290        saltClimFile    = ' '        saltClimFile    = ' '
291        EmPmRfile       = ' '        EmPmRfile       = ' '
292          saltFluxFile    = ' '
293        surfQfile       = ' '        surfQfile       = ' '
294        surfQnetFile    = ' '        surfQnetFile    = ' '
295        surfQswFile     = ' '        surfQswFile     = ' '
# Line 261  C--   Input files Line 298  C--   Input files
298        pSurfInitFile   = ' '        pSurfInitFile   = ' '
299        dQdTFile        = ' '        dQdTFile        = ' '
300        ploadFile       = ' '        ploadFile       = ' '
301          eddyTauxFile    = ' '
302          eddyTauyFile    = ' '
303          lambdaThetaFile = ' '
304          lambdaSaltFile  = ' '
305        mdsioLocalDir   = ' '        mdsioLocalDir   = ' '
306    
307  C  C
Legend:
Removed from v.1.70  
changed lines
  Added in v.1.103
  ViewVC Help
Powered by ViewVC 1.1.22