[MITgcm-support] advection scheme

David Ferreira dfer at mit.edu
Thu Sep 6 10:43:07 EDT 2007


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
> 





More information about the MITgcm-support mailing list