[MITgcm-support] Query regarding overlapping in output files

Martin Losch Martin.Losch at awi.de
Sat Jun 11 11:20:03 EDT 2011


That's also my impression. Matlab's index convention is based on matrix algebra (so i=row, j=column). In Fortran (and in the MITgcm) the first index (i) is the x-direction and the second (j) is the y-direction, so if you plot any output file like Depth.data you have to transpose:
b=rdmds('Depth');
contour(b')
to get the right plot and you also need to save the input in the correct way as indicated by Serge.

Another thing to be aware of: The model bathymetry will be different from the one you prescribed, because there are some restrictions on the minimum cell thickness (hFacMin=0.2 sets these restrictions), but these differences are small.

Martin

On Jun 11, 2011, at 6:42 AM, fancer fancer wrote:

> Hi Himansu Pradhan,
> 
> You should be careful with initial data. You set inverted bathymetry data. Correct a bathy_binary.м on:
> ieee='b';accuracy='real*8';
> t = load('bathy2DIM.txt');
> fid=fopen('bathyBOB.bin','wb',ieee);
> fwrite(fid,t',accuracy);
> fclose(fid);
> here t' - inverted matrix.
> 
> Yours sincerely,
> M.S. student of NNSTU named after Alekseev R.E.
> Nizhny Novgorod, Russia
> Serge Semin
> 
> On Sat, Jun 11, 2011 at 3:41 PM, himansu pradhan <oceancalling at gmail.com> wrote:
> Hello MITgcm users, 
> I am a new user to MITgcm.
> I am trying to configure MITgcm for observing Internal waves propagation for my region (Bay of Bengal,Indian Ocean).
> 
> 
> 
> Initially, my input files are (a) bathymetry (360*270 matrix)(b) 3 levels temperature file (c) Wind files. 
> My /input/data file looks like : 
> ====================
> # | Model parameters |
> 
> # ====================
> 
> 
> #
> 
> # Continuous equation parameters
>  &PARM01
>  tRef= 26.29,24.29,22.29,
>  sRef= 3*35.,
>  viscAz=1.E-3,
>  viscAh=1.E-2,
>  no_slip_sides=.TRUE.,
>  no_slip_bottom=.TRUE.,
>  diffKhT=1.E-2,
>  diffKzT=1.E-3,
> 
> 
>  diffKhS=1.E-2,
> 
>  diffKzS=1.E-5,
>  f0=0,
>  beta=0,
>  gravity=9.81,
>  rigidLid=.FALSE.,
>  implicitFreeSurface=.TRUE.,
>  exactConserv=.TRUE.,
>  eosType='LINEAR',
>  hFacMin=0.2,
>  nonHydrostatic=.FALSE.,
> 
> 
>  readBinaryPrec=64,
> 
>  writeBinaryPrec=64,
>  globalFiles=.TRUE.,
>  &
> 
> # Elliptic solver parameters
>  &PARM02
>  cg2dMaxIters=1000,
>  cg2dTargetResidual=1.E-13,
>  cg3dMaxIters=50,
>  cg3dTargetResidual=1.E-13,
> 
> 
>  &
> 
> 
> # Time stepping parameters
>  &PARM03
> startTime=0,
>  endTime=1080,
>  deltaTmom=10,
>  deltaTtracer=10,
>  abEps=0.02,
>  pChkptFreq=4000000.0,
>  chkptFreq=0.0,
>  dumpFreq=20,
>  monitorFreq=20,
> 
> 
>  &
> 
> 
> # Gridding parameters
>  &PARM04
>  usingCartesianGrid=.TRUE.,
>  usingSphericalPolarGrid=.FALSE.,
>  dXspacing=3600,
>  dYspacing=3600,
>  delZ=3*50,
>  &
> 
> # Input datasets
>  &PARM05
> 
> 
>  bathyFile='bathyBOB.bin',
> 
>  hydrogThetaFile='temBOB3lvl.bin',
>  zonalWindFile='WINuBOB.bin',
>  meridWindFile='WINvBOB.bin',
>  checkIniTemp=.false.,
>  &
> 
> My query is: the output depth files (Depth.data and Depth.meta) plots is seen different to that of input bathymetry file(bathyBOB.bin,bathy2DIM.txt). I am confused 
> 
> 
> : Is the model overwritting/overlapping the output file or anything else......?????? Do help.
> 
> 
> The images of model output and input bathymetry is in the attached tar file (himansu_project.tar.gz). The other required informations ( code and input directory)
> 
> 
> is also in the tar file.
> 
> 
> Thanking you in advance
> 
> Himansu , PhD student
> Indian Institute of Technology,Delhi
> India.
> 
> _______________________________________________
> 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




More information about the MITgcm-support mailing list