--- MITgcm_contrib/osse/codemod/set_defaults.F	2004/06/22 19:44:41	1.1
+++ MITgcm_contrib/osse/codemod/set_defaults.F	2004/06/24 17:52:39	1.3
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/osse/codemod/set_defaults.F,v 1.1 2004/06/22 19:44:41 afe Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/osse/codemod/set_defaults.F,v 1.3 2004/06/24 17:52:39 afe Exp $
 C $Name:  $
 
 #include "CPP_OPTIONS.h"
@@ -65,7 +65,7 @@
 
 C     Horizontal gridding
 C     In cartesian coords distances are in metres
-      usingCartesianGrid = .TRUE.
+      usingCartesianGrid = .FALSE.
       DO i=1,Nx
        delX(i) = UNSET_RL
       ENDDO
@@ -80,20 +80,27 @@
 C     General curvilinear coordinate system
       usingCurvilinearGrid = .FALSE.
 C     General cylindrical coordinate system
-      bUseCylindricalGrid = .FALSE.
+      usingCylindricalGrid = .FALSE.
+
+
 C--   Set default "physical" parameters
       DO K =1,Nr
        tRef(K) = 30.D0 - FLOAT( K )
 Cml       sRef(K) = 35.D0
       ENDDO
-      gravitySign         = 1.D0
       gravity             = 9.81D0
       rhoNil              = 999.8D0
-      rhoConstFresh       = 999.8D0
+C-- jmc : the default is to set rhoConstFresh to rhoConst (=rhoNil by default)
+C         (so that the default produces same results as before)
+c     rhoConstFresh       = 999.8D0
       f0                  = 1.D-4
       beta                = 1.D-11
-      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)
+      rotationPeriod      = 86400. _d 0
       viscAh              = 0.D3
+      viscAhGrid          = 0.D0
+      viscAhMax           = 1.D21
+      viscC2leith         = 0.D0
       viscAstrain         = 0.D3
       viscAtension        = 0.D3
       diffKhT             = 0.D3
@@ -103,7 +110,14 @@
       no_slip_bottom      = .TRUE.
       diffKrTDefault      = 0.D-3
       diffKrSDefault      = 0.D-3
+      diffKrBL79surf      = 0.D0
+      diffKrBL79deep      = 0.D0
+      diffKrBL79scl       = 200.D0
+      diffKrBL79Ho        = -2000.D0
       viscA4              = 0.D11
+      viscA4Grid          = 0.D0
+      viscA4Max           = 1.D21
+      viscC4leith         = 0.D0
       diffK4T             = 0.D11
       diffK4S             = 0.D11
       cosPower            = 0.
@@ -130,13 +144,18 @@
       saltForcing         = .TRUE.
       tr1Stepping         = .FALSE.
       metricTerms         = .TRUE.
-      useNHMTerms         = .TRUE.
+      useNHMTerms         = .FALSE.
       implicitDiffusion   = .FALSE.
       implicitViscosity   = .FALSE.
+      momImplVertAdv      = .FALSE.
+      tempImplVertAdv     = .FALSE.
+      saltImplVertAdv     = .FALSE.
       nonHydrostatic      = .FALSE.
       quasiHydrostatic    = .FALSE.
       globalFiles         = .FALSE.
+      useSingleCpuIO      = .FALSE.
       allowFreezing       = .FALSE.
+      useOldFreezing      = .FALSE.
       ivdc_kappa          = 0.D0
       groundAtK1          = .FALSE.
       bottomDragLinear    = 0.
@@ -149,8 +168,15 @@
       saltAdvScheme       = 2
       tracerAdvScheme     = 2
       multiDimAdvection   = .TRUE.
+      useCDscheme         = .FALSE.
       useEnergyConservingCoriolis = .FALSE.
       useJamartWetPoints  = .FALSE.
+      useJamartMomAdv     = .FALSE.
+      SadournyCoriolis    = .FALSE.
+      upwindVorticity     = .FALSE.
+      highOrderVorticity  = .FALSE.
+      useAbsVorticity     = .FALSE.
+      debugLevel          = debLevA
 
 C--   Set (free)surface-related parameters
       implicitFreeSurface = .TRUE.
@@ -160,18 +186,21 @@
       exactConserv        = .FALSE.
       uniformLin_PhiSurf  = .TRUE.
       nonlinFreeSurf      = 0
-      hFacInf             = 1.D0
-      hFacSup             = 1.D0
+      hFacInf             = 0.2 _d 0
+      hFacSup             = 2.0 _d 0
+      select_rStar        = 0
       useRealFreshWaterFlux = .FALSE.
       temp_EvPrRn = UNSET_RL
       salt_EvPrRn = 0.
       trac_EvPrRn = UNSET_RL
 
 C--   Atmospheric physical parameters (e.g.: EOS)
-      atm_Po =  1.D5
-      atm_Cp =  1004.D0
-      atm_Rd =  UNSET_RL
-      atm_kappa = 2.D0 / 7.D0
+      celsius2K = 273.16 _d 0
+      atm_Po =  1. _d 5
+      atm_Cp = 1004. _d 0
+      atm_Rd = UNSET_RL
+      atm_kappa = 2. _d 0 / 7. _d 0
+      atm_Rq = 0. _d 0
       integr_GeoPot = 2
       selectFindRoSurf = 0
 
@@ -189,6 +218,7 @@
       deltaT            = 0.
       nIter0            = 0
       startTime         = deltaT*float(nIter0)
+      pickupSuff        = ' '
       nTimeSteps        = 0
       nEndIter          = nIter0+nTimeSteps
       endTime           = deltaT*float(nEndIter)
@@ -197,6 +227,7 @@
       pchkPtFreq        = deltaT*0
       chkPtFreq         = deltaT*0
       dumpFreq          = deltaT*0
+      adjDumpFreq       = deltaT*0
       diagFreq          = deltaT*0
       monitorFreq       = -1.
       taveFreq          = deltaT*0
@@ -215,7 +246,8 @@
       periodicExternalForcing = .FALSE.
       externForcingPeriod     = 0.
       externForcingCycle      = 0.
-      tCyl		= 0.
+      tCyl              = 0.
+
 
 C--   Input files
       bathyFile       = ' '
@@ -234,6 +266,7 @@
       pSurfInitFile   = ' '
       dQdTFile        = ' '
       ploadFile       = ' '
+      mdsioLocalDir   = ' '
 
 C
       RETURN