/[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.43 by jmc, Tue Dec 10 02:55:47 2002 UTC revision 1.61 by adcroft, Wed Jun 2 13:23:55 2004 UTC
# Line 65  C     Vertical gridding Line 65  C     Vertical gridding
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 85  C--   Set default "physical" parameters Line 85  C--   Set default "physical" parameters
85         tRef(K) = 30.D0 - FLOAT( K )         tRef(K) = 30.D0 - FLOAT( K )
86  Cml       sRef(K) = 35.D0  Cml       sRef(K) = 35.D0
87        ENDDO        ENDDO
       gravitySign         = 1.D0  
88        gravity             = 9.81D0        gravity             = 9.81D0
89        rhoNil              = 999.8D0        rhoNil              = 999.8D0
90        rhoConstFresh       = 999.8D0  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        omega               = 2.D0 * PI / ( 3600.D0 * 24.D0 )  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        viscAstrain         = 0.D3
102        viscAtension        = 0.D3        viscAtension        = 0.D3
103        diffKhT             = 0.D3        diffKhT             = 0.D3
# Line 103  Cml       sRef(K) = 35.D0 Line 108  Cml       sRef(K) = 35.D0
108        diffKrTDefault      = 0.D-3        diffKrTDefault      = 0.D-3
109        diffKrSDefault      = 0.D-3        diffKrSDefault      = 0.D-3
110        viscA4              = 0.D11        viscA4              = 0.D11
111          viscA4Grid          = 0.D0
112          viscA4Max           = 1.D21
113          viscC4leith         = 0.D0
114        diffK4T             = 0.D11        diffK4T             = 0.D11
115        diffK4S             = 0.D11        diffK4S             = 0.D11
116        cosPower            = 0.        cosPower            = 0.
# Line 129  Cml      sBeta               = 7.4D-4 Line 137  Cml      sBeta               = 7.4D-4
137        saltForcing         = .TRUE.        saltForcing         = .TRUE.
138        tr1Stepping         = .FALSE.        tr1Stepping         = .FALSE.
139        metricTerms         = .TRUE.        metricTerms         = .TRUE.
140        useNHMTerms         = .TRUE.        useNHMTerms         = .FALSE.
141        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
142        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
143          momImplVertAdv      = .FALSE.
144          tempImplVertAdv     = .FALSE.
145          saltImplVertAdv     = .FALSE.
146        nonHydrostatic      = .FALSE.        nonHydrostatic      = .FALSE.
147        quasiHydrostatic    = .FALSE.        quasiHydrostatic    = .FALSE.
148        globalFiles         = .FALSE.        globalFiles         = .FALSE.
149          useSingleCpuIO      = .FALSE.
150        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
151          useOldFreezing      = .FALSE.
152        ivdc_kappa          = 0.D0        ivdc_kappa          = 0.D0
153        groundAtK1          = .FALSE.        groundAtK1          = .FALSE.
154        bottomDragLinear    = 0.        bottomDragLinear    = 0.
# Line 148  Cml      sBeta               = 7.4D-4 Line 161  Cml      sBeta               = 7.4D-4
161        saltAdvScheme       = 2        saltAdvScheme       = 2
162        tracerAdvScheme     = 2        tracerAdvScheme     = 2
163        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
164          useCDscheme         = .FALSE.
165        useEnergyConservingCoriolis = .FALSE.        useEnergyConservingCoriolis = .FALSE.
166        useJamartWetPoints  = .FALSE.        useJamartWetPoints  = .FALSE.
167          SadournyCoriolis    = .FALSE.
168          upwindVorticity     = .FALSE.
169          highOrderVorticity  = .FALSE.
170          useAbsVorticity     = .FALSE.
171          debugLevel          = debLevA
172    
173  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
174        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
# Line 159  C--   Set (free)surface-related paramete Line 178  C--   Set (free)surface-related paramete
178        exactConserv        = .FALSE.        exactConserv        = .FALSE.
179        uniformLin_PhiSurf  = .TRUE.        uniformLin_PhiSurf  = .TRUE.
180        nonlinFreeSurf      = 0        nonlinFreeSurf      = 0
181        hFacInf             = 1.D0        hFacInf             = 0.2 _d 0
182        hFacSup             = 1.D0        hFacSup             = 2.0 _d 0
183          select_rStar        = 0
184        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
185        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
186        salt_EvPrRn = 0.        salt_EvPrRn = 0.
187        trac_EvPrRn = UNSET_RL        trac_EvPrRn = UNSET_RL
188    
189  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
190        atm_Po =  1.D5        celsius2K = 273.16 _d 0
191        atm_Cp =  1004.D0        atm_Po =  1. _d 5
192        atm_Rd =  UNSET_RL        atm_Cp = 1004. _d 0
193        atm_kappa = 2.D0 / 7.D0        atm_Rd = UNSET_RL
194          atm_kappa = 2. _d 0 / 7. _d 0
195          atm_Rq = 0. _d 0
196        integr_GeoPot = 2        integr_GeoPot = 2
197        selectFindRoSurf = 0        selectFindRoSurf = 0
198    
# Line 188  C--   Time stepping parameters Line 210  C--   Time stepping parameters
210        deltaT            = 0.        deltaT            = 0.
211        nIter0            = 0        nIter0            = 0
212        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
213          pickupSuff        = ' '
214        nTimeSteps        = 0        nTimeSteps        = 0
215        nEndIter          = nIter0+nTimeSteps        nEndIter          = nIter0+nTimeSteps
216        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
# Line 196  C--   Time stepping parameters Line 219  C--   Time stepping parameters
219        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
220        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
221        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
222          adjDumpFreq       = deltaT*0
223        diagFreq          = deltaT*0        diagFreq          = deltaT*0
224        monitorFreq       = -1.        monitorFreq       = -1.
225        taveFreq          = deltaT*0        taveFreq          = deltaT*0
# Line 232  C--   Input files Line 256  C--   Input files
256        pSurfInitFile   = ' '        pSurfInitFile   = ' '
257        dQdTFile        = ' '        dQdTFile        = ' '
258        ploadFile       = ' '        ploadFile       = ' '
259          mdsioLocalDir   = ' '
260    
261  C  C
262        RETURN        RETURN
Legend:
Removed from v.1.43  
changed lines
  Added in v.1.61
  ViewVC Help
Powered by ViewVC 1.1.22