/[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.42 by adcroft, Thu Nov 7 21:51:15 2002 UTC revision 1.60 by adcroft, Wed May 26 14:50:10 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          debugLevel          = debLevA
168    
169  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
170        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
# Line 159  C--   Set (free)surface-related paramete Line 174  C--   Set (free)surface-related paramete
174        exactConserv        = .FALSE.        exactConserv        = .FALSE.
175        uniformLin_PhiSurf  = .TRUE.        uniformLin_PhiSurf  = .TRUE.
176        nonlinFreeSurf      = 0        nonlinFreeSurf      = 0
177        hFacInf             = 1.D0        hFacInf             = 0.2 _d 0
178        hFacSup             = 1.D0        hFacSup             = 2.0 _d 0
179          select_rStar        = 0
180        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
181        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
182        salt_EvPrRn = 0.        salt_EvPrRn = 0.
183        trac_EvPrRn = UNSET_RL        trac_EvPrRn = UNSET_RL
184    
185  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
186        atm_po =  1.D5        celsius2K = 273.16 _d 0
187        atm_cp =  1004.D0        atm_Po =  1. _d 5
188        atm_kappa = 2.D0 / 7.D0        atm_Cp = 1004. _d 0
189        Integr_GeoPot = 2        atm_Rd = UNSET_RL
190          atm_kappa = 2. _d 0 / 7. _d 0
191          atm_Rq = 0. _d 0
192          integr_GeoPot = 2
193          selectFindRoSurf = 0
194    
195  C--   Elliptic solver parameters  C--   Elliptic solver parameters
196        cg2dMaxIters       = 150        cg2dMaxIters       = 150
# Line 186  C--   Time stepping parameters Line 206  C--   Time stepping parameters
206        deltaT            = 0.        deltaT            = 0.
207        nIter0            = 0        nIter0            = 0
208        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
209          pickupSuff        = ' '
210        nTimeSteps        = 0        nTimeSteps        = 0
211        nEndIter          = nIter0+nTimeSteps        nEndIter          = nIter0+nTimeSteps
212        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
# Line 194  C--   Time stepping parameters Line 215  C--   Time stepping parameters
215        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
216        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
217        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
218          adjDumpFreq       = deltaT*0
219        diagFreq          = deltaT*0        diagFreq          = deltaT*0
220        monitorFreq       = -1.        monitorFreq       = -1.
221        taveFreq          = deltaT*0        taveFreq          = deltaT*0
# Line 230  C--   Input files Line 252  C--   Input files
252        pSurfInitFile   = ' '        pSurfInitFile   = ' '
253        dQdTFile        = ' '        dQdTFile        = ' '
254        ploadFile       = ' '        ploadFile       = ' '
255          mdsioLocalDir   = ' '
256    
257  C  C
258        RETURN        RETURN
Legend:
Removed from v.1.42  
changed lines
  Added in v.1.60
  ViewVC Help
Powered by ViewVC 1.1.22