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

[Patrick Heimbach]

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

-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1588 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20131210/c3a6249c/attachment.p7s>