[MITgcm-support] advection scheme
Matthew Mazloff
mmazloff at MIT.EDU
Thu Sep 6 11:34:44 EDT 2007
Hi Julie,
Make sure you don't have
autodiff
in packages.conf
-Matt
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
More information about the MITgcm-support
mailing list