[MITgcm-devel] overlap question for diagnostic field oceSPtnd

[An Nguyen]

Hello again, 

To answer my own question, I think I have found a bug in the saltplume code in:
salt_plume_apply.F, lines 148-149, where instead of looping through the range 1-OL[x,y],sN[x,y]+OL[x,y], it is currently doing the loop only through:

         DO j=1,OLy
          DO i=1,OLx

I'll debug some more, and will send updated code through github soonish.
Cheers,
An

> On Dec 31, 2018, at 04:57, An Nguyen <antnguyen13 at gmail.com> wrote:
> 
> Hello MITgcm gurus,
> 
> I am using advection scheme 7 with overlapping Olx=Oly=8, and I have been getting the 3D diagnostics oceSPtnd looking like this (I hope the attachment work, if not i'm including also a link to a figure) where the checkerboard squares are exactly 8x8 in size.  This is a plot for the top level (k=1).. The problem exists for other depths too. 
> 
> This diagnostic is filled inside salt_plume_tendency_apply_s.F 
>        CALL DIAGNOSTICS_FILL (
>      &      plumetend,'oceSPtnd',k,1,2,bi,bj,myThid )
> 
> with arguments specifying that we're passing in a tile with overlaping regions.
> 
> I can't seem to find where the problem might be.  There is an exact diagnostics, which is filled outside of that subroutine under the diagnostic name SPtd_SLT, and it has the same problem.  
> 
> Could you suggest where I can get started to diagnose the problem? it looks like there is a missing exchange somewhere ?
> 
> Thanks,
> An
> 
> link to figure in case it's not seen with this message:
> https://utexas.box.com/s/pg8s0zq9gv47cg6xtht755shndk7s9or
> 
> <oceSPtnd.jpg>