[MITgcm-devel] prather and os7mp

[Jean-Michel Campin]

Hi Dimitris,

I don't remember paying attention to OBCs in Prather implementation.
Couldt it be the reason why it's not working ?
Jean-Michel

On Sun, May 20, 2007 at 08:17:53AM -0700, Dimitris Menemenlis wrote:
> Jean-Michel, I am trying os7mp and prather in the Arctic domain.  os7mp is 
> running but prather crashed after 18-24 hour integration.  Data file is 
> attached.  Anything specific I need to watch out for in getting prather to 
> work?
> 
> Dimitris

>  &PARM01
>  tRef               = 3*23.,3*22.,21.,2*20.,19.,2*18.,17.,2*16.,15.,14.,13.,
>                       12.,11.,2*9.,8.,7.,2*6.,2*5.,3*4.,3*3.,4*2.,12*1.,
>  sRef               = 50*34.5,
>  viscAr= 6.0826e-04
>  no_slip_sides=.FALSE.,
>  no_slip_bottom=.TRUE.,
>  diffKhT=0.,
>  diffK4T=0.,
>  diffKrT=1.4022e-05
>  diffKhS=0.,
>  diffK4S=0.,
>  diffKrS=1.4022e-05
>  tAlpha=2.E-4,
>  sBeta =7.4E-4,
>  gravity=9.81,
>  eosType='JMD95Z',
>  readBinaryPrec=32,
>  writeBinaryPrec=32,
>  useSingleCPUio=.TRUE.
>  globalfiles=.true.,
>  hFacMinDr=50.,
>  hFacMin=0.3,
>  hFacInf=0.1,
>  hFacSup=5.,
>  implicitDiffusion=.true.,
>  implicitViscosity=.true.,
>  allowFreezing=.false.,
>  viscAh=0.,
>  viscAhGrid=0.,
>  viscC4Leith=1.5,
>  viscC4Leithd=1.5,
>  viscA4GridMax=0.5,
>  viscC2Leith=0.,
>  viscA4=0.,
>  useAreaViscLength=.TRUE.,
>  sideDragFactor=0.,
>  SadournyCoriolis=.FALSE.,
>  useAbsVorticity=.FALSE.,
>  highOrderVorticity  = .TRUE.
>  bottomDragQuadratic = 0.0012,
>  tempAdvScheme=80,
>  saltAdvScheme=80,
>  StaggerTimeStep=.true.,
>  multiDimAdvection=.true.,
>  vectorInvariantMomentum=.TRUE.,
>  rigidLid=.FALSE.,
>  implicitFreeSurface=.TRUE.,
>  exactConserv=.TRUE.,
>  debuglevel=-1
>  &
> 
>  &PARM02
>  cg2dMaxIters=100,
>  cg2dTargetResidual=1.E-5,
>  &
> 
>  &PARM03
>  niter0=0
>  endtime=347155200.,
>  abEps=0.1,
>  deltaT=1200,
>  cAdjFreq=0.,
>  pChkptFreq=31536000,
>  chkptFreq=2592000
>  tauSaltClimRelax = 8722203.,
>  monitorFreq=21600
>  forcing_In_AB = .FALSE.
>  dumpInitAndLast=.FALSE.
>  &
> 
>  &PARM04
>  usingCartesianGrid=.FALSE.,
>  usingSphericalPolarGrid=.FALSE.,
>  usingCurvilinearGrid=.TRUE.,
>  delR   = 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.00, 10.01,
>  10.03, 10.11, 10.32, 10.80, 11.76, 13.42, 16.04 , 19.82, 24.85, 
>  31.10, 38.42, 46.50, 55.00, 63.50, 71.58, 78.90, 85.15, 90.18, 
>  93.96, 96.58, 98.25, 99.25,100.01,101.33,104.56,111.33,122.83,
>  139.09,158.94,180.83,203.55,226.50,249.50,272.50,295.50,318.50,
>  341.50,364.50,387.50,410.50,433.50,456.50,
>  &
> 
> # Input datasets
>  &PARM05
>  bathyFile      ='ETOPO2_420x384_arctic',
>  hydrogThetaFile='WGHC_T_420x384x50_arctic',
>  hydrogSaltFile ='WGHC_S_420x384x50_arctic',
>  &

> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel