| 7 | 
  tRef               = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13., | 
  tRef               = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13., | 
| 8 | 
                       12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1., | 
                       12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1., | 
| 9 | 
  sRef               = 50*34.5, | 
  sRef               = 50*34.5, | 
| 10 | 
  viscAr= 1.77e-5 | 
  viscAr= 6.0826e-04, | 
| 11 | 
  no_slip_sides=.TRUE., | 
  no_slip_sides=.FALSE., | 
| 12 | 
  no_slip_bottom=.TRUE., | 
  no_slip_bottom=.TRUE., | 
| 13 | 
  diffKhT=0., | 
  diffKhT=0., | 
| 14 | 
  diffK4T=0., | 
  diffK4T=0., | 
| 15 | 
  diffKrT=1.51e-5 | 
  diffKrT=0., | 
| 16 | 
  diffKhS=0., | 
  diffKhS=0., | 
| 17 | 
  diffK4S=0., | 
  diffK4S=0., | 
| 18 | 
  diffKrS=1.51e-5 | 
  diffKrS=0., | 
| 19 | 
  | 
 # Bryan and Lewis 1979 vertical diffusivity parameters | 
| 20 | 
  | 
 # with latitudinal dependence | 
| 21 | 
  | 
  diffKrBL79surf      = 3.0e-5, | 
| 22 | 
  | 
  diffKrBL79deep      = 1.3e-4, | 
| 23 | 
  | 
  diffKrBL79scl       = 200., | 
| 24 | 
  | 
  diffKrBL79Ho        = -2000., | 
| 25 | 
  | 
  BL79LatVary         = 30., | 
| 26 | 
  | 
  diffKrBLEQsurf      = 8.0e-6, | 
| 27 | 
  | 
  diffKrBLEQdeep      = 1.3e-4, | 
| 28 | 
  | 
  diffKrBLEQscl       = 200., | 
| 29 | 
  | 
  diffKrBLEQHo        = -2000., | 
| 30 | 
  | 
 # | 
| 31 | 
  tAlpha=2.E-4, | 
  tAlpha=2.E-4, | 
| 32 | 
  sBeta =7.4E-4, | 
  sBeta =7.4E-4, | 
| 33 | 
  gravity=9.81, | 
  gravity=9.81, | 
 | 
  rigidLid=.FALSE., | 
  | 
 | 
  implicitFreeSurface=.TRUE., | 
  | 
| 34 | 
  eosType='JMD95Z', | 
  eosType='JMD95Z', | 
| 35 | 
  readBinaryPrec=32, | 
  readBinaryPrec=32, | 
| 36 | 
  writeBinaryPrec=32, | 
  writeBinaryPrec=32, | 
| 37 | 
  useSingleCPUio=.TRUE. | 
  useSingleCPUio=.TRUE., | 
| 38 | 
  hFacMinDr=50., | 
  hFacMinDr=50., | 
| 39 | 
  hFacMin=0.1, | 
  hFacMin=0.3, | 
| 40 | 
  hFacInf=0.1, | 
  hFacInf=0.1, | 
| 41 | 
  hFacSup=5., | 
  hFacSup=5., | 
| 42 | 
 # z* - begin | 
 # z* - begin | 
| 46 | 
  implicitDiffusion=.true., | 
  implicitDiffusion=.true., | 
| 47 | 
  implicitViscosity=.true., | 
  implicitViscosity=.true., | 
| 48 | 
  allowFreezing=.false., | 
  allowFreezing=.false., | 
| 49 | 
  viscA4     = 0., | 
  viscAh=0., | 
| 50 | 
  viscA4Grid = 4e-3, | 
  viscAhGrid=0., | 
| 51 | 
  viscAh     = 0., | 
 # Following change took place starting in checkpoint57v_post | 
| 52 | 
  | 
 # Redefine viscC?Leith*.  Change the coefficients of Leith viscosity  | 
| 53 | 
  | 
 # to be proportional to enstrophy dissipation length.  Thus, | 
| 54 | 
  | 
 #  viscC?leith?_new=pi*(viscC?leith?_old)^(1/3) | 
| 55 | 
  | 
 # This is now consistent with Smagorinsky and with the manual. | 
| 56 | 
  | 
 # viscC4Leith=1.0, | 
| 57 | 
  | 
 # viscC4Leithd=1.0, | 
| 58 | 
  | 
 # viscC4Leith=3.1416, | 
| 59 | 
  | 
 # viscC4Leithd=3.1416, | 
| 60 | 
  | 
  viscC4Leith=1.5, | 
| 61 | 
  | 
  viscC4Leithd=1.5, | 
| 62 | 
  | 
  viscA4GridMax=0.5, | 
| 63 | 
  | 
 ########################################################## | 
| 64 | 
  | 
  viscC2Leith=0., | 
| 65 | 
  | 
  viscA4=0., | 
| 66 | 
  | 
  useAreaViscLength=.TRUE., | 
| 67 | 
  | 
  sideDragFactor=0., | 
| 68 | 
  | 
  SadournyCoriolis=.FALSE., | 
| 69 | 
  | 
  useAbsVorticity=.FALSE., | 
| 70 | 
  | 
  highOrderVorticity  = .TRUE., | 
| 71 | 
  bottomDragQuadratic = 0.002, | 
  bottomDragQuadratic = 0.002, | 
| 72 | 
  tempAdvScheme=33 | 
  tempAdvScheme=33, | 
| 73 | 
  saltAdvScheme=33 | 
  saltAdvScheme=33, | 
| 74 | 
  StaggerTimeStep=.true. | 
  StaggerTimeStep=.true., | 
| 75 | 
  multiDimAdvection=.true. | 
  multiDimAdvection=.true., | 
| 76 | 
  vectorInvariantMomentum=.TRUE., | 
  vectorInvariantMomentum=.TRUE., | 
| 77 | 
  rigidLid=.FALSE., | 
  rigidLid=.FALSE., | 
| 78 | 
  implicitFreeSurface=.TRUE., | 
  implicitFreeSurface=.TRUE., | 
| 79 | 
  exactConserv=.TRUE., | 
  exactConserv=.TRUE., | 
| 80 | 
  debuglevel=0 | 
  debuglevel=-1, | 
| 81 | 
  | 
  useRealFreshWaterFlux=.TRUE., | 
| 82 | 
  & | 
  & | 
| 83 | 
  | 
  | 
| 84 | 
 # Elliptic solver parameters | 
 # Elliptic solver parameters | 
| 89 | 
  | 
  | 
| 90 | 
 # Time stepping parameters | 
 # Time stepping parameters | 
| 91 | 
  &PARM03 | 
  &PARM03 | 
| 92 | 
  niter0=0, | 
  niter0=216, | 
| 93 | 
  endTime=259200 | 
 # endtime for December 31, 2002 | 
| 94 | 
  | 
  endtime=347155200., | 
| 95 | 
  abEps=0.1, | 
  abEps=0.1, | 
| 96 | 
  deltaT=300., | 
  deltaT=1200, | 
| 97 | 
  cAdjFreq=0., | 
  cAdjFreq=0., | 
| 98 | 
  pChkptFreq=259200 | 
  pChkptFreq=31536000., | 
| 99 | 
  tavefreq=259200 | 
  chkptFreq=2592000., | 
| 100 | 
  tauSaltClimRelax = 3842380., | 
 # blending: 12/24/06 | 
| 101 | 
  monitorFreq=1 | 
   tauSaltClimRelax = 8722203., | 
| 102 | 
  | 
 # tauSaltClimRelax = 3842380., | 
| 103 | 
  | 
  monitorFreq=86400., | 
| 104 | 
  | 
  forcing_In_AB = .FALSE., | 
| 105 | 
  | 
  dumpInitAndLast=.FALSE., | 
| 106 | 
  & | 
  & | 
| 107 | 
  | 
  | 
| 108 | 
 # Gridding parameters | 
 # Gridding parameters | 
| 120 | 
  | 
  | 
| 121 | 
 # Input datasets | 
 # Input datasets | 
| 122 | 
  &PARM05 | 
  &PARM05 | 
| 123 | 
  bathyFile      ='ETOPO2_18360x85_ver5e_filled.bin', | 
  bathyFile      ='ETOPO2_510x6x510_ver5e_filled.bin', | 
| 124 | 
  hydrogThetaFile='LEVT01_JAN_18360x85x50.bin', | 
  hydrogThetaFile='WGHC_T_510x6x510.bin', | 
| 125 | 
  hydrogSaltFile ='LEVS01_JAN_18360x85x50.bin', | 
  hydrogSaltFile ='WGHC_S_510x6x510.bin', | 
| 126 | 
  & | 
  & |