/[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.23 by heimbach, Fri Jul 13 14:26:57 2001 UTC revision 1.36 by mlosch, Wed Jul 31 16:38:31 2002 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  #include "GRID.h"  #include "GRID.h"
35    
36    C     !INPUT/OUTPUT PARAMETERS:
37  C     === Routine arguments ===  C     === Routine arguments ===
38  C     myThid - Number of this instance of INI_PARMS  C     myThid - Number of this instance of INI_PARMS
39        INTEGER myThid        INTEGER myThid
# Line 36  C     myThid - Number of this instance o Line 44  C     myThid - Number of this instance o
44        _RL delRDefault(Nr)        _RL delRDefault(Nr)
45        _RS rkFacDefault        _RS rkFacDefault
46    
47    C     !LOCAL VARIABLES:
48  C     === Local variables ===  C     === Local variables ===
49  C     K, I, J - Loop counters  C     K, I, J - Loop counters
50        INTEGER K, I, J        INTEGER K, I, J
51    CEOP
52    
53  C--   Grid parameters  C--   Grid parameters
54  C     Vertical gridding  C     Vertical gridding
# Line 48  C     Vertical gridding Line 58  C     Vertical gridding
58        DO k=1,Nr        DO k=1,Nr
59         delRdefault(k) = 0.         delRdefault(k) = 0.
60        ENDDO        ENDDO
61          DO k=1,Nr+1
62           delRc(k) = UNSET_RL
63          ENDDO
64    
65  C     Horizontal gridding  C     Horizontal gridding
66  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
# Line 76  C--   Set default "physical" parameters Line 89  C--   Set default "physical" parameters
89        beta                = 1.D-11        beta                = 1.D-11
90        omega               = 2.D0 * PI / ( 3600.D0 * 24.D0 )        omega               = 2.D0 * PI / ( 3600.D0 * 24.D0 )
91        viscAh              = 0.D3        viscAh              = 0.D3
92          viscAstrain         = 0.D3
93          viscAtension        = 0.D3
94        diffKhT             = 0.D3        diffKhT             = 0.D3
95        diffKhS             = 0.D3        diffKhS             = 0.D3
96        viscArDefault       = 0.D-3        viscArDefault       = 0.D-3
# Line 92  C--   Set default "physical" parameters Line 107  C--   Set default "physical" parameters
107        sBeta               = 7.4D-4        sBeta               = 7.4D-4
108        eosType             = 'LINEAR'        eosType             = 'LINEAR'
109        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
       implicitFreeSurface = .TRUE.  
       rigidLid            = .FALSE.  
       implicSurfPress     = 1.D0  
       implicDiv2DFlow     = 1.D0  
110        hFacMin             = 1.D0        hFacMin             = 1.D0
111        hFacMinDrDefault    = 0.D0        hFacMinDrDefault    = 0.D0
       exactConserv        = .FALSE.  
       uniformLin_PhiSurf  = .TRUE.  
       nonlinFreeSurf      = 0  
       hFacInf             = 1.D0  
       hFacSup             = 1.D0  
112        staggerTimeStep     = .FALSE.        staggerTimeStep     = .FALSE.
113        momViscosity        = .TRUE.        momViscosity        = .TRUE.
114        momAdvection        = .TRUE.        momAdvection        = .TRUE.
# Line 110  C--   Set default "physical" parameters Line 116  C--   Set default "physical" parameters
116        useCoriolis         = .TRUE.        useCoriolis         = .TRUE.
117        momPressureForcing  = .TRUE.        momPressureForcing  = .TRUE.
118        momStepping         = .TRUE.        momStepping         = .TRUE.
119          vectorInvariantMomentum = .FALSE.
120        tempStepping        = .TRUE.        tempStepping        = .TRUE.
121          tempAdvection       = .TRUE.
122          tempForcing         = .TRUE.
123        saltStepping        = .TRUE.        saltStepping        = .TRUE.
124          saltAdvection       = .TRUE.
125          saltForcing         = .TRUE.
126        tr1Stepping         = .FALSE.        tr1Stepping         = .FALSE.
127        metricTerms         = .TRUE.        metricTerms         = .TRUE.
128        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
# Line 121  C--   Set default "physical" parameters Line 132  C--   Set default "physical" parameters
132        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
133        ivdc_kappa          = 0.D0        ivdc_kappa          = 0.D0
134        groundAtK1          = .FALSE.        groundAtK1          = .FALSE.
       zonal_filt_lat      = 90.  
       zonal_filt_sinpow   = 2  
       zonal_filt_cospow   = 2  
135        bottomDragLinear    = 0.        bottomDragLinear    = 0.
136        bottomDragQuadratic = 0.        bottomDragQuadratic = 0.
137        usePickupBeforeC35    = .FALSE.        usePickupBeforeC35    = .FALSE.
138        debugMode             = .FALSE.        debugMode             = .FALSE.
139        readPickupWithTracer  = .FALSE.        readPickupWithTracer  = .FALSE.
140        writePickupWithTracer = .FALSE.        writePickupWithTracer = .FALSE.
141          tempAdvScheme       = 2
142          saltAdvScheme       = 2
143          tracerAdvScheme     = 2
144          multiDimAdvection   = .TRUE.
145          useEnergyConservingCoriolis = .FALSE.
146          useJamartWetPoints  = .FALSE.
147    
148    C--   Set (free)surface-related parameters
149          implicitFreeSurface = .TRUE.
150          rigidLid            = .FALSE.
151          implicSurfPress     = 1.D0
152          implicDiv2DFlow     = 1.D0
153          exactConserv        = .FALSE.
154          uniformLin_PhiSurf  = .TRUE.
155          nonlinFreeSurf      = 0
156          hFacInf             = 1.D0
157          hFacSup             = 1.D0
158          useRealFreshWaterFlux = .FALSE.
159          temp_EvPrRn = UNSET_RL
160          salt_EvPrRn = 0.
161          trac_EvPrRn = UNSET_RL
162    
163  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
164        atm_po =  1.D5        atm_po =  1.D5
# Line 154  C--   Time stepping parameters Line 183  C--   Time stepping parameters
183        nTimeSteps        = 0        nTimeSteps        = 0
184        nEndIter          = nIter0+nTimeSteps        nEndIter          = nIter0+nTimeSteps
185        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
186          forcing_In_AB     = .TRUE.
187        abEps             = 0.01        abEps             = 0.01
188        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
189        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
190        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
191          diagFreq          = deltaT*0
192        monitorFreq       = -1.        monitorFreq       = -1.
193        taveFreq          = deltaT*0        taveFreq          = deltaT*0
194        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
# Line 191  C--   Input files Line 222  C--   Input files
222        vVelInitFile    = ' '        vVelInitFile    = ' '
223        pSurfInitFile   = ' '        pSurfInitFile   = ' '
224        dQdTFile        = ' '        dQdTFile        = ' '
225          ploadFile       = ' '
226    
227  C  C
228        RETURN        RETURN
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.36
  ViewVC Help
Powered by ViewVC 1.1.22