[MITgcm-support] Problem of loading bathmetry file
何映晖
heyinghui1984 at gmail.com
Tue May 18 12:48:44 EDT 2010
Dear Martin:
Thanks to your hint, I halved the grid numbers and then found the point.
I got something wrong in file *data*. The model is running correctly when I
replaced *dyspacing* by *delY.*
* *Best wishes.
Yinghui He
2010/5/15 Martin Losch <Martin.Losch at awi.de>
> Try to reproduce the problem with a much smaller domain (order 40x20x15)
> and then put the entire configuration ("code" directory and "input"
> directory) somewhere to download.
>
> Martin
>
> On May 14, 2010, at 6:50 PM, 何映晖 wrote:
>
> > Dear all,
> > I am using MITgcm to run a global ocean model with varible
> grid-space. In order to fit the configuration i want, i have modified the
> tutorial example "tutorial_global_latlon". after I checked the results, I
> found that only a part of grids have been assigned. There were 663*558
> grids, but only the former part(663*170 grids) have been correctly assigned
> in the depth output file, while the depth on rest grids are zero.
> > The "data" file is in attachment. And the following is "SIZE.h" file:
> >
> > C $Header:
> /u/gcmpack/MITgcm/verification/tutorial_global_oce_latlon/code/SIZE.h,v 1.1
> 2006/07/14 20:35:22 jmc Exp $
> > C $Name: $
> > C
> > C /==========================================================\
> > C | SIZE.h Declare size of underlying computational grid. |
> > C |==========================================================|
> > C | The design here support a three-dimensional model grid |
> > C | with indices I,J and K. The three-dimensional domain |
> > C | is comprised of nPx*nSx blocks of size sNx along one axis|
> > C | nPy*nSy blocks of size sNy along another axis and one |
> > C | block of size Nz along the final axis. |
> > C | Blocks have overlap regions of size OLx and OLy along the|
> > C | dimensions that are subdivided. |
> > C \==========================================================/
> > C Voodoo numbers controlling data layout.
> > C sNx - No. X points in sub-grid.
> > C sNy - No. Y points in sub-grid.
> > C OLx - Overlap extent in X.
> > C OLy - Overlat extent in Y.
> > C nSx - No. sub-grids in X.
> > C nSy - No. sub-grids in Y.
> > C nPx - No. of processes to use in X.
> > C nPy - No. of processes to use in Y.
> > C Nx - No. points in X for the total domain.
> > C Ny - No. points in Y for the total domain.
> > C Nr - No. points in Z for full process domain.
> > INTEGER sNx
> > INTEGER sNy
> > INTEGER OLx
> > INTEGER OLy
> > INTEGER nSx
> > INTEGER nSy
> > INTEGER nPx
> > INTEGER nPy
> > INTEGER Nx
> > INTEGER Ny
> > INTEGER Nr
> > PARAMETER (
> > & sNx = 221,
> > & sNy = 279,
> > & OLx = 3,
> > & OLy = 3,
> > & nSx = 1,
> > & nSy = 1,
> > & nPx = 3,
> > & nPy = 2,
> > & Nx = sNx*nSx*nPx,
> > & Ny = sNy*nSy*nPy,
> > & Nr = 25)
> > C MAX_OLX - Set to the maximum overlap region size of any array
> > C MAX_OLY that will be exchanged. Controls the sizing of exch
> > C routine buufers.
> > INTEGER MAX_OLX
> > INTEGER MAX_OLY
> > PARAMETER ( MAX_OLX = OLx,
> > & MAX_OLY = OLy )
> > Any helps will be appreciate!
> >
> > --
> > Yinghui He
> > South China Sea Institute of Oceanology,
> > China Academic of Science
> > <data>_______________________________________________
> > 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
>
--
Yinghui He
South China Sea Institute of Oceanology,
China Academic of Science
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20100519/e536ec33/attachment.htm>
More information about the MITgcm-support
mailing list