[MITgcm-support] format of input files wrong in regional model with open boundary conditions?

Jonny Williams Jonny.Williams at bristol.ac.uk
Tue Dec 10 10:12:40 EST 2013


Hi Patrick

Thank you very much for your email.

What I am actually doing now is changing the model so that I can use it for
regional studies in paleoclimate.

I believe I have created all the necessary input files that I will need to
replace those supplied with the Arctic regional model. I am firstly trying
to use an un-varying forcing just to get the model running.

Given what you have just said , would it therefore be reasonable to assume
that the dimensions of the input files I have created are wrong?

Many thanks!

Jonny






On 10 December 2013 15:03, Patrick Heimbach <heimbach at mit.edu> wrote:

> Hi Jonny,
>
> if I remember correctly you are running the regional 36 km (or 18 km)
> Arctic, right?
> If so, you are probably running longer than the available open boundary
> conditions provided with the setup
> (e.g., there are 15 years of monthly OB's in the file, but you are trying
> to run 20 years).
> Eventually, the model is trying to access a non-existing record.
>
> p.
>
> On Dec 10, 2013, at 9:55 AM, Jonny Williams <Jonny.Williams at bristol.ac.uk>
> wrote:
>
> > Hi there
> >
> > I am trying to get a regional model version of the MITgcm running but am
> having a problem which I have been stuck on for the last few days.
> >
> > I can seemingly compile my model successfully but at run time I get the
> following error...
> >
> > PGFIO-F-253/unformatted read/unit=9/attempt to read non-existent record
> (direct access).
> >  File name = OBNu.bin    unformatted, direct access   record = 556
> >  In source file mdsio_read_section.f, at line number 1851
> >
> > I do not understand this because my input file OBNu.bin is definitely
> present in the run directory.
> >
> > Line 1851 of mdsio_read_section.f is illustrated in bold, underlined
> font in the snippet below...
> >
> > Does anyone have any ideas about what the issue could be with this?
> >
> > Many thanks!
> >
> >         IF (fileIsOpen) THEN
> >          DO k=1,kSize
> >            IF (globalFile) THEN
> >             iG = myXGlobalLo-1 + (bi-1)*sNx
> >             jG = (myYGlobalLo-1)/sNy + (bj-1)
> >             irec=1 + INT(iG/sNx) + nSx*nPx*(k-1)
> >      &           + nSx*nPx*kSize*(irecord-1)
> >            ELSE
> >             iG = 0
> >             jG = 0
> >             irec=k + kSize*(irecord-1)
> >            ENDIF
> >            IF (filePrec .EQ. precFloat32) THEN
> >             READ(dUnit,rec=irec) r4seg
> >             IF (arrType .EQ. 'RS') THEN
> >              CALL MDS_SEG4toRS_2D( sNx,oLx,kSize,bi,bj,k,.TRUE.,
> >      &                             r4seg,fldRS )
> >             ELSEIF (arrType .EQ. 'RL') THEN
> >              CALL MDS_SEG4toRL_2D( sNx,oLx,kSize,bi,bj,k,.TRUE.,
> >      &                             r4seg,fldRL )
> >             ELSE
> >              WRITE(msgBuf,'(A)')
> >      &         ' MDS_READ_SEC_XZ: illegal value for arrType'
> >              CALL PRINT_ERROR( msgBuf, myThid )
> >              STOP 'ABNORMAL END: S/R MDS_READ_SEC_XZ'
> >             ENDIF
> >            ELSEIF (filePrec .EQ. precFloat64) THEN
> >             READ(dUnit,rec=irec) r8seg
> >             IF (arrType .EQ. 'RS') THEN
> >              CALL MDS_SEG8toRS_2D(sNx,oLx,kSize,bi,bj,k,.TRUE.,
> >      &                             r8seg,fldRS )
> >             ELSEIF (arrType .EQ. 'RL') THEN
> >              CALL MDS_SEG8toRL_2D(sNx,oLx,kSize,bi,bj,k,.TRUE.,
> >      &                             r8seg,fldRL )
> >             ELSE
> >              WRITE(msgBuf,'(A)')
> >      &         ' MDS_READ_SEC_XZ: illegal value for arrType'
> >              CALL PRINT_ERROR( msgBuf, myThid )
> >              STOP 'ABNORMAL END: S/R MDS_READ_SEC_XZ'
> >             ENDIF
> >            ELSE
> >             WRITE(msgBuf,'(A)')
> >      &        ' MDS_READ_SEC_XZ: illegal value for filePrec'
> >             CALL PRINT_ERROR( msgBuf, myThid )
> >             STOP 'ABNORMAL END: S/R MDS_READ_SEC_XZ'
> >            ENDIF
> > C End of k loop
> >          ENDDO
> >
> >
> > --
> > Dr Jonny Williams
> > School of Geographical Sciences
> > University of Bristol
> > University Road
> > BS8 1SS
> >
> > +44 (0)117 3318352
> > jonny.williams at bristol.ac.uk
> > bit.ly/jonnywilliams
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> ---
> Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
> MIT | EAPS 54-1420 | 77 Massachusetts Ave | Cambridge MA 02139 USA
> FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>


-- 
Dr Jonny Williams
School of Geographical Sciences
University of Bristol
University Road
BS8 1SS

+44 (0)117 3318352
jonny.williams at bristol.ac.uk
bit.ly/jonnywilliams
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20131210/9a13500c/attachment-0001.htm>


More information about the MITgcm-support mailing list