/[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.150 by jmc, Fri Nov 12 03:16:16 2010 UTC revision 1.172 by gforget, Sat Jun 28 22:34:16 2014 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
# Line 20  C     !USES: Line 21  C     !USES:
21  #include "SIZE.h"  #include "SIZE.h"
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
23  #include "PARAMS.h"  #include "PARAMS.h"
24  Cml#include "EOS.h"  #ifdef ALLOW_EXCH2
25  c#include "GRID.h"  # include "W2_EXCH2_SIZE.h"
26    #endif /* ALLOW_EXCH2 */
27    #include "SET_GRID.h"
28    
29  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
30  C     myThid :: my Thread Id. Number  C     myThid :: my Thread Id. Number
# Line 55  C-    vertical profile Line 58  C-    vertical profile
58        tRefFile            = ' '        tRefFile            = ' '
59        sRefFile            = ' '        sRefFile            = ' '
60        rhoRefFile          = ' '        rhoRefFile          = ' '
61          thetaConst          = UNSET_RL
62          DO k=1,Nr
63           tRef(k)            = UNSET_RL
64           sRef(k)            = UNSET_RL
65          ENDDO
66    
67  C-    Horizontal gridding  C-    Horizontal gridding
68        delXFile            = ' '        delXFile            = ' '
# Line 63  C-    Horizontal gridding Line 71  C-    Horizontal gridding
71        deepAtmosphere      = .FALSE.        deepAtmosphere      = .FALSE.
72        xgOrigin            = UNSET_RL        xgOrigin            = UNSET_RL
73        ygOrigin            = UNSET_RL        ygOrigin            = UNSET_RL
74        DO i=1,Nx        DO i=1,grid_maxNx
75         delX(i)            = UNSET_RL         delX(i)            = UNSET_RL
76        ENDDO        ENDDO
77        DO j=1,Ny        DO j=1,grid_maxNy
78         delY(j)            = UNSET_RL         delY(j)            = UNSET_RL
79        ENDDO        ENDDO
80  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
81        usingCartesianGrid  = .FALSE.        usingCartesianGrid  = .FALSE.
82  C     In spherical polar distances are in degrees  C     In spherical polar distances are in degrees
83        usingSphericalPolarGrid = .FALSE.        usingSphericalPolarGrid = .FALSE.
84        rSphere             = 6370. _d 3        rSphere             = UNSET_RL
85  C     General curvilinear coordinate system  C     General curvilinear coordinate system
86        usingCurvilinearGrid= .FALSE.        usingCurvilinearGrid= .FALSE.
87          radius_fromHorizGrid= UNSET_RL
88          hasWetCSCorners     = .FALSE.
89  C     General cylindrical coordinate system  C     General cylindrical coordinate system
90        usingCylindricalGrid= .FALSE.        usingCylindricalGrid= .FALSE.
91  C     Coriolis map:  C     Coriolis map:
# Line 90  C     grid rotation Line 100  C     grid rotation
100  C--   Set default "physical" parameters  C--   Set default "physical" parameters
101        nh_Am2              = 1. _d 0        nh_Am2              = 1. _d 0
102        gravity             = 9.81 _d 0        gravity             = 9.81 _d 0
103          gBaro               = UNSET_RL
104        rhoNil              = 999.8 _d 0        rhoNil              = 999.8 _d 0
105          rhoConst            = UNSET_RL
106  C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)  C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)
107  C         (so that the default produces same results as before)  C         (so that the default produces same results as before)
108  c     rhoConstFresh       = 999.8 _d 0  c     rhoConstFresh       = 999.8 _d 0
109          rhoConstFresh       = UNSET_RL
110          convertFW2Salt      = UNSET_RL
111        f0                  = 1. _d -4        f0                  = 1. _d -4
112        beta                = 1. _d -11        beta                = 1. _d -11
113        fPrime              = 0. _d 0        fPrime              = 0. _d 0
114  C-    Earth rotation period is 86400*365.25/366.25 (use to be 1.day)  C-    Earth rotation period is 86400*365.25/366.25 (use to be 1.day)
115        rotationPeriod      = 86164. _d 0        rotationPeriod      = 86164. _d 0
116          omega               = UNSET_RL
117    C-    viscosity and diffusivity default value:
118        viscAh              = 0. _d 3        viscAh              = 0. _d 3
119          smag3D_coeff        = 1. _d -2
120        viscAhGrid          = 0. _d 0        viscAhGrid          = 0. _d 0
121        viscAhGridMin       = 0. _d 0        viscAhGridMin       = 0. _d 0
122        viscAhGridMax       = 1. _d 21        viscAhGridMax       = 1. _d 21
# Line 108  C-    Earth rotation period is 86400*365 Line 125  C-    Earth rotation period is 86400*365
125        viscC2leith         = 0. _d 0        viscC2leith         = 0. _d 0
126        viscC2leithD        = 0. _d 0        viscC2leithD        = 0. _d 0
127        viscC2smag          = 0. _d 0        viscC2smag          = 0. _d 0
       diffKhT             = 0. _d 3  
       diffKhS             = 0. _d 3  
128        viscArDefault       = 0. _d -3        viscArDefault       = 0. _d -3
129          viscFacAdj          = 1. _d 0
130        no_slip_sides       = .TRUE.        no_slip_sides       = .TRUE.
131        no_slip_bottom      = .TRUE.        no_slip_bottom      = .TRUE.
132        sideDragFactor      = 2. _d 0        sideDragFactor      = 2. _d 0
133        bottomDragLinear    = 0.        bottomDragLinear    = 0.
134        bottomDragQuadratic = 0.        bottomDragQuadratic = 0.
       smoothAbsFuncRange  = 0. _d 0  
       diffKrTDefault      = 0. _d -3  
       diffKrSDefault      = 0. _d -3  
       diffKrBL79surf      = 0. _d 0  
       diffKrBL79deep      = 0. _d 0  
       diffKrBL79scl       = 200. _d 0  
       diffKrBL79Ho        = -2000. _d 0  
       BL79LatVary         = 30  
       diffKrBLEQsurf      = UNSET_RL  
       diffKrBLEQdeep      = UNSET_RL  
       diffKrBLEQscl       = UNSET_RL  
       diffKrBLEQHo        = UNSET_RL  
135        viscA4              = 0. _d 11        viscA4              = 0. _d 11
136        viscA4Grid          = 0. _d 0        viscA4Grid          = 0. _d 0
137        viscA4GridMax       = 1. _d 21        viscA4GridMax       = 1. _d 21
# Line 137  C-    Earth rotation period is 86400*365 Line 141  C-    Earth rotation period is 86400*365
141        viscC4leith         = 0. _d 0        viscC4leith         = 0. _d 0
142        viscC4leithD        = 0. _d 0        viscC4leithD        = 0. _d 0
143        viscC4smag          = 0. _d 0        viscC4smag          = 0. _d 0
144          DO k=1,Nr
145           viscArNr(k)        = UNSET_RL
146          ENDDO
147          cosPower            = 0.
148          diffKhT             = 0. _d 3
149          diffKhS             = 0. _d 3
150        diffK4T             = 0. _d 11        diffK4T             = 0. _d 11
151        diffK4S             = 0. _d 11        diffK4S             = 0. _d 11
152        cosPower            = 0.        diffKrTDefault      = 0. _d -3
153          diffKrSDefault      = 0. _d -3
154          diffKrBL79surf      = 0. _d 0
155          diffKrBL79deep      = 0. _d 0
156          diffKrBL79scl       = 200. _d 0
157          diffKrBL79Ho        = -2000. _d 0
158          BL79LatVary         = 30.
159          diffKrBLEQsurf      = UNSET_RL
160          diffKrBLEQdeep      = UNSET_RL
161          diffKrBLEQscl       = UNSET_RL
162          diffKrBLEQHo        = UNSET_RL
163          DO k=1,Nr
164           diffKrNrT(k)       = UNSET_RL
165           diffKrNrS(k)       = UNSET_RL
166           diffKr4T(k)        = 0.
167           diffKr4S(k)        = 0.
168          ENDDO
169        HeatCapacity_Cp     = 3994. _d 0        HeatCapacity_Cp     = 3994. _d 0
 Cml      tAlpha              = 2. _d -4  
 Cml      sBeta               = 7.4 _d -4  
170        eosType             = 'LINEAR'        eosType             = 'LINEAR'
171        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
172          smoothAbsFuncRange  = 0. _d 0
173        hFacMin             = 1. _d 0        hFacMin             = 1. _d 0
174        hFacMinDrDefault    = 0. _d 0        hFacMinDrDefault    = 0. _d 0
175        implicitIntGravWave = .FALSE.        implicitIntGravWave = .FALSE.
176        staggerTimeStep     = .FALSE.        staggerTimeStep     = .FALSE.
177          doResetHFactors     = .FALSE.
178        momViscosity        = .TRUE.        momViscosity        = .TRUE.
179        momAdvection        = .TRUE.        momAdvection        = .TRUE.
180        momForcing          = .TRUE.        momForcing          = .TRUE.
# Line 162  Cml      sBeta               = 7.4 _d -4 Line 188  Cml      sBeta               = 7.4 _d -4
188        saltStepping        = .TRUE.        saltStepping        = .TRUE.
189        saltAdvection       = .TRUE.        saltAdvection       = .TRUE.
190        saltForcing         = .TRUE.        saltForcing         = .TRUE.
191          addFrictionHeating  = .FALSE.
192        metricTerms         = .TRUE.        metricTerms         = .TRUE.
193        useNHMTerms         = .FALSE.        useNHMTerms         = .FALSE.
194          useSmag3D           = .FALSE.
195        useFullLeith        = .FALSE.        useFullLeith        = .FALSE.
196        useAreaViscLength   = .FALSE.        useAreaViscLength   = .FALSE.
197        useStrainTensionVisc= .FALSE.        useStrainTensionVisc= .FALSE.
# Line 176  Cml      sBeta               = 7.4 _d -4 Line 204  Cml      sBeta               = 7.4 _d -4
204        quasiHydrostatic    = .FALSE.        quasiHydrostatic    = .FALSE.
205        globalFiles         = .FALSE.        globalFiles         = .FALSE.
206        useSingleCpuIO      = .FALSE.        useSingleCpuIO      = .FALSE.
207          useSingleCpuInput   = .FALSE.
208        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
       useOldFreezing      = .FALSE.  
209        ivdc_kappa          = 0. _d 0        ivdc_kappa          = 0. _d 0
210        hMixCriteria        = -.8 _d 0        hMixCriteria        = -.8 _d 0
211        dRhoSmall           = 1. _d -6        dRhoSmall           = 1. _d -6
212        hMixSmooth          = 0. _d 0        hMixSmooth          = 0. _d 0
213        usePickupBeforeC54    = .FALSE.        usePickupBeforeC54    = .FALSE.
       debugMode             = .FALSE.  
214        tempAdvScheme       = 2        tempAdvScheme       = 2
215        saltAdvScheme       = 2        saltAdvScheme       = 2
216        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
# Line 198  Cml      sBeta               = 7.4 _d -4 Line 225  Cml      sBeta               = 7.4 _d -4
225        useAbsVorticity     = .FALSE.        useAbsVorticity     = .FALSE.
226        upwindShear         = .FALSE.        upwindShear         = .FALSE.
227        selectKEscheme      = 0        selectKEscheme      = 0
228        debugLevel          = debLevA        IF ( debugMode ) THEN
229        inAdMode            = .FALSE.          debugLevel        = debLevD
230        inAdExact           = .TRUE.        ELSE
231            debugLevel        = debLevB
232    #ifdef ALLOW_AUTODIFF
233            debugLevel        = debLevA
234    #endif
235          ENDIF
236    
237  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
238        implicitFreeSurface = .FALSE.        implicitFreeSurface = .FALSE.
# Line 224  C--   Set (free)surface-related paramete Line 256  C--   Set (free)surface-related paramete
256        balanceEmPmR        = .FALSE.        balanceEmPmR        = .FALSE.
257        balanceQnet         = .FALSE.        balanceQnet         = .FALSE.
258        balancePrintMean    = .FALSE.        balancePrintMean    = .FALSE.
259          balanceThetaClimRelax = .FALSE.
260          balanceSaltClimRelax  = .FALSE.
261    
262  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
263        celsius2K = 273.16 _d 0        celsius2K = 273.15 _d 0
264        atm_Po =  1. _d 5        atm_Po =  1. _d 5
265        atm_Cp = 1004. _d 0        atm_Cp = 1004. _d 0
266        atm_Rd = UNSET_RL        atm_Rd = UNSET_RL
# Line 240  C--   Elliptic solver parameters Line 274  C--   Elliptic solver parameters
274        cg2dTargetResidual = 1. _d -7        cg2dTargetResidual = 1. _d -7
275        cg2dTargetResWunit = -1.        cg2dTargetResWunit = -1.
276        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
277          cg2dUseMinResSol   = UNSET_I
278        cg2dpcOffDFac      = 0.51 _d 0        cg2dpcOffDFac      = 0.51 _d 0
279        cg2dPreCondFreq    = 1        cg2dPreCondFreq    = 1
280        cg3dMaxIters       = 150        cg3dMaxIters       = 150
# Line 249  C--   Elliptic solver parameters Line 284  C--   Elliptic solver parameters
284    
285  C--   Time stepping parameters  C--   Time stepping parameters
286        deltaT            = 0. _d 0        deltaT            = 0. _d 0
287        deltaTmom         = 0. _d 0        deltaTMom         = 0. _d 0
288        deltaTfreesurf    = 0. _d 0        deltaTFreeSurf    = 0. _d 0
289        DO k=1,Nr        DO k=1,Nr
290          dTtracerLev(k)  = 0. _d 0          dTtracerLev(k)  = 0. _d 0
291        ENDDO        ENDDO
292        baseTime          = 0. _d 0        baseTime          = 0. _d 0
293        nIter0            = 0        nIter0            = -1
294        startTime         = deltaT*float(nIter0)        startTime         = UNSET_RL
295        pickupSuff        = ' '        pickupSuff        = ' '
296        pickupStrictlyMatch = .TRUE.        pickupStrictlyMatch = .TRUE.
297        nTimeSteps        = 0        nTimeSteps        = 0
298        nEndIter          = nIter0+nTimeSteps        nEndIter          = 0
299        endTime           = deltaT*float(nEndIter)        endTime           = 0. _d 0
300        momForcingOutAB   = UNSET_I        momForcingOutAB   = UNSET_I
301        tracForcingOutAB  = UNSET_I        tracForcingOutAB  = UNSET_I
302        momDissip_In_AB   = .TRUE.        momDissip_In_AB   = .TRUE.
# Line 311  C--   Time stepping parameters Line 346  C--   Time stepping parameters
346  C--   Input files  C--   Input files
347        bathyFile       = ' '        bathyFile       = ' '
348        topoFile        = ' '        topoFile        = ' '
349          addWwallFile    = ' '
350          addSwallFile    = ' '
351        hydrogSaltFile  = ' '        hydrogSaltFile  = ' '
352        hydrogThetaFile = ' '        hydrogThetaFile = ' '
353        maskIniTemp     = .TRUE.        maskIniTemp     = .TRUE.
# Line 334  C--   Input files Line 371  C--   Input files
371        uVelInitFile    = ' '        uVelInitFile    = ' '
372        vVelInitFile    = ' '        vVelInitFile    = ' '
373        pSurfInitFile   = ' '        pSurfInitFile   = ' '
374        dQdTFile        = ' '        pLoadFile       = ' '
375        ploadFile       = ' '        addMassFile     = ' '
376        eddyPsiXFile    = ' '        eddyPsiXFile    = ' '
377        eddyPsiYFile    = ' '        eddyPsiYFile    = ' '
378        lambdaThetaFile = ' '        lambdaThetaFile = ' '
379        lambdaSaltFile  = ' '        lambdaSaltFile  = ' '
380        mdsioLocalDir   = ' '        mdsioLocalDir   = ' '
381          adTapeDir       = ' '
382        the_run_name    = ' '        the_run_name    = ' '
383    
384        RETURN        RETURN
Legend:
Removed from v.1.150  
changed lines
  Added in v.1.172
  ViewVC Help
Powered by ViewVC 1.1.22