/[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.34 by jmc, Thu Mar 7 14:09:02 2002 UTC revision 1.61 by adcroft, Wed Jun 2 13:23:55 2004 UTC
# Line 31  C     === Global variables === Line 31  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"
36    
37  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
# Line 64  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 82  C     General curvilinear coordinate sys Line 83  C     General curvilinear coordinate sys
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        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 99  C--   Set default "physical" parameters Line 108  C--   Set default "physical" parameters
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.
117        HeatCapacity_Cp     = 3994.D0        HeatCapacity_Cp     = 3994.D0
118        tAlpha              = 2.D-4  Cml      tAlpha              = 2.D-4
119        sBeta               = 7.4D-4  Cml      sBeta               = 7.4D-4
120        eosType             = 'LINEAR'        eosType             = 'LINEAR'
121        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
122        hFacMin             = 1.D0        hFacMin             = 1.D0
# Line 118  C--   Set default "physical" parameters Line 130  C--   Set default "physical" parameters
130        momStepping         = .TRUE.        momStepping         = .TRUE.
131        vectorInvariantMomentum = .FALSE.        vectorInvariantMomentum = .FALSE.
132        tempStepping        = .TRUE.        tempStepping        = .TRUE.
133          tempAdvection       = .TRUE.
134          tempForcing         = .TRUE.
135        saltStepping        = .TRUE.        saltStepping        = .TRUE.
136          saltAdvection       = .TRUE.
137          saltForcing         = .TRUE.
138        tr1Stepping         = .FALSE.        tr1Stepping         = .FALSE.
139        metricTerms         = .TRUE.        metricTerms         = .TRUE.
140          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.
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 138  C--   Set default "physical" parameters Line 161  C--   Set default "physical" parameters
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 149  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_kappa = 2.D0 / 7.D0        atm_Cp = 1004. _d 0
193        Integr_GeoPot = 2        atm_Rd = UNSET_RL
194          atm_kappa = 2. _d 0 / 7. _d 0
195          atm_Rq = 0. _d 0
196          integr_GeoPot = 2
197          selectFindRoSurf = 0
198    
199  C--   Elliptic solver parameters  C--   Elliptic solver parameters
200        cg2dMaxIters       = 150        cg2dMaxIters       = 150
# Line 176  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)
217          forcing_In_AB     = .TRUE.
218        abEps             = 0.01        abEps             = 0.01
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
226          tave_lastIter     = 0.5 _d 0
227        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
228        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
229        readBinaryPrec    = precFloat32        readBinaryPrec    = precFloat32
# Line 217  C--   Input files Line 255  C--   Input files
255        vVelInitFile    = ' '        vVelInitFile    = ' '
256        pSurfInitFile   = ' '        pSurfInitFile   = ' '
257        dQdTFile        = ' '        dQdTFile        = ' '
258          ploadFile       = ' '
259          mdsioLocalDir   = ' '
260    
261  C  C
262        RETURN        RETURN
263        END        END
   

Legend:
Removed from v.1.34  
changed lines
  Added in v.1.61

  ViewVC Help
Powered by ViewVC 1.1.22