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

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

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


Revision 1.79 - (hide annotations) (download)
Mon Jan 3 23:16:32 2005 UTC (19 years, 6 months ago) by jmc
Branch: MAIN
Changes since 1.78: +2 -4 lines
debugMode & debugLevel: go back to what it was before.
  version 1.137 (ini_parms.F) & 1.77 (set_defaults.F)

1 jmc 1.79 C $Header: /u/gcmpack/MITgcm/model/src/set_defaults.F,v 1.78 2004/12/13 17:24:55 dimitri Exp $
2 heimbach 1.18 C $Name: $
3 adcroft 1.1
4     #include "CPP_OPTIONS.h"
5    
6 edhill 1.70 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7 cnh 1.29 CBOP
8     C !ROUTINE: SET_DEFAULTS
9     C !INTERFACE:
10 adcroft 1.1 SUBROUTINE SET_DEFAULTS(
11     O viscArDefault, diffKrTDefault, diffKrSDefault,
12 adcroft 1.2 O hFacMinDrDefault, delRdefault, rkFacDefault,
13 adcroft 1.1 I myThid )
14 cnh 1.29
15 edhill 1.70 C !DESCRIPTION:
16     C Routine to set model "parameter defaults".
17 cnh 1.29
18     C !USES:
19 adcroft 1.1 IMPLICIT NONE
20     #include "SIZE.h"
21     #include "EEPARAMS.h"
22     #include "PARAMS.h"
23 mlosch 1.37 Cml#include "EOS.h"
24 adcroft 1.1 #include "GRID.h"
25    
26 cnh 1.29 C !INPUT/OUTPUT PARAMETERS:
27 adcroft 1.1 C myThid - Number of this instance of INI_PARMS
28     INTEGER myThid
29     _RL viscArDefault
30     _RL diffKrTDefault
31     _RL diffKrSDefault
32     _RL hFacMinDrDefault
33     _RL delRDefault(Nr)
34 adcroft 1.2 _RS rkFacDefault
35 adcroft 1.1
36 cnh 1.29 C !LOCAL VARIABLES:
37 adcroft 1.1 C K, I, J - Loop counters
38     INTEGER K, I, J
39 cnh 1.29 CEOP
40 adcroft 1.1
41     C-- Grid parameters
42     C Vertical gridding
43 adcroft 1.2 rkFacDefault = 1.D0
44     horiVertRatio = 1.D0
45 adcroft 1.12 Ro_SeaLevel = 0.
46 adcroft 1.1 DO k=1,Nr
47 adcroft 1.2 delRdefault(k) = 0.
48 jmc 1.34 ENDDO
49     DO k=1,Nr+1
50     delRc(k) = UNSET_RL
51 adcroft 1.1 ENDDO
52    
53     C Horizontal gridding
54     C In cartesian coords distances are in metres
55 jmc 1.56 usingCartesianGrid = .FALSE.
56 adcroft 1.1 DO i=1,Nx
57 adcroft 1.8 delX(i) = UNSET_RL
58 adcroft 1.1 ENDDO
59     DO j=1,Ny
60 adcroft 1.8 delY(j) = UNSET_RL
61 adcroft 1.1 ENDDO
62     C In spherical polar distances are in degrees
63     usingSphericalPolarGrid = .FALSE.
64 adcroft 1.2 phiMin = 0.0
65     thetaMin = 0.
66     rSphere = 6370. * 1.D3
67 adcroft 1.19 C General curvilinear coordinate system
68     usingCurvilinearGrid = .FALSE.
69 afe 1.64 C General cylindrical coordinate system
70     usingCylindricalGrid = .FALSE.
71    
72 adcroft 1.1
73     C-- Set default "physical" parameters
74     DO K =1,Nr
75     tRef(K) = 30.D0 - FLOAT( K )
76 mlosch 1.39 Cml sRef(K) = 35.D0
77 adcroft 1.1 ENDDO
78 adcroft 1.2 gravity = 9.81D0
79     rhoNil = 999.8D0
80 jmc 1.55 C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)
81     C (so that the default produces same results as before)
82     c rhoConstFresh = 999.8D0
83 adcroft 1.2 f0 = 1.D-4
84     beta = 1.D-11
85 jmc 1.46 C- Always use 1 day in the past but should be 86164 (=86400*365.25/366.25)
86     rotationPeriod = 86400. _d 0
87 adcroft 1.2 viscAh = 0.D3
88 adcroft 1.58 viscAhGrid = 0.D0
89     viscAhMax = 1.D21
90 adcroft 1.60 viscC2leith = 0.D0
91 adcroft 1.27 viscAstrain = 0.D3
92     viscAtension = 0.D3
93 adcroft 1.2 diffKhT = 0.D3
94     diffKhS = 0.D3
95     viscArDefault = 0.D-3
96     no_slip_sides = .TRUE.
97     no_slip_bottom = .TRUE.
98     diffKrTDefault = 0.D-3
99     diffKrSDefault = 0.D-3
100 adcroft 1.62 diffKrBL79surf = 0.D0
101     diffKrBL79deep = 0.D0
102     diffKrBL79scl = 200.D0
103     diffKrBL79Ho = -2000.D0
104 adcroft 1.2 viscA4 = 0.D11
105 adcroft 1.58 viscA4Grid = 0.D0
106 dimitri 1.75 viscA4GridMax = 0.D0
107     viscA4GridMin = 0.D0
108 adcroft 1.58 viscA4Max = 1.D21
109 adcroft 1.60 viscC4leith = 0.D0
110 adcroft 1.2 diffK4T = 0.D11
111     diffK4S = 0.D11
112 adcroft 1.4 cosPower = 0.
113 adcroft 1.2 HeatCapacity_Cp = 3994.D0
114 mlosch 1.37 Cml tAlpha = 2.D-4
115     Cml sBeta = 7.4D-4
116 adcroft 1.2 eosType = 'LINEAR'
117 adcroft 1.1 buoyancyRelation = 'OCEANIC'
118     hFacMin = 1.D0
119     hFacMinDrDefault = 0.D0
120 adcroft 1.12 staggerTimeStep = .FALSE.
121 adcroft 1.1 momViscosity = .TRUE.
122     momAdvection = .TRUE.
123     momForcing = .TRUE.
124     useCoriolis = .TRUE.
125     momPressureForcing = .TRUE.
126     momStepping = .TRUE.
127 adcroft 1.25 vectorInvariantMomentum = .FALSE.
128 adcroft 1.1 tempStepping = .TRUE.
129 jmc 1.35 tempAdvection = .TRUE.
130     tempForcing = .TRUE.
131 adcroft 1.1 saltStepping = .TRUE.
132 jmc 1.35 saltAdvection = .TRUE.
133     saltForcing = .TRUE.
134 jmc 1.16 metricTerms = .TRUE.
135 jmc 1.53 useNHMTerms = .FALSE.
136 adcroft 1.1 implicitDiffusion = .FALSE.
137 adcroft 1.5 implicitViscosity = .FALSE.
138 jmc 1.57 momImplVertAdv = .FALSE.
139     tempImplVertAdv = .FALSE.
140     saltImplVertAdv = .FALSE.
141 adcroft 1.2 nonHydrostatic = .FALSE.
142 adcroft 1.42 quasiHydrostatic = .FALSE.
143 adcroft 1.3 globalFiles = .FALSE.
144 dimitri 1.45 useSingleCpuIO = .FALSE.
145 adcroft 1.3 allowFreezing = .FALSE.
146 jmc 1.54 useOldFreezing = .FALSE.
147 adcroft 1.6 ivdc_kappa = 0.D0
148 adcroft 1.12 groundAtK1 = .FALSE.
149     bottomDragLinear = 0.
150     bottomDragQuadratic = 0.
151 heimbach 1.23 usePickupBeforeC35 = .FALSE.
152 jmc 1.65 usePickupBeforeC54 = .FALSE.
153 heimbach 1.23 debugMode = .FALSE.
154     readPickupWithTracer = .FALSE.
155     writePickupWithTracer = .FALSE.
156 adcroft 1.24 tempAdvScheme = 2
157     saltAdvScheme = 2
158 adcroft 1.28 multiDimAdvection = .TRUE.
159 jmc 1.47 useCDscheme = .FALSE.
160 adcroft 1.30 useEnergyConservingCoriolis = .FALSE.
161 adcroft 1.31 useJamartWetPoints = .FALSE.
162 adcroft 1.63 useJamartMomAdv = .FALSE.
163 adcroft 1.61 SadournyCoriolis = .FALSE.
164     upwindVorticity = .FALSE.
165     highOrderVorticity = .FALSE.
166     useAbsVorticity = .FALSE.
167 jmc 1.79 debugLevel = debLevA
168 jmc 1.33
169     C-- Set (free)surface-related parameters
170     implicitFreeSurface = .TRUE.
171     rigidLid = .FALSE.
172     implicSurfPress = 1.D0
173     implicDiv2DFlow = 1.D0
174     exactConserv = .FALSE.
175     uniformLin_PhiSurf = .TRUE.
176     nonlinFreeSurf = 0
177 jmc 1.46 hFacInf = 0.2 _d 0
178     hFacSup = 2.0 _d 0
179 jmc 1.44 select_rStar = 0
180 jmc 1.33 useRealFreshWaterFlux = .FALSE.
181     temp_EvPrRn = UNSET_RL
182     salt_EvPrRn = 0.
183 adcroft 1.1
184 jmc 1.22 C-- Atmospheric physical parameters (e.g.: EOS)
185 jmc 1.49 celsius2K = 273.16 _d 0
186     atm_Po = 1. _d 5
187     atm_Cp = 1004. _d 0
188     atm_Rd = UNSET_RL
189     atm_kappa = 2. _d 0 / 7. _d 0
190 jmc 1.51 atm_Rq = 0. _d 0
191 jmc 1.43 integr_GeoPot = 2
192     selectFindRoSurf = 0
193 jmc 1.22
194 adcroft 1.1 C-- Elliptic solver parameters
195     cg2dMaxIters = 150
196     cg2dTargetResidual = 1.D-7
197 adcroft 1.19 cg2dTargetResWunit = -1.
198 adcroft 1.1 cg2dChkResFreq = 1
199 adcroft 1.2 cg2dpcOffDFac = 0.51D0
200 jmc 1.66 cg2dPreCondFreq = 1
201 adcroft 1.1 cg3dMaxIters = 150
202     cg3dTargetResidual = 1.D-7
203     cg3dChkResFreq = 1
204    
205     C-- Time stepping parameters
206 jmc 1.77 deltaT = 0. _d 0
207     deltaTmom = 0. _d 0
208     deltaTfreesurf = 0. _d 0
209     DO k=1,Nr
210     dTtracerLev(k) = 0. _d 0
211     ENDDO
212 adcroft 1.1 nIter0 = 0
213     startTime = deltaT*float(nIter0)
214 adcroft 1.48 pickupSuff = ' '
215 adcroft 1.2 nTimeSteps = 0
216 adcroft 1.6 nEndIter = nIter0+nTimeSteps
217     endTime = deltaT*float(nEndIter)
218 jmc 1.35 forcing_In_AB = .TRUE.
219 adcroft 1.1 abEps = 0.01
220     pchkPtFreq = deltaT*0
221     chkPtFreq = deltaT*0
222 edhill 1.74 outputTypesInclusive = .FALSE.
223 edhill 1.70 pickup_read_mdsio = .TRUE.
224     pickup_write_mdsio= .TRUE.
225     pickup_write_immed= .FALSE.
226 adcroft 1.2 dumpFreq = deltaT*0
227 heimbach 1.59 adjDumpFreq = deltaT*0
228 adcroft 1.26 diagFreq = deltaT*0
229 edhill 1.70 snapshot_mdsio = .TRUE.
230 adcroft 1.20 monitorFreq = -1.
231 edhill 1.71 monitor_stdio = .TRUE.
232 adcroft 1.1 taveFreq = deltaT*0
233 edhill 1.70 timeave_mdsio = .TRUE.
234 jmc 1.40 tave_lastIter = 0.5 _d 0
235 adcroft 1.1 writeStatePrec = precFloat64
236     writeBinaryPrec = precFloat32
237     readBinaryPrec = precFloat32
238     nCheckLev = 1
239     checkPtSuff(1) = 'ckptA'
240     checkPtSuff(2) = 'ckptB'
241 adcroft 1.6 cAdjFreq = 0.D0
242 adcroft 1.1 tauCD = 0.D0
243     tauThetaClimRelax = 0.D0
244     tauSaltClimRelax = 0.D0
245 heimbach 1.23 tauTr1ClimRelax = 0.D0
246 adcroft 1.1 periodicExternalForcing = .FALSE.
247     externForcingPeriod = 0.
248     externForcingCycle = 0.
249 afe 1.64 tCyl = 0.
250    
251 adcroft 1.1 C-- Input files
252     bathyFile = ' '
253 jmc 1.22 topoFile = ' '
254 adcroft 1.1 hydrogSaltFile = ' '
255     hydrogThetaFile = ' '
256     zonalWindFile = ' '
257     meridWindFile = ' '
258     thetaClimFile = ' '
259     saltClimFile = ' '
260 adcroft 1.2 EmPmRfile = ' '
261     surfQfile = ' '
262 jmc 1.67 surfQnetFile = ' '
263     surfQswFile = ' '
264 adcroft 1.7 uVelInitFile = ' '
265     vVelInitFile = ' '
266 adcroft 1.9 pSurfInitFile = ' '
267 heimbach 1.17 dQdTFile = ' '
268 mlosch 1.36 ploadFile = ' '
269 heimbach 1.50 mdsioLocalDir = ' '
270 adcroft 1.1
271     C
272     RETURN
273     END

  ViewVC Help
Powered by ViewVC 1.1.22