--- MITgcm/verification/global_ocean.cs32x15/input.thsice/data	2004/09/10 12:19:34	1.5
+++ MITgcm/verification/global_ocean.cs32x15/input.thsice/data	2017/08/21 18:22:51	1.12
@@ -8,16 +8,14 @@
  sRef=15*35.,
  viscAh =3.E5,
  viscAr =1.E-3,
- diffKhT=0., 
- diffK4T=0., 
  diffKrT=3.E-5,
- diffKhS=0.,  
- diffK4S=0., 
  diffKrS=3.E-5,
  ivdc_kappa=10.,
  implicitDiffusion=.TRUE.,
+ implicitViscosity=.TRUE.,
+ selectImplicitDrag=2,
  gravity=9.81,
- rhonil=1035.,
+ rhoConst=1035.,
  rhoConstFresh=1000.,
  eosType='JMD95Z',
  staggerTimeStep=.TRUE.,
@@ -33,6 +31,7 @@
  hFacMinDr=20.,
  readBinaryPrec=64,
 #writeBinaryPrec=64,
+ plotLevel=0,
  &
 
 # Elliptic solver parameters
@@ -46,12 +45,13 @@
  &PARM03
  nIter0=36000,
  nTimeSteps=20,
- deltaTmom   =1200.,
+ deltaTMom   =1200.,
  deltaTtracer=86400.,
- deltaTfreesurf=86400.,
+ deltaTFreeSurf=86400.,
  deltaTClock =86400.,
  abEps = 0.1,
- forcing_In_AB=.FALSE.,
+#forcing_In_AB=.FALSE.,
+ tracForcingOutAB=1,
  pChkptFreq  =31104000.,
  taveFreq    =31104000.,
  dumpFreq    =31104000.,
@@ -66,14 +66,13 @@
  tauSaltClimRelax = 15552000.,
  latBandClimRelax=50.,
  monitorFreq =1.,
- monitor_ioinc=.TRUE.,
  &
 
 # Gridding parameters
  &PARM04
- usingCartesianGrid=.FALSE.,
- usingSphericalPolarGrid=.FALSE.,
  usingCurvilinearGrid=.TRUE.,
+ horizGridFile='grid_cs32',
+ radius_fromHorizGrid=6370.E3,
  delR= 50., 70., 100., 140., 190.,
        240., 290., 340., 390., 440.,
        490., 540., 590., 640., 690.,