/[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.3 by adcroft, Wed May 5 18:32:35 1999 UTC revision 1.13 by cnh, Sun Feb 4 14:38:48 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 44  C--   Grid parameters Line 45  C--   Grid parameters
45  C     Vertical gridding  C     Vertical gridding
46        rkFacDefault         = 1.D0        rkFacDefault         = 1.D0
47        horiVertRatio        = 1.D0        horiVertRatio        = 1.D0
48          Ro_SeaLevel = 0.
49        DO k=1,Nr        DO k=1,Nr
50         delRdefault(k) = 0.         delRdefault(k) = 0.
51        ENDDO        ENDDO
# Line 52  C     Horizontal gridding Line 54  C     Horizontal gridding
54  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
55        usingCartesianGrid = .TRUE.        usingCartesianGrid = .TRUE.
56        DO i=1,Nx        DO i=1,Nx
57         delX(i) = 0.         delX(i) = UNSET_RL
58        ENDDO        ENDDO
59        DO j=1,Ny        DO j=1,Ny
60         delY(j) = 0.         delY(j) = UNSET_RL
61        ENDDO        ENDDO
62  C     In spherical polar distances are in degrees  C     In spherical polar distances are in degrees
63        usingSphericalPolarGrid = .FALSE.        usingSphericalPolarGrid = .FALSE.
# Line 82  C--   Set default "physical" parameters Line 84  C--   Set default "physical" parameters
84        viscA4              = 0.D11        viscA4              = 0.D11
85        diffK4T             = 0.D11        diffK4T             = 0.D11
86        diffK4S             = 0.D11        diffK4S             = 0.D11
87        GMmaxslope          = 1.D-2        cosPower            = 0.
       GMlength            = 200.D3  
       GMalpha             = 0.D0  
       GMdepth             = 1000.D0  
       GMkbackground       = 0.D0  
       GMmaxval            = 2500.D0  
88        HeatCapacity_Cp     = 3994.D0        HeatCapacity_Cp     = 3994.D0
89        tAlpha              = 2.D-4        tAlpha              = 2.D-4
90        sBeta               = 7.4D-4        sBeta               = 7.4D-4
# Line 97  C--   Set default "physical" parameters Line 94  C--   Set default "physical" parameters
94        rigidLid            = .FALSE.        rigidLid            = .FALSE.
95        hFacMin             = 1.D0        hFacMin             = 1.D0
96        hFacMinDrDefault    = 0.D0        hFacMinDrDefault    = 0.D0
97          staggerTimeStep     = .FALSE.
98        momViscosity        = .TRUE.        momViscosity        = .TRUE.
99        momAdvection        = .TRUE.        momAdvection        = .TRUE.
100        momForcing          = .TRUE.        momForcing          = .TRUE.
# Line 107  C--   Set default "physical" parameters Line 105  C--   Set default "physical" parameters
105        saltStepping        = .TRUE.        saltStepping        = .TRUE.
106        metricTerms         = .FALSE.        metricTerms         = .FALSE.
107        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
108        openBoundaries      = .FALSE.        implicitViscosity   = .FALSE.
109        nonHydrostatic      = .FALSE.        nonHydrostatic      = .FALSE.
110        globalFiles         = .FALSE.        globalFiles         = .FALSE.
111        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
112          ivdc_kappa          = 0.D0
113          groundAtK1          = .FALSE.
114          nShap               = 0
115          zonal_filt_lat      = 90.
116          zonal_filt_sinpow   = 2
117          zonal_filt_cospow   = 4
118          bottomDragLinear    = 0.
119          bottomDragQuadratic = 0.
120    
121  C--   Elliptic solver parameters  C--   Elliptic solver parameters
122        cg2dMaxIters       = 150        cg2dMaxIters       = 150
# Line 126  C--   Time stepping parameters Line 132  C--   Time stepping parameters
132        nIter0            = 0        nIter0            = 0
133        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
134        nTimeSteps        = 0        nTimeSteps        = 0
135        endTime           = deltaT*float(nTimeSteps)        nEndIter          = nIter0+nTimeSteps
136          endTime           = deltaT*float(nEndIter)
137        abEps             = 0.01        abEps             = 0.01
138        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
139        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
# Line 138  C--   Time stepping parameters Line 145  C--   Time stepping parameters
145        nCheckLev         = 1        nCheckLev         = 1
146        checkPtSuff(1)    = 'ckptA'        checkPtSuff(1)    = 'ckptA'
147        checkPtSuff(2)    = 'ckptB'        checkPtSuff(2)    = 'ckptB'
148        cAdjFreq          = -1.D0        cAdjFreq          =  0.D0
149        tauCD             =  0.D0        tauCD             =  0.D0
150        tauThetaClimRelax =  0.D0        tauThetaClimRelax =  0.D0
151        tauSaltClimRelax  =  0.D0        tauSaltClimRelax  =  0.D0
# Line 156  C--   Input files Line 163  C--   Input files
163        saltClimFile    = ' '        saltClimFile    = ' '
164        EmPmRfile       = ' '        EmPmRfile       = ' '
165        surfQfile       = ' '        surfQfile       = ' '
166          surfQswfile     = ' '
167  C--   OBCS        uVelInitFile    = ' '
168        DO I=1,Nx        vVelInitFile    = ' '
169         OB_Jnorth(I)=0        pSurfInitFile   = ' '
        OB_Jsouth(I)=0  
       ENDDO  
       DO J=1,Ny  
        OB_Ieast(J)=0  
        OB_Iwest(J)=0  
       ENDDO  
170    
171  C  C
172        RETURN        RETURN
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22