[MITgcm-support] ploblem with a variable-grid global ocean circulation model
Martin Losch
Martin.Losch at awi.de
Tue May 8 03:50:19 EDT 2007
Hi lianzhan,
do I interpret this correctly: you use delx/dely of 0.1667 to 2 deg?
If that's the case then your horizontal viscosity is way too hight
> viscAh=1.E4,
> viscAhGrid=2,
When you look at mom_calc_visc.F you'll see that both viscAh and
viscAhGrid contribute to the overall viscosity. viscAh=1e4 is
suitable for a 2deg grid and to that you add even more viscosity, so
you may be running into a cfl-problem with the viscosities. Courtesy
of mom_calc_visc.F:
C harmonic viscosity=
C viscAh (or viscAhD on div pts and viscAhZ on zeta pts)
C +0.25*L**2*viscAhGrid/deltaT
So in your case 1e4 + 2 * .25*(2*pi/180*6371000)^2/60 = 4.1215e+08 to
2 * .25*(.1667*pi/180*6371000)^2/60 + 1e4 = 2.8733e+06
set viscAh=0 and viscAhGrid to a much smaller number so that your cfl
will not blow:
deltaT * (viscAh/L^2) = deltaT * viscAhGrid < 1
Or have a look at
http://mitgcm.org/r2_web_testing/latest/online_documents/node85.html
and subsections to figure out what you really need.
Martin
On 8 May 2007, at 08:35, lianzhan at fio.org.cn wrote:
> Hi mitgcm,
>
> I met a ploblem with a variable-grid global ocean circulation
> model. There are 546 grids in meridian direction,and 521 grids in
> latitudinal direction, as following. I seem to get noise in all the
> domain.
>
> Is the grid suitable for the MITgcm? What do I need to modify?
>
> That's all.Thanks for your help.
>
>
> The delx is
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 1.9961 1.9805 1.9495
>>
>> 1.9037 1.8439 1.7710 1.6864 1.5915 1.4879
>> 1.3774 1.2618 1.1432 1.0235
>>
>> 0.9049 0.7893 0.6788 0.5752 0.4802 0.3956
>> 0.3228 0.2630 0.2172 0.1862
>>
>> 0.1706 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1706
>>
>> 0.1862 0.2172 0.2630 0.3228 0.3956 0.4802
>> 0.5752 0.6788 0.7893 0.9049
>>
>> 1.0235 1.1432 1.2618 1.3774 1.4879 1.5915
>> 1.6864 1.7710 1.8439 1.9037
>>
>> 1.9495 1.9805 1.9961 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000
>>
>> The dely is
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 2.0000 2.0000
>>
>> 2.0000 2.0000 2.0000 2.0000 2.0000 2.0000
>> 2.0000 2.0000 1.9961 1.9805
>>
>> 1.9495 1.9037 1.8439 1.7710 1.6864 1.5915
>> 1.4879 1.3774 1.2618 1.1432
>>
>> 1.0235 0.9049 0.7893 0.6788 0.5752 0.4802
>> 0.3956 0.3228 0.2630 0.2172
>>
>> 0.1862 0.1706 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1667 0.1667 0.1667 0.1667
>> 0.1667 0.1667 0.1667 0.1667
>>
>> 0.1667 0.1667 0.1706 0.1862 0.2172 0.2630
>> 0.3228 0.3956 0.4802 0.5752
>>
>> 0.6788 0.7893 0.9049 1.0235 1.1432 1.2618
>> 1.3774 1.4879 1.5915 1.6864
>>
>> 1.7710 1.8439 1.9037 1.9495 1.9805 1.9961
>> 1.9961
>>
>>
>>
>> # ====================
>>
>> # | Model parameters |
>>
>> # ====================
>>
>> #
>>
>> # Continuous equation parameters
>>
>> &PARM01
>>
>> tRef= 16.0 , 15.2 , 14.5 , 13.9 , 13.3 ,
>>
>> 12.4 , 11.3 , 9.9 , 8.4 , 6.7 ,
>>
>> 5.2 , 3.8 , 2.9 , 2.3 , 1.8 ,
>>
>> 1.5 , 1.1 , 0.8 ,
>>
>> sRef= 34.65, 34.75, 34.82, 34.87, 34.90,
>>
>> 34.90, 34.86, 34.78, 34.69, 34.60,
>>
>> 34.58, 34.62, 34.68, 34.72, 34.73,
>>
>> 34.74, 34.73, 34.73,
>>
>> viscAz=1.E-4,
>>
>> viscAh=1.E4,
>>
>> viscAhGrid=2
>>
>> no_slip_sides=.FALSE.,
>>
>> no_slip_bottom=.TRUE.,
>>
>> diffKhT=1.E3,
>>
>> diffKzT=2.E-5,
>>
>> diffKhS=1.E3,
>>
>> diffKzS=2.E-5,
>>
>> tAlpha=2.E-4,
>>
>> sBeta =7.4E-4,
>>
>> gravity=9.81,
>>
>> rigidLid=.FALSE.,
>>
>> implicitFreeSurface=.TRUE.,
>>
>> useCDscheme=.TRUE.,
>>
>> useNHMTerms=.TRUE.,
>>
>> eosType='LINEAR',
>>
>> readBinaryPrec=32,
>>
>> tempStepping=.FALSE.,
>>
>> saltStepping=.FALSE.,
>>
>> implicitDiffusion=.TRUE.,
>>
>> implicitViscosity=.TRUE.,
>>
>> &
>>
>>
>>
>> # Elliptic solver parameters
>>
>> &PARM02
>>
>> cg2dMaxIters=1000,
>>
>> cg2dTargetResidual=1.E-13,
>>
>> &
>>
>> # Time stepping parameters
>>
>> &PARM03
>>
>> startTime=0,
>>
>>
>>
>> endTime=311040000,
>>
>> # endTime=864000.,
>>
>> deltaTmom=60.0,
>>
>> tauCD=321428.,
>>
>> deltaTtracer=3600.0,
>>
>> deltaTClock =3600.0,
>>
>> cAdjFreq=0.,
>>
>> abEps=0.1,
>>
>> pChkptFreq=311040000.,
>>
>> chkptFreq=311040000.,
>>
>> dumpFreq=3600.0,
>>
>> tauThetaClimRelax=2592000.0,
>>
>> tauSaltClimRelax=2592000.0,
>>
>> monitorFreq=1.,
>>
>> &
>>
>>
>>
>> # Gridding parameters
>>
>> &PARM04
>>
>> usingCartesianGrid=.FALSE.,
>>
>> usingSphericalPolarGrid=.TRUE.,
>>
>> delZ=4.,6.,10.,15.,23.,32.,48.,68.94,
>>
>> 88.1964,132.3541,216.3472,331.9537,
>>
>> 467.8573, 610.7548, 746.6584, 862.2648,
>>
>> 946.2580, 990.4156,
>>
>> phiMin=-80.,
>>
>> # delY=160*1,
>>
>> # delX=360*1,
>>
>> delXfile='dx.bin',
>>
>> delYfile='dy.bin',
>>
>> &
>>
>>
>>
>> # Input datasets
>>
>> &PARM05
>>
>> hydrogThetaFile='temp.bin',
>>
>> hydrogSaltFile='salt.bin',
>>
>> bathyFile='biangrid.bin',
>>
>> hydrogSaltFile='salt.bin',
>>
>> bathyFile='biangrid.bin',
>>
>> zonalWindFile='fx.bin',
>>
>> meridWindFile='fy.bin',
>>
>> thetaClimFile='sst.bin',
>>
>> saltClimFile='sss.bin',
>>
>> &
>>
>>
>>
>>
>>
>> lianzhan
>>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list