[MITgcm-support] bathy file with Spherical Polar Grid

[Lequan Chi]

Hi, Martin

Thank you for your reply. Here is the related part of the "data" file.  The 
Luzon_3d.topo is the bathy file derived from ETOP1 
(http://www.ngdc.noaa.gov/mgg/gdas/gd_designagrid.html?dbase=grdet1), whose 
value is negative in sea area.

The Luzon_3d.topo.mat ,which is more convenient to read by matlab, is the 
.mat version of Luzon_3d.topo.

========================================
# Gridding parameters
&PARM04
usingCartesianGrid=.FALSE.,
usingSphericalPolarGrid=.TRUE.,
delZ=50*10.,20*50.,20*150.,
ygOrigin=-14.649090,
# xgOrigin= 63.898265194,
zCoordInputData=.TRUE.,
delXfile='dlon.init',
# delX=3424*250.,
delYfile='dlat.init',
# delY=158*1000.,
&

# Input datasets
&PARM05
bathyFile='Luzon_3d.topo',
========================================
And another param you may need:
readBinaryPrec=64,

By the way, I have tried to use dXspaing and dYspacing instead of delXfile 
and delYfile on account of the nonuniform of the grid. But it is still 
useless.

Lequan Chi

-----原始邮件----- 
From: Martin Losch
Sent: Thursday, October 27, 2011 9:28 PM
To: mitgcm-support at mitgcm.org
Subject: Re: [MITgcm-support] bathy file with Spherical Polar Grid

Hi,

the bathymetry file needs to have value < 0 to be recognized as ocean depth.
You need to specify the grid spacing in both cartesian (in meters) and 
spherical grid (in deg), please send your file "data".

Martin

On Oct 27, 2011, at 3:07 PM, Lequan Chi wrote:

> Hi All,
>
> I am trying to use spherical ploar grid in my case. To do this, I only 
> turn‘usingSphericalPolarGrid’on ( usingCartesianGrid=.FALSE. at the same 
> time ), and give a value to ygOrigin. Nothing else has been done. The 
> MITgcm just runs normally, no err. But I find the values in 
> “Depth.***.***.data”are all set to 0. Then I set the 
> ‘usingCartesianGrid’to .TRUE. ( usingSphericalPolarGrid=.FALSE. ) and 
> then restart the model. This time the values in “Depth.***.***.data”files 
> are right.
>
> I use the same barthy file during the two experiments above. I can’t 
> understand why values in Depth.***.***.data are all 0 when I use Spherical 
> Polar Grid.
>
> Anything else is needed in the spherical polar grid?
>
>
> Wish best wishes
>
> Lequan Chi,
> _______________________________________________
> 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 --------------
A non-text attachment was scrubbed...
Name: dlat.init
Type: application/octet-stream
Size: 1536 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20111028/965582f7/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dlon.init
Type: application/octet-stream
Size: 25472 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20111028/965582f7/attachment-0005.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Luzon_3d.topo
Type: application/octet-stream
Size: 4890624 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20111028/965582f7/attachment-0006.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Luzon_3d.topo.mat
Type: application/octet-stream
Size: 4358043 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20111028/965582f7/attachment-0007.obj>