[MITgcm-devel] bug in autodiff_restore

[Martin Losch]

Hi there,

I just found a bug in autodiff_restore.F:

lines like
           xx_obcsn0(I,K,i,j,1)   = StoreOBCSN(I,K,bi,bj,11)
should really be i,j-> bi,bj
           xx_obcsn0(I,K,bi,bj,1)   = StoreOBCSN(I,K,bi,bj,11)

I am not sure how much of an effect this will have on any existing   
simulation, but this is clearly wrong and caused serious memory issues  
with gfortran. I will check the fix in, unless someone objects ...

M.