[MITgcm-support] ploblem with a variable-grid global ocean circulation model

lianzhan at fio.org.cn lianzhan at fio.org.cn
Tue May 8 02:35:49 EDT 2007


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
>




More information about the MITgcm-support mailing list