<div dir="ltr"><div>Thank you Martin and Estanislao,</div><div><br></div><div>in my setup I only have a full Northern OB, Southern closed boundary and West and East partial land, partial sea. <br></div><div>My problem was that I saw values from the W boundary repeated on the Eastern boundary for eta, which is not possible to specify as OBetaFile with my options. So I was wondering if there was this periodicity at work. <br></div><div>Moreover I experienced an error at initialization of OBCS, a mismatch between the OB which I specified and the mask, where I think the East OB was checked against the land/sea mask on the West side. Setting either insideOBmaskFile, or specify also the OB files on West side seem to have solved the issue.</div><div><br></div><div>Kind regards,</div><div><br></div><div>Enrico<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mer 13 nov 2019 alle ore 10:25 Martin Losch <<a href="mailto:Martin.Losch@awi.de">Martin.Losch@awi.de</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Enrico,<br>
<br>
I am not sure what you are trying to do. The “standard” way of supressing the periodicity is to impose one wall at either end of the domain to arrive at a channel or, if you do it for both horizontal directions, a box with solid walls.<br>
<br>
If you want to use open boundaries (i.e. the pkg/obcs), then the periodicity is still there as the default behaviour, but it is overwritten by the obcs values. The “halos” may still contain values based on the periodicity, but they are never used.<br>
<br>
Martin<br>
<br>
> On 13. Nov 2019, at 07:34, ESTANISLAO GAVILAN PASCUAL-AHU <<a href="mailto:e.gavilan@hhu.edu.cn" target="_blank">e.gavilan@hhu.edu.cn</a>> wrote:<br>
> <br>
> <br>
> Hi Pochini,<br>
> <br>
> I am no sure about that option, but if you do not want to have a periodic open boundary (i.e. time dependent). You could just set the forcing period to zero. From the data file you could write this:<br>
> periodicExternalForcing= .FALSE.,<br>
> externForcingPeriod= 0.,<br>
> externForcingCycle = 120000.,<br>
> or from the the data.exf:<br>
> &EXF_NML_OBCS<br>
> <br>
> obcsNstartdate1=19781216,<br>
> obcsNstartdate2=180000,<br>
> obcsNperiod=0,<br>
> obcsNrepCycle=120000,<br>
> obcsSstartdate1=19781216,<br>
> obcsSstartdate2=180000,<br>
> obcsSperiod=0,<br>
> obcsSrepCycle=120000,<br>
> <br>
> Kind regards,<br>
> <br>
> Estanislao<br>
> <br>
>> Hi all,<br>
>> <br>
>> I am trying to deactivate the default behaviour of periodic boundaries. I<br>
>> tried setting in eedata notUsingXPeriodicity=.TRUE.,<br>
>> notUsingYPeriodicity=.TRUE., but the option is not recognized. Setting<br>
>> forcefully in eeparms_set.F, where it appears in COMMON block EEPARMS,<br>
>> cause the following error to appear:<br>
>> <br>
>> *** An error occurred in MPI_Cart_rank<br>
>>> *** reported by process [3771400193,0]<br>
>>> *** on communicator MPI_COMMUNICATOR 3<br>
>>> *** MPI_ERR_ARG: invalid argument of some other kind<br>
>>> *** MPI_ERRORS_ARE_FATAL (processes in this communicator will now abort,<br>
>>> ***    and potentially your MPI job)<br>
>>> 6 more processes have sent help message help-mpi-errors.txt /<br>
>>> mpi_errors_are_fatal<br>
>>> Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error<br>
>>> messages<br>
>>> <br>
>> <br>
>> The two lines<br>
>> #define CAN_PREVENT_X_PERIODICITY<br>
>> #define CAN_PREVENT_Y_PERIODICITY<br>
>> are defined in CPP_EEOPTIONS.h<br>
>> <br>
>> Does the periodic domain definition interact with OB specifications?<br>
>> Is there a way to solve it without setting bathymetry=0 at margins?<br>
>> <br>
>> Thanks,<br>
>> <br>
>> Enrico<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>