<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">I don’t know for sure, but I suspect if you set<div class="">externForcingPeriod = 0.</div><div class="">And give the model two time records (ny*nz*2) to be safe all will be fine</div><div class=""><br class=""></div><div class="">Matt</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Jun 26, 2019, at 8:37 AM, Kevin Ha <<a href="mailto:kevin.ha.pro@gmail.com" class="">kevin.ha.pro@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class="">ah now it's working perfectly (the mistake was that East OBCS was not activated on the OBCS_OPTIONS.h). My last question would be, regarding the answer of Martin : do we need to necesarily set an externForcingPeriod ? In fact my boundary condition at x=0 is a steady flow (U = cst) so is it possible to define a ny*nz file instead of a ny*nz*mytime file to specify my OBCS on the left side of my domain? Below my data file and gendata.m . Thank you in advance<br class=""><br class=""></div>My current data file :<br class=""># ====================<br class=""># | Model parameters |<br class=""># ====================<br class="">#<br class=""># Continuous equation parameters<br class=""> &PARM01<br class=""><br class=""> tRef = 350*0., <br class=""> sRef = 350*35.,<br class=""><br class=""> no_slip_sides = .TRUE.,<br class=""> no_slip_bottom = .TRUE.,<br class=""><br class=""> viscAh = 2.E-3,<br class=""> viscAz = 2.E-3,<br class=""><br class=""> diffKhT = 1.E-4,<br class=""> diffKzT = 1.E-4,<br class=""><br class=""> f0 = 0.53E-4,<br class=""> beta = 0.E-11,<br class=""><br class=""> tAlpha = 2.E-4,<br class=""> sBeta = 0.E-4,<br class=""><br class=""> gravity = 9.81,<br class=""> gBaro = 9.81,<br class=""><br class=""> rigidLid = .FALSE.,<br class=""> implicitFreeSurface=.TRUE.,<br class=""> exactConserv = .TRUE.<br class=""><br class=""> nonHydrostatic = .TRUE.,<br class=""> useSingleCpuIO = .TRUE.,<br class=""><br class=""> readBinaryPrec = 64,<br class=""> writeBinaryPrec = 64,<br class=""> writeStatePrec = 64,<br class=""><br class=""> staggerTimeStep=.TRUE.<br class=""><br class=""> saltStepping = .FALSE.,<br class=""><br class=""> tempAdvScheme = 7, !7th Order One Step method with Monotonicity Preserving Limiter ???? GAD.h<br class=""><br class=""> hFacMin=0.1,<br class=""><br class=""> &<br class=""><br class=""># Elliptic solver parameters<br class=""> &PARM02<br class=""> cg2dMaxIters = 10000,<br class=""> cg2dTargetResidual = 1.E-14,<br class=""><br class=""> cg3dMaxIters = 400,<br class=""> cg3dTargetResidual = 1.E-14,<br class=""> &<br class=""><br class=""># Time stepping parameters<br class=""> &PARM03<br class=""><br class=""> niter0 = 0,<br class=""><br class=""> nTimeSteps = 28800, ! 10j avec deltaT=30<br class=""> deltaT = 30.,<br class=""> dumpFreq = 1800.,<br class=""> monitorFreq = 43200.,<br class=""><br class=""> abEps = 0.1,<br class=""> pChkptFreq = 86400.,<br class=""> chkptFreq = 0.0,<br class=""><br class=""> &<br class=""><br class=""># Gridding parameters<br class=""> &PARM04<br class=""> usingCartesianGrid=.TRUE.,<br class=""> usingSphericalPolarGrid=.FALSE.,<br class=""><br class=""> delY = 1*30., <br class=""> delRFile = 'delZ.init',<br class=""> delX = 2000*30., ! 2000 pts * 30m =60km = L<br class=""> &<br class=""><br class=""># Input datasets<br class=""> &PARM05<br class=""><br class=""> uVelInitFile = 'U.init',<br class=""> vVelInitFile = 'V.init',<br class=""> hydrogThetaFile = 'T.init',<br class=""> bathyFile = 'topog.init',<br class=""> &<br class=""><br class=""></div><div class="">My current gendata.m :</div><div class=""><br class=""></div><div class=""><br class="">clear<br class=""><br class="">%--------<br class="">nx = 2000;<br class="">ny = 1;<br class="">nz = 350;<br class=""><br class="">dx = 30;<br class="">dz = 10;<br class=""><br class="">dz = dz * ones(nz,1);<br class="">dz(251:350) = dz(250) * 1.02.^(1:100);<br class="">%--------<br class=""><br class="">%--------<br class="">Lx = nx*dx;<br class="">Hz = sum(dz);<br class=""><br class="">x = (dx/2:dx:Lx-dx/2);<br class="">z = - Hz + cumsum(dz) - dz/2;<br class=""><br class="">[x,z] = meshgrid(x,z);<br class="">%--------<br class=""><br class="">%--------<br class="">xbot = x(1,:);<br class="">%Hbot = 110*(sin(2*pi*5/(60*10^3).*xbot(1:end))+1)- Hz;<br class="">Hbot = [110*(zeros(1,577))- Hz 110*(sin(10^(-3).*xbot(578:1413))+1)- Hz 110*(zeros(1,2000-1413))- Hz];<br class="">%--------<br class=""><br class=""><br class="">%--------<br class="">U = zeros(nz,nx);<br class="">V = zeros(nz,nx);<br class="">T = (1e-3)^2 *(z+Hz) /9.81/2e-4;<br class="">R = T(:,1);<br class="">%--------<br class=""><br class="">%--------<br class="">U = flipud(U)';<br class="">V = flipud(V)';<br class="">T = flipud(T)';<br class="">R = flipud(R);<br class="">dz = flipud(dz);<br class=""><br class="">fid = fopen('U.init','w','b'); <br class="">fwrite(fid,U,'real*8'); <br class="">fclose(fid);<br class=""><br class="">fid = fopen('V.init','w','b'); <br class="">fwrite(fid,V,'real*8'); <br class="">fclose(fid);<br class=""><br class="">fid = fopen('T.init','w','b'); <br class="">fwrite(fid,T,'real*8'); <br class="">fclose(fid);<br class=""><br class="">save('R.init','R','-ascii');<br class=""><br class="">fid = fopen('topog.init','w','b'); <br class="">fwrite(fid,Hbot(1,:)','real*8'); <br class="">fclose(fid);<br class=""><br class="">fid=fopen('delZ.init','w','b'); <br class="">fwrite(fid,dz,'real*8'); <br class="">fclose(fid);<br class=""><br class="">U = flipud(U');<br class="">V = flipud(V');<br class="">T = flipud(T');<br class="">R = flipud(R);<br class="">dz = flipud(dz);<br class="">%--------<br class=""><br class="">%--------OBCS<br class="">Uobcs = 0.01*ones(nz,ny);<br class="">Uobcsref = flipud(Uobcs');<br class="">Uobcsnew = flipud(Uobcs');<br class=""><br class="">%tfin = 28800*30/21600;<br class="">tfin=1000;<br class="">for i=1:tfin<br class=""> Uobcsnew=cat(3,Uobcsref,Uobcsnew);<br class="">end<br class="">Uobcs = Uobcsnew ;<br class="">Vobcs = Uobcs * 0.53e-4/1.4e-4 ;<br class="">Wobcs = zeros(ny,nz,tfin+1) ;<br class=""><br class="">fid = fopen('u.box','w','b'); <br class="">fwrite(fid,Uobcs,'real*8'); <br class="">fclose(fid);<br class=""><br class="">fid = fopen('v.box','w','b'); <br class="">fwrite(fid,Vobcs,'real*8'); <br class="">fclose(fid);<br class=""><br class="">fid = fopen('w.box','w','b'); <br class="">fwrite(fid,Wobcs,'real*8'); <br class="">fclose(fid);</div><div class=""><br class=""></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le mar. 25 juin 2019 à 19:29, Matthew Mazloff <<a href="mailto:mmazloff@ucsd.edu" class="">mmazloff@ucsd.edu</a>> a écrit :<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;" class="">You have something set wrong. But this is not enough information for me to determine what. Look in other outputs (e.g. STDOUT) to see if any other information is given.<div class=""><br class=""></div><div class="">Matt</div><div class=""><br class=""></div><div class=""><br class=""><div class=""><br class=""><blockquote type="cite" class=""><div class="">On Jun 25, 2019, at 4:44 AM, Kevin Ha <<a href="mailto:kevin.ha.pro@gmail.com" target="_blank" class="">kevin.ha.pro@gmail.com</a>> wrote:</div><br class="gmail-m_-411519106908035758Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="">Ok so now after putting a [(nx * ny) * nTimeSteps * deltaT/externForcingPeriod]+1, I got this message (see below). Can someone explain ?<br class=""></div><div class=""><br class=""></div><div class="">STOP ABNORMAL END: S/R OBCS_CHECK</div>STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class="">STOP ABNORMAL END: S/R OBCS_CHECK<br class=""></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le mar. 25 juin 2019 à 10:51, Martin Losch <<a href="mailto:Martin.Losch@awi.de" target="_blank" class="">Martin.Losch@awi.de</a>> a écrit :<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Kevin,<br class="">
<br class="">
for 2D forcing you need to have [(nx * ny) * nTimeSteps * deltaT/externForcingPeriod]+1, if you use externForcingPeriod in data. Check out some of the examples, e.g. global_ocean.90x40x15, exp4, tutorial_global_oce_latlon ...<br class="">
<br class="">
If you are using pkg/exf to read the forcing, you do something similar: based on your ${var}period and your number of timesteps (i.e. model integration time span nTimeSteps * deltaT) you determine how many times you need to read data. In exf it is slightly more complicated, because the time for which the forcing is meant, can be specified explicitly. Compare example global_with_exf to global_ocean.90x40x15. They do the same, one with exf, one without.<br class="">
<br class="">
Another common mistake is to prepare real*4 data and try to read it as real*8. Then you'll have only half the records that you need (plus they are completely garbled).<br class="">
<br class="">
M.<br class="">
<br class="">
> On 25. Jun 2019, at 09:05, Kevin Ha <<a href="mailto:kevin.ha.pro@gmail.com" target="_blank" class="">kevin.ha.pro@gmail.com</a>> wrote:<br class="">
> <br class="">
> Hello Matt,<br class="">
> <br class="">
> Thank you for your quick reply. My time parameters are the following (see below). I have 480 time records in my u.box file. How many should I have ?<br class="">
> <br class="">
> # Time stepping parameters<br class="">
> &PARM03<br class="">
> <br class="">
> niter0 = 0,<br class="">
> <br class="">
> nTimeSteps = 28800, ! 10j avec deltaT=30<br class="">
> deltaT = 30.,<br class="">
> dumpFreq = 1800.,<br class="">
> monitorFreq = 43200.,<br class="">
> <br class="">
> abEps = 0.1,<br class="">
> pChkptFreq = 86400.,<br class="">
> chkptFreq = 0.0,<br class="">
> &<br class="">
> <br class="">
> Le lun. 24 juin 2019 à 20:32, Matthew Mazloff <<a href="mailto:mmazloff@ucsd.edu" target="_blank" class="">mmazloff@ucsd.edu</a>> a écrit :<br class="">
> Hi Kevin<br class="">
> <br class="">
> You don’t have enough time records in your u.box file. The model wants more. <br class="">
> <br class="">
> If you think that you do have enough time records, then I would need more info to diagnose the issue.<br class="">
> <br class="">
> Matt<br class="">
> <br class="">
> <br class="">
> > On Jun 24, 2019, at 7:38 AM, Kevin Ha <<a href="mailto:kevin.ha.pro@gmail.com" target="_blank" class="">kevin.ha.pro@gmail.com</a>> wrote:<br class="">
> > <br class="">
> > Dear community,<br class="">
> > <br class="">
> > I am having an error message when implementing obcs package to my configuration. After the run I get the message written below. Can someone explain what is wrong ? The u.box file is my specified boundary condition for velocity u at the north face of the domain (the file is generated using matlab code and has dimension Ny*Nz*Nt).<br class="">
> > <br class="">
> > Thank you in advance<br class="">
> > <br class="">
> > Best<br class="">
> > <br class="">
> > Kevin<br class="">
> > <br class="">
> > Fortran runtime error: Non-existing record number<br class="">
> > <br class="">
> > Error termination. Backtrace:<br class="">
> > At line 1873 of file mdsio_read_section.f (unit = 9, file = 'u.box')<br class="">
> > At line 1873 of file mdsio_read_section.f (unit = 9, file = 'u.box')<br class="">
> > Fortran runtime error: Non-existing record number<br class="">
> > <br class="">
> > Error termination. Backtrace:<br class="">
> > Fortran runtime error: Non-existing record number<br class="">
> > <br class="">
> > Error termination. Backtrace:<br class="">
> > #0 0x7fb72c0b7607 in ???<br class="">
> > #1 0x7fb72c0b8115 in ???<br class="">
> > #2 0x7fb72c0b8869 in ???<br class="">
> > #3 0x7fb72c17edc9 in ???<br class="">
> > #0 0x7ff59cbfc607 in ???<br class="">
> > #1 0x7ff59cbfd115 in ???<br class="">
> > #2 0x7ff59cbfd869 in ???<br class="">
> > #3 0x7ff59ccc3dc9 in ???<br class="">
> > #4 0x43d946 in ???<br class="">
> > #5 0x4b76f4 in ???<br class="">
> > #6 0x4a1086 in ???<br class="">
> > #7 0x49a19e in ???<br class="">
> > #8 0x4a60a5 in ???<br class="">
> > #9 0x5bdbb7 in ???<br class="">
> > #10 0x5b0924 in ???<br class="">
> > #11 0x5db6fe in ???<br class="">
> > #12 0x5db8c4 in ???<br class="">
> > #13 0x511c35 in ???<br class="">
> > #14 0x402366 in ???<br class="">
> > #15 0x7ff59bee9724 in ???<br class="">
> > #16 0x402398 in ???<br class="">
> > at ../sysdeps/x86_64/start.S:118<br class="">
> > #17 0xffffffffffffffff in ???<br class="">
> > #4 0x43d946 in ???<br class="">
> > #5 0x4b76f4 in ???<br class="">
> > #6 0x4a1086 in ???<br class="">
> > #7 0x49a19e in ???<br class="">
> > #8 0x4a60a5 in ???<br class="">
> > #9 0x5bdbb7 in ???<br class="">
> > #10 0x5b0924 in ???<br class="">
> > #11 0x5db6fe in ???<br class="">
> > -------------------------------------------------------<br class="">
> > Primary job terminated normally, but 1 process returned<br class="">
> > a non-zero exit code.. Per user-direction, the job has been aborted.<br class="">
> > -------------------------------------------------------<br class="">
> > --------------------------------------------------------------------------<br class="">
> > mpirun detected that one or more processes exited with non-zero status, thus causing<br class="">
> > the job to be terminated. The first process to do so was:<br class="">
> > <br class="">
> > Process name: [[9039,1],27]<br class="">
> > Exit code: 2<br class="">
> > --------------------------------------------------------------------------<br class="">
> > _______________________________________________<br class="">
> > MITgcm-support mailing list<br class="">
> > <a href="mailto:MITgcm-support@mitgcm.org" target="_blank" class="">MITgcm-support@mitgcm.org</a><br class="">
> > <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">
> <br class="">
> _______________________________________________<br class="">
> MITgcm-support mailing list<br class="">
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank" class="">MITgcm-support@mitgcm.org</a><br class="">
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">
> _______________________________________________<br class="">
> MITgcm-support mailing list<br class="">
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank" class="">MITgcm-support@mitgcm.org</a><br class="">
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">
<br class="">
_______________________________________________<br class="">
MITgcm-support mailing list<br class="">
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank" class="">MITgcm-support@mitgcm.org</a><br class="">
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">
</blockquote></div>
_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" target="_blank" class="">MITgcm-support@mitgcm.org</a><br class=""><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" target="_blank" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class=""></div></blockquote></div><br class=""></div></div>_______________________________________________<br class="">
MITgcm-support mailing list<br class="">
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank" class="">MITgcm-support@mitgcm.org</a><br class="">
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">
</blockquote></div>
_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br class=""></div></blockquote></div><br class=""></div></body></html>