/[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.11 by heimbach, Mon Sep 11 20:51:48 2000 UTC revision 1.63 by adcroft, Mon Jun 14 17:48:17 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6    CBOP
7    C     !ROUTINE: SET_DEFAULTS
8    C     !INTERFACE:
9        SUBROUTINE SET_DEFAULTS(        SUBROUTINE SET_DEFAULTS(
10       O   viscArDefault, diffKrTDefault, diffKrSDefault,       O   viscArDefault, diffKrTDefault, diffKrSDefault,
11       O   hFacMinDrDefault, delRdefault, rkFacDefault,       O   hFacMinDrDefault, delRdefault, rkFacDefault,
12       I   myThid )       I   myThid )
 C     /==========================================================\  
 C     | SUBROUTINE SET_DEFAULTS                                  |  
 C     | o Routine to set model "parameters"                      |  
 C     |==========================================================|  
 C     | Notes:                                                   |  
 C     | ======                                                   |  
 C     | The present version of this routine is a place-holder.   |  
 C     | A production version needs to handle parameters from an  |  
 C     | external file and possibly reading in some initial field |  
 C     | values.                                                  |  
 C     \==========================================================/  
       IMPLICIT NONE  
13    
14    C     !DESCRIPTION: \bv
15    C     *==========================================================*
16    C     | SUBROUTINE SET_DEFAULTS                                  
17    C     | o Routine to set model "parameters"                      
18    C     *==========================================================*
19    C     | Notes:                                                    
20    C     | ======                                                    
21    C     | The present version of this routine is a place-holder.    
22    C     | A production version needs to handle parameters from an  
23    C     | external file and possibly reading in some initial field  
24    C     | values.                                                  
25    C     *==========================================================*
26    C     \ev
27    
28    C     !USES:
29          IMPLICIT NONE
30  C     === Global variables ===  C     === Global variables ===
31  #include "SIZE.h"  #include "SIZE.h"
32  #include "EEPARAMS.h"  #include "EEPARAMS.h"
33  #include "PARAMS.h"  #include "PARAMS.h"
34    Cml#include "EOS.h"
35  #include "GRID.h"  #include "GRID.h"
 #include "CG2D.h"  
36    
37    C     !INPUT/OUTPUT PARAMETERS:
38  C     === Routine arguments ===  C     === Routine arguments ===
39  C     myThid - Number of this instance of INI_PARMS  C     myThid - Number of this instance of INI_PARMS
40        INTEGER myThid        INTEGER myThid
# Line 36  C     myThid - Number of this instance o Line 45  C     myThid - Number of this instance o
45        _RL delRDefault(Nr)        _RL delRDefault(Nr)
46        _RS rkFacDefault        _RS rkFacDefault
47    
48    C     !LOCAL VARIABLES:
49  C     === Local variables ===  C     === Local variables ===
50  C     K, I, J - Loop counters  C     K, I, J - Loop counters
51        INTEGER K, I, J        INTEGER K, I, J
52    CEOP
53    
54  C--   Grid parameters  C--   Grid parameters
55  C     Vertical gridding  C     Vertical gridding
56        rkFacDefault         = 1.D0        rkFacDefault         = 1.D0
57        horiVertRatio        = 1.D0        horiVertRatio        = 1.D0
58          Ro_SeaLevel = 0.
59        DO k=1,Nr        DO k=1,Nr
60         delRdefault(k) = 0.         delRdefault(k) = 0.
61        ENDDO        ENDDO
62          DO k=1,Nr+1
63           delRc(k) = UNSET_RL
64          ENDDO
65    
66  C     Horizontal gridding  C     Horizontal gridding
67  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
68        usingCartesianGrid = .TRUE.        usingCartesianGrid = .FALSE.
69        DO i=1,Nx        DO i=1,Nx
70         delX(i) = UNSET_RL         delX(i) = UNSET_RL
71        ENDDO        ENDDO
# Line 62  C     In spherical polar distances are i Line 77  C     In spherical polar distances are i
77        phiMin               = 0.0        phiMin               = 0.0
78        thetaMin             = 0.        thetaMin             = 0.
79        rSphere              = 6370. * 1.D3        rSphere              = 6370. * 1.D3
80    C     General curvilinear coordinate system
81          usingCurvilinearGrid = .FALSE.
82    
83  C--   Set default "physical" parameters  C--   Set default "physical" parameters
84        DO K =1,Nr        DO K =1,Nr
85         tRef(K) = 30.D0 - FLOAT( K )         tRef(K) = 30.D0 - FLOAT( K )
86    Cml       sRef(K) = 35.D0
87        ENDDO        ENDDO
88        gravity             = 9.81D0        gravity             = 9.81D0
89        rhoNil              = 999.8D0        rhoNil              = 999.8D0
90    C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)
91    C         (so that the default produces same results as before)
92    c     rhoConstFresh       = 999.8D0
93        f0                  = 1.D-4        f0                  = 1.D-4
94        beta                = 1.D-11        beta                = 1.D-11
95    C-    Always use 1 day in the past but should be 86164 (=86400*365.25/366.25)
96          rotationPeriod      = 86400. _d 0
97        viscAh              = 0.D3        viscAh              = 0.D3
98          viscAhGrid          = 0.D0
99          viscAhMax           = 1.D21
100          viscC2leith         = 0.D0
101          viscAstrain         = 0.D3
102          viscAtension        = 0.D3
103        diffKhT             = 0.D3        diffKhT             = 0.D3
104        diffKhS             = 0.D3        diffKhS             = 0.D3
105        viscArDefault       = 0.D-3        viscArDefault       = 0.D-3
# Line 79  C--   Set default "physical" parameters Line 107  C--   Set default "physical" parameters
107        no_slip_bottom      = .TRUE.        no_slip_bottom      = .TRUE.
108        diffKrTDefault      = 0.D-3        diffKrTDefault      = 0.D-3
109        diffKrSDefault      = 0.D-3        diffKrSDefault      = 0.D-3
110          diffKrBL79surf      = 0.D0
111          diffKrBL79deep      = 0.D0
112          diffKrBL79scl       = 200.D0
113          diffKrBL79Ho        = -2000.D0
114        viscA4              = 0.D11        viscA4              = 0.D11
115          viscA4Grid          = 0.D0
116          viscA4Max           = 1.D21
117          viscC4leith         = 0.D0
118        diffK4T             = 0.D11        diffK4T             = 0.D11
119        diffK4S             = 0.D11        diffK4S             = 0.D11
120        cosPower            = 0.        cosPower            = 0.
121        HeatCapacity_Cp     = 3994.D0        HeatCapacity_Cp     = 3994.D0
122        tAlpha              = 2.D-4  Cml      tAlpha              = 2.D-4
123        sBeta               = 7.4D-4  Cml      sBeta               = 7.4D-4
124        eosType             = 'LINEAR'        eosType             = 'LINEAR'
125        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
       implicitFreeSurface = .TRUE.  
       rigidLid            = .FALSE.  
126        hFacMin             = 1.D0        hFacMin             = 1.D0
127        hFacMinDrDefault    = 0.D0        hFacMinDrDefault    = 0.D0
128          staggerTimeStep     = .FALSE.
129        momViscosity        = .TRUE.        momViscosity        = .TRUE.
130        momAdvection        = .TRUE.        momAdvection        = .TRUE.
131        momForcing          = .TRUE.        momForcing          = .TRUE.
132        useCoriolis         = .TRUE.        useCoriolis         = .TRUE.
133        momPressureForcing  = .TRUE.        momPressureForcing  = .TRUE.
134        momStepping         = .TRUE.        momStepping         = .TRUE.
135          vectorInvariantMomentum = .FALSE.
136        tempStepping        = .TRUE.        tempStepping        = .TRUE.
137          tempAdvection       = .TRUE.
138          tempForcing         = .TRUE.
139        saltStepping        = .TRUE.        saltStepping        = .TRUE.
140        metricTerms         = .FALSE.        saltAdvection       = .TRUE.
141          saltForcing         = .TRUE.
142          tr1Stepping         = .FALSE.
143          metricTerms         = .TRUE.
144          useNHMTerms         = .FALSE.
145        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
146        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
147        openBoundaries      = .FALSE.        momImplVertAdv      = .FALSE.
148          tempImplVertAdv     = .FALSE.
149          saltImplVertAdv     = .FALSE.
150        nonHydrostatic      = .FALSE.        nonHydrostatic      = .FALSE.
151          quasiHydrostatic    = .FALSE.
152        globalFiles         = .FALSE.        globalFiles         = .FALSE.
153          useSingleCpuIO      = .FALSE.
154        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
155          useOldFreezing      = .FALSE.
156        ivdc_kappa          = 0.D0        ivdc_kappa          = 0.D0
157          groundAtK1          = .FALSE.
158          bottomDragLinear    = 0.
159          bottomDragQuadratic = 0.
160          usePickupBeforeC35    = .FALSE.
161          debugMode             = .FALSE.
162          readPickupWithTracer  = .FALSE.
163          writePickupWithTracer = .FALSE.
164          tempAdvScheme       = 2
165          saltAdvScheme       = 2
166          tracerAdvScheme     = 2
167          multiDimAdvection   = .TRUE.
168          useCDscheme         = .FALSE.
169          useEnergyConservingCoriolis = .FALSE.
170          useJamartWetPoints  = .FALSE.
171          useJamartMomAdv     = .FALSE.
172          SadournyCoriolis    = .FALSE.
173          upwindVorticity     = .FALSE.
174          highOrderVorticity  = .FALSE.
175          useAbsVorticity     = .FALSE.
176          debugLevel          = debLevA
177    
178    C--   Set (free)surface-related parameters
179          implicitFreeSurface = .TRUE.
180          rigidLid            = .FALSE.
181          implicSurfPress     = 1.D0
182          implicDiv2DFlow     = 1.D0
183          exactConserv        = .FALSE.
184          uniformLin_PhiSurf  = .TRUE.
185          nonlinFreeSurf      = 0
186          hFacInf             = 0.2 _d 0
187          hFacSup             = 2.0 _d 0
188          select_rStar        = 0
189          useRealFreshWaterFlux = .FALSE.
190          temp_EvPrRn = UNSET_RL
191          salt_EvPrRn = 0.
192          trac_EvPrRn = UNSET_RL
193    
194    C--   Atmospheric physical parameters (e.g.: EOS)
195          celsius2K = 273.16 _d 0
196          atm_Po =  1. _d 5
197          atm_Cp = 1004. _d 0
198          atm_Rd = UNSET_RL
199          atm_kappa = 2. _d 0 / 7. _d 0
200          atm_Rq = 0. _d 0
201          integr_GeoPot = 2
202          selectFindRoSurf = 0
203    
204  C--   Elliptic solver parameters  C--   Elliptic solver parameters
205        cg2dMaxIters       = 150        cg2dMaxIters       = 150
206        cg2dTargetResidual = 1.D-7        cg2dTargetResidual = 1.D-7
207          cg2dTargetResWunit = -1.
208        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
209        cg2dpcOffDFac      = 0.51D0        cg2dpcOffDFac      = 0.51D0
210        cg3dMaxIters       = 150        cg3dMaxIters       = 150
# Line 122  C--   Time stepping parameters Line 215  C--   Time stepping parameters
215        deltaT            = 0.        deltaT            = 0.
216        nIter0            = 0        nIter0            = 0
217        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
218          pickupSuff        = ' '
219        nTimeSteps        = 0        nTimeSteps        = 0
220        nEndIter          = nIter0+nTimeSteps        nEndIter          = nIter0+nTimeSteps
221        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
222          forcing_In_AB     = .TRUE.
223        abEps             = 0.01        abEps             = 0.01
224        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
225        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
226        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
227          adjDumpFreq       = deltaT*0
228          diagFreq          = deltaT*0
229          monitorFreq       = -1.
230        taveFreq          = deltaT*0        taveFreq          = deltaT*0
231          tave_lastIter     = 0.5 _d 0
232        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
233        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
234        readBinaryPrec    = precFloat32        readBinaryPrec    = precFloat32
# Line 140  C--   Time stepping parameters Line 239  C--   Time stepping parameters
239        tauCD             =  0.D0        tauCD             =  0.D0
240        tauThetaClimRelax =  0.D0        tauThetaClimRelax =  0.D0
241        tauSaltClimRelax  =  0.D0        tauSaltClimRelax  =  0.D0
242          tauTr1ClimRelax   =  0.D0
243        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
244        externForcingPeriod     = 0.        externForcingPeriod     = 0.
245        externForcingCycle      = 0.        externForcingCycle      = 0.
246    
247  C--   Input files  C--   Input files
248        bathyFile       = ' '        bathyFile       = ' '
249          topoFile        = ' '
250        hydrogSaltFile  = ' '        hydrogSaltFile  = ' '
251        hydrogThetaFile = ' '        hydrogThetaFile = ' '
252        zonalWindFile   = ' '        zonalWindFile   = ' '
# Line 158  C--   Input files Line 259  C--   Input files
259        uVelInitFile    = ' '        uVelInitFile    = ' '
260        vVelInitFile    = ' '        vVelInitFile    = ' '
261        pSurfInitFile   = ' '        pSurfInitFile   = ' '
262          dQdTFile        = ' '
263  C--   OBCS        ploadFile       = ' '
264        DO I=1,Nx        mdsioLocalDir   = ' '
        OB_Jnorth(I)=0  
        OB_Jsouth(I)=0  
       ENDDO  
       DO J=1,Ny  
        OB_Ieast(J)=0  
        OB_Iwest(J)=0  
       ENDDO  
265    
266  C  C
267        RETURN        RETURN
268        END        END
   
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.63
  ViewVC Help
Powered by ViewVC 1.1.22