[MITgcm-devel] exch2 problem?

[Martin Losch]

Hi there,

Yossi is running a zonally averaged model with Nx=1, Ny = 84, Nr=120 on a lat-lon grid. I attached a simple configuration which starts from uniform salinity/temperature and has some heat flux to make the model move.

We find that it works with one tile, but not with more than one (the attached SIZE.h uses 2 tiles on 1processor, but 2processore with 1 tile each gives the same). It explodes once the cg2d solver has something to do (basically cg2d no longer converges). Yossi finds that it also works with more than one tile, if he increases: sNx=2.

The problem goes away when I recompile without exch2, so it seems to be related to exch2. While I don't think that exch2 has to work with such a configuation (and exch2 is only in there because Yossi went from a cs32 configuration to this lon-lat configuration), I do think it would be good, if we could either fix the problem, or add a stop when (sNx=1||sNy=1)&&exch2. What do you think? [Maybe I overlooked something and there is problem with the configuration, but not with exch2?]

Martin


-------------- next part --------------
A non-text attachment was scrubbed...
Name: forYossi.tgz
Type: application/octet-stream
Size: 34205 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20120324/5fc105ca/attachment-0001.obj>