<div dir="ltr"><div>Thank you Gus. Martin, rightly pointed out, the error was caused in the NOBC boundary file as I had 1 open point (1x48x2), but it was not in the required 660x48x2 format. I added zeros for every other point and it worked! Thanks a lot Martin.</div><div><br></div><div>I am trying to keep 19 open points in NOBC boundary through data.obcs:</div><div> OB_Jnorth(550)=-100,<br> OB_Jnorth(551)=-100,<br> OB_Jnorth(552)=-99,<br> OB_Jnorth(553)=-98,<br> OB_Jnorth(554)=-98,<br> OB_Jnorth(555)=-97,<br> OB_Jnorth(556)=-97,<br> OB_Jnorth(557)=-97,<br> OB_Jnorth(558)=-97,<br> OB_Jnorth(559)=-97,<br> OB_Jnorth(560)=-97,<br> OB_Jnorth(561)=-96,<br> OB_Jnorth(562)=-96,<br> OB_Jnorth(563)=-95,<br> OB_Jnorth(564)=-95,<br> OB_Jnorth(565)=-95,<br> OB_Jnorth(566)=-95,<br> OB_Jnorth(567)=-94,<br> OB_Jnorth(568)=-94,<br> OB_Jsouth=660*1,<br> OB_Iwest=480*1,<br> OB_Ieast=480*-1,</div><div><br></div><div>Is there a simpler way to feed it? I referred to the OBCS section in the manual and tried a few things but it didn't work out. Any suggestions?</div><div><br></div><div>Regards,</div><div>Kunal<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 19, 2021 at 7:32 AM Martin Losch <<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Kunal,<br>
<br>
the error message (as usual) means that there are not enough records in your input file (‘OBNu.seaice_obcs’)<br>
you domain dimensions are (according to STDOUT.0000):<br>
(PID.TID 0000.0001) nPx = 10 ; /* No. processes in X */<br>
(PID.TID 0000.0001) nPy = 6 ; /* No. processes in Y */<br>
(PID.TID 0000.0001) nSx = 1 ; /* No. tiles in X per process */<br>
(PID.TID 0000.0001) nSy = 1 ; /* No. tiles in Y per process */<br>
(PID.TID 0000.0001) sNx = 66 ; /* Tile size in X */<br>
(PID.TID 0000.0001) sNy = 80 ; /* Tile size in Y */<br>
(PID.TID 0000.0001) OLx = 4 ; /* Tile overlap distance in X */<br>
(PID.TID 0000.0001) OLy = 4 ; /* Tile overlap distance in Y */<br>
(PID.TID 0000.0001) nTx = 1 ; /* No. threads in X per process */<br>
(PID.TID 0000.0001) nTy = 1 ; /* No. threads in Y per process */<br>
(PID.TID 0000.0001) Nr = 48 ; /* No. levels in the vertical */<br>
(PID.TID 0000.0001) Nx = 660 ; /* Total domain size in X ( = nPx*nSx*sNx ) */<br>
(PID.TID 0000.0001) Ny = 480 ; /* Total domain size in Y ( = nPy*nSy*sNy ) */<br>
<br>
So your input file for the OBCS pkg for the northern or southern boundary for a 3D field like temperature needs to have a field with the dimensions (t,z,x) = (nrec,48,660) and for the eastern/western boundaries (nrec,48,480). nrec needs to be at least 1 (for constant obcs forcing). This is independent of the number of actual boundary points. I.e. if you have only one open boundary point (OB_Jnorth(550)=-100,) your file still need the entire array, it can be zero everywhere except for the point i=550 (for i starting at 1, for i starting a 0 as in python, it should be at 549).<br>
<br>
For a 2D field like sea ice velocity you drop the z-dimension (but I guess you only give you files these names, and don’t use sea ice).<br>
<br>
does that make sense?<br>
<br>
Martin<br>
> On 18. Aug 2021, at 15:44, kunal madkaiker <<a href="mailto:kunal.madkaiker02@gmail.com" target="_blank">kunal.madkaiker02@gmail.com</a>> wrote:<br>
> <br>
> Dear All,<br>
> <br>
> I am using MITgcm for regional analysis. I have kept a single OB point on the northern boundary. But I keep getting backtrace error and model blows.<br>
> <br>
> Error termination. Backtrace:<br>
> At line 1946 of file mdsio_read_section.f (unit = 9, file = 'OBNu.seaice_obcs')<br>
> Fortran runtime error: Non-existing record number<br>
> <br>
> I have gone through similar asked questions on the support group and checked accordingly. The dimensions for my NOBCS pt is 1x48x2 (nx x nz x nt). I have prepared the T,S,U and V data and written as real*8. But still error persists.<br>
> Attaching my data, data.exf, data.obcs files. Kindly assist.<br>
> <br>
> Regards,<br>
> Kunal<br>
> <data><data.obcs><data.exf><STDOUT.0000>_______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>