/[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.1 by adcroft, Tue Dec 15 00:20:35 1998 UTC revision 1.25 by adcroft, Thu Aug 16 17:12:24 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    
6        SUBROUTINE SET_DEFAULTS(        SUBROUTINE SET_DEFAULTS(
7       O   viscArDefault, diffKrTDefault, diffKrSDefault,       O   viscArDefault, diffKrTDefault, diffKrSDefault,
8       O   hFacMinDrDefault, delRdefault,       O   hFacMinDrDefault, delRdefault, rkFacDefault,
9       I   myThid )       I   myThid )
10  C     /==========================================================\  C     /==========================================================\
11  C     | SUBROUTINE SET_DEFAULTS                                  |  C     | SUBROUTINE SET_DEFAULTS                                  |
# 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 34  C     myThid - Number of this instance o Line 34  C     myThid - Number of this instance o
34        _RL diffKrSDefault        _RL diffKrSDefault
35        _RL hFacMinDrDefault        _RL hFacMinDrDefault
36        _RL delRDefault(Nr)        _RL delRDefault(Nr)
37          _RS rkFacDefault
38    
39  C     === Local variables ===  C     === Local variables ===
40  C     K, I, J - Loop counters  C     K, I, J - Loop counters
# Line 41  C     K, I, J - Loop counters Line 42  C     K, I, J - Loop counters
42    
43  C--   Grid parameters  C--   Grid parameters
44  C     Vertical gridding  C     Vertical gridding
45        rkFac = 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) = 1.D2         delRdefault(k) = 0.
50        ENDDO        ENDDO
51    
52  C     Horizontal gridding  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) = 20.D0 * 1000.D0         delX(i) = UNSET_RL
57        ENDDO        ENDDO
58        DO j=1,Ny        DO j=1,Ny
59         delY(j) = 20.D0 * 1000.D0         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    = -5.0        phiMin               = 0.0
64        thetaMin  = 0.        thetaMin             = 0.
65        rSphere   = 6370. * 1.D3        rSphere              = 6370. * 1.D3
66        IF ( usingSphericalPolarGrid ) THEN  C     General curvilinear coordinate system
67         DO i=1,Nx        usingCurvilinearGrid = .FALSE.
         delX(i) = 1.D0  
        ENDDO  
        DO j=1,Ny  
         delY(j) = 1.D0  
        ENDDO  
       ENDIF  
 C  
68    
69  C--   Set default "physical" parameters  C--   Set default "physical" parameters
70        DO K =1,Nr        DO K =1,Nr
71         tRef(K) = 30.D0 - FLOAT( K )         tRef(K) = 30.D0 - FLOAT( K )
72        ENDDO        ENDDO
73        gravity  = 9.81D0        gravity             = 9.81D0
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        viscAh   = 1.D3        omega               = 2.D0 * PI / ( 3600.D0 * 24.D0 )
78        diffKhT  = 1.D3        viscAh              = 0.D3
79        diffKhS  = 1.D3        diffKhT             = 0.D3
80        viscArDefault = 1.D-3        diffKhS             = 0.D3
81        no_slip_sides = .TRUE.        viscArDefault       = 0.D-3
82        no_slip_bottom = .TRUE.        no_slip_sides       = .TRUE.
83        diffKrTDefault = 1.D-5        no_slip_bottom      = .TRUE.
84        diffKrSDefault = 1.D-5        diffKrTDefault      = 0.D-3
85        viscA4   = 0.        diffKrSDefault      = 0.D-3
86        diffK4T  = 0.        viscA4              = 0.D11
87        diffK4S  = 0.        diffK4T             = 0.D11
88        GMmaxslope   =   1.D-2        diffK4S             = 0.D11
89        GMlength     = 200.D3        cosPower            = 0.
90        GMalpha      = 0.D0        HeatCapacity_Cp     = 3994.D0
91        GMdepth      = 1000.D0        tAlpha              = 2.D-4
92        GMkbackground= 0.D0        sBeta               = 7.4D-4
93        GMmaxval     = 2500.D0        eosType             = 'LINEAR'
       HeatCapacity_Cp = 3994.D0  
       tAlpha       = 2.D-4  
       sBeta        = 7.4D-4  
       eosType      = 'LINEAR'  
94        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
95        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
96        rigidLid            = .FALSE.        rigidLid            = .FALSE.
97          implicSurfPress     = 1.D0
98          implicDiv2DFlow     = 1.D0
99        hFacMin             = 1.D0        hFacMin             = 1.D0
100        hFacMinDrDefault    = 0.D0        hFacMinDrDefault    = 0.D0
101          exactConserv        = .FALSE.
102          uniformLin_PhiSurf  = .TRUE.
103          nonlinFreeSurf      = 0
104          hFacInf             = 1.D0
105          hFacSup             = 1.D0
106          staggerTimeStep     = .FALSE.
107        momViscosity        = .TRUE.        momViscosity        = .TRUE.
108        momAdvection        = .TRUE.        momAdvection        = .TRUE.
109        momForcing          = .TRUE.        momForcing          = .TRUE.
110        useCoriolis         = .TRUE.        useCoriolis         = .TRUE.
111        momPressureForcing  = .TRUE.        momPressureForcing  = .TRUE.
112        momStepping         = .TRUE.        momStepping         = .TRUE.
113          vectorInvariantMomentum = .FALSE.
114        tempStepping        = .TRUE.        tempStepping        = .TRUE.
115        saltStepping        = .TRUE.        saltStepping        = .TRUE.
116        metricTerms         = .FALSE.        tr1Stepping         = .FALSE.
117          metricTerms         = .TRUE.
118        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
119        openBoundaries      = .FALSE.        implicitViscosity   = .FALSE.
120          nonHydrostatic      = .FALSE.
121          globalFiles         = .FALSE.
122          allowFreezing       = .FALSE.
123          ivdc_kappa          = 0.D0
124          groundAtK1          = .FALSE.
125          zonal_filt_lat      = 90.
126          zonal_filt_sinpow   = 2
127          zonal_filt_cospow   = 2
128          bottomDragLinear    = 0.
129          bottomDragQuadratic = 0.
130          usePickupBeforeC35    = .FALSE.
131          debugMode             = .FALSE.
132          readPickupWithTracer  = .FALSE.
133          writePickupWithTracer = .FALSE.
134          tempAdvScheme       = 2
135          saltAdvScheme       = 2
136          tracerAdvScheme     = 2
137    
138    C--   Atmospheric physical parameters (e.g.: EOS)
139          atm_po =  1.D5
140          atm_cp =  1004.D0
141          atm_kappa = 2.D0 / 7.D0
142          Integr_GeoPot = 2
143    
144  C--   Elliptic solver parameters  C--   Elliptic solver parameters
145        cg2dMaxIters       = 150        cg2dMaxIters       = 150
146        cg2dTargetResidual = 1.D-7        cg2dTargetResidual = 1.D-7
147          cg2dTargetResWunit = -1.
148        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
149          cg2dpcOffDFac      = 0.51D0
150        cg3dMaxIters       = 150        cg3dMaxIters       = 150
151        cg3dTargetResidual = 1.D-7        cg3dTargetResidual = 1.D-7
152        cg3dChkResFreq     = 1        cg3dChkResFreq     = 1
       cg2dpcOffDFac      = 0.51D0  
153    
154  C--   Time stepping parameters  C--   Time stepping parameters
155        deltaT            = 0.        deltaT            = 0.
156        nIter0            = 0        nIter0            = 0
157        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
158        nTimeSteps        = 100        nTimeSteps        = 0
159        endTime           = deltaT*float(nTimeSteps)        nEndIter          = nIter0+nTimeSteps
160          endTime           = deltaT*float(nEndIter)
161        abEps             = 0.01        abEps             = 0.01
162        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
163        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
164        dumpFreq          = deltaT*50        dumpFreq          = deltaT*0
165          monitorFreq       = -1.
166        taveFreq          = deltaT*0        taveFreq          = deltaT*0
167        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
168        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
# Line 143  C--   Time stepping parameters Line 170  C--   Time stepping parameters
170        nCheckLev         = 1        nCheckLev         = 1
171        checkPtSuff(1)    = 'ckptA'        checkPtSuff(1)    = 'ckptA'
172        checkPtSuff(2)    = 'ckptB'        checkPtSuff(2)    = 'ckptB'
173        cAdjFreq          = -1.D0        cAdjFreq          =  0.D0
174        tauCD             =  0.D0        tauCD             =  0.D0
175        tauThetaClimRelax =  0.D0        tauThetaClimRelax =  0.D0
176        tauSaltClimRelax  =  0.D0        tauSaltClimRelax  =  0.D0
177          tauTr1ClimRelax   =  0.D0
178        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
179        externForcingPeriod     = 0.        externForcingPeriod     = 0.
180        externForcingCycle      = 0.        externForcingCycle      = 0.
181    
182  C--   Input files  C--   Input files
183        bathyFile       = ' '        bathyFile       = ' '
184          topoFile        = ' '
185        hydrogSaltFile  = ' '        hydrogSaltFile  = ' '
186        hydrogThetaFile = ' '        hydrogThetaFile = ' '
187        zonalWindFile   = ' '        zonalWindFile   = ' '
188        meridWindFile   = ' '        meridWindFile   = ' '
189        thetaClimFile   = ' '        thetaClimFile   = ' '
190        saltClimFile    = ' '        saltClimFile    = ' '
191          EmPmRfile       = ' '
192  C--   OBCS        surfQfile       = ' '
193        DO I=1,Nx        surfQswfile     = ' '
194         OB_Jnorth(I)=0        uVelInitFile    = ' '
195         OB_Jsouth(I)=0        vVelInitFile    = ' '
196        ENDDO        pSurfInitFile   = ' '
197        DO J=1,Ny        dQdTFile        = ' '
        OB_Ieast(J)=0  
        OB_Iwest(J)=0  
       ENDDO  
198    
199  C  C
200        RETURN        RETURN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22