/[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.5 by adcroft, Wed Jul 28 16:32:14 1999 UTC revision 1.27 by adcroft, Thu Sep 6 14:23:57 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
# Line 24  C     === Global variables === Line 25  C     === Global variables ===
25  #include "EEPARAMS.h"  #include "EEPARAMS.h"
26  #include "PARAMS.h"  #include "PARAMS.h"
27  #include "GRID.h"  #include "GRID.h"
 #include "CG2D.h"  
28    
29  C     === Routine arguments ===  C     === Routine arguments ===
30  C     myThid - Number of this instance of INI_PARMS  C     myThid - Number of this instance of INI_PARMS
# Line 44  C--   Grid parameters Line 44  C--   Grid parameters
44  C     Vertical gridding  C     Vertical gridding
45        rkFacDefault         = 1.D0        rkFacDefault         = 1.D0
46        horiVertRatio        = 1.D0        horiVertRatio        = 1.D0
47          Ro_SeaLevel = 0.
48        DO k=1,Nr        DO k=1,Nr
49         delRdefault(k) = 0.         delRdefault(k) = 0.
50        ENDDO        ENDDO
# Line 52  C     Horizontal gridding Line 53  C     Horizontal gridding
53  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
54        usingCartesianGrid = .TRUE.        usingCartesianGrid = .TRUE.
55        DO i=1,Nx        DO i=1,Nx
56         delX(i) = 0.         delX(i) = UNSET_RL
57        ENDDO        ENDDO
58        DO j=1,Ny        DO j=1,Ny
59         delY(j) = 0.         delY(j) = UNSET_RL
60        ENDDO        ENDDO
61  C     In spherical polar distances are in degrees  C     In spherical polar distances are in degrees
62        usingSphericalPolarGrid = .FALSE.        usingSphericalPolarGrid = .FALSE.
63        phiMin               = 0.0        phiMin               = 0.0
64        thetaMin             = 0.        thetaMin             = 0.
65        rSphere              = 6370. * 1.D3        rSphere              = 6370. * 1.D3
66    C     General curvilinear coordinate system
67          usingCurvilinearGrid = .FALSE.
68    
69  C--   Set default "physical" parameters  C--   Set default "physical" parameters
70        DO K =1,Nr        DO K =1,Nr
# Line 71  C--   Set default "physical" parameters Line 74  C--   Set default "physical" parameters
74        rhoNil              = 999.8D0        rhoNil              = 999.8D0
75        f0                  = 1.D-4        f0                  = 1.D-4
76        beta                = 1.D-11        beta                = 1.D-11
77          omega               = 2.D0 * PI / ( 3600.D0 * 24.D0 )
78        viscAh              = 0.D3        viscAh              = 0.D3
79          viscAstrain         = 0.D3
80          viscAtension        = 0.D3
81        diffKhT             = 0.D3        diffKhT             = 0.D3
82        diffKhS             = 0.D3        diffKhS             = 0.D3
83        viscArDefault       = 0.D-3        viscArDefault       = 0.D-3
# Line 83  C--   Set default "physical" parameters Line 89  C--   Set default "physical" parameters
89        diffK4T             = 0.D11        diffK4T             = 0.D11
90        diffK4S             = 0.D11        diffK4S             = 0.D11
91        cosPower            = 0.        cosPower            = 0.
       GMmaxslope          = 1.D-2  
       GMlength            = 200.D3  
       GMalpha             = 0.D0  
       GMdepth             = 1000.D0  
       GMkbackground       = 0.D0  
       GMmaxval            = 2500.D0  
92        HeatCapacity_Cp     = 3994.D0        HeatCapacity_Cp     = 3994.D0
93        tAlpha              = 2.D-4        tAlpha              = 2.D-4
94        sBeta               = 7.4D-4        sBeta               = 7.4D-4
# Line 96  C--   Set default "physical" parameters Line 96  C--   Set default "physical" parameters
96        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
97        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
98        rigidLid            = .FALSE.        rigidLid            = .FALSE.
99          implicSurfPress     = 1.D0
100          implicDiv2DFlow     = 1.D0
101        hFacMin             = 1.D0        hFacMin             = 1.D0
102        hFacMinDrDefault    = 0.D0        hFacMinDrDefault    = 0.D0
103          exactConserv        = .FALSE.
104          uniformLin_PhiSurf  = .TRUE.
105          nonlinFreeSurf      = 0
106          hFacInf             = 1.D0
107          hFacSup             = 1.D0
108          staggerTimeStep     = .FALSE.
109        momViscosity        = .TRUE.        momViscosity        = .TRUE.
110        momAdvection        = .TRUE.        momAdvection        = .TRUE.
111        momForcing          = .TRUE.        momForcing          = .TRUE.
112        useCoriolis         = .TRUE.        useCoriolis         = .TRUE.
113        momPressureForcing  = .TRUE.        momPressureForcing  = .TRUE.
114        momStepping         = .TRUE.        momStepping         = .TRUE.
115          vectorInvariantMomentum = .FALSE.
116        tempStepping        = .TRUE.        tempStepping        = .TRUE.
117        saltStepping        = .TRUE.        saltStepping        = .TRUE.
118        metricTerms         = .FALSE.        tr1Stepping         = .FALSE.
119          metricTerms         = .TRUE.
120        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
121        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
       openBoundaries      = .FALSE.  
122        nonHydrostatic      = .FALSE.        nonHydrostatic      = .FALSE.
123        globalFiles         = .FALSE.        globalFiles         = .FALSE.
124        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
125          ivdc_kappa          = 0.D0
126          groundAtK1          = .FALSE.
127          zonal_filt_lat      = 90.
128          zonal_filt_sinpow   = 2
129          zonal_filt_cospow   = 2
130          bottomDragLinear    = 0.
131          bottomDragQuadratic = 0.
132          usePickupBeforeC35    = .FALSE.
133          debugMode             = .FALSE.
134          readPickupWithTracer  = .FALSE.
135          writePickupWithTracer = .FALSE.
136          tempAdvScheme       = 2
137          saltAdvScheme       = 2
138          tracerAdvScheme     = 2
139    
140    C--   Atmospheric physical parameters (e.g.: EOS)
141          atm_po =  1.D5
142          atm_cp =  1004.D0
143          atm_kappa = 2.D0 / 7.D0
144          Integr_GeoPot = 2
145    
146  C--   Elliptic solver parameters  C--   Elliptic solver parameters
147        cg2dMaxIters       = 150        cg2dMaxIters       = 150
148        cg2dTargetResidual = 1.D-7        cg2dTargetResidual = 1.D-7
149          cg2dTargetResWunit = -1.
150        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
151        cg2dpcOffDFac      = 0.51D0        cg2dpcOffDFac      = 0.51D0
152        cg3dMaxIters       = 150        cg3dMaxIters       = 150
# Line 128  C--   Time stepping parameters Line 158  C--   Time stepping parameters
158        nIter0            = 0        nIter0            = 0
159        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
160        nTimeSteps        = 0        nTimeSteps        = 0
161        endTime           = deltaT*float(nTimeSteps)        nEndIter          = nIter0+nTimeSteps
162          endTime           = deltaT*float(nEndIter)
163        abEps             = 0.01        abEps             = 0.01
164        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
165        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
166        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
167          diagFreq          = deltaT*0
168          monitorFreq       = -1.
169        taveFreq          = deltaT*0        taveFreq          = deltaT*0
170        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
171        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
# Line 140  C--   Time stepping parameters Line 173  C--   Time stepping parameters
173        nCheckLev         = 1        nCheckLev         = 1
174        checkPtSuff(1)    = 'ckptA'        checkPtSuff(1)    = 'ckptA'
175        checkPtSuff(2)    = 'ckptB'        checkPtSuff(2)    = 'ckptB'
176        cAdjFreq          = -1.D0        cAdjFreq          =  0.D0
177        tauCD             =  0.D0        tauCD             =  0.D0
178        tauThetaClimRelax =  0.D0        tauThetaClimRelax =  0.D0
179        tauSaltClimRelax  =  0.D0        tauSaltClimRelax  =  0.D0
180          tauTr1ClimRelax   =  0.D0
181        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
182        externForcingPeriod     = 0.        externForcingPeriod     = 0.
183        externForcingCycle      = 0.        externForcingCycle      = 0.
184    
185  C--   Input files  C--   Input files
186        bathyFile       = ' '        bathyFile       = ' '
187          topoFile        = ' '
188        hydrogSaltFile  = ' '        hydrogSaltFile  = ' '
189        hydrogThetaFile = ' '        hydrogThetaFile = ' '
190        zonalWindFile   = ' '        zonalWindFile   = ' '
# Line 158  C--   Input files Line 193  C--   Input files
193        saltClimFile    = ' '        saltClimFile    = ' '
194        EmPmRfile       = ' '        EmPmRfile       = ' '
195        surfQfile       = ' '        surfQfile       = ' '
196          surfQswfile     = ' '
197  C--   OBCS        uVelInitFile    = ' '
198        DO I=1,Nx        vVelInitFile    = ' '
199         OB_Jnorth(I)=0        pSurfInitFile   = ' '
200         OB_Jsouth(I)=0        dQdTFile        = ' '
       ENDDO  
       DO J=1,Ny  
        OB_Ieast(J)=0  
        OB_Iwest(J)=0  
       ENDDO  
201    
202  C  C
203        RETURN        RETURN
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22