[MITgcm-support] advection scheme

[Patrick Heimbach]

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