[MITgcm-support] advection scheme
Matthew Mazloff
mmazloff at MIT.EDU
Thu Sep 6 13:34:05 EDT 2007
Hi Julie,
when in the directory you want to update you can use
cvs update -d -P -A
for the latest or
cvs update -d -P -r release1_p5
for a specific version
see: http://mitgcm.org/using_cvs.html
for more info.
Caution though, updating one package often leads to version
incompatibilities...it is not so straightforward...so before updating
a package make a backup
or, if you can, just update the whole model....that usually is fine
-Matt
On Sep 6, 2007, at 12:47 PM, Julie Deshayes wrote:
> 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
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20070906/f39fce93/attachment.htm>
More information about the MITgcm-support
mailing list