/[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.109 by heimbach, Tue Dec 13 19:45:19 2005 UTC revision 1.129 by jmc, Thu Aug 23 19:08:25 2007 UTC
# Line 21  C     !USES: Line 21  C     !USES:
21  #include "EEPARAMS.h"  #include "EEPARAMS.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
23  Cml#include "EOS.h"  Cml#include "EOS.h"
24  #include "GRID.h"  c#include "GRID.h"
25    
26  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
27  C     myThid - Number of this instance of INI_PARMS  C     myThid - Number of this instance of INI_PARMS
# Line 39  C     K, I, J - Loop counters Line 39  C     K, I, J - Loop counters
39  CEOP  CEOP
40    
41  C--   Grid parameters  C--   Grid parameters
42  C     Vertical gridding  C-    Vertical gridding
43        rkFacDefault         = 1.D0        delRFile            = ' '
44        horiVertRatio        = 1.D0        delRcFile           = ' '
45        Ro_SeaLevel = 0.        rkFacDefault        = 1. _d 0
46          Ro_SeaLevel         = 0.
47        DO k=1,Nr        DO k=1,Nr
48         delRdefault(k) = 0.         delRdefault(k)     = 0.
49        ENDDO        ENDDO
50        DO k=1,Nr+1        DO k=1,Nr+1
51         delRc(k) = UNSET_RL         delRc(k)           = UNSET_RL
52        ENDDO        ENDDO
53    C-    vertical profile
54  C     Horizontal gridding        tRefFile            = ' '
55        delXFile        = ' '        sRefFile            = ' '
56        delYFile        = ' '        rhoRefFile          = ' '
57        horizGridFile   = ' '  
58    C-    Horizontal gridding
59          delXFile            = ' '
60          delYFile            = ' '
61          horizGridFile       = ' '
62          deepAtmosphere      = .FALSE.
63  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
64        usingCartesianGrid = .FALSE.        usingCartesianGrid  = .FALSE.
65        DO i=1,Nx        DO i=1,Nx
66         delX(i) = UNSET_RL         delX(i)            = UNSET_RL
67        ENDDO        ENDDO
68        DO j=1,Ny        DO j=1,Ny
69         delY(j) = UNSET_RL         delY(j)            = UNSET_RL
70        ENDDO        ENDDO
71  C     In spherical polar distances are in degrees  C     In spherical polar distances are in degrees
72        usingSphericalPolarGrid = .FALSE.        usingSphericalPolarGrid = .FALSE.
73        phiMin               = 0.0        phiMin              = 0.0
74        thetaMin             = 0.        thetaMin            = 0.
75        rSphere              = 6370. * 1.D3        rSphere             = 6370. _d 3
76  C     General curvilinear coordinate system  C     General curvilinear coordinate system
77        usingCurvilinearGrid = .FALSE.        usingCurvilinearGrid= .FALSE.
78  C     General cylindrical coordinate system  C     General cylindrical coordinate system
79        usingCylindricalGrid = .FALSE.        usingCylindricalGrid= .FALSE.
80  C     Coriolis map:  C     Coriolis map:
81        useConstantF  = .FALSE.        useConstantF        = .FALSE.
82        useBetaPlaneF = .FALSE.        useBetaPlaneF       = .FALSE.
83        useSphereF    = .FALSE.        useSphereF          = .TRUE.
84          use3dCoriolis       = .TRUE.
85    
86  C--   Set default "physical" parameters  C--   Set default "physical" parameters
87        nh_Am2              = 1.D0        nh_Am2              = 1. _d 0
88        gravity             = 9.81D0        gravity             = 9.81 _d 0
89        rhoNil              = 999.8D0        rhoNil              = 999.8 _d 0
90  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)
91  C         (so that the default produces same results as before)  C         (so that the default produces same results as before)
92  c     rhoConstFresh       = 999.8D0  c     rhoConstFresh       = 999.8 _d 0
93        f0                  = 1.D-4        f0                  = 1. _d -4
94        beta                = 1.D-11        beta                = 1. _d -11
95  C-    Always use 1 day in the past but should be 86164 (=86400*365.25/366.25)  C-    Earth rotation period is 86400*365.25/366.25 (use to be 1.day)
96        rotationPeriod      = 86400. _d 0        rotationPeriod      = 86164. _d 0
97        viscAh              = 0.D3        viscAh              = 0. _d 3
98        viscAhGrid          = 0.D0        viscAhGrid          = 0. _d 0
99        viscAhGridMin       = 0.D0        viscAhGridMin       = 0. _d 0
100        viscAhGridMax       = 1.D21        viscAhGridMax       = 1. _d 21
101        viscAhMax           = 1.D21        viscAhMax           = 1. _d 21
102        viscAhReMax         = 0.D0        viscAhReMax         = 0. _d 0
103        viscC2leith         = 0.D0        viscC2leith         = 0. _d 0
104        viscC2leithD        = 0.D0        viscC2leithD        = 0. _d 0
105        viscC2smag          = 0.D0        viscC2smag          = 0. _d 0
106        diffKhT             = 0.D3        diffKhT             = 0. _d 3
107        diffKhS             = 0.D3        diffKhS             = 0. _d 3
108        viscArDefault       = 0.D-3        viscArDefault       = 0. _d -3
109        no_slip_sides       = .TRUE.        no_slip_sides       = .TRUE.
110        no_slip_bottom      = .TRUE.        no_slip_bottom      = .TRUE.
111        sideDragFactor      = 2. _d 0        sideDragFactor      = 2. _d 0
112        bottomDragLinear    = 0.        bottomDragLinear    = 0.
113        bottomDragQuadratic = 0.        bottomDragQuadratic = 0.
114        diffKrTDefault      = 0.D-3        smoothAbsFuncRange  = 0. _d 0
115        diffKrSDefault      = 0.D-3        diffKrTDefault      = 0. _d -3
116        diffKrBL79surf      = 0.D0        diffKrSDefault      = 0. _d -3
117        diffKrBL79deep      = 0.D0        diffKrBL79surf      = 0. _d 0
118        diffKrBL79scl       = 200.D0        diffKrBL79deep      = 0. _d 0
119        diffKrBL79Ho        = -2000.D0        diffKrBL79scl       = 200. _d 0
120        viscA4              = 0.D11        diffKrBL79Ho        = -2000. _d 0
121        viscA4Grid          = 0.D0        BL79LatVary         = 30
122        viscA4GridMax       = 1.D21        diffKrBLEQsurf      = UNSET_RL
123        viscA4GridMin       = 0.D0        diffKrBLEQdeep      = UNSET_RL
124        viscA4Max           = 1.D21        diffKrBLEQscl       = UNSET_RL
125        viscA4ReMax         = 0.D0        diffKrBLEQHo        = UNSET_RL
126        viscC4leith         = 0.D0        viscA4              = 0. _d 11
127        viscC4leithD        = 0.D0        viscA4Grid          = 0. _d 0
128        viscC4smag          = 0.D0        viscA4GridMax       = 1. _d 21
129        diffK4T             = 0.D11        viscA4GridMin       = 0. _d 0
130        diffK4S             = 0.D11        viscA4Max           = 1. _d 21
131          viscA4ReMax         = 0. _d 0
132          viscC4leith         = 0. _d 0
133          viscC4leithD        = 0. _d 0
134          viscC4smag          = 0. _d 0
135          diffK4T             = 0. _d 11
136          diffK4S             = 0. _d 11
137        cosPower            = 0.        cosPower            = 0.
138        HeatCapacity_Cp     = 3994.D0        HeatCapacity_Cp     = 3994. _d 0
139  Cml      tAlpha              = 2.D-4  Cml      tAlpha              = 2. _d -4
140  Cml      sBeta               = 7.4D-4  Cml      sBeta               = 7.4 _d -4
141        eosType             = 'LINEAR'        eosType             = 'LINEAR'
142        buoyancyRelation    = 'OCEANIC'        buoyancyRelation    = 'OCEANIC'
143        hFacMin             = 1.D0        hFacMin             = 1. _d 0
144        hFacMinDrDefault    = 0.D0        hFacMinDrDefault    = 0. _d 0
145        implicitIntGravWave = .FALSE.        implicitIntGravWave = .FALSE.
146        staggerTimeStep     = .FALSE.        staggerTimeStep     = .FALSE.
147        momViscosity        = .TRUE.        momViscosity        = .TRUE.
# Line 149  Cml      sBeta               = 7.4D-4 Line 161  Cml      sBeta               = 7.4D-4
161        useNHMTerms         = .FALSE.        useNHMTerms         = .FALSE.
162        useFullLeith        = .FALSE.        useFullLeith        = .FALSE.
163        useAreaViscLength   = .FALSE.        useAreaViscLength   = .FALSE.
164        useStrainTensionVisc = .FALSE.        useStrainTensionVisc= .FALSE.
165        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
166        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
167        momImplVertAdv      = .FALSE.        momImplVertAdv      = .FALSE.
# Line 161  Cml      sBeta               = 7.4D-4 Line 173  Cml      sBeta               = 7.4D-4
173        useSingleCpuIO      = .FALSE.        useSingleCpuIO      = .FALSE.
174        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
175        useOldFreezing      = .FALSE.        useOldFreezing      = .FALSE.
176        ivdc_kappa          = 0.D0        ivdc_kappa          = 0. _d 0
177        usePickupBeforeC35    = .FALSE.        hMixCriteria        = -.8 _d 0
178        usePickupBeforeC54    = .FALSE.        usePickupBeforeC54    = .FALSE.
179        debugMode             = .FALSE.        debugMode             = .FALSE.
180        tempAdvScheme       = 2        tempAdvScheme       = 2
181        saltAdvScheme       = 2        saltAdvScheme       = 2
182        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
183          useMultiDimAdvec    = .FALSE.
184        useCDscheme         = .FALSE.        useCDscheme         = .FALSE.
185        useEnergyConservingCoriolis = .FALSE.        useEnergyConservingCoriolis = .FALSE.
186        useJamartWetPoints  = .FALSE.        useJamartWetPoints  = .FALSE.
# Line 185  Cml      sBeta               = 7.4D-4 Line 198  Cml      sBeta               = 7.4D-4
198  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
199        implicitFreeSurface = .FALSE.        implicitFreeSurface = .FALSE.
200        rigidLid            = .FALSE.        rigidLid            = .FALSE.
201        implicSurfPress     = 1.D0        implicSurfPress     = 1. _d 0
202        implicDiv2DFlow     = 1.D0        implicDiv2DFlow     = 1. _d 0
203        exactConserv        = .FALSE.        exactConserv        = .FALSE.
204          linFSConserveTr     = .FALSE.
205        uniformLin_PhiSurf  = .TRUE.        uniformLin_PhiSurf  = .TRUE.
206        nonlinFreeSurf      = 0        nonlinFreeSurf      = 0
207        hFacInf             = 0.2 _d 0        hFacInf             = 0.2 _d 0
# Line 196  C--   Set (free)surface-related paramete Line 210  C--   Set (free)surface-related paramete
210        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
211        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
212        salt_EvPrRn = 0.        salt_EvPrRn = 0.
213          balanceEmPmR        = .FALSE.
214          balanceQnet         = .FALSE.
215          balancePrintMean    = .FALSE.
216    
217  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
218        celsius2K = 273.16 _d 0        celsius2K = 273.16 _d 0
# Line 209  C--   Atmospheric physical parameters (e Line 226  C--   Atmospheric physical parameters (e
226    
227  C--   Elliptic solver parameters  C--   Elliptic solver parameters
228        cg2dMaxIters       = 150        cg2dMaxIters       = 150
229        cg2dTargetResidual = 1.D-7        cg2dTargetResidual = 1. _d -7
230        cg2dTargetResWunit = -1.        cg2dTargetResWunit = -1.
231        cg2dChkResFreq     = 1        cg2dChkResFreq     = 1
232        cg2dpcOffDFac      = 0.51D0        cg2dpcOffDFac      = 0.51 _d 0
233        cg2dPreCondFreq    = 1        cg2dPreCondFreq    = 1
234        cg3dMaxIters       = 150        cg3dMaxIters       = 150
235        cg3dTargetResidual = 1.D-7        cg3dTargetResidual = 1. _d -7
236        cg3dChkResFreq     = 1        cg3dChkResFreq     = 1
237    
238  C--   Time stepping parameters  C--   Time stepping parameters
# Line 232  C--   Time stepping parameters Line 249  C--   Time stepping parameters
249        nTimeSteps        = 0        nTimeSteps        = 0
250        nEndIter          = nIter0+nTimeSteps        nEndIter          = nIter0+nTimeSteps
251        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
252        forcing_In_AB     = .TRUE.        momForcingOutAB   = UNSET_I
253          tracForcingOutAB  = UNSET_I
254          momDissip_In_AB   = .TRUE.
255        doAB_onGtGs       = .TRUE.        doAB_onGtGs       = .TRUE.
256        abEps             = 0.01        abEps             = 0.01 _d 0
257  #ifdef ALLOW_ADAMSBASHFORTH_3  #ifdef ALLOW_ADAMSBASHFORTH_3
258        alph_AB           = 0.5 _d 0        alph_AB           = 0.5 _d 0
259        beta_AB           = 5. _d 0 / 12. _d 0        beta_AB           = 5. _d 0 / 12. _d 0
# Line 244  C--   Time stepping parameters Line 263  C--   Time stepping parameters
263        beta_AB           = UNSET_RL        beta_AB           = UNSET_RL
264        startFromPickupAB2= .TRUE.        startFromPickupAB2= .TRUE.
265  #endif  #endif
266        pchkPtFreq        = deltaT*0        pChkPtFreq        = deltaT*0
267        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
268        outputTypesInclusive = .FALSE.        outputTypesInclusive = .FALSE.
269        pickup_read_mdsio = .TRUE.        pickup_read_mdsio = .TRUE.
# Line 267  C--   Time stepping parameters Line 286  C--   Time stepping parameters
286        nCheckLev         = 1        nCheckLev         = 1
287        checkPtSuff(1)    = 'ckptA'        checkPtSuff(1)    = 'ckptA'
288        checkPtSuff(2)    = 'ckptB'        checkPtSuff(2)    = 'ckptB'
289        cAdjFreq          =  0.D0        cAdjFreq          =  0. _d 0
290        tauCD             =  0.D0        tauCD             =  0. _d 0
291        tauThetaClimRelax =  0.D0        tauThetaClimRelax =  0. _d 0
292        tauSaltClimRelax  =  0.D0        tauSaltClimRelax  =  0. _d 0
293        tauThetaClimRelax3Dim =  0.D0        tauTr1ClimRelax   =  0. _d 0
       tauSaltClimRelax3Dim  =  0.D0  
       tauTr1ClimRelax   =  0.D0  
294        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
295        externForcingPeriod     = 0.        externForcingPeriod     = 0.
296        externForcingCycle      = 0.        externForcingCycle      = 0.
297        tCylIn             = 0.        tCylIn             = 0.
298        tCylOut            = 20.        tCylOut            = 20.
       calendarDumps     = .FALSE.  
299    
300  C--   Input files  C--   Input files
301        bathyFile       = ' '        bathyFile       = ' '
302        topoFile        = ' '        topoFile        = ' '
303          shelfIceFile    = ' '
304        hydrogSaltFile  = ' '        hydrogSaltFile  = ' '
305        hydrogThetaFile = ' '        hydrogThetaFile = ' '
306          diffKrFile      = ' '
307        zonalWindFile   = ' '        zonalWindFile   = ' '
308        meridWindFile   = ' '        meridWindFile   = ' '
309        thetaClimFile   = ' '        thetaClimFile   = ' '
# Line 305  C--   Input files Line 323  C--   Input files
323        lambdaThetaFile = ' '        lambdaThetaFile = ' '
324        lambdaSaltFile  = ' '        lambdaSaltFile  = ' '
325        mdsioLocalDir   = ' '        mdsioLocalDir   = ' '
326          the_run_name    = ' '
327    
 C  
328        RETURN        RETURN
329        END        END
Legend:
Removed from v.1.109  
changed lines
  Added in v.1.129
  ViewVC Help
Powered by ViewVC 1.1.22