/[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.48 by adcroft, Thu Jun 5 16:03:05 2003 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 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          gravitySign         = 1.D0
89        gravity             = 9.81D0        gravity             = 9.81D0
90        rhoNil              = 999.8D0        rhoNil              = 999.8D0
91          rhoConstFresh       = 999.8D0
92        f0                  = 1.D-4        f0                  = 1.D-4
93        beta                = 1.D-11        beta                = 1.D-11
94        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)
95          rotationPeriod      = 86400. _d 0
96        viscAh              = 0.D3        viscAh              = 0.D3
97        viscAstrain         = 0.D3        viscAstrain         = 0.D3
98        viscAtension        = 0.D3        viscAtension        = 0.D3
# Line 103  C--   Set default "physical" parameters Line 108  C--   Set default "physical" parameters
108        diffK4S             = 0.D11        diffK4S             = 0.D11
109        cosPower            = 0.        cosPower            = 0.
110        HeatCapacity_Cp     = 3994.D0        HeatCapacity_Cp     = 3994.D0
111        tAlpha              = 2.D-4  Cml      tAlpha              = 2.D-4
112        sBeta               = 7.4D-4  Cml      sBeta               = 7.4D-4
113        eosType             = 'LINEAR'        eosType             = 'LINEAR'
114        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
115        hFacMin             = 1.D0        hFacMin             = 1.D0
# Line 118  C--   Set default "physical" parameters Line 123  C--   Set default "physical" parameters
123        momStepping         = .TRUE.        momStepping         = .TRUE.
124        vectorInvariantMomentum = .FALSE.        vectorInvariantMomentum = .FALSE.
125        tempStepping        = .TRUE.        tempStepping        = .TRUE.
126          tempAdvection       = .TRUE.
127          tempForcing         = .TRUE.
128        saltStepping        = .TRUE.        saltStepping        = .TRUE.
129          saltAdvection       = .TRUE.
130          saltForcing         = .TRUE.
131        tr1Stepping         = .FALSE.        tr1Stepping         = .FALSE.
132        metricTerms         = .TRUE.        metricTerms         = .TRUE.
133          useNHMTerms         = .TRUE.
134        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
135        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
136        nonHydrostatic      = .FALSE.        nonHydrostatic      = .FALSE.
137          quasiHydrostatic    = .FALSE.
138        globalFiles         = .FALSE.        globalFiles         = .FALSE.
139          useSingleCpuIO      = .FALSE.
140        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
141        ivdc_kappa          = 0.D0        ivdc_kappa          = 0.D0
142        groundAtK1          = .FALSE.        groundAtK1          = .FALSE.
# Line 138  C--   Set default "physical" parameters Line 150  C--   Set default "physical" parameters
150        saltAdvScheme       = 2        saltAdvScheme       = 2
151        tracerAdvScheme     = 2        tracerAdvScheme     = 2
152        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
153          useCDscheme         = .FALSE.
154        useEnergyConservingCoriolis = .FALSE.        useEnergyConservingCoriolis = .FALSE.
155        useJamartWetPoints  = .FALSE.        useJamartWetPoints  = .FALSE.
156    
# Line 149  C--   Set (free)surface-related paramete Line 162  C--   Set (free)surface-related paramete
162        exactConserv        = .FALSE.        exactConserv        = .FALSE.
163        uniformLin_PhiSurf  = .TRUE.        uniformLin_PhiSurf  = .TRUE.
164        nonlinFreeSurf      = 0        nonlinFreeSurf      = 0
165        hFacInf             = 1.D0        hFacInf             = 0.2 _d 0
166        hFacSup             = 1.D0        hFacSup             = 2.0 _d 0
167          select_rStar        = 0
168        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
169        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
170        salt_EvPrRn = 0.        salt_EvPrRn = 0.
171        trac_EvPrRn = UNSET_RL        trac_EvPrRn = UNSET_RL
172    
173  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
174        atm_po =  1.D5        atm_Po =  1.D5
175        atm_cp =  1004.D0        atm_Cp =  1004.D0
176          atm_Rd =  UNSET_RL
177        atm_kappa = 2.D0 / 7.D0        atm_kappa = 2.D0 / 7.D0
178        Integr_GeoPot = 2        integr_GeoPot = 2
179          selectFindRoSurf = 0
180    
181  C--   Elliptic solver parameters  C--   Elliptic solver parameters
182        cg2dMaxIters       = 150        cg2dMaxIters       = 150
# Line 176  C--   Time stepping parameters Line 192  C--   Time stepping parameters
192        deltaT            = 0.        deltaT            = 0.
193        nIter0            = 0        nIter0            = 0
194        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
195          pickupSuff        = ' '
196        nTimeSteps        = 0        nTimeSteps        = 0
197        nEndIter          = nIter0+nTimeSteps        nEndIter          = nIter0+nTimeSteps
198        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
199          forcing_In_AB     = .TRUE.
200        abEps             = 0.01        abEps             = 0.01
201        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
202        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
# Line 186  C--   Time stepping parameters Line 204  C--   Time stepping parameters
204        diagFreq          = deltaT*0        diagFreq          = deltaT*0
205        monitorFreq       = -1.        monitorFreq       = -1.
206        taveFreq          = deltaT*0        taveFreq          = deltaT*0
207          tave_lastIter     = 0.5 _d 0
208        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
209        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
210        readBinaryPrec    = precFloat32        readBinaryPrec    = precFloat32
# Line 217  C--   Input files Line 236  C--   Input files
236        vVelInitFile    = ' '        vVelInitFile    = ' '
237        pSurfInitFile   = ' '        pSurfInitFile   = ' '
238        dQdTFile        = ' '        dQdTFile        = ' '
239          ploadFile       = ' '
240    
241  C  C
242        RETURN        RETURN
243        END        END
   
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.48
  ViewVC Help
Powered by ViewVC 1.1.22