[MITgcm-support] Troubleshooting Gridding for Near N-Pole Simulation

[Benjamin Ocampo]

Hi Jean-Michel:

The message from using the debugLevel is in this form:

forrtl: severe (17): syntax error in NAMELIST input, unit 11, file
/data2/PBStmp/10708346.jasper-usradm/fortsx3FcS, line 26, position 2
Image              PC                Routine            Line        Source
libintlc.so.5      00002B0DBF0A9A1E  Unknown               Unknown  Unknown
libintlc.so.5      00002B0DBF0A84B6  Unknown               Unknown  Unknown
libifcore.so.5     00002B0DBE563E1E  Unknown               Unknown  Unknown
libifcore.so.5     00002B0DBE4D3A5F  Unknown               Unknown  Unknown
libifcore.so.5     00002B0DBE4D2F1D  Unknown               Unknown  Unknown
libifcore.so.5     00002B0DBE51622B  Unknown               Unknown  Unknown
mitgcmuv           00000000006B0CC5  Unknown               Unknown  Unknown
mitgcmuv           00000000006C7666  Unknown               Unknown  Unknown
mitgcmuv           00000000006FE0CF  Unknown               Unknown  Unknown
mitgcmuv           000000000061B0D5  Unknown               Unknown  Unknown
mitgcmuv           000000000040D8DC  Unknown               Unknown  Unknown
libc.so.6          000000318B41ED5D  Unknown               Unknown  Unknown
mitgcmuv           000000000040D7D9  Unknown               Unknown  Unknown

for:

ygOrigin = 64.0
delY = 50.0*0.5

Is there a sample data file I can check out for a spherical polar case that
analyzes GFD conditions encompassing the pole (i.e. theta_max = 90^oN)?
I seem to be having a hard time finding one to understand why this program
behaves this way.

Cheers,
Benjamin


On Fri, Apr 8, 2016 at 3:13 PM, Jean-Michel Campin <jmc at mit.edu> wrote:

> Hi Benjamin,
>
> I would recommand not to use:
> >  debugLevel=-1,
> but rather just the default value for parameter debugLevel.
> debugLevel=-1 prevent you to get some useful message/warning.
>
> Cheers,
> Jean-Michel
>
> On Fri, Apr 08, 2016 at 02:50:59PM -0600, Benjamin Ocampo wrote:
> > Hi All:
> >
> > I am a little confused about using MITgcm to model fluid behaviour near
> the
> > pole,
> >
> > When I set the following variables in the data file:
> >
> > ygOrigin = 66.0
> > delY = 50.0*0.5
> >
> > The program does not crash (i.e. the temperature average does not diverge
> > to NaN with increasing time). However, when I set the variables to:
> >
> > ygOrigin = 65.0
> > delY = 50.0*0.5
> >
> > or
> >
> > ygOrigin = 64.0
> > delY = 50.0*0.5
> >
> > the temperature average of the fluid diverges to NaN (typically to +
> > infinity at roughly 40 time steps). I am wondering if anyone knows why
> this
> > is the case.
> >
> > This is the current data file for the case that did work:
> >
> > # Model parameters
> > # Continuous equation parameters
> >  &PARM01
> >  tRef=40*25.,
> >  sRef=40*25.,
> >  viscAz=7.5E-2,
> >  viscAh=7.5E-2,
> >  no_slip_sides=.FALSE.,
> >  no_slip_bottom=.FALSE.,
> >  diffKhT=7.5E-2,
> >  diffKzT=7.5E-2,
> >  rotationPeriod=6283185,
> >  beta=0.0,
> >  tAlpha=1.3106E-12,
> >  sBeta =0.,
> >  gravity=9.81,
> >  rhoConst=1000.0,
> >  rhoNil=1000.0,
> >  saltStepping=.FALSE.,
> >  saltAdvection=.FALSE.,
> >  saltForcing=.FALSE.,
> >  rigidLid=.FALSE.,
> >  implicitFreeSurface=.TRUE.,
> >  nonHydrostatic=.TRUE.,
> >  eosType='LINEAR',
> >  readBinaryPrec=64,
> >  useSingleCpuIO=.TRUE.,
> >  debugLevel=-1,
> >  debugMode=.FALSE.,
> >  &
> > # Elliptic solver parameters
> >  &PARM02
> >  cg2dMaxIters=1000,
> >  cg2dTargetResidual=1.E-13,
> >  &
> > # Time stepping parameters
> >  &PARM03
> >  startTime=0.0,
> >  endTime=30000000.0,
> > # startTime=30000000.0,
> > # endTime=120000000.0,
> >  deltaT=750.0,
> > # deltaTmom=1200.0,
> > # deltaTtracer=1200.0,
> >  abEps=0.1,
> >  pChkptFreq=30000000.0,
> >  chkptFreq=0.0,
> >  dumpFreq=300000.0,
> >  monitorFreq=30000000.,
> >  &
> > # Gridding parameters
> >  &PARM04
> >  usingSphericalPolarGrid=.TRUE.,
> >  delZ = 14.6712, 29.2561, 43.6692, 57.8262, 71.6437,
> >         85.0408, 97.9389, 110.262, 121.939, 132.899,
> >         143.080, 152.421, 160.868, 168.370, 174.885,
> >         180.373, 184.803, 188.148, 190.389, 191.513,
> >         191.513, 190.389, 188.148, 184.803, 180.373,
> >         174.885, 168.370, 160.868, 152.421, 143.080,
> >         132.899, 121.939, 110.262, 97.9389, 85.0408,
> >         71.6437, 57.8262, 43.6692, 29.2561, 14.6712,
> >  ygOrigin=66.0,
> >  rSphere=6.37E+6,
> >  delX=360*0.5,
> >  delY=50*0.5,
> >  &
> >  &PARM05
> > # bathyFile='topog.box',
> >  hydrogThetaFile='ThetaR0.bin',
> >  uVelInitFile =  'uvelinit.bin',
> >  vVelInitFile =  'vvelinit.bin',
> >
> > Cheers,
> > Benjamin
>
> > _______________________________________________
> > 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 --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160408/d912a35b/attachment-0001.htm>