/[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.10 by adcroft, Wed Jun 21 19:18:59 2000 UTC revision 1.21 by adcroft, Wed Jun 6 14:55:45 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 24  C     === Global variables === Line 25  C     === Global variables ===
25  #include "EEPARAMS.h"  #include "EEPARAMS.h"
26  #include "PARAMS.h"  #include "PARAMS.h"
27  #include "GRID.h"  #include "GRID.h"
 #include "CG2D.h"  
28    
29  C     === Routine arguments ===  C     === Routine arguments ===
30  C     myThid - Number of this instance of INI_PARMS  C     myThid - Number of this instance of INI_PARMS
# Line 44  C--   Grid parameters Line 44  C--   Grid parameters
44  C     Vertical gridding  C     Vertical gridding
45        rkFacDefault         = 1.D0        rkFacDefault         = 1.D0
46        horiVertRatio        = 1.D0        horiVertRatio        = 1.D0
47          Ro_SeaLevel = 0.
48        DO k=1,Nr        DO k=1,Nr
49         delRdefault(k) = 0.         delRdefault(k) = 0.
50        ENDDO        ENDDO
# Line 62  C     In spherical polar distances are i Line 63  C     In spherical polar distances are i
63        phiMin               = 0.0        phiMin               = 0.0
64        thetaMin             = 0.        thetaMin             = 0.
65        rSphere              = 6370. * 1.D3        rSphere              = 6370. * 1.D3
66    C     General curvilinear coordinate system
67          usingCurvilinearGrid = .FALSE.
68    
69  C--   Set default "physical" parameters  C--   Set default "physical" parameters
70        DO K =1,Nr        DO K =1,Nr
# Line 71  C--   Set default "physical" parameters Line 74  C--   Set default "physical" parameters
74        rhoNil              = 999.8D0        rhoNil              = 999.8D0
75        f0                  = 1.D-4        f0                  = 1.D-4
76        beta                = 1.D-11        beta                = 1.D-11
77          omega               = 2.D0 * PI / ( 3600.D0 * 24.D0 )
78        viscAh              = 0.D3        viscAh              = 0.D3
79        diffKhT             = 0.D3        diffKhT             = 0.D3
80        diffKhS             = 0.D3        diffKhS             = 0.D3
# Line 90  C--   Set default "physical" parameters Line 94  C--   Set default "physical" parameters
94        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
95        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
96        rigidLid            = .FALSE.        rigidLid            = .FALSE.
97          implicSurfPress     = 1.D0
98          implicDiv2DFlow     = 1.D0
99        hFacMin             = 1.D0        hFacMin             = 1.D0
100        hFacMinDrDefault    = 0.D0        hFacMinDrDefault    = 0.D0
101          exactConserv        = .FALSE.
102          nonlinFreeSurf      = .FALSE.
103          hFacInf             = 1.D0
104          hFacSup             = 1.D0
105          staggerTimeStep     = .FALSE.
106        momViscosity        = .TRUE.        momViscosity        = .TRUE.
107        momAdvection        = .TRUE.        momAdvection        = .TRUE.
108        momForcing          = .TRUE.        momForcing          = .TRUE.
# Line 100  C--   Set default "physical" parameters Line 111  C--   Set default "physical" parameters
111        momStepping         = .TRUE.        momStepping         = .TRUE.
112        tempStepping        = .TRUE.        tempStepping        = .TRUE.
113        saltStepping        = .TRUE.        saltStepping        = .TRUE.
114        metricTerms         = .FALSE.        metricTerms         = .TRUE.
115        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
116        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
       openBoundaries      = .FALSE.  
117        nonHydrostatic      = .FALSE.        nonHydrostatic      = .FALSE.
118        globalFiles         = .FALSE.        globalFiles         = .FALSE.
119        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
120        ivdc_kappa          = 0.D0        ivdc_kappa          = 0.D0
121          groundAtK1          = .FALSE.
122          zonal_filt_lat      = 90.
123          zonal_filt_sinpow   = 2
124          zonal_filt_cospow   = 2
125          bottomDragLinear    = 0.
126          bottomDragQuadratic = 0.
127          usePickupBeforeC35  = .FALSE.
128          debugMode           = .FALSE.
129    
130  C--   Elliptic solver parameters  C--   Elliptic solver parameters
131        cg2dMaxIters       = 150        cg2dMaxIters       = 150
132        cg2dTargetResidual = 1.D-7        cg2dTargetResidual = 1.D-7
133          cg2dTargetResWunit = -1.
134        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
135        cg2dpcOffDFac      = 0.51D0        cg2dpcOffDFac      = 0.51D0
136        cg3dMaxIters       = 150        cg3dMaxIters       = 150
# Line 129  C--   Time stepping parameters Line 148  C--   Time stepping parameters
148        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
149        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
150        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
151          monitorFreq       = -1.
152        taveFreq          = deltaT*0        taveFreq          = deltaT*0
153        writeStatePrec    = precFloat64        writeStatePrec    = precFloat64
154        writeBinaryPrec   = precFloat32        writeBinaryPrec   = precFloat32
# Line 154  C--   Input files Line 174  C--   Input files
174        saltClimFile    = ' '        saltClimFile    = ' '
175        EmPmRfile       = ' '        EmPmRfile       = ' '
176        surfQfile       = ' '        surfQfile       = ' '
177          surfQswfile     = ' '
178        uVelInitFile    = ' '        uVelInitFile    = ' '
179        vVelInitFile    = ' '        vVelInitFile    = ' '
180        pSurfInitFile   = ' '        pSurfInitFile   = ' '
181          dQdTFile        = ' '
 C--   OBCS  
       DO I=1,Nx  
        OB_Jnorth(I)=0  
        OB_Jsouth(I)=0  
       ENDDO  
       DO J=1,Ny  
        OB_Ieast(J)=0  
        OB_Iwest(J)=0  
       ENDDO  
182    
183  C  C
184        RETURN        RETURN
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22