/[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.27 by adcroft, Thu Sep 6 14:23:57 2001 UTC revision 1.48 by adcroft, Thu Jun 5 16:03:05 2003 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6    CBOP
7    C     !ROUTINE: SET_DEFAULTS
8    C     !INTERFACE:
9        SUBROUTINE SET_DEFAULTS(        SUBROUTINE SET_DEFAULTS(
10       O   viscArDefault, diffKrTDefault, diffKrSDefault,       O   viscArDefault, diffKrTDefault, diffKrSDefault,
11       O   hFacMinDrDefault, delRdefault, rkFacDefault,       O   hFacMinDrDefault, delRdefault, rkFacDefault,
12       I   myThid )       I   myThid )
 C     /==========================================================\  
 C     | SUBROUTINE SET_DEFAULTS                                  |  
 C     | o Routine to set model "parameters"                      |  
 C     |==========================================================|  
 C     | Notes:                                                   |  
 C     | ======                                                   |  
 C     | The present version of this routine is a place-holder.   |  
 C     | A production version needs to handle parameters from an  |  
 C     | external file and possibly reading in some initial field |  
 C     | values.                                                  |  
 C     \==========================================================/  
       IMPLICIT NONE  
13    
14    C     !DESCRIPTION: \bv
15    C     *==========================================================*
16    C     | SUBROUTINE SET_DEFAULTS                                  
17    C     | o Routine to set model "parameters"                      
18    C     *==========================================================*
19    C     | Notes:                                                    
20    C     | ======                                                    
21    C     | The present version of this routine is a place-holder.    
22    C     | A production version needs to handle parameters from an  
23    C     | external file and possibly reading in some initial field  
24    C     | values.                                                  
25    C     *==========================================================*
26    C     \ev
27    
28    C     !USES:
29          IMPLICIT NONE
30  C     === Global variables ===  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:
38  C     === Routine arguments ===  C     === Routine arguments ===
39  C     myThid - Number of this instance of INI_PARMS  C     myThid - Number of this instance of INI_PARMS
40        INTEGER myThid        INTEGER myThid
# Line 36  C     myThid - Number of this instance o Line 45  C     myThid - Number of this instance o
45        _RL delRDefault(Nr)        _RL delRDefault(Nr)
46        _RS rkFacDefault        _RS rkFacDefault
47    
48    C     !LOCAL VARIABLES:
49  C     === Local variables ===  C     === Local variables ===
50  C     K, I, J - Loop counters  C     K, I, J - Loop counters
51        INTEGER K, I, J        INTEGER K, I, J
52    CEOP
53    
54  C--   Grid parameters  C--   Grid parameters
55  C     Vertical gridding  C     Vertical gridding
# Line 48  C     Vertical gridding Line 59  C     Vertical gridding
59        DO k=1,Nr        DO k=1,Nr
60         delRdefault(k) = 0.         delRdefault(k) = 0.
61        ENDDO        ENDDO
62          DO k=1,Nr+1
63           delRc(k) = UNSET_RL
64          ENDDO
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
# Line 69  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 90  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'
       implicitFreeSurface = .TRUE.  
       rigidLid            = .FALSE.  
       implicSurfPress     = 1.D0  
       implicDiv2DFlow     = 1.D0  
115        hFacMin             = 1.D0        hFacMin             = 1.D0
116        hFacMinDrDefault    = 0.D0        hFacMinDrDefault    = 0.D0
       exactConserv        = .FALSE.  
       uniformLin_PhiSurf  = .TRUE.  
       nonlinFreeSurf      = 0  
       hFacInf             = 1.D0  
       hFacSup             = 1.D0  
117        staggerTimeStep     = .FALSE.        staggerTimeStep     = .FALSE.
118        momViscosity        = .TRUE.        momViscosity        = .TRUE.
119        momAdvection        = .TRUE.        momAdvection        = .TRUE.
# Line 114  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.
       zonal_filt_lat      = 90.  
       zonal_filt_sinpow   = 2  
       zonal_filt_cospow   = 2  
143        bottomDragLinear    = 0.        bottomDragLinear    = 0.
144        bottomDragQuadratic = 0.        bottomDragQuadratic = 0.
145        usePickupBeforeC35    = .FALSE.        usePickupBeforeC35    = .FALSE.
# Line 136  C--   Set default "physical" parameters Line 149  C--   Set default "physical" parameters
149        tempAdvScheme       = 2        tempAdvScheme       = 2
150        saltAdvScheme       = 2        saltAdvScheme       = 2
151        tracerAdvScheme     = 2        tracerAdvScheme     = 2
152          multiDimAdvection   = .TRUE.
153          useCDscheme         = .FALSE.
154          useEnergyConservingCoriolis = .FALSE.
155          useJamartWetPoints  = .FALSE.
156    
157    C--   Set (free)surface-related parameters
158          implicitFreeSurface = .TRUE.
159          rigidLid            = .FALSE.
160          implicSurfPress     = 1.D0
161          implicDiv2DFlow     = 1.D0
162          exactConserv        = .FALSE.
163          uniformLin_PhiSurf  = .TRUE.
164          nonlinFreeSurf      = 0
165          hFacInf             = 0.2 _d 0
166          hFacSup             = 2.0 _d 0
167          select_rStar        = 0
168          useRealFreshWaterFlux = .FALSE.
169          temp_EvPrRn = UNSET_RL
170          salt_EvPrRn = 0.
171          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 157  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 167  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 198  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.27  
changed lines
  Added in v.1.48
  ViewVC Help
Powered by ViewVC 1.1.22