model/src/set_defaults.F

C $Header: /u/gcmpack/MITgcm/model/src/set_defaults.F,v 1.153 2011/04/20 13:45:14 dimitri Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP
C     !ROUTINE: SET_DEFAULTS
C     !INTERFACE:
      SUBROUTINE SET_DEFAULTS(
     O   viscArDefault, diffKrTDefault, diffKrSDefault,
     O   hFacMinDrDefault, delRdefault,
     I   myThid )

C     !DESCRIPTION:
C     Routine to set model "parameter defaults".

C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
Cml#include "EOS.h"
c#include "GRID.h"

C     !INPUT/OUTPUT PARAMETERS:
C     myThid :: my Thread Id. Number
      INTEGER myThid
      _RL viscArDefault
      _RL diffKrTDefault
      _RL diffKrSDefault
      _RL hFacMinDrDefault
      _RL delRDefault(Nr)

C     !LOCAL VARIABLES:
C     i, j, k :: Loop counters
      INTEGER i, j, k
CEOP

C--   Grid parameters
C-    Vertical gridding
      delRFile            = ' '
      delRcFile           = ' '
      hybSigmFile         = ' '
      Ro_SeaLevel         = 0.
      rSigmaBnd           = UNSET_RL
      selectSigmaCoord    = 0
      DO k=1,Nr
       delRdefault(k)     = 0.
      ENDDO
      DO k=1,Nr+1
       delRc(k)           = UNSET_RL
      ENDDO
C-    vertical profile
      tRefFile            = ' '
      sRefFile            = ' '
      rhoRefFile          = ' '

C-    Horizontal gridding
      delXFile            = ' '
      delYFile            = ' '
      horizGridFile       = ' '
      deepAtmosphere      = .FALSE.
      xgOrigin            = UNSET_RL
      ygOrigin            = UNSET_RL
      DO i=1,Nx
       delX(i)            = UNSET_RL
      ENDDO
      DO j=1,Ny
       delY(j)            = UNSET_RL
      ENDDO
C     In cartesian coords distances are in metres
      usingCartesianGrid  = .FALSE.
C     In spherical polar distances are in degrees
      usingSphericalPolarGrid = .FALSE.
      rSphere             = UNSET_RL
C     General curvilinear coordinate system
      usingCurvilinearGrid= .FALSE.
      radius_fromHorizGrid= UNSET_RL
C     General cylindrical coordinate system
      usingCylindricalGrid= .FALSE.
C     Coriolis map:
      selectCoriMap       = -1
      use3dCoriolis       = .TRUE.
C     grid rotation
      rotateGrid          = .FALSE.
      phiEuler            = 0. _d 0
      thetaEuler          = 0. _d 0
      psiEuler            = 0. _d 0

C--   Set default "physical" parameters
      nh_Am2              = 1. _d 0
      gravity             = 9.81 _d 0
      rhoNil              = 999.8 _d 0
C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)
C         (so that the default produces same results as before)
c     rhoConstFresh       = 999.8 _d 0
      f0                  = 1. _d -4
      beta                = 1. _d -11
      fPrime              = 0. _d 0
C-    Earth rotation period is 86400*365.25/366.25 (use to be 1.day)
      rotationPeriod      = 86164. _d 0
      viscAh              = 0. _d 3
      viscAhGrid          = 0. _d 0
      viscAhGridMin       = 0. _d 0
      viscAhGridMax       = 1. _d 21
      viscAhMax           = 1. _d 21
      viscAhReMax         = 0. _d 0
      viscC2leith         = 0. _d 0
      viscC2leithD        = 0. _d 0
      viscC2smag          = 0. _d 0
      diffKhT             = 0. _d 3
      diffKhS             = 0. _d 3
      viscArDefault       = 0. _d -3
      no_slip_sides       = .TRUE.
      no_slip_bottom      = .TRUE.
      sideDragFactor      = 2. _d 0
      bottomDragLinear    = 0.
      bottomDragQuadratic = 0.
      smoothAbsFuncRange  = 0. _d 0
      diffKrTDefault      = 0. _d -3
      diffKrSDefault      = 0. _d -3
      diffKrBL79surf      = 0. _d 0
      diffKrBL79deep      = 0. _d 0
      diffKrBL79scl       = 200. _d 0
      diffKrBL79Ho        = -2000. _d 0
      BL79LatVary         = 30
      diffKrBLEQsurf      = UNSET_RL
      diffKrBLEQdeep      = UNSET_RL
      diffKrBLEQscl       = UNSET_RL
      diffKrBLEQHo        = UNSET_RL
      viscA4              = 0. _d 11
      viscA4Grid          = 0. _d 0
      viscA4GridMax       = 1. _d 21
      viscA4GridMin       = 0. _d 0
      viscA4Max           = 1. _d 21
      viscA4ReMax         = 0. _d 0
      viscC4leith         = 0. _d 0
      viscC4leithD        = 0. _d 0
      viscC4smag          = 0. _d 0
      diffK4T             = 0. _d 11
      diffK4S             = 0. _d 11
      cosPower            = 0.
      HeatCapacity_Cp     = 3994. _d 0
Cml      tAlpha              = 2. _d -4
Cml      sBeta               = 7.4 _d -4
      eosType             = 'LINEAR'
      buoyancyRelation    = 'OCEANIC'
      hFacMin             = 1. _d 0
      hFacMinDrDefault    = 0. _d 0
      implicitIntGravWave = .FALSE.
      staggerTimeStep     = .FALSE.
      momViscosity        = .TRUE.
      momAdvection        = .TRUE.
      momForcing          = .TRUE.
      useCoriolis         = .TRUE.
      momPressureForcing  = .TRUE.
      momStepping         = .TRUE.
      vectorInvariantMomentum = .FALSE.
      tempStepping        = .TRUE.
      tempAdvection       = .TRUE.
      tempForcing         = .TRUE.
      saltStepping        = .TRUE.
      saltAdvection       = .TRUE.
      saltForcing         = .TRUE.
      metricTerms         = .TRUE.
      useNHMTerms         = .FALSE.
      useFullLeith        = .FALSE.
      useAreaViscLength   = .FALSE.
      useStrainTensionVisc= .FALSE.
      implicitDiffusion   = .FALSE.
      implicitViscosity   = .FALSE.
      momImplVertAdv      = .FALSE.
      tempImplVertAdv     = .FALSE.
      saltImplVertAdv     = .FALSE.
      nonHydrostatic      = .FALSE.
      quasiHydrostatic    = .FALSE.
      globalFiles         = .FALSE.
      useSingleCpuIO      = .FALSE.
      allowFreezing       = .FALSE.
      allowInteriorFreezing = .FALSE.
      useOldFreezing      = .FALSE.
      ivdc_kappa          = 0. _d 0
      hMixCriteria        = -.8 _d 0
      dRhoSmall           = 1. _d -6
      hMixSmooth          = 0. _d 0
      usePickupBeforeC54    = .FALSE.
      debugMode             = .FALSE.
      tempAdvScheme       = 2
      saltAdvScheme       = 2
      multiDimAdvection   = .TRUE.
      useMultiDimAdvec    = .FALSE.
      useCDscheme         = .FALSE.
      useEnergyConservingCoriolis = .FALSE.
      useJamartWetPoints  = .FALSE.
      useJamartMomAdv     = .FALSE.
      selectVortScheme    = UNSET_I
      upwindVorticity     = .FALSE.
      highOrderVorticity  = .FALSE.
      useAbsVorticity     = .FALSE.
      upwindShear         = .FALSE.
      selectKEscheme      = 0
      debugLevel          = debLevA
      inAdMode            = .FALSE.
      inAdExact           = .TRUE.

C--   Set (free)surface-related parameters
      implicitFreeSurface = .FALSE.
      rigidLid            = .FALSE.
      implicSurfPress     = 1. _d 0
      implicDiv2Dflow     = 1. _d 0
      exactConserv        = .FALSE.
      linFSConserveTr     = .FALSE.
      uniformLin_PhiSurf  = .TRUE.
      nonlinFreeSurf      = 0
      hFacInf             = 0.2 _d 0
      hFacSup             = 2.0 _d 0
      select_rStar        = 0
      selectNHfreeSurf    = 0
      selectAddFluid      = 0
      useRealFreshWaterFlux = .FALSE.
      temp_EvPrRn = UNSET_RL
      salt_EvPrRn = 0.
      temp_addMass = UNSET_RL
      salt_addMass = UNSET_RL
      balanceEmPmR        = .FALSE.
      balanceQnet         = .FALSE.
      balancePrintMean    = .FALSE.

C--   Atmospheric physical parameters (e.g.: EOS)
      celsius2K = 273.16 _d 0
      atm_Po =  1. _d 5
      atm_Cp = 1004. _d 0
      atm_Rd = UNSET_RL
      atm_kappa = 2. _d 0 / 7. _d 0
      atm_Rq = 0. _d 0
      integr_GeoPot = 2
      selectFindRoSurf = 0

C--   Elliptic solver parameters
      cg2dMaxIters       = 150
      cg2dTargetResidual = 1. _d -7
      cg2dTargetResWunit = -1.
      cg2dChkResFreq     = 1
      cg2dpcOffDFac      = 0.51 _d 0
      cg2dPreCondFreq    = 1
      cg3dMaxIters       = 150
      cg3dTargetResidual = 1. _d -7
      cg3dChkResFreq     = 1
      useSRCGSolver      = .FALSE.

C--   Time stepping parameters
      deltaT            = 0. _d 0
      deltaTmom         = 0. _d 0
      deltaTfreesurf    = 0. _d 0
      DO k=1,Nr
        dTtracerLev(k)  = 0. _d 0
      ENDDO
      baseTime          = 0. _d 0
      nIter0            = -1
      startTime         = UNSET_RL
      pickupSuff        = ' '
      pickupStrictlyMatch = .TRUE.
      nTimeSteps        = 0
      nEndIter          = 0
      endTime           = 0. _d 0
      momForcingOutAB   = UNSET_I
      tracForcingOutAB  = UNSET_I
      momDissip_In_AB   = .TRUE.
      doAB_onGtGs       = .TRUE.
      abEps             = 0.01 _d 0
#ifdef ALLOW_ADAMSBASHFORTH_3
      alph_AB           = 0.5 _d 0
      beta_AB           = 5. _d 0 / 12. _d 0
      startFromPickupAB2= .FALSE.
#else
      alph_AB           = UNSET_RL
      beta_AB           = UNSET_RL
      startFromPickupAB2= .TRUE.
#endif
      pChkPtFreq        = deltaT*0
      chkPtFreq         = deltaT*0
      outputTypesInclusive = .FALSE.
      pickup_read_mdsio = .TRUE.
      pickup_write_mdsio= .TRUE.
      pickup_write_immed= .FALSE.
      writePickupAtEnd  = .TRUE.
      dumpFreq          = deltaT*0
      adjDumpFreq       = deltaT*0
      diagFreq          = deltaT*0
      dumpInitAndLast   = .TRUE.
      snapshot_mdsio    = .TRUE.
      monitorFreq       = -1.
      adjMonitorFreq    = 0.
      monitorSelect     = UNSET_I
      monitor_stdio     = .TRUE.
      taveFreq          = deltaT*0
      timeave_mdsio     = .TRUE.
      tave_lastIter     = 0.5 _d 0
      writeStatePrec    = precFloat64
      writeBinaryPrec   = precFloat32
      readBinaryPrec    = precFloat32
      cAdjFreq          =  0. _d 0
      tauCD             =  0. _d 0
      tauThetaClimRelax =  0. _d 0
      tauSaltClimRelax  =  0. _d 0
      periodicExternalForcing = .FALSE.
      externForcingPeriod     = 0.
      externForcingCycle      = 0.
      tCylIn             = 0.
      tCylOut            = 20.

C--   Input files
      bathyFile       = ' '
      topoFile        = ' '
      hydrogSaltFile  = ' '
      hydrogThetaFile = ' '
      maskIniTemp     = .TRUE.
      maskIniSalt     = .TRUE.
      checkIniTemp    = .TRUE.
      checkIniSalt    = .TRUE.
      diffKrFile      = ' '
      viscAhDfile     = ' '
      viscAhZfile     = ' '
      viscA4Dfile     = ' '
      viscA4Zfile     = ' '
      zonalWindFile   = ' '
      meridWindFile   = ' '
      thetaClimFile   = ' '
      saltClimFile    = ' '
      EmPmRfile       = ' '
      saltFluxFile    = ' '
      surfQfile       = ' '
      surfQnetFile    = ' '
      surfQswFile     = ' '
      uVelInitFile    = ' '
      vVelInitFile    = ' '
      pSurfInitFile   = ' '
      dQdTFile        = ' '
      ploadFile       = ' '
      eddyPsiXFile    = ' '
      eddyPsiYFile    = ' '
      lambdaThetaFile = ' '
      lambdaSaltFile  = ' '
      mdsioLocalDir   = ' '
      adTapeDir       = ' '
      the_run_name    = ' '

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/model/src/set_defaults.F,v 1.153 2011/04/20 13:45:14 dimitri Exp $
2	C $Name: $
3
4	#include "CPP_OPTIONS.h"
5
6	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7	CBOP
8	C !ROUTINE: SET_DEFAULTS
9	C !INTERFACE:
10	SUBROUTINE SET_DEFAULTS(
11	O viscArDefault, diffKrTDefault, diffKrSDefault,
12	O hFacMinDrDefault, delRdefault,
13	I myThid )
14
15	C !DESCRIPTION:
16	C Routine to set model "parameter defaults".
17
18	C !USES:
19	IMPLICIT NONE
20	#include "SIZE.h"
21	#include "EEPARAMS.h"
22	#include "PARAMS.h"
23	Cml#include "EOS.h"
24	c#include "GRID.h"
25
26	C !INPUT/OUTPUT PARAMETERS:
27	C myThid :: my Thread Id. Number
28	INTEGER myThid
29	_RL viscArDefault
30	_RL diffKrTDefault
31	_RL diffKrSDefault
32	_RL hFacMinDrDefault
33	_RL delRDefault(Nr)
34
35	C !LOCAL VARIABLES:
36	C i, j, k :: Loop counters
37	INTEGER i, j, k
38	CEOP
39
40	C-- Grid parameters
41	C- Vertical gridding
42	delRFile = ' '
43	delRcFile = ' '
44	hybSigmFile = ' '
45	Ro_SeaLevel = 0.
46	rSigmaBnd = UNSET_RL
47	selectSigmaCoord = 0
48	DO k=1,Nr
49	delRdefault(k) = 0.
50	ENDDO
51	DO k=1,Nr+1
52	delRc(k) = UNSET_RL
53	ENDDO
54	C- vertical profile
55	tRefFile = ' '
56	sRefFile = ' '
57	rhoRefFile = ' '
58
59	C- Horizontal gridding
60	delXFile = ' '
61	delYFile = ' '
62	horizGridFile = ' '
63	deepAtmosphere = .FALSE.
64	xgOrigin = UNSET_RL
65	ygOrigin = UNSET_RL
66	DO i=1,Nx
67	delX(i) = UNSET_RL
68	ENDDO
69	DO j=1,Ny
70	delY(j) = UNSET_RL
71	ENDDO
72	C In cartesian coords distances are in metres
73	usingCartesianGrid = .FALSE.
74	C In spherical polar distances are in degrees
75	usingSphericalPolarGrid = .FALSE.
76	rSphere = UNSET_RL
77	C General curvilinear coordinate system
78	usingCurvilinearGrid= .FALSE.
79	radius_fromHorizGrid= UNSET_RL
80	C General cylindrical coordinate system
81	usingCylindricalGrid= .FALSE.
82	C Coriolis map:
83	selectCoriMap = -1
84	use3dCoriolis = .TRUE.
85	C grid rotation
86	rotateGrid = .FALSE.
87	phiEuler = 0. _d 0
88	thetaEuler = 0. _d 0
89	psiEuler = 0. _d 0
90
91	C-- Set default "physical" parameters
92	nh_Am2 = 1. _d 0
93	gravity = 9.81 _d 0
94	rhoNil = 999.8 _d 0
95	C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)
96	C (so that the default produces same results as before)
97	c rhoConstFresh = 999.8 _d 0
98	f0 = 1. _d -4
99	beta = 1. _d -11
100	fPrime = 0. _d 0
101	C- Earth rotation period is 86400*365.25/366.25 (use to be 1.day)
102	rotationPeriod = 86164. _d 0
103	viscAh = 0. _d 3
104	viscAhGrid = 0. _d 0
105	viscAhGridMin = 0. _d 0
106	viscAhGridMax = 1. _d 21
107	viscAhMax = 1. _d 21
108	viscAhReMax = 0. _d 0
109	viscC2leith = 0. _d 0
110	viscC2leithD = 0. _d 0
111	viscC2smag = 0. _d 0
112	diffKhT = 0. _d 3
113	diffKhS = 0. _d 3
114	viscArDefault = 0. _d -3
115	no_slip_sides = .TRUE.
116	no_slip_bottom = .TRUE.
117	sideDragFactor = 2. _d 0
118	bottomDragLinear = 0.
119	bottomDragQuadratic = 0.
120	smoothAbsFuncRange = 0. _d 0
121	diffKrTDefault = 0. _d -3
122	diffKrSDefault = 0. _d -3
123	diffKrBL79surf = 0. _d 0
124	diffKrBL79deep = 0. _d 0
125	diffKrBL79scl = 200. _d 0
126	diffKrBL79Ho = -2000. _d 0
127	BL79LatVary = 30
128	diffKrBLEQsurf = UNSET_RL
129	diffKrBLEQdeep = UNSET_RL
130	diffKrBLEQscl = UNSET_RL
131	diffKrBLEQHo = UNSET_RL
132	viscA4 = 0. _d 11
133	viscA4Grid = 0. _d 0
134	viscA4GridMax = 1. _d 21
135	viscA4GridMin = 0. _d 0
136	viscA4Max = 1. _d 21
137	viscA4ReMax = 0. _d 0
138	viscC4leith = 0. _d 0
139	viscC4leithD = 0. _d 0
140	viscC4smag = 0. _d 0
141	diffK4T = 0. _d 11
142	diffK4S = 0. _d 11
143	cosPower = 0.
144	HeatCapacity_Cp = 3994. _d 0
145	Cml tAlpha = 2. _d -4
146	Cml sBeta = 7.4 _d -4
147	eosType = 'LINEAR'
148	buoyancyRelation = 'OCEANIC'
149	hFacMin = 1. _d 0
150	hFacMinDrDefault = 0. _d 0
151	implicitIntGravWave = .FALSE.
152	staggerTimeStep = .FALSE.
153	momViscosity = .TRUE.
154	momAdvection = .TRUE.
155	momForcing = .TRUE.
156	useCoriolis = .TRUE.
157	momPressureForcing = .TRUE.
158	momStepping = .TRUE.
159	vectorInvariantMomentum = .FALSE.
160	tempStepping = .TRUE.
161	tempAdvection = .TRUE.
162	tempForcing = .TRUE.
163	saltStepping = .TRUE.
164	saltAdvection = .TRUE.
165	saltForcing = .TRUE.
166	metricTerms = .TRUE.
167	useNHMTerms = .FALSE.
168	useFullLeith = .FALSE.
169	useAreaViscLength = .FALSE.
170	useStrainTensionVisc= .FALSE.
171	implicitDiffusion = .FALSE.
172	implicitViscosity = .FALSE.
173	momImplVertAdv = .FALSE.
174	tempImplVertAdv = .FALSE.
175	saltImplVertAdv = .FALSE.
176	nonHydrostatic = .FALSE.
177	quasiHydrostatic = .FALSE.
178	globalFiles = .FALSE.
179	useSingleCpuIO = .FALSE.
180	allowFreezing = .FALSE.
181	allowInteriorFreezing = .FALSE.
182	useOldFreezing = .FALSE.
183	ivdc_kappa = 0. _d 0
184	hMixCriteria = -.8 _d 0
185	dRhoSmall = 1. _d -6
186	hMixSmooth = 0. _d 0
187	usePickupBeforeC54 = .FALSE.
188	debugMode = .FALSE.
189	tempAdvScheme = 2
190	saltAdvScheme = 2
191	multiDimAdvection = .TRUE.
192	useMultiDimAdvec = .FALSE.
193	useCDscheme = .FALSE.
194	useEnergyConservingCoriolis = .FALSE.
195	useJamartWetPoints = .FALSE.
196	useJamartMomAdv = .FALSE.
197	selectVortScheme = UNSET_I
198	upwindVorticity = .FALSE.
199	highOrderVorticity = .FALSE.
200	useAbsVorticity = .FALSE.
201	upwindShear = .FALSE.
202	selectKEscheme = 0
203	debugLevel = debLevA
204	inAdMode = .FALSE.
205	inAdExact = .TRUE.
206
207	C-- Set (free)surface-related parameters
208	implicitFreeSurface = .FALSE.
209	rigidLid = .FALSE.
210	implicSurfPress = 1. _d 0
211	implicDiv2Dflow = 1. _d 0
212	exactConserv = .FALSE.
213	linFSConserveTr = .FALSE.
214	uniformLin_PhiSurf = .TRUE.
215	nonlinFreeSurf = 0
216	hFacInf = 0.2 _d 0
217	hFacSup = 2.0 _d 0
218	select_rStar = 0
219	selectNHfreeSurf = 0
220	selectAddFluid = 0
221	useRealFreshWaterFlux = .FALSE.
222	temp_EvPrRn = UNSET_RL
223	salt_EvPrRn = 0.
224	temp_addMass = UNSET_RL
225	salt_addMass = UNSET_RL
226	balanceEmPmR = .FALSE.
227	balanceQnet = .FALSE.
228	balancePrintMean = .FALSE.
229
230	C-- Atmospheric physical parameters (e.g.: EOS)
231	celsius2K = 273.16 _d 0
232	atm_Po = 1. _d 5
233	atm_Cp = 1004. _d 0
234	atm_Rd = UNSET_RL
235	atm_kappa = 2. _d 0 / 7. _d 0
236	atm_Rq = 0. _d 0
237	integr_GeoPot = 2
238	selectFindRoSurf = 0
239
240	C-- Elliptic solver parameters
241	cg2dMaxIters = 150
242	cg2dTargetResidual = 1. _d -7
243	cg2dTargetResWunit = -1.
244	cg2dChkResFreq = 1
245	cg2dpcOffDFac = 0.51 _d 0
246	cg2dPreCondFreq = 1
247	cg3dMaxIters = 150
248	cg3dTargetResidual = 1. _d -7
249	cg3dChkResFreq = 1
250	useSRCGSolver = .FALSE.
251
252	C-- Time stepping parameters
253	deltaT = 0. _d 0
254	deltaTmom = 0. _d 0
255	deltaTfreesurf = 0. _d 0
256	DO k=1,Nr
257	dTtracerLev(k) = 0. _d 0
258	ENDDO
259	baseTime = 0. _d 0
260	nIter0 = -1
261	startTime = UNSET_RL
262	pickupSuff = ' '
263	pickupStrictlyMatch = .TRUE.
264	nTimeSteps = 0
265	nEndIter = 0
266	endTime = 0. _d 0
267	momForcingOutAB = UNSET_I
268	tracForcingOutAB = UNSET_I
269	momDissip_In_AB = .TRUE.
270	doAB_onGtGs = .TRUE.
271	abEps = 0.01 _d 0
272	#ifdef ALLOW_ADAMSBASHFORTH_3
273	alph_AB = 0.5 _d 0
274	beta_AB = 5. _d 0 / 12. _d 0
275	startFromPickupAB2= .FALSE.
276	#else
277	alph_AB = UNSET_RL
278	beta_AB = UNSET_RL
279	startFromPickupAB2= .TRUE.
280	#endif
281	pChkPtFreq = deltaT*0
282	chkPtFreq = deltaT*0
283	outputTypesInclusive = .FALSE.
284	pickup_read_mdsio = .TRUE.
285	pickup_write_mdsio= .TRUE.
286	pickup_write_immed= .FALSE.
287	writePickupAtEnd = .TRUE.
288	dumpFreq = deltaT*0
289	adjDumpFreq = deltaT*0
290	diagFreq = deltaT*0
291	dumpInitAndLast = .TRUE.
292	snapshot_mdsio = .TRUE.
293	monitorFreq = -1.
294	adjMonitorFreq = 0.
295	monitorSelect = UNSET_I
296	monitor_stdio = .TRUE.
297	taveFreq = deltaT*0
298	timeave_mdsio = .TRUE.
299	tave_lastIter = 0.5 _d 0
300	writeStatePrec = precFloat64
301	writeBinaryPrec = precFloat32
302	readBinaryPrec = precFloat32
303	cAdjFreq = 0. _d 0
304	tauCD = 0. _d 0
305	tauThetaClimRelax = 0. _d 0
306	tauSaltClimRelax = 0. _d 0
307	periodicExternalForcing = .FALSE.
308	externForcingPeriod = 0.
309	externForcingCycle = 0.
310	tCylIn = 0.
311	tCylOut = 20.
312
313	C-- Input files
314	bathyFile = ' '
315	topoFile = ' '
316	hydrogSaltFile = ' '
317	hydrogThetaFile = ' '
318	maskIniTemp = .TRUE.
319	maskIniSalt = .TRUE.
320	checkIniTemp = .TRUE.
321	checkIniSalt = .TRUE.
322	diffKrFile = ' '
323	viscAhDfile = ' '
324	viscAhZfile = ' '
325	viscA4Dfile = ' '
326	viscA4Zfile = ' '
327	zonalWindFile = ' '
328	meridWindFile = ' '
329	thetaClimFile = ' '
330	saltClimFile = ' '
331	EmPmRfile = ' '
332	saltFluxFile = ' '
333	surfQfile = ' '
334	surfQnetFile = ' '
335	surfQswFile = ' '
336	uVelInitFile = ' '
337	vVelInitFile = ' '
338	pSurfInitFile = ' '
339	dQdTFile = ' '
340	ploadFile = ' '
341	eddyPsiXFile = ' '
342	eddyPsiYFile = ' '
343	lambdaThetaFile = ' '
344	lambdaSaltFile = ' '
345	mdsioLocalDir = ' '
346	adTapeDir = ' '
347	the_run_name = ' '
348
349	RETURN
350	END