/[MITgcm]/MITgcm_contrib/osse/codemod/set_defaults.F
ViewVC logotype

Diff of /MITgcm_contrib/osse/codemod/set_defaults.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by afe, Tue Jun 22 19:44:41 2004 UTC revision 1.2 by afe, Wed Jun 23 19:56:48 2004 UTC
# Line 65  C     Vertical gridding Line 65  C     Vertical gridding
65    
66  C     Horizontal gridding  C     Horizontal gridding
67  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
68        usingCartesianGrid = .TRUE.        usingCartesianGrid = .FALSE.
69        DO i=1,Nx        DO i=1,Nx
70         delX(i) = UNSET_RL         delX(i) = UNSET_RL
71        ENDDO        ENDDO
# Line 81  C     General curvilinear coordinate sys Line 81  C     General curvilinear coordinate sys
81        usingCurvilinearGrid = .FALSE.        usingCurvilinearGrid = .FALSE.
82  C     General cylindrical coordinate system  C     General cylindrical coordinate system
83        bUseCylindricalGrid = .FALSE.        bUseCylindricalGrid = .FALSE.
84    
85    
86  C--   Set default "physical" parameters  C--   Set default "physical" parameters
87        DO K =1,Nr        DO K =1,Nr
88         tRef(K) = 30.D0 - FLOAT( K )         tRef(K) = 30.D0 - FLOAT( K )
89  Cml       sRef(K) = 35.D0  Cml       sRef(K) = 35.D0
90        ENDDO        ENDDO
       gravitySign         = 1.D0  
91        gravity             = 9.81D0        gravity             = 9.81D0
92        rhoNil              = 999.8D0        rhoNil              = 999.8D0
93        rhoConstFresh       = 999.8D0  C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)
94    C         (so that the default produces same results as before)
95    c     rhoConstFresh       = 999.8D0
96        f0                  = 1.D-4        f0                  = 1.D-4
97        beta                = 1.D-11        beta                = 1.D-11
98        omega               = 2.D0 * PI / ( 3600.D0 * 24.D0 )  C-    Always use 1 day in the past but should be 86164 (=86400*365.25/366.25)
99          rotationPeriod      = 86400. _d 0
100        viscAh              = 0.D3        viscAh              = 0.D3
101          viscAhGrid          = 0.D0
102          viscAhMax           = 1.D21
103          viscC2leith         = 0.D0
104        viscAstrain         = 0.D3        viscAstrain         = 0.D3
105        viscAtension        = 0.D3        viscAtension        = 0.D3
106        diffKhT             = 0.D3        diffKhT             = 0.D3
# Line 103  Cml       sRef(K) = 35.D0 Line 110  Cml       sRef(K) = 35.D0
110        no_slip_bottom      = .TRUE.        no_slip_bottom      = .TRUE.
111        diffKrTDefault      = 0.D-3        diffKrTDefault      = 0.D-3
112        diffKrSDefault      = 0.D-3        diffKrSDefault      = 0.D-3
113          diffKrBL79surf      = 0.D0
114          diffKrBL79deep      = 0.D0
115          diffKrBL79scl       = 200.D0
116          diffKrBL79Ho        = -2000.D0
117        viscA4              = 0.D11        viscA4              = 0.D11
118          viscA4Grid          = 0.D0
119          viscA4Max           = 1.D21
120          viscC4leith         = 0.D0
121        diffK4T             = 0.D11        diffK4T             = 0.D11
122        diffK4S             = 0.D11        diffK4S             = 0.D11
123        cosPower            = 0.        cosPower            = 0.
# Line 130  Cml      sBeta               = 7.4D-4 Line 144  Cml      sBeta               = 7.4D-4
144        saltForcing         = .TRUE.        saltForcing         = .TRUE.
145        tr1Stepping         = .FALSE.        tr1Stepping         = .FALSE.
146        metricTerms         = .TRUE.        metricTerms         = .TRUE.
147        useNHMTerms         = .TRUE.        useNHMTerms         = .FALSE.
148        implicitDiffusion   = .FALSE.        implicitDiffusion   = .FALSE.
149        implicitViscosity   = .FALSE.        implicitViscosity   = .FALSE.
150          momImplVertAdv      = .FALSE.
151          tempImplVertAdv     = .FALSE.
152          saltImplVertAdv     = .FALSE.
153        nonHydrostatic      = .FALSE.        nonHydrostatic      = .FALSE.
154        quasiHydrostatic    = .FALSE.        quasiHydrostatic    = .FALSE.
155        globalFiles         = .FALSE.        globalFiles         = .FALSE.
156          useSingleCpuIO      = .FALSE.
157        allowFreezing       = .FALSE.        allowFreezing       = .FALSE.
158          useOldFreezing      = .FALSE.
159        ivdc_kappa          = 0.D0        ivdc_kappa          = 0.D0
160        groundAtK1          = .FALSE.        groundAtK1          = .FALSE.
161        bottomDragLinear    = 0.        bottomDragLinear    = 0.
# Line 149  Cml      sBeta               = 7.4D-4 Line 168  Cml      sBeta               = 7.4D-4
168        saltAdvScheme       = 2        saltAdvScheme       = 2
169        tracerAdvScheme     = 2        tracerAdvScheme     = 2
170        multiDimAdvection   = .TRUE.        multiDimAdvection   = .TRUE.
171          useCDscheme         = .FALSE.
172        useEnergyConservingCoriolis = .FALSE.        useEnergyConservingCoriolis = .FALSE.
173        useJamartWetPoints  = .FALSE.        useJamartWetPoints  = .FALSE.
174          useJamartMomAdv     = .FALSE.
175          SadournyCoriolis    = .FALSE.
176          upwindVorticity     = .FALSE.
177          highOrderVorticity  = .FALSE.
178          useAbsVorticity     = .FALSE.
179          debugLevel          = debLevA
180    
181  C--   Set (free)surface-related parameters  C--   Set (free)surface-related parameters
182        implicitFreeSurface = .TRUE.        implicitFreeSurface = .TRUE.
# Line 160  C--   Set (free)surface-related paramete Line 186  C--   Set (free)surface-related paramete
186        exactConserv        = .FALSE.        exactConserv        = .FALSE.
187        uniformLin_PhiSurf  = .TRUE.        uniformLin_PhiSurf  = .TRUE.
188        nonlinFreeSurf      = 0        nonlinFreeSurf      = 0
189        hFacInf             = 1.D0        hFacInf             = 0.2 _d 0
190        hFacSup             = 1.D0        hFacSup             = 2.0 _d 0
191          select_rStar        = 0
192        useRealFreshWaterFlux = .FALSE.        useRealFreshWaterFlux = .FALSE.
193        temp_EvPrRn = UNSET_RL        temp_EvPrRn = UNSET_RL
194        salt_EvPrRn = 0.        salt_EvPrRn = 0.
195        trac_EvPrRn = UNSET_RL        trac_EvPrRn = UNSET_RL
196    
197  C--   Atmospheric physical parameters (e.g.: EOS)  C--   Atmospheric physical parameters (e.g.: EOS)
198        atm_Po =  1.D5        celsius2K = 273.16 _d 0
199        atm_Cp =  1004.D0        atm_Po =  1. _d 5
200        atm_Rd =  UNSET_RL        atm_Cp = 1004. _d 0
201        atm_kappa = 2.D0 / 7.D0        atm_Rd = UNSET_RL
202          atm_kappa = 2. _d 0 / 7. _d 0
203          atm_Rq = 0. _d 0
204        integr_GeoPot = 2        integr_GeoPot = 2
205        selectFindRoSurf = 0        selectFindRoSurf = 0
206    
# Line 189  C--   Time stepping parameters Line 218  C--   Time stepping parameters
218        deltaT            = 0.        deltaT            = 0.
219        nIter0            = 0        nIter0            = 0
220        startTime         = deltaT*float(nIter0)        startTime         = deltaT*float(nIter0)
221          pickupSuff        = ' '
222        nTimeSteps        = 0        nTimeSteps        = 0
223        nEndIter          = nIter0+nTimeSteps        nEndIter          = nIter0+nTimeSteps
224        endTime           = deltaT*float(nEndIter)        endTime           = deltaT*float(nEndIter)
# Line 197  C--   Time stepping parameters Line 227  C--   Time stepping parameters
227        pchkPtFreq        = deltaT*0        pchkPtFreq        = deltaT*0
228        chkPtFreq         = deltaT*0        chkPtFreq         = deltaT*0
229        dumpFreq          = deltaT*0        dumpFreq          = deltaT*0
230          adjDumpFreq       = deltaT*0
231        diagFreq          = deltaT*0        diagFreq          = deltaT*0
232        monitorFreq       = -1.        monitorFreq       = -1.
233        taveFreq          = deltaT*0        taveFreq          = deltaT*0
# Line 215  C--   Time stepping parameters Line 246  C--   Time stepping parameters
246        periodicExternalForcing = .FALSE.        periodicExternalForcing = .FALSE.
247        externForcingPeriod     = 0.        externForcingPeriod     = 0.
248        externForcingCycle      = 0.        externForcingCycle      = 0.
249        tCyl              = 0.        tCyl              = 0.
250    
251    
252  C--   Input files  C--   Input files
253        bathyFile       = ' '        bathyFile       = ' '
# Line 234  C--   Input files Line 266  C--   Input files
266        pSurfInitFile   = ' '        pSurfInitFile   = ' '
267        dQdTFile        = ' '        dQdTFile        = ' '
268        ploadFile       = ' '        ploadFile       = ' '
269          mdsioLocalDir   = ' '
270    
271  C  C
272        RETURN        RETURN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22