[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 09:55:19 EST 2013


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<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=jonny.williams@bristol.ac.uk>
bit.ly/jonnywilliams
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20131210/a31c5b97/attachment-0002.htm>


More information about the MITgcm-support mailing list