/[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.115 by jmc, Mon Mar 20 15:11:18 2006 UTC
# Line 40  CEOP Line 40  CEOP
40    
41  C--   Grid parameters  C--   Grid parameters
42  C     Vertical gridding  C     Vertical gridding
43        rkFacDefault         = 1.D0        rkFacDefault        = 1. _d 0
44        horiVertRatio        = 1.D0        horiVertRatio       = 1. _d 0
45        Ro_SeaLevel = 0.        Ro_SeaLevel         = 0.
46        DO k=1,Nr        DO k=1,Nr
47         delRdefault(k) = 0.         delRdefault(k)     = 0.
48        ENDDO        ENDDO
49        DO k=1,Nr+1        DO k=1,Nr+1
50         delRc(k) = UNSET_RL         delRc(k)           = UNSET_RL
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
60         delX(i) = UNSET_RL         delX(i)            = UNSET_RL
61        ENDDO        ENDDO
62        DO j=1,Ny        DO j=1,Ny
63         delY(j) = UNSET_RL         delY(j)            = UNSET_RL
64        ENDDO        ENDDO
65  C     In spherical polar distances are in degrees  C     In spherical polar distances are in degrees
66        usingSphericalPolarGrid = .FALSE.        usingSphericalPolarGrid = .FALSE.
67        phiMin               = 0.0        phiMin              = 0.0
68        thetaMin             = 0.        thetaMin            = 0.
69        rSphere              = 6370. * 1.D3        rSphere             = 6370. _d 3
70  C     General curvilinear coordinate system  C     General curvilinear coordinate system
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
81        DO K =1,Nr        nh_Am2              = 1. _d 0
82         tRef(K) = 30.D0 - FLOAT( K )        gravity             = 9.81 _d 0
83  Cml       sRef(K) = 35.D0        rhoNil              = 999.8 _d 0
       ENDDO  
       gravity             = 9.81D0  
       rhoNil              = 999.8D0  
84  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)
85  C         (so that the default produces same results as before)  C         (so that the default produces same results as before)
86  c     rhoConstFresh       = 999.8D0  c     rhoConstFresh       = 999.8 _d 0
87        f0                  = 1.D-4        f0                  = 1. _d -4
88        beta                = 1.D-11        beta                = 1. _d -11
89  C-    Always use 1 day in the past but should be 86164 (=86400*365.25/366.25)  C-    Earth rotation period is 86400*365.25/366.25 (use to be 1.day)
90        rotationPeriod      = 86400. _d 0        rotationPeriod      = 86164. _d 0
91        viscAh              = 0.D3        viscAh              = 0. _d 3
92        viscAhGrid          = 0.D0        viscAhGrid          = 0. _d 0
93        viscAhMax           = 1.D21        viscAhGridMin       = 0. _d 0
94        viscC2leith         = 0.D0        viscAhGridMax       = 1. _d 21
95        viscAstrain         = 0.D3        viscAhMax           = 1. _d 21
96        viscAtension        = 0.D3        viscAhReMax         = 0. _d 0
97        diffKhT             = 0.D3        viscC2leith         = 0. _d 0
98        diffKhS             = 0.D3        viscC2leithD        = 0. _d 0
99        viscArDefault       = 0.D-3        viscC2smag          = 0. _d 0
100          diffKhT             = 0. _d 3
101          diffKhS             = 0. _d 3
102          viscArDefault       = 0. _d -3
103        no_slip_sides       = .TRUE.        no_slip_sides       = .TRUE.
104        no_slip_bottom      = .TRUE.        no_slip_bottom      = .TRUE.
105        diffKrTDefault      = 0.D-3        sideDragFactor      = 2. _d 0
106        diffKrSDefault      = 0.D-3        bottomDragLinear    = 0.
107        diffKrBL79surf      = 0.D0        bottomDragQuadratic = 0.
108        diffKrBL79deep      = 0.D0        diffKrTDefault      = 0. _d -3
109        diffKrBL79scl       = 200.D0        diffKrSDefault      = 0. _d -3
110        diffKrBL79Ho        = -2000.D0        diffKrBL79surf      = 0. _d 0
111        viscA4              = 0.D11        diffKrBL79deep      = 0. _d 0
112        viscA4Grid          = 0.D0        diffKrBL79scl       = 200. _d 0
113        viscA4Max           = 1.D21        diffKrBL79Ho        = -2000. _d 0
114        viscC4leith         = 0.D0        viscA4              = 0. _d 11
115        diffK4T             = 0.D11        viscA4Grid          = 0. _d 0
116        diffK4S             = 0.D11        viscA4GridMax       = 1. _d 21
117          viscA4GridMin       = 0. _d 0
118          viscA4Max           = 1. _d 21
119          viscA4ReMax         = 0. _d 0
120          viscC4leith         = 0. _d 0
121          viscC4leithD        = 0. _d 0
122          viscC4smag          = 0. _d 0
123          diffK4T             = 0. _d 11
124          diffK4S             = 0. _d 11
125        cosPower            = 0.        cosPower            = 0.
126        HeatCapacity_Cp     = 3994.D0        HeatCapacity_Cp     = 3994. _d 0
127  Cml      tAlpha              = 2.D-4  Cml      tAlpha              = 2. _d -4
128  Cml      sBeta               = 7.4D-4  Cml      sBeta               = 7.4 _d -4
129        eosType             = 'LINEAR'        eosType             = 'LINEAR'
130        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
131        hFacMin             = 1.D0        hFacMin             = 1. _d 0
132        hFacMinDrDefault    = 0.D0        hFacMinDrDefault    = 0. _d 0
133          implicitIntGravWave = .FALSE.
134        staggerTimeStep     = .FALSE.        staggerTimeStep     = .FALSE.
135        momViscosity        = .TRUE.        momViscosity        = .TRUE.
136        momAdvection        = .TRUE.        momAdvection        = .TRUE.
# 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 143  Cml      sBeta               = 7.4D-4 Line 161  Cml      sBeta               = 7.4D-4
161        useSingleCpuIO      = .FALSE.        useSingleCpuIO      = .FALSE.
162        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
163        useOldFreezing      = .FALSE.        useOldFreezing      = .FALSE.
164        ivdc_kappa          = 0.D0        ivdc_kappa          = 0. _d 0
       groundAtK1          = .FALSE.  
       bottomDragLinear    = 0.  
       bottomDragQuadratic = 0.  
165        usePickupBeforeC35    = .FALSE.        usePickupBeforeC35    = .FALSE.
166        usePickupBeforeC54    = .FALSE.        usePickupBeforeC54    = .FALSE.
167        debugMode             = .FALSE.        debugMode             = .FALSE.
       readPickupWithTracer  = .FALSE.  
       writePickupWithTracer = .FALSE.  
168        tempAdvScheme       = 2        tempAdvScheme       = 2
169        saltAdvScheme       = 2        saltAdvScheme       = 2
       tracerAdvScheme     = 2  
170        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
171        useCDscheme         = .FALSE.        useCDscheme         = .FALSE.
172        useEnergyConservingCoriolis = .FALSE.        useEnergyConservingCoriolis = .FALSE.
# Line 164  Cml      sBeta               = 7.4D-4 Line 176  Cml      sBeta               = 7.4D-4
176        upwindVorticity     = .FALSE.        upwindVorticity     = .FALSE.
177        highOrderVorticity  = .FALSE.        highOrderVorticity  = .FALSE.
178        useAbsVorticity     = .FALSE.        useAbsVorticity     = .FALSE.
179          upwindShear         = .FALSE.
180          selectKEscheme      = 0
181        debugLevel          = debLevA        debugLevel          = debLevA
182          inAdMode            = .FALSE.
183          inAdExact           = .TRUE.
184    
185  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
186        implicitFreeSurface = .TRUE.        implicitFreeSurface = .FALSE.
187        rigidLid            = .FALSE.        rigidLid            = .FALSE.
188        implicSurfPress     = 1.D0        implicSurfPress     = 1. _d 0
189        implicDiv2DFlow     = 1.D0        implicDiv2DFlow     = 1. _d 0
190        exactConserv        = .FALSE.        exactConserv        = .FALSE.
191        uniformLin_PhiSurf  = .TRUE.        uniformLin_PhiSurf  = .TRUE.
192        nonlinFreeSurf      = 0        nonlinFreeSurf      = 0
# Line 180  C--   Set (free)surface-related paramete Line 196  C--   Set (free)surface-related paramete
196        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
197        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
198        salt_EvPrRn = 0.        salt_EvPrRn = 0.
199        trac_EvPrRn = UNSET_RL        balanceEmPmR        = .FALSE.
200          balanceQnet         = .FALSE.
201          balancePrintMean    = .FALSE.
202    
203  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
204        celsius2K = 273.16 _d 0        celsius2K = 273.16 _d 0
# Line 194  C--   Atmospheric physical parameters (e Line 212  C--   Atmospheric physical parameters (e
212    
213  C--   Elliptic solver parameters  C--   Elliptic solver parameters
214        cg2dMaxIters       = 150        cg2dMaxIters       = 150
215        cg2dTargetResidual = 1.D-7        cg2dTargetResidual = 1. _d -7
216        cg2dTargetResWunit = -1.        cg2dTargetResWunit = -1.
217        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
218        cg2dpcOffDFac      = 0.51D0        cg2dpcOffDFac      = 0.51 _d 0
219        cg2dPreCondFreq    = 1        cg2dPreCondFreq    = 1
220        cg3dMaxIters       = 150        cg3dMaxIters       = 150
221        cg3dTargetResidual = 1.D-7        cg3dTargetResidual = 1. _d -7
222        cg3dChkResFreq     = 1        cg3dChkResFreq     = 1
223    
224  C--   Time stepping parameters  C--   Time stepping parameters
225        deltaT            = 0.        deltaT            = 0. _d 0
226          deltaTmom         = 0. _d 0
227          deltaTfreesurf    = 0. _d 0
228          DO k=1,Nr
229            dTtracerLev(k)  = 0. _d 0
230          ENDDO
231          baseTime          = 0. _d 0
232        nIter0            = 0        nIter0            = 0
233        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
234        pickupSuff        = ' '        pickupSuff        = ' '
235        nTimeSteps        = 0        nTimeSteps        = 0
236        nEndIter          = nIter0+nTimeSteps        nEndIter          = nIter0+nTimeSteps
237        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
238        forcing_In_AB     = .TRUE.        momForcingOutAB   = UNSET_I
239        abEps             = 0.01        tracForcingOutAB  = UNSET_I
240        pchkPtFreq        = deltaT*0        momDissip_In_AB   = .TRUE.
241          doAB_onGtGs       = .TRUE.
242          abEps             = 0.01 _d 0
243    #ifdef ALLOW_ADAMSBASHFORTH_3
244          alph_AB           = 0.5 _d 0
245          beta_AB           = 5. _d 0 / 12. _d 0
246          startFromPickupAB2= .FALSE.
247    #else
248          alph_AB           = UNSET_RL
249          beta_AB           = UNSET_RL
250          startFromPickupAB2= .TRUE.
251    #endif
252          pChkPtFreq        = deltaT*0
253        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
254          outputTypesInclusive = .FALSE.
255        pickup_read_mdsio = .TRUE.        pickup_read_mdsio = .TRUE.
256        pickup_write_mdsio= .TRUE.        pickup_write_mdsio= .TRUE.
257        pickup_write_immed= .FALSE.        pickup_write_immed= .FALSE.
258        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
259        adjDumpFreq       = deltaT*0        adjDumpFreq       = deltaT*0
260        diagFreq          = deltaT*0        diagFreq          = deltaT*0
261          dumpInitAndLast   = .TRUE.
262        snapshot_mdsio    = .TRUE.        snapshot_mdsio    = .TRUE.
263        monitorFreq       = -1.        monitorFreq       = -1.
264        monitor_mdsio     = .TRUE.        adjMonitorFreq    = 0.
265          monitor_stdio     = .TRUE.
266        taveFreq          = deltaT*0        taveFreq          = deltaT*0
267        timeave_mdsio     = .TRUE.        timeave_mdsio     = .TRUE.
268        tave_lastIter     = 0.5 _d 0        tave_lastIter     = 0.5 _d 0
# Line 233  C--   Time stepping parameters Line 272  C--   Time stepping parameters
272        nCheckLev         = 1        nCheckLev         = 1
273        checkPtSuff(1)    = 'ckptA'        checkPtSuff(1)    = 'ckptA'
274        checkPtSuff(2)    = 'ckptB'        checkPtSuff(2)    = 'ckptB'
275        cAdjFreq          =  0.D0        cAdjFreq          =  0. _d 0
276        tauCD             =  0.D0        tauCD             =  0. _d 0
277        tauThetaClimRelax =  0.D0        tauThetaClimRelax =  0. _d 0
278        tauSaltClimRelax  =  0.D0        tauSaltClimRelax  =  0. _d 0
279        tauTr1ClimRelax   =  0.D0        tauTr1ClimRelax   =  0. _d 0
280        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
281        externForcingPeriod     = 0.        externForcingPeriod     = 0.
282        externForcingCycle      = 0.        externForcingCycle      = 0.
283        tCyl              = 0.        tCylIn             = 0.
284          tCylOut            = 20.
285    
286  C--   Input files  C--   Input files
287        bathyFile       = ' '        bathyFile       = ' '
288        topoFile        = ' '        topoFile        = ' '
289          shelfIceFile    = ' '
290        hydrogSaltFile  = ' '        hydrogSaltFile  = ' '
291        hydrogThetaFile = ' '        hydrogThetaFile = ' '
292        zonalWindFile   = ' '        zonalWindFile   = ' '
# Line 253  C--   Input files Line 294  C--   Input files
294        thetaClimFile   = ' '        thetaClimFile   = ' '
295        saltClimFile    = ' '        saltClimFile    = ' '
296        EmPmRfile       = ' '        EmPmRfile       = ' '
297          saltFluxFile    = ' '
298        surfQfile       = ' '        surfQfile       = ' '
299        surfQnetFile    = ' '        surfQnetFile    = ' '
300        surfQswFile     = ' '        surfQswFile     = ' '
# Line 261  C--   Input files Line 303  C--   Input files
303        pSurfInitFile   = ' '        pSurfInitFile   = ' '
304        dQdTFile        = ' '        dQdTFile        = ' '
305        ploadFile       = ' '        ploadFile       = ' '
306          eddyTauxFile    = ' '
307          eddyTauyFile    = ' '
308          lambdaThetaFile = ' '
309          lambdaSaltFile  = ' '
310        mdsioLocalDir   = ' '        mdsioLocalDir   = ' '
311    
 C  
312        RETURN        RETURN
313        END        END
Legend:
Removed from v.1.70  
changed lines
  Added in v.1.115
  ViewVC Help
Powered by ViewVC 1.1.22