/[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.2 by adcroft, Mon Mar 22 15:54:04 1999 UTC
# Line 4  C $Header$ Line 4  C $Header$
4    
5        SUBROUTINE SET_DEFAULTS(        SUBROUTINE SET_DEFAULTS(
6       O   viscArDefault, diffKrTDefault, diffKrSDefault,       O   viscArDefault, diffKrTDefault, diffKrSDefault,
7       O   hFacMinDrDefault, delRdefault,       O   hFacMinDrDefault, delRdefault, rkFacDefault,
8       I   myThid )       I   myThid )
9  C     /==========================================================\  C     /==========================================================\
10  C     | SUBROUTINE SET_DEFAULTS                                  |  C     | SUBROUTINE SET_DEFAULTS                                  |
# 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        DO k=1,Nr        DO k=1,Nr
48         delRdefault(k) = 1.D2         delRdefault(k) = 0.
49        ENDDO        ENDDO
50    
51  C     Horizontal gridding  C     Horizontal gridding
52  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
53        usingCartesianGrid = .TRUE.        usingCartesianGrid = .TRUE.
54        DO i=1,Nx        DO i=1,Nx
55         delX(i) = 20.D0 * 1000.D0         delX(i) = 0.
56        ENDDO        ENDDO
57        DO j=1,Ny        DO j=1,Ny
58         delY(j) = 20.D0 * 1000.D0         delY(j) = 0.
59        ENDDO        ENDDO
60  C     In spherical polar distances are in degrees  C     In spherical polar distances are in degrees
61        usingSphericalPolarGrid = .FALSE.        usingSphericalPolarGrid = .FALSE.
62        phiMin    = -5.0        phiMin               = 0.0
63        thetaMin  = 0.        thetaMin             = 0.
64        rSphere   = 6370. * 1.D3        rSphere              = 6370. * 1.D3
       IF ( usingSphericalPolarGrid ) THEN  
        DO i=1,Nx  
         delX(i) = 1.D0  
        ENDDO  
        DO j=1,Ny  
         delY(j) = 1.D0  
        ENDDO  
       ENDIF  
 C  
65    
66  C--   Set default "physical" parameters  C--   Set default "physical" parameters
67        DO K =1,Nr        DO K =1,Nr
68         tRef(K) = 30.D0 - FLOAT( K )         tRef(K) = 30.D0 - FLOAT( K )
69        ENDDO        ENDDO
70        gravity  = 9.81D0        gravity             = 9.81D0
71        rhoNil   = 999.8D0        rhoNil              = 999.8D0
72        f0       = 1.D-4        f0                  = 1.D-4
73        beta     = 1.D-11        beta                = 1.D-11
74        viscAh   = 1.D3        viscAh              = 0.D3
75        diffKhT  = 1.D3        diffKhT             = 0.D3
76        diffKhS  = 1.D3        diffKhS             = 0.D3
77        viscArDefault = 1.D-3        viscArDefault       = 0.D-3
78        no_slip_sides = .TRUE.        no_slip_sides       = .TRUE.
79        no_slip_bottom = .TRUE.        no_slip_bottom      = .TRUE.
80        diffKrTDefault = 1.D-5        diffKrTDefault      = 0.D-3
81        diffKrSDefault = 1.D-5        diffKrSDefault      = 0.D-3
82        viscA4   = 0.        viscA4              = 0.D11
83        diffK4T  = 0.        diffK4T             = 0.D11
84        diffK4S  = 0.        diffK4S             = 0.D11
85        GMmaxslope   =   1.D-2        GMmaxslope          = 1.D-2
86        GMlength     = 200.D3        GMlength            = 200.D3
87        GMalpha      = 0.D0        GMalpha             = 0.D0
88        GMdepth      = 1000.D0        GMdepth             = 1000.D0
89        GMkbackground= 0.D0        GMkbackground       = 0.D0
90        GMmaxval     = 2500.D0        GMmaxval            = 2500.D0
91        HeatCapacity_Cp = 3994.D0        HeatCapacity_Cp     = 3994.D0
92        tAlpha       = 2.D-4        tAlpha              = 2.D-4
93        sBeta        = 7.4D-4        sBeta               = 7.4D-4
94        eosType      = 'LINEAR'        eosType             = 'LINEAR'
95        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
96        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
97        rigidLid            = .FALSE.        rigidLid            = .FALSE.
# Line 116  C--   Set default "physical" parameters Line 108  C--   Set default "physical" parameters
108        metricTerms         = .FALSE.        metricTerms         = .FALSE.
109        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
110        openBoundaries      = .FALSE.        openBoundaries      = .FALSE.
111          nonHydrostatic      = .FALSE.
112    
113  C--   Elliptic solver parameters  C--   Elliptic solver parameters
114        cg2dMaxIters       = 150        cg2dMaxIters       = 150
115        cg2dTargetResidual = 1.D-7        cg2dTargetResidual = 1.D-7
116        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
117          cg2dpcOffDFac      = 0.51D0
118        cg3dMaxIters       = 150        cg3dMaxIters       = 150
119        cg3dTargetResidual = 1.D-7        cg3dTargetResidual = 1.D-7
120        cg3dChkResFreq     = 1        cg3dChkResFreq     = 1
       cg2dpcOffDFac      = 0.51D0  
121    
122  C--   Time stepping parameters  C--   Time stepping parameters
123        deltaT            = 0.        deltaT            = 0.
124        nIter0            = 0        nIter0            = 0
125        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
126        nTimeSteps        = 100        nTimeSteps        = 0
127        endTime           = deltaT*float(nTimeSteps)        endTime           = deltaT*float(nTimeSteps)
128        abEps             = 0.01        abEps             = 0.01
129        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
130        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
131        dumpFreq          = deltaT*50        dumpFreq          = deltaT*0
132        taveFreq          = deltaT*0        taveFreq          = deltaT*0
133        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
134        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
# Line 159  C--   Input files Line 152  C--   Input files
152        meridWindFile   = ' '        meridWindFile   = ' '
153        thetaClimFile   = ' '        thetaClimFile   = ' '
154        saltClimFile    = ' '        saltClimFile    = ' '
155          EmPmRfile       = ' '
156          surfQfile       = ' '
157    
158  C--   OBCS  C--   OBCS
159        DO I=1,Nx        DO I=1,Nx
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22