[MITgcm-support] advection scheme

[Jean-Michel Campin]

Hi,

I would recommand a full update, since lot of little things
have changed (e.g., rhoFacC, deepFac2C ... were not there at all
in Jul 06, but needed now).
And recommand to try first:
1) cvs -n update (don't do anything, just list the files that need
to be updated + possible conflict)
2) cvs update -P -d 
3) if you have conflict that are hard to solve,
  cvs update -P -d -A (force to get the current version)

Jean-Michel

On Thu, Sep 06, 2007 at 12:47:31PM -0400, Julie Deshayes wrote:
> indeed ...
> C $Header: /u/gcmpack/MITgcm/pkg/generic_advdiff/gad_advection.F,v 1.42 
> 2006/07/18 03:23:50 heimbach Exp $
> 
> what is the easiest way to update the package ?
> 
> Julie
> 
> Patrick Heimbach wrote:
> >
> >Hi there,
> >
> >ah colleagues, don't always blame the adjoint ;o)
> >(just kidding).
> >
> >All parameter files look ok.
> >Here's my suspicion.
> >Julie, you have code that pre-dates
> >Alistair's addition of the OS7MP code.
> >That would have been January 20th, 2007:
> >http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm/pkg/generic_advdiff/gad_advection.F 
> >
> >Could you check?
> >
> >-Patrick
> >
> >
> >
> >On Sep 6, 2007, at 10:43 AM, David Ferreira wrote:
> >
> >>Julie,
> >>From a quick look at gad_advection.F, it seems that to obtain
> >>the first error you described, ALLOW_AUTODIFF_TAMC would need to
> >>be defined (so that OS7MP is skipped and you end up in the STOP).
> >>Are you using the adjoint ?
> >>david
> >>
> >>
> >>Julie Deshayes wrote:
> >>>Hi Dimitris and other MITgcm users,
> >>>Following your advice, I tried using tempAdvScheme=7. At the first 
> >>>time step of calculation, I obtained the following error message :
> >>>STOP GAD_ADVECTION: adv. scheme incompatibale with multi-dim 
> >>>statement executed
> >>>when adding "multiDimAdvection=.FALSE.," in my data file, I then have 
> >>>the following error message (always at first time step of the 
> >>>calculation) :
> >>>STOP GAD_CALC_RHS: Bad advectionScheme (X) statement executed
> >>>There seem to be an incompatibility in my configuration which impedes 
> >>>me from using tempAdvScheme=7. Do you have an idea which it may be ? 
> >>>at the end of this email is my full data file, data.pkg and 
> >>>packages.conf files.
> >>>thanks for your help
> >>>Julie
> >>>Dimitris Menemenlis wrote:
> >>>>Julie, you can also try Alistai's OS7MP scheme:
> >>>>tempAdvScheme=7,
> >>>>saltAdvScheme=7,
> >>>>It is much less diffusive than 33.
> >>>>
> >>>>D.
> >>>>_______________________________________________
> >>>>MITgcm-support mailing list
> >>>>MITgcm-support at mitgcm.org
> >>>>http://mitgcm.org/mailman/listinfo/mitgcm-support
> >>>>
> >>>************************** data file
> >>># Model parameters
> >>>#
> >>># Continuous equation parameters
> >>> &PARM01
> >>> readBinaryPrec=64,
> >>> useSingleCpuIO=.TRUE.,
> >>> tRef=  5.35, 5.20, 5.05, 4.9, 4.75, 4.6, 4.45, 4.3, 4.15, 4.0,
> >>> sRef=10*10.,
> >>> rigidLid=.FALSE.,
> >>> implicitFreeSurface=.TRUE.,
> >>> hFacMin=0.2,
> >>> no_slip_sides=.TRUE.,
> >>> no_slip_bottom=.TRUE.,
> >>> bottomDragLinear=1.e-6,
> >>> eosType='LINEAR',
> >>> tAlpha=2.E-4,
> >>> sBeta =0.,
> >>> saltStepping=.FALSE.,
> >>> tempAdvScheme=7,
> >>> saltAdvScheme=7,
> >>> multiDimAdvection=.FALSE.,
> >>> diffKhT=0.,
> >>> diffK4T=1.E10,
> >>> diffKzT=1.E-5,
> >>> f0=1.0e-4,
> >>> beta=0.E-11,
> >>># useCDscheme=.FALSE.,
> >>> useNHMTerms=.FALSE.,
> >>> implicitDiffusion=.TRUE.,
> >>> ivdc_kappa=1000.0,
> >>> viscAz=2.E-4,
> >>> viscAh=100.,
> >>> &
> >>>#
> >>># Elliptic solver parameters
> >>> &PARM02
> >>> cg2dMaxIters=300,
> >>> cg2dTargetResidual=1.E-8,
> >>> &
> >>>#
> >>># Time stepping parameters
> >>> &PARM03
> >>>startTime=622080000.,
> >>> endTime=653184000.,,
> >>> deltaT=900.0,
> >>> abEps=0.1,
> >>> pChkptFreq=0.0,
> >>> chkptFreq=311040000.,
> >>> dumpFreq=311040000.,
> >>># taveFreq=86400.,
> >>> cAdjFreq=0,
> >>> &
> >>>#
> >>># Gridding parameters
> >>> &PARM04
> >>> usingCartesianGrid=.TRUE.,
> >>> usingSphericalPolarGrid=.FALSE.,
> >>> dXspacing=10.e+3,
> >>> dYspacing=10.e+3,
> >>> delZ=10*100.,
> >>> &
> >>>#
> >>># topography, bathymetry and initial fields
> >>> &PARM05
> >>> bathyFile='topog_10km_obcs.bin',
> >>> hydrogThetaFile='',
> >>> hydrogSaltFile=,
> >>> zonalWindFile=,
> >>> meridWindFile=,
> >>> uVelInitFile='',
> >>> vVelInitFile='',
> >>> surfQfile='Qnet_10km_obcs.bin',
> >>> &
> >>>#
> >>># heat flux forcing
> >>> &PARM06
> >>> HfluxP=360.,
> >>> HfluxW=60.,
> >>> qmax=1.2,
> >>> &
> >>>************************** data.pkg
> >>># Packages
> >>> &PACKAGES
> >>> useOBCS=.TRUE.,
> >>> usediagnostics=.TRUE.,
> >>> &
> >>>************************** packages.conf
> >>>generic_advdiff
> >>>obcs
> >>>diagnostics
> >>>mdsio
> >>>rw
> >>>mom_fluxform
> >>>_______________________________________________
> >>>MITgcm-support mailing list
> >>>MITgcm-support at mitgcm.org
> >>>http://mitgcm.org/mailman/listinfo/mitgcm-support
> >>
> >>
> >>_______________________________________________
> >>MITgcm-support mailing list
> >>MITgcm-support at mitgcm.org
> >>http://mitgcm.org/mailman/listinfo/mitgcm-support
> >
> >---
> >Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
> >MIT | EAPS 54-1518 | 77 Massachusetts Ave | Cambridge MA 02139 USA
> >FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach
> >
> >
> >_______________________________________________
> >MITgcm-support mailing list
> >MITgcm-support at mitgcm.org
> >http://mitgcm.org/mailman/listinfo/mitgcm-support
> >
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support