[MITgcm-support] advection scheme

Julie Deshayes jdeshayes at whoi.edu
Thu Sep 6 12:47:31 EDT 2007


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
> 




More information about the MITgcm-support mailing list