[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