<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi Estanislao<div class=""><br class=""></div><div class="">Glad you got it all sorted out.</div><div class=""><br class=""></div><div class="">One other suggestion; try and avoid odd-sized dimensions like 43. Some optimizers have issues with odd-numbered dimensions. I had a heck of a problem trying to debug a model run that kept blowing up near tile boundaries until I figured out that th problem was that my grid size was 61, and changing to 62 made it all run fine. </div><div class=""><br class=""></div><div class="">Cheers, Jody<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Aug 26, 2019, at 1:26 AM, Estanislao Gavilan Pascual-Ahuir <<a href="mailto:dramauh@hotmail.com" class="">dramauh@hotmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Helvetica, sans-serif; font-size: 12pt;" class="">Hi Michael and Martin</div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Helvetica, sans-serif; font-size: 12pt;" class="">I figured out the problem, Michael pointed me in the correct way. I did not notice that the forcing cycle cannot be bigger than the forcing period times the number of time steps of your forcing files (ie.1200x100). I thought that if I set a big forcing cycle the forcing files were going to keep constant once the model gets to the last time step. In this experiment, I was trying to ramp up a source-sink circulation during the first 100 time steps. Then, I wanted to keep it constant (That is why I set a huge period cycle). It seems I will need to used restart files to do this. Thanks so much for your time and all your help. I would not be able to figure out this by myself</div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Helvetica, sans-serif; font-size: 12pt;" class="">Kind regards,</div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Helvetica, sans-serif; font-size: 12pt;" class="">Estanislao</div><div id="appendonsend" style="caret-color: rgb(0, 0, 0); font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div class="BodyFragment" style="caret-color: rgb(0, 0, 0); font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><font size="2" class=""><span style="font-size: 11pt;" class=""><div class="PlainText"><br class=""><br class=""><br class="">----------------------------------------------------------------------<br class=""><br class="">Message: 1<br class="">Date: Sun, 25 Aug 2019 21:04:55 -0500<br class="">From: Michael Schaferkotter <<a href="mailto:schaferk@bellsouth.net" class="">schaferk@bellsouth.net</a>><br class="">To:<span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a><br class="">Subject: Re: [MITgcm-support] Reading errors<br class="">Message-ID: <<a href="mailto:0A6D5C87-53FC-4FF9-B7E7-34CBB169C096@bellsouth.net" class="">0A6D5C87-53FC-4FF9-B7E7-34CBB169C096@bellsouth.net</a>><br class="">Content-Type: text/plain; charset="us-ascii"<br class=""><br class="">0. i misread the original post and had the idea that the simulation was 2d.<br class="">that/s why there was the comment about setting delY.<br class=""><br class="">1. the parameters given in SIZE.h imply that you should use 43 in the MATLAB scripts as stated by another.<br class=""><br class="">2.<br class="">> & OLx = 2,<br class="">> & OLy = 2,<br class=""><br class=""><br class="">i never found success with using these values (2), so i always set them to 3.<br class=""><br class="">3.<br class="">> > > > OB_Jnorth(16:28)=13*41,<br class="">> > > > > OB_Jsouth(16:28)=13*3,<br class=""><br class="">namelists are not standard and can cause problems for various reasons.<br class=""><br class="">so i would always explicitly list the values using only -1,0,1 (-1 on north and east)<br class=""><br class="">OB_Jnorth=15*0,13*-1,15*0<br class="">OB_Jsouth=15*0,13*1,15*0<br class=""><br class="">4.<br class="">i never cared for using reentrant methods and would force the boundaries with values such that the transport was corrected for incoming signals. For outgoing signals and to reduce boundary reflections, a method by Lavelle and Thacker was incorporated into the OBCS package.<br class=""><br class="">5. the periodic forcing parameters must be consistent. (i never used periodic forcing)<br class=""><br class=""><br class="">> On Aug 24, 2019, at 21:11, Estanislao Gavilan Pascual-Ahuir <<a href="mailto:dramauh@hotmail.com" class="">dramauh@hotmail.com</a>> wrote:<br class="">><span class="Apple-converted-space"> </span><br class="">> Hi Michael,<br class="">><span class="Apple-converted-space"> </span><br class="">> My SIZE.h file is this one<span class="Apple-converted-space"> </span><br class="">><span class="Apple-converted-space"> </span><br class="">> INTEGER sNx<br class="">> INTEGER sNy<br class="">> INTEGER OLx<br class="">> INTEGER OLy<br class="">> INTEGER nSx<br class="">> INTEGER nSy<br class="">> INTEGER nPx<br class="">> INTEGER nPy<br class="">> INTEGER Nx<br class="">> INTEGER Ny<br class="">> INTEGER Nr<br class="">> PARAMETER (<br class="">> & sNx = 43,<br class="">> & sNy = 43,<br class="">> & OLx = 2,<br class="">> & OLy = 2,<br class="">> & nSx = 1,<br class="">> & nSy = 1,<br class="">> & nPx = 1,<br class="">> & nPy = 1,<br class="">> & Nx = sNx*nSx*nPx,<br class="">> & Ny = sNy*nSy*nPy,<br class="">> & Nr = 2)<br class="">> C MAX_OLX :: Set to the maximum overlap region size of any array<br class="">> C MAX_OLY that will be exchanged. Controls the sizing of exch<br class="">> C routine buffers.<br class="">> INTEGER MAX_OLX<br class="">> INTEGER MAX_OLY<br class="">> PARAMETER ( MAX_OLX = OLx,<br class="">> & MAX_OLY = OLy )<br class="">><span class="Apple-converted-space"> </span><br class="">> About the second comment. I set it just in case because I found in the example. About the OBCS, the indexes can be other than 1,0 or -1. My guess, if you have a big domain and you want to open a boundary in the last row it is easier to put -1 than 5984. For my experiments, I have designed a circular domain and I decided to leave the two first and last rows and columns as land. I do not think this is causing the problem because if I run the model without forcing period the model works perfectly fine. He reads only the first time step of my forcing bin file but it works. It is just when I activate the forcing period in the data namelist.<span class="Apple-converted-space"> </span><br class="">><span class="Apple-converted-space"> </span><br class="">><span class="Apple-converted-space"> </span><br class="">> Kind regards,<br class="">><span class="Apple-converted-space"> </span><br class="">> Estanislao<br class="">><span class="Apple-converted-space"> </span><br class="">> PS:thank you in advance for all your help. I am really stuck with this matter.<br class="">><span class="Apple-converted-space"> </span><br class="">> ----------------------------------------------------------------------<br class="">><span class="Apple-converted-space"> </span><br class="">> Message: 1<br class="">> Date: Sat, 24 Aug 2019 09:49:09 -0500<br class="">> From: Michael Schaferkotter <<a href="mailto:schaferk@bellsouth.net" class="">schaferk@bellsouth.net</a>><br class="">> To:<span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a><br class="">> Subject: Re: [MITgcm-support] MITgcm-support Digest, Vol 194, Issue 16<br class="">> Message-ID: <<a href="mailto:ACF3B949-4FA8-403F-A0F4-FCFFEB5B7616@bellsouth.net" class="">ACF3B949-4FA8-403F-A0F4-FCFFEB5B7616@bellsouth.net</a>><br class="">> Content-Type: text/plain; charset="utf-8"<br class="">><span class="Apple-converted-space"> </span><br class="">> 1. please post the SIZE.h file.<br class="">> unless the behaviors induced by the paramters has changed, the info in thisfile will permit straightening out ofthese extents issues.<br class="">><span class="Apple-converted-space"> </span><br class="">> 2. below in parm04<br class="">> delY is set to 43*1 (why set it? i would not assume that there is no side effect from this, maybe there isn/t, i don/t know)<br class="">><span class="Apple-converted-space"> </span><br class="">> below in OBCS_PARMO1<br class="">><span class="Apple-converted-space"> </span><br class="">> > > OB_Jnorth(16:28)=13*41,<br class="">> > > > > OB_Jsouth(16:28)=13*3,<br class="">><span class="Apple-converted-space"> </span><br class="">> iirc, the values of these arrays are -1,0,1<br class="">> it/s been a while since i setup and ran a simulation<br class="">><span class="Apple-converted-space"> </span><br class="">><span class="Apple-converted-space"> </span><br class="">> Sent from Here3.<br class="">><span class="Apple-converted-space"> </span><br class="">> > On Aug 24, 2019, at 09:11, Estanislao Gavilan Pascual-Ahuir <<a href="mailto:dramauh@hotmail.com" class="">dramauh@hotmail.com</a>> wrote:<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > Hi Martin,<span class="Apple-converted-space"> </span><br class="">> ><span class="Apple-converted-space"> </span><br class="">> > I tried increasing the nx from 13 to 43 and I have the same error...This is a nightmare. Do you think I might need to make a bigger boundary?Maybe 13 is cells is too small.. Thank you so much in advance for your help.<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > Kind regards,<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > Estanislao<br class="">> > De: MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" class="">mitgcm-support-bounces@mitgcm.org</a>> en nombre de<span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support-request@mitgcm.org" class="">mitgcm-support-request@mitgcm.org</a><span class="Apple-converted-space"> </span><<a href="mailto:mitgcm-support-request@mitgcm.org" class="">mitgcm-support-request@mitgcm.org</a>><br class="">> > Enviado: s?bado, 24 de agosto de 2019 8:45<br class="">> > Para:<span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a><span class="Apple-converted-space"> </span><<a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a>><br class="">> > Asunto: MITgcm-support Digest, Vol 194, Issue 16<br class="">> > <span class="Apple-converted-space"> </span><br class="">> > Send MITgcm-support mailing list submissions to<br class="">> > <span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a><br class="">> ><span class="Apple-converted-space"> </span><br class="">> > To subscribe or unsubscribe via the World Wide Web, visit<br class="">> > <span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">> > or, via email, send a message with subject or body 'help' to<br class="">> > <span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support-request@mitgcm.org" class="">mitgcm-support-request@mitgcm.org</a><br class="">> ><span class="Apple-converted-space"> </span><br class="">> > You can reach the person managing the list at<br class="">> > <span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support-owner@mitgcm.org" class="">mitgcm-support-owner@mitgcm.org</a><br class="">> ><span class="Apple-converted-space"> </span><br class="">> > When replying, please edit your Subject line so it is more specific<br class="">> > than "Re: Contents of MITgcm-support digest..."<br class="">> ><span class="Apple-converted-space"> </span><br class="">> ><span class="Apple-converted-space"> </span><br class="">> > Today's Topics:<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > 1. Re: Reading errors (Martin Losch)<br class="">> ><span class="Apple-converted-space"> </span><br class="">> ><span class="Apple-converted-space"> </span><br class="">> > ----------------------------------------------------------------------<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > Message: 1<br class="">> > Date: Sat, 24 Aug 2019 08:45:29 +0200<br class="">> > From: Martin Losch <<a href="mailto:Martin.Losch@awi.de" class="">Martin.Losch@awi.de</a>><br class="">> > To: <<a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a>><br class="">> > Subject: Re: [MITgcm-support] Reading errors<br class="">> > Message-ID: <<a href="mailto:C4098FFF-029D-4EEC-AC58-566F448B660A@awi.de" class="">C4098FFF-029D-4EEC-AC58-566F448B660A@awi.de</a>><br class="">> > Content-Type: text/plain; charset="us-ascii"<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > Hi, I do think, that this may be your problem, because your fields are actually too small! You use nx=13, when it should be 43. Even if only part of your Northern or Southern boundary is really "open", you need to provide the full field. Please refer to the documentation for more details. Unfortunately, we don't have any verification experiments with exactly your configuration. Most of the time all boundaries are declared as open, even when only a few are actually "wet", e.g. the verification/dome experiment, where the open point along the northern bounary are only a few.<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > So, I think you can fix your problem by changing 13 to 43 in your matlab scripts.<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > Martin<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > > On 24. Aug 2019, at 03:07, Estanislao Gavilan Pascual-Ahuir <<a href="mailto:dramauh@hotmail.com" class="">dramauh@hotmail.com</a>> wrote:<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > Hi Martin<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > I am not sure if that is the problem because I used readbin.m to check the input files and the matrix were correct. This is the MATLAB code I used.<span class="Apple-converted-space"> </span><br class="">> > > % creat N time slabs for testing<br class="">> > > ieee='b';<br class="">> > > prec='real*4';<br class="">> > > v0=0.0041;<br class="">> > > N=100<br class="">> > > vvel=ones(13,2,N);<span class="Apple-converted-space"> </span><br class="">> > > uMerid=zeros(13,2,N);<span class="Apple-converted-space"> </span><br class="">> > > for i=1:N<br class="">> > > vZonal(:,:,i)=v0*vvel(:,:,i)*i/N;<br class="">> > > end<br class="">> > > fid=fopen('OBmeridU.bin','w',ieee); fwrite(fid,uMerid,prec);fclose(fid);<br class="">> > > fid=fopen('OBzonalV.bin','w',ieee); fwrite(fid,vZonal,prec); fclose(fid);<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > I have even tried to use your code, but the files do not change the size. They are always 11KB<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > % creat N time slabs for testing<br class="">> > > v0=0.0041;<br class="">> > > N=100<br class="">> > > vvel=ones(13,2,N);<span class="Apple-converted-space"> </span><br class="">> > > uMerid=zeros(13,2,N);<span class="Apple-converted-space"> </span><br class="">> > > for i=1:N<br class="">> > > vZonal(:,:,i)=v0*vvel(:,:,i)*i/N;<br class="">> > > end<br class="">> > > fid=fopen('OBmeridU.bin','w',ieee);<span class="Apple-converted-space"> </span><br class="">> > > for i=1:N<br class="">> > > fwrite(fid,squeeze(uMerid(:,:,i)),prec);<br class="">> > > end<br class="">> > > fclose(fid);<br class="">> > > fid=fopen('OBzonalV.bin','w',ieee);<span class="Apple-converted-space"> </span><br class="">> > > for i=1:N<br class="">> > > fwrite(fid,squeeze(vZonal(:,:,i)),prec);<br class="">> > > end<br class="">> > > fclose(fid);<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > Kind regards,<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > Estanislao<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > De: MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" class="">mitgcm-support-bounces@mitgcm.org</a>> en nombre de<span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support-request@mitgcm.org" class="">mitgcm-support-request@mitgcm.org</a><span class="Apple-converted-space"> </span><<a href="mailto:mitgcm-support-request@mitgcm.org" class="">mitgcm-support-request@mitgcm.org</a>><br class="">> > > Enviado: viernes, 23 de agosto de 2019 18:00<br class="">> > > Para:<span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a><span class="Apple-converted-space"> </span><<a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a>><br class="">> > > Asunto: MITgcm-support Digest, Vol 194, Issue 14<br class="">> > > <span class="Apple-converted-space"> </span><br class="">> > > Send MITgcm-support mailing list submissions to<br class="">> > > <span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > To subscribe or unsubscribe via the World Wide Web, visit<br class="">> > > <span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">> > > or, via email, send a message with subject or body 'help' to<br class="">> > > <span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support-request@mitgcm.org" class="">mitgcm-support-request@mitgcm.org</a><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > You can reach the person managing the list at<br class="">> > > <span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support-owner@mitgcm.org" class="">mitgcm-support-owner@mitgcm.org</a><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > When replying, please edit your Subject line so it is more specific<br class="">> > > than "Re: Contents of MITgcm-support digest..."<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > Today's Topics:<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > 1. Re: Reading errors (Martin Losch)<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > ----------------------------------------------------------------------<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > Message: 1<br class="">> > > Date: Fri, 23 Aug 2019 10:02:49 +0200<br class="">> > > From: Martin Losch <<a href="mailto:Martin.Losch@awi.de" class="">Martin.Losch@awi.de</a>><br class="">> > > To: MITgcm Support <<a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a>><br class="">> > > Subject: Re: [MITgcm-support] Reading errors<br class="">> > > Message-ID: <<a href="mailto:1AB158B6-A07E-43EB-A3F6-6530BDE345C9@awi.de" class="">1AB158B6-A07E-43EB-A3F6-6530BDE345C9@awi.de</a>><br class="">> > > Content-Type: text/plain; charset="utf-8"<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > Hi Estanis,<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > this most likely means that your OBzonalV.bin is ?too short?.<br class="">> > > Guess from you parameters you'll need fields with dimension (nt,nz,nx) = (100,2,43). In you data file you do not specify the ?readBinaryPrec?, so that will default to ?real*4?, single precision, ?float32?, but check in your stdout, so in python this would mean writing a file like this:<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > import numpy as np<br class="">> > > import sys<br class="">> > > nt=100<br class="">> > > nx=43<br class="">> > > nz=2<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > data = np.zeros((nt,nz,nx),dtype=?float32') # or some other values, but dtype is important<br class="">> > > if sys.byteorder == 'little': data.byteswap(True)<br class="">> > > #<br class="">> > > fid = open(?OBzonalV.bin',"wb")<br class="">> > > data.tofile(fid)<br class="">> > > fid.close()<br class="">> > > # switch back to machine format<br class="">> > > if sys.byteorder == 'little': data.byteswap(True)<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > With this I get a file of this size:<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > >>> du -h ./OBzonalV.bin<br class="">> > > 36K ./OBzonalV.bin<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > >>> ls -l OBzonalV.bin<br class="">> > > -rw-r--r-- 1 mlosch DMAWI\KLI04 34400 Aug 23 09:53 OBzonalV.bin<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > # in matlab it would be similar, but note the different order of dimensions:<br class="">> > > data = zeros(nx,nz,nt);<br class="">> > > accuracy = ?real*4?;<br class="">> > > [fid message] = fopen(?OBzonalV.bin','w','ieee-be');<br class="">> > > count = fwrite(fid,data,accuracy);<br class="">> > > for k=1:nt<br class="">> > > fwrite(fid,squeeze(data(:,:,k)),accuracy);<br class="">> > > end<br class="">> > > fclose(fid);<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > About the NETCDF_ROOT and other INCLUDES: this file is meant to be as generic as possible, i.e. the code at the end of the build options file will look for paths in the usual places. In 95% of the cases this works. If you have a special installation, you might want to hard-wire the paths, as you have done, then you don?t need the rest anymore. Alternatively you could set the environment varialbes NETCDF_ROOT and MPI_INC_DIR, e.g. like this<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > export NETCDF_ROOT=/WORK/app/netcdf/4.3.2/01-CF-14/<br class="">> > > export MPI_INC_DIR=/usr/local/mpi3-dynamic/include<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > ***before*** you run genmake2, then it should work, too.<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > Martin<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > > On 23. Aug 2019, at 04:37, Estanislao Gavilan Pascual-Ahuir <<a href="mailto:dramauh@hotmail.com" class="">dramauh@hotmail.com</a>> wrote:<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Hi Martin and Jean-Michel<br class="">> > > > I still have the error forrtl: severe (36): attempt to access non-existent record, unit 16, file /bblablabla/run/OBzonalV.bin. In addition, I saw this message after doing make depend f90mkdepend: no source file found for module this. I am no sure if this is important.<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > I followed your advice. I start doing a simulation in serial mode without obcs and the model ran perfectly fine. Once I opened the obcs (i.e useOBCS= TRUE), I had that error. Following the architecture of my cluster change a little bit the linux_ia64_ifort. By the way, do I need the last set of conditions of NETCDF_ROOT and netcdf test? I was thinking to remove them.<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Kind regards,<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Estanis<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > #!/bin/bash<br class="">> > > > #<br class="">> > > > # Tested on<span class="Apple-converted-space"> </span><a href="http://uv100.awi.de/" class="">uv100.awi.de</a><span class="Apple-converted-space"> </span>(SGI UV 100, details:<br class="">> > > > # <span class="Apple-converted-space"> </span><a href="http://www.sgi.com/products/servers/uv/specs.html" class="">http://www.sgi.com/products/servers/uv/specs.html</a>)<br class="">> > > > # a) For more speed, provided your data size does not exceed 2GB you can<br class="">> > > > # remove -fPIC which carries a performance penalty of 2-6%.<br class="">> > > > # b) You can replace -fPIC with '-mcmodel=medium -shared-intel' which may<br class="">> > > > # perform faster than -fPIC and still support data sizes over 2GB per<br class="">> > > > # process but all the libraries you link to must be compiled with<br class="">> > > > # -fPIC or -mcmodel=medium<br class="">> > > > # c) flags adjusted for ifort 12.1.0<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > FC=ifort<br class="">> > > > F90C=ifort<br class="">> > > > CC=icc<br class="">> > > > # requires that all static libraries are available:<br class="">> > > > #LINK='ifort -static'<br class="">> > > > LINK='ifort'<br class="">> > > > # for adjoint runs the default makedepend often cannot handle enough files<br class="">> > > > #MAKEDEPEND=tools_xmakedepend<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > DEFINES='-DWORDLENGTH=4'<br class="">> > > > CPP='cpp -traditional -P'<br class="">> > > > F90FIXEDFORMAT='-fixed -Tf'<br class="">> > > > EXTENDED_SRC_FLAG='-132'<br class="">> > > > #OMPFLAG='-openmp'<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > NOOPTFLAGS="-O0 -g -m64"<br class="">> > > > NOOPTFILES=''<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > MCMODEL='-fPIC'<br class="">> > > > # for large memory requirements uncomment this line<br class="">> > > > #MCMODEL='-mcmodel=medium -shared-intel'<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > FFLAGS="$FFLAGS -W0 -WB -convert big_endian -assume byterecl $MCMODEL"<br class="">> > > > #- might want to use '-r8' for fizhi pkg:<br class="">> > > > #FFLAGS="$FFLAGS -r8"<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > if test "x$IEEE" = x ; then #- with optimisation:<br class="">> > > > FOPTIM='-O3 -align'<br class="">> > > > # does not work when -static does not work<br class="">> > > > # FOPTIM='-fast -align'<br class="">> > > > # instead you can use<br class="">> > > > # FOPTIM='-O3 -ipo -align'<br class="">> > > > else<br class="">> > > > if test "x$DEVEL" = x ; then #- no optimisation + IEEE :<br class="">> > > > FOPTIM='-O0 -noalign -fp-model precise'<br class="">> > > > # -fltconsistency<br class="">> > > > else #- development/check options:<br class="">> > > > FOPTIM='-O0 -noalign -fp-model precise'<br class="">> > > > FOPTIM="$FOPTIM -g -check all -fpe0 -traceback -ftrapuv -fp-model except -warn all"<br class="">> > > > fi<br class="">> > > > fi<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > F90FLAGS=$FFLAGS<br class="">> > > > F90OPTIM=$FOPTIM<br class="">> > > > CFLAGS="-O0 -ip $MCMODEL"<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > INCLUDEDIRS=''<br class="">> > > > INCLUDES='-I. -I$NETCDF/include -I/WORK/app/netcdf/4.3.2/01-CF-14/include -I/usr/local/mpi3-dynamic/include'<br class="">> > > > LIBS='-L$NETCDF/lib -lnetcdff -lnetcdf -I/WORK/app/netcdf/4.3.2/01-CF-14/lib'<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > if [ "x$NETCDF_ROOT" != x ] ; then<br class="">> > > > INCLUDEDIRS="${NETCDF_ROOT}/include"<br class="">> > > > INCLUDES="-I${NETCDF_ROOT}/include"<br class="">> > > > LIBS="-L${NETCDF_ROOT}/lib"<br class="">> > > > elif [ "x$NETCDF_HOME" != x ]; then<br class="">> > > > INCLUDEDIRS="${NETCDF_HOME}/include"<br class="">> > > > INCLUDES="-I${NETCDF_HOME}/include"<br class="">> > > > LIBS="-L${NETCDF_HOME}/lib"<br class="">> > > > elif [ "x$NETCDF_INC" != x -a"x$NETCDF_LIB" != x ]; then<br class="">> > > > NETCDF_INC=`echo $NETCDF_INC | sed 's/-I//g'`<br class="">> > > > NETCDF_LIB=`echo $NETCDF_LIB | sed 's/-L//g'`<br class="">> > > > INCLUDEDIRS="${NETCDF_INC}"<br class="">> > > > INCLUDES="-I${NETCDF_INC}"<br class="">> > > > LIBS="-L${NETCDF_LIB}"<br class="">> > > > elif [ "x$NETCDF_INCDIR" != x -a "x$NETCDF_LIBDIR" != x ]; then<br class="">> > > > INCLUDEDIRS="${NETCDF_INCDIR}"<br class="">> > > > INCLUDES="-I${NETCDF_INCDIR}"<br class="">> > > > LIBS="-L${NETCDF_LIBDIR}"<br class="">> > > > elif test -d /usr/include/netcdf-3 ; then<br class="">> > > > INCLUDEDIRS='/usr/include/netcdf-3'<br class="">> > > > INCLUDES='-I/usr/include/netcdf-3'<br class="">> > > > LIBS='-L/usr/lib/netcdf-3 -L/usr/lib64/netcdf-3'<br class="">> > > > elif test -d /usr/include/netcdf ; then<br class="">> > > > INCLUDEDIRS='/usr/include/netcdf'<br class="">> > > > INCLUDES='-I/usr/include/netcdf'<br class="">> > > > elif test -d /usr/local/netcdf ; then<br class="">> > > > INCLUDEDIRS='/usr/include/netcdf/include'<br class="">> > > > INCLUDES='-I/usr/local/netcdf/include'<br class="">> > > > LIBS='-L/usr/local/netcdf/lib'<br class="">> > > > elif test -d /usr/local/include/netcdf.inc ; then<br class="">> > > > INCLUDEDIRS='/usr/local/include'<br class="">> > > > INCLUDES='-I/usr/local/include'<br class="">> > > > LIBS='-L/usr/local/lib64'<br class="">> > > > fi<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > if [ -n "$MPI_INC_DIR" -a "x$MPI" = xtrue ] ; then<br class="">> > > > LIBS="$LIBS -lmpi"<br class="">> > > > INCLUDES="$INCLUDES -I$MPI_INC_DIR"<br class="">> > > > INCLUDEDIRS="$INCLUDEDIRS $MPI_INC_DIR"<br class="">> > > > #- used for parallel (MPI) DIVA<br class="">> > > > MPIINCLUDEDIR="$MPI_INC_DIR"<br class="">> > > > fi<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > ----------------------------------------------------------------------<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Message: 1<br class="">> > > > Date: Thu, 22 Aug 2019 13:10:55 +0200<br class="">> > > > From: Martin Losch <<a href="mailto:Martin.Losch@awi.de" class="">Martin.Losch@awi.de</a>><br class="">> > > > To: MITgcm Support <<a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a>><br class="">> > > > Subject: Re: [MITgcm-support] Reading errors<br class="">> > > > Message-ID: <<a href="mailto:325AC4E8-1648-4C6A-BFED-7722921C733E@awi.de" class="">325AC4E8-1648-4C6A-BFED-7722921C733E@awi.de</a>><br class="">> > > > Content-Type: text/plain; charset="utf-8"<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Hi Estanis,<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > thanks for the details. This is what I would do:<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > - At the compile level use a standard build options file, with an intel compiler on a linux system I would start with MITgcm/tools/build_options/linux_amd64_ifort, or linux_ia64_ifort (depending on the output of uname -a, in fact, genmake2 is probably able to pick the correct file if you don?t specify it), and since your domain is small I would first try without MPI, ie. like this:<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > ${somepath}/tools/genmake2 -of ${somepath}/tools/build_options/linux_amd64_ifort -mods ../code<br class="">> > > > make CLEAN && make depend && make<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > - With this non-MPI configuration I would try to run the model. First with useOBCS=.FALSE. (just a few timesteps), and then with .TRUE.<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > - once this works, you can recompile with MPI (if you really need it), like this:<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > ${somepath}/tools/genmake2 -of ${somepath}/tools/build_options/linux_amd64_ifort -mods ../code -mpi<br class="">> > > > make CLEAN && make depend && make<br class="">> > > > (note that the extra flag ?-mpi" is enough)<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > and check if you get the same. For further help, you should record the potential error messages after each step.<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Martin<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > PS. Some comments about your namelist below:<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > > On 22. Aug 2019, at 12:39, Estanislao Gavilan Pascual-Ahuir <<a href="mailto:dramauh@hotmail.com" class="">dramauh@hotmail.com</a>> wrote:<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > Hi Martin ,<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > Before anything thank you so much for your help. I will try to answer all you questions.<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > what is the platform, the compiler?<br class="">> > > > > The platform in linux 2.6.32-431.TH.x86_64 GNU/Linux. Red Hat Enterprise Linux Server release 6.5. I am using intel compilers wrapped them in mpi. The version of the compiler is 14.0.2<br class="">> > > > > details of the configuration (content of code-directory and namelist files)<br class="">> > > > > I am running a simple simulation with open boundaries. I load the packages gfd, obcs, mnc and diagnostics using the packagesconfig. The frequency of the open boundaries is stated in the data file. This is the data file<br class="">> > > > > Model parameters<br class="">> > > > > # Continuous equation parameters<br class="">> > > > > &PARM01<br class="">> > > > > tRef=23.,23.,<br class="">> > > > > sRef=35.,35.,<br class="">> > > > > selectCoriMap=4,<br class="">> > > > > viscAh=4.E2,<br class="">> > > > with your grid choice (sphericalPolarGrid), the coriolis parameter is computed and these values are not used.<br class="">> > > > > f0=1.E-4,<br class="">> > > > > beta=1.E-11,<br class="">> > > > > rhoNil=1000.,<br class="">> > > > > gBaro=9.81,<br class="">> > > > > rigidLid=.FALSE.,<br class="">> > > > > implicitFreeSurface=.TRUE.,<br class="">> > > > > # momAdvection=.FALSE.,<br class="">> > > > > tempStepping=.FALSE.,<br class="">> > > > > saltStepping=.FALSE.,<br class="">> > > > > &<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > # Elliptic solver parameters<br class="">> > > > > &PARM02<br class="">> > > > > cg2dTargetResidual=1.E-7,<br class="">> > > > > cg2dMaxIters=1000,<br class="">> > > > > &<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > # Time stepping parameters<br class="">> > > > > &PARM03<br class="">> > > > > nIter0=0,<br class="">> > > > > nTimeSteps=100,<br class="">> > > > > deltaT=1200.0,<br class="">> > > > > pChkptFreq=31104000.0,<br class="">> > > > > chkptFreq=15552000.0,<br class="">> > > > > dumpFreq=15552000.0,<br class="">> > > > # this will give you monitor output every timestep (which is what you want while debugging), later I would set it to something like 20-50 * deltaT<br class="">> > > > > monitorFreq=1200.,<br class="">> > > > > monitorSelect=2,<br class="">> > > > > periodicExternalForcing=.TRUE.,<br class="">> > > > # this means that you will read data each time step. Is that what you want?<br class="">> > > > > externForcingPeriod= 1200.,<br class="">> > > > # with your choice of externForcingPeriod, this requires that you have 1000. records in the file(s)<br class="">> > > > > externForcingCycle = 12000000.,<br class="">> > > > > &<br class="">> > > > > # Gridding parameters<br class="">> > > > > &PARM04<br class="">> > > > > usingSphericalPolarGrid=.TRUE.,<br class="">> > > > # alternatively you can say dxSpacing = 1., dySpacing = 1.,<br class="">> > > > > delX=43*1.,<br class="">> > > > > delY=43*1.,<br class="">> > > > > xgOrigin=-21.,<br class="">> > > > > ygOrigin=-21.,<br class="">> > > > > delR=2*500.,<br class="">> > > > > &<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > # Input datasets<br class="">> > > > > &PARM05<br class="">> > > > > bathyFile='bathy_cir.bin'<br class="">> > > > > meridWindFile=,<br class="">> > > > > &<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > This is the data.obcs<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > # Open-boundaries<br class="">> > > > > &OBCS_PARM01<br class="">> > > > > OBCSfixTopo=.FALSE.,<br class="">> > > > # if I understand the configuration correctly, you have a zonally reentrant channel with wall in the north and the south (python notation: bathy[0+2,:] = 0, and bathy[ny-1,:] = 0, except where you have the open boundaries)? you could actually save two grid rows (have 40 instead of 43 point in j-direction and set bathy[0,:]=0, bathy[ny,:]=0)<br class="">> > > > > OB_Ieast=0,<br class="">> > > > > OB_Iwest=0,<br class="">> > > > > OB_Jnorth(16:28)=13*41,<br class="">> > > > > OB_Jsouth(16:28)=13*3,<br class="">> > > > > useOBCSprescribe = .TRUE.,<br class="">> > > > These files should be found, if they are in the same directory where you run your model. They should each contain (according to you dimensions and time parameters) for 100 timesteps 100 fields of dimension (nx,nz). For anything above 1000 timesteps, they should have 1000 fields (because after the 1000ths record, the model starts from the beginning again, according to you externForcingCycle)<br class="">> > > > > OBNvFile = 'OBzonalV.bin',<br class="">> > > > > OBSvFile = 'OBzonalV.bin',<br class="">> > > > > OBNuFile = 'OBmeridU.bin',<br class="">> > > > > OBSuFile = 'OBmeridU.bin?,<br class="">> > > > # same as before this will give you a lot of output. You may want to comment out this line, becasue OBCS_monitorFreq = monitorFreq by default<br class="">> > > > > OBCS_monitorFreq=1200.00,<br class="">> > > > > OBCS_monSelect = 1,<br class="">> > > > > &<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > &OBCS_PARM02<br class="">> > > > > &<br class="">> > > > > are you using latest code (some of the flags in the build-option look very outdated ?)?<br class="">> > > > > Yes, it is the latest code (version MITgcm_c67k). About the flags in my build version, I did not make my own one. I used one that I found in our research group.<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > Kind regards,<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > Estanislao<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > ------------------------------<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Message: 2<br class="">> > > > Date: Thu, 22 Aug 2019 09:33:18 -0400<br class="">> > > > From: Jean-Michel Campin <<a href="mailto:jmc@mit.edu" class="">jmc@mit.edu</a>><br class="">> > > > To:<span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" class="">mitgcm-support@mitgcm.org</a><br class="">> > > > Subject: Re: [MITgcm-support] Reading errors<br class="">> > > > Message-ID: <<a href="mailto:20190822133318.GA13562@ocean.mit.edu" class="">20190822133318.GA13562@ocean.mit.edu</a>><br class="">> > > > Content-Type: text/plain; charset=us-ascii<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Hi Estanis,<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Just a small adjustment:<span class="Apple-converted-space"> </span><br class="">> > > > the standard optfile for intel compiler (version 11 and newer) is:<br class="">> > > > linux_amd64_ifort11<br class="">> > > > in MITgcm/tools/build_options<br class="">> > > > The optfile "linux_amd64_ifort" is for older version (10 and older).<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > However, if you are compiling with intel MPI (recent version of the compiler),<br class="">> > > > then you need to use: linux_amd64_ifort+impi<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Cheers,<br class="">> > > > Jean-Michel<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > On Thu, Aug 22, 2019 at 01:10:55PM +0200, Martin Losch wrote:<br class="">> > > > > Hi Estanis,<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > thanks for the details. This is what I would do:<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > - At the compile level use a standard build options file, with an intel compiler on a linux system I would start with MITgcm/tools/build_options/linux_amd64_ifort, or linux_ia64_ifort (depending on the output of uname -a, in fact, genmake2 is probably able to pick the correct file if you don???t specify it), and since your domain is small I would first try without MPI, ie. like this:<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > ${somepath}/tools/genmake2 -of ${somepath}/tools/build_options/linux_amd64_ifort -mods ../code<br class="">> > > > > make CLEAN && make depend && make<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > - With this non-MPI configuration I would try to run the model. First with useOBCS=.FALSE. (just a few timesteps), and then with .TRUE.<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > - once this works, you can recompile with MPI (if you really need it), like this:<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > ${somepath}/tools/genmake2 -of ${somepath}/tools/build_options/linux_amd64_ifort -mods ../code -mpi<br class="">> > > > > make CLEAN && make depend && make<br class="">> > > > > (note that the extra flag ???-mpi" is enough)<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > and check if you get the same. For further help, you should record the potential error messages after each step.<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > Martin<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > PS. Some comments about your namelist below:<br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > On 22. Aug 2019, at 12:39, Estanislao Gavilan Pascual-Ahuir <<a href="mailto:dramauh@hotmail.com" class="">dramauh@hotmail.com</a>> wrote:<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > Hi Martin ,<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > Before anything thank you so much for your help. I will try to answer all you questions.<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > what is the platform, the compiler?<br class="">> > > > > > The platform in linux 2.6.32-431.TH.x86_64 GNU/Linux. Red Hat Enterprise Linux Server release 6.5. I am using intel compilers wrapped them in mpi. The version of the compiler is 14.0.2<br class="">> > > > > > details of the configuration (content of code-directory and namelist files)<br class="">> > > > > > I am running a simple simulation with open boundaries. I load the packages gfd, obcs, mnc and diagnostics using the packagesconfig. The frequency of the open boundaries is stated in the data file. This is the data file<br class="">> > > > > > Model parameters<br class="">> > > > > > # Continuous equation parameters<br class="">> > > > > > &PARM01<br class="">> > > > > > tRef=23.,23.,<br class="">> > > > > > sRef=35.,35.,<br class="">> > > > > > selectCoriMap=4,<br class="">> > > > > > viscAh=4.E2,<br class="">> > > > > with your grid choice (sphericalPolarGrid), the coriolis parameter is computed and these values are not used.<br class="">> > > > > > f0=1.E-4,<br class="">> > > > > > beta=1.E-11,<br class="">> > > > > > rhoNil=1000.,<br class="">> > > > > > gBaro=9.81,<br class="">> > > > > > rigidLid=.FALSE.,<br class="">> > > > > > implicitFreeSurface=.TRUE.,<br class="">> > > > > > # momAdvection=.FALSE.,<br class="">> > > > > > tempStepping=.FALSE.,<br class="">> > > > > > saltStepping=.FALSE.,<br class="">> > > > > > &<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > # Elliptic solver parameters<br class="">> > > > > > &PARM02<br class="">> > > > > > cg2dTargetResidual=1.E-7,<br class="">> > > > > > cg2dMaxIters=1000,<br class="">> > > > > > &<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > # Time stepping parameters<br class="">> > > > > > &PARM03<br class="">> > > > > > nIter0=0,<br class="">> > > > > > nTimeSteps=100,<br class="">> > > > > > deltaT=1200.0,<br class="">> > > > > > pChkptFreq=31104000.0,<br class="">> > > > > > chkptFreq=15552000.0,<br class="">> > > > > > dumpFreq=15552000.0,<br class="">> > > > > # this will give you monitor output every timestep (which is what you want while debugging), later I would set it to something like 20-50 * deltaT<br class="">> > > > > > monitorFreq=1200.,<br class="">> > > > > > monitorSelect=2,<br class="">> > > > > > periodicExternalForcing=.TRUE.,<br class="">> > > > > # this means that you will read data each time step. Is that what you want?<br class="">> > > > > > externForcingPeriod= 1200.,<br class="">> > > > > # with your choice of externForcingPeriod, this requires that you have 1000. records in the file(s)<br class="">> > > > > > externForcingCycle = 12000000.,<br class="">> > > > > > &<br class="">> > > > > > # Gridding parameters<br class="">> > > > > > &PARM04<br class="">> > > > > > usingSphericalPolarGrid=.TRUE.,<br class="">> > > > > # alternatively you can say dxSpacing = 1., dySpacing = 1.,<br class="">> > > > > > delX=43*1.,<br class="">> > > > > > delY=43*1.,<br class="">> > > > > > xgOrigin=-21.,<br class="">> > > > > > ygOrigin=-21.,<br class="">> > > > > > delR=2*500.,<br class="">> > > > > > &<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > # Input datasets<br class="">> > > > > > &PARM05<br class="">> > > > > > bathyFile='bathy_cir.bin'<br class="">> > > > > > meridWindFile=,<br class="">> > > > > > &<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > This is the data.obcs<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > # Open-boundaries<br class="">> > > > > > &OBCS_PARM01<br class="">> > > > > > OBCSfixTopo=.FALSE.,<br class="">> > > > > # if I understand the configuration correctly, you have a zonally reentrant channel with wall in the north and the south (python notation: bathy[0+2,:] = 0, and bathy[ny-1,:] = 0, except where you have the open boundaries)? you could actually save two grid rows (have 40 instead of 43 point in j-direction and set bathy[0,:]=0, bathy[ny,:]=0)<br class="">> > > > > > OB_Ieast=0,<br class="">> > > > > > OB_Iwest=0,<br class="">> > > > > > OB_Jnorth(16:28)=13*41,<br class="">> > > > > > OB_Jsouth(16:28)=13*3,<br class="">> > > > > > useOBCSprescribe = .TRUE.,<br class="">> > > > > These files should be found, if they are in the same directory where you run your model. They should each contain (according to you dimensions and time parameters) for 100 timesteps 100 fields of dimension (nx,nz). For anything above 1000 timesteps, they should have 1000 fields (because after the 1000ths record, the model starts from the beginning again, according to you externForcingCycle)<br class="">> > > > > > OBNvFile = 'OBzonalV.bin',<br class="">> > > > > > OBSvFile = 'OBzonalV.bin',<br class="">> > > > > > OBNuFile = 'OBmeridU.bin',<br class="">> > > > > > OBSuFile = 'OBmeridU.bin???,<br class="">> > > > > # same as before this will give you a lot of output. You may want to comment out this line, becasue OBCS_monitorFreq = monitorFreq by default<br class="">> > > > > > OBCS_monitorFreq=1200.00,<br class="">> > > > > > OBCS_monSelect = 1,<br class="">> > > > > > &<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > &OBCS_PARM02<br class="">> > > > > > &<br class="">> > > > > > are you using latest code (some of the flags in the build-option look very outdated ?)?<br class="">> > > > > > Yes, it is the latest code (version MITgcm_c67k). About the flags in my build version, I did not make my own one. I used one that I found in our research group.<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > Kind regards,<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > > Estanislao<br class="">> > > > > ><span class="Apple-converted-space"> </span><br class="">> > > > ><span class="Apple-converted-space"> </span><br class="">> > > > > _______________________________________________<br class="">> > > > > MITgcm-support mailing list<br class="">> > > > ><span class="Apple-converted-space"> </span><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">> > > > ><span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > ------------------------------<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > Subject: Digest Footer<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > _______________________________________________<br class="">> > > > MITgcm-support mailing list<br class="">> > > ><span class="Apple-converted-space"> </span><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">> > > ><span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > ------------------------------<br class="">> > > ><span class="Apple-converted-space"> </span><br class="">> > > > End of MITgcm-support Digest, Vol 194, Issue 12<br class="">> > > > ***********************************************<br class="">> > > > _______________________________________________<br class="">> > > > MITgcm-support mailing list<br class="">> > > ><span class="Apple-converted-space"> </span><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">> > > ><span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > ------------------------------<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > Subject: Digest Footer<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > _______________________________________________<br class="">> > > MITgcm-support mailing list<br class="">> > ><span class="Apple-converted-space"> </span><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">> > ><span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > ------------------------------<br class="">> > ><span class="Apple-converted-space"> </span><br class="">> > > End of MITgcm-support Digest, Vol 194, Issue 14<br class="">> > > ***********************************************<br class="">> > > _______________________________________________<br class="">> > > MITgcm-support mailing list<br class="">> > ><span class="Apple-converted-space"> </span><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">> > ><span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">> > -------------- next part --------------<br class="">> > An HTML attachment was scrubbed...<br class="">> > URL: <<a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190824/b7756f1c/attachment.html" class="">http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190824/b7756f1c/attachment.html</a>><br class="">> ><span class="Apple-converted-space"> </span><br class="">> > ------------------------------<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > Subject: Digest Footer<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > _______________________________________________<br class="">> > MITgcm-support mailing list<br class="">> ><span class="Apple-converted-space"> </span><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">> ><span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">> ><span class="Apple-converted-space"> </span><br class="">> ><span class="Apple-converted-space"> </span><br class="">> > ------------------------------<br class="">> ><span class="Apple-converted-space"> </span><br class="">> > End of MITgcm-support Digest, Vol 194, Issue 16<br class="">> > ***********************************************<br class="">> > _______________________________________________<br class="">> > MITgcm-support mailing list<br class="">> ><span class="Apple-converted-space"> </span><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">> ><span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">> -------------- next part --------------<br class="">> An HTML attachment was scrubbed...<br class="">> URL: <<a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190824/e57f7f6a/attachment.html" class="">http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190824/e57f7f6a/attachment.html</a>><br class="">><span class="Apple-converted-space"> </span><br class="">> ------------------------------<br class="">><span class="Apple-converted-space"> </span><br class="">> Subject: Digest Footer<br class="">><span class="Apple-converted-space"> </span><br class="">> _______________________________________________<br class="">> MITgcm-support mailing list<br class="">><span class="Apple-converted-space"> </span><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">><span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">><span class="Apple-converted-space"> </span><br class="">><span class="Apple-converted-space"> </span><br class="">> ------------------------------<br class="">><span class="Apple-converted-space"> </span><br class="">> End of MITgcm-support Digest, Vol 194, Issue 19<br class="">> ***********************************************<br class="">><span class="Apple-converted-space"> </span><br class="">> Virus-free.<span class="Apple-converted-space"> </span><a href="http://www.avast.com/" class="">www.avast.com</a><br class="">> _______________________________________________<br class="">> MITgcm-support mailing list<br class="">><span class="Apple-converted-space"> </span><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">><span class="Apple-converted-space"> </span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">-------------- next part --------------<br class="">An HTML attachment was scrubbed...<br class="">URL: <<a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190825/4b07aa90/attachment.html" class="">http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190825/4b07aa90/attachment.html</a>><br class=""><br class="">------------------------------<br class=""><br class="">Subject: Digest Footer<br class=""><br class="">_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class=""><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class=""><br class=""><br class="">------------------------------<br class=""><br class="">End of MITgcm-support Digest, Vol 194, Issue 21<br class="">***********************************************<br class=""></div></span></font></div><div id="DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" style="caret-color: rgb(0, 0, 0); font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><br class=""><table style="border-top-width: 1px; border-top-style: solid; border-top-color: rgb(211, 212, 222);" class=""><tbody class=""><tr class=""><td style="width: 55px; padding-top: 13px;" class=""><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank" class=""><img alt="" width="46" height="29" src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" style="width: 46px; height: 29px;" class=""></a></td><td style="width: 470px; padding-top: 12px; color: rgb(65, 66, 78); font-size: 13px; font-family: Arial, Helvetica, sans-serif; line-height: 18px;" class="">Virus-free.<span class="Apple-converted-space"> </span><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank" style="color: rgb(68, 83, 234);" class="">www.avast.com</a></td></tr></tbody></table><a href="x-msg://5/#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1" class=""></a></div><span style="caret-color: rgb(0, 0, 0); font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><span style="caret-color: rgb(0, 0, 0); font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">MITgcm-support mailing list</span><br style="caret-color: rgb(0, 0, 0); font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><a href="mailto:MITgcm-support@mitgcm.org" style="font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">MITgcm-support@mitgcm.org</a><br style="caret-color: rgb(0, 0, 0); font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" style="font-family: LucidaSans-Typewriter; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a></div></blockquote></div><br class=""></div></body></html>