[MITgcm-support] Gridding problem in the c64f release

[François Lefeuvre]

Hi all,

I am a student in Master2 who is recently working on the MITgcm. I was
trying to change the size of my domain (make me change the overlap values)
and when I run my simulation, I got an error which is:
 S/R EXCH1_RS: myOLn .GT. MAX_OLY_EXCH

So, I went for this problem in the code in the file exch1_rs.F
(MITgcm_c64f/eesupp/src) and I found the place where the problem is. Here
is two line from this file:
 IF ( myOLn      .GT. MAX_*OL**X*_EXCH )
     &  STOP ' S/R EXCH1_RS: myOLn .GT. MAX_*OLY*_EXCH'

If you look well to these two lines, you see a contradiction between the
error message and the condition in the if (I put it in bold above). My
feeling, when I watched to the others "error checks" which are in the same
part of the file, was that the error is contained in the if condition. The
"_OLX_" should be an "_OLY_" (as the error message suggests it). Indeed,
The comparison between the maximum of Olx (zonal direction) with the north
overlap region size doesn't appear necessary to me. And the others "error
checks" make me think that I am right.
This error is minor because we often use kind of "symmetric" grid but I
want to know if I am right or if there is some subtlety that I miss. I have
also tried to look at this problem in the newiest release of the code but
this file doesn't exist. It seems that all the "exch_" files were
reorganized. So I don't know if this error has been propagated in new
releases.

My apologies for the english and thank you!

François Lefeuvre.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130808/5d4f9bf1/attachment.htm>