[MITgcm-support] pickup problem

[Tom Haine]

Dear Jean-Michel, 

Thanks for your help. My setup has an implicit free surface (I'm unsure
if this is the NLFS or not), and simple obcs (Orlanski is switched off,
sponge is switched on in data.obcs). In addition, I'm using pkg/seaice,
pkg/kpp, pkg/ecco (ecco is just to handle exf/cal/seaice properly). My
code mods are very small and I'm running version: checkpoint57e_post.
I'm using MPI on the GFDL lsc (optfile = irix64_ip35+mpi_gfdl).

My data file is attached below. I hope this helps you.

Thanks again. Tom.

-- 
---------------------------------------------
 Dr. Thomas W. N. Haine,
 Associate Professor of Physical Oceanography

 Department of Earth & Planetary Sciences,
 329 Olin Hall, The Johns Hopkins University,
 Baltimore, MD 21218, USA.
 Tel : 410 516 7048,   Fax : 410 516 7933
---------------------------------------------
-------------- next part --------------
# ====================
# | Model parameters |
# ====================
#
# DJL (23/4/02)
# Modified by twnh 27Jul03.
#
# Continuous equation parameters
 &PARM01
 tref = 97*4.0,
 sref = 97*35.0,
 viscAz=1.E-3,
 viscA4Grid=0.,
 viscC4Leith=1.0,
 viscA4GridMax=0.005,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
 diffK4T=0.
 diffKzT=0.
 diffK4S=0.
 diffKzS=0.
 bottomDragQuadratic = 0.001,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 eosType='JMD95Z',
 readBinaryPrec=64,
 writeBinaryPrec=64,
 globalFiles=.TRUE.,
 useSingleCpuIo=.TRUE.,
 hFacMin=0.1,            
 hFacMinDz=5.,
 allowFreezing=.FALSE.,
 HeatCapacity_Cp = 3986.D0
 gravity         = 9.8156,
 gBaro           = 9.8156,
 rhoNil          = 1027.D0,
 implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
 staggerTimeStep=.TRUE.,
 tempAdvScheme=33,
 saltAdvScheme=33,
#this was moved somewhere else in last quarter of 2004 tracerAdvScheme=33,
 multiDimAdvection=.TRUE.,
# z* - begin
  select_rStar=2,
  nonlinFreeSurf=4,
  hFacInf=0.01,
  hFacSup=22.5,
  exactConserv=.TRUE.,
# z* - end
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-8,
 &

# Time stepping parameters
# 8hrs with 81 lsc cpus can manage 2880 steps (2days at 60s steps).
# 4.75hrs with 36 lsc cpus can manage 2880 steps (10days at 300s steps).
 &PARM03
 nIter0=4320,
 endTime=2160000.,
 deltaT=300.,
 tauCD=0.,
 cAdjFreq=0.,
 abEps=0.1,
 forcing_In_AB=.FALSE.,
 pChkptFreq=1296000.,
 dumpFreq=7200.,
 diagFreq=7200.,
 monitorFreq=7200.0,
 tauSaltClimRelax=0.,
 tauThetaClimRelax=0.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.TRUE.,
 delZ= 2,   3,   4,   5,   6,   7,   8,   9,  10,  11, 12,
      13,  14,  15,  16,  17,  18,  19,  20,  20,  20,
      20,  20,  20,  20,  20,  20,  20,  20,  20,  20,
      20,  20,  20,  20,  20,  20,  20,  20,  20,  20,
      20,  20,  20,  20,  20,  20,  20,  20,  20,  20,
      20,  20,  20,  20,  20,  20,  20,  20,  20,  20,
      20,  20,  20,  20,  20,  20,  20,  20,  20,  20,
      20,  20,  20,  20,  20,  20,  20,  20,  30,  40,
      50,  60,  70,  80,  90, 100, 110, 120, 130,
     140, 150, 160, 170, 180, 190, 200,
 phiMin=63.5, 
 delY=288*0.01735000,
 delX=432*0.04166666,
 rSphere         = 6371.D3
 &

# Input datasets: data.seaice has initial ice thickness field.
# Start from a previous iteration with both data and meta files available. Just do 'T.0000000720' (etc).
 &PARM05
 bathyFile      ='bathymetry.bin',
 hydrogThetaFile='tinit_9km_model.bin',
 hydrogSaltFile ='sinit_9km_model.bin',
 uVelInitFile   ='uinit_9km_model.bin',
 vVelInitFile   ='vinit_9km_model.bin',
 &