[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