Bug in mod_boundary.F

[annsp]

I found a bug in mod_boundary causing segmentation fault when ADJUST_BOUNDARY is defined.
It is only a problem when the southern boundary is closed while the northern boundary is open.
Changing:

Code: Select all

      IF (tl_LBC(isouth,isUbar,ng)%acquire) THEN
        allocate ( BOUNDARY(ng) % tl_ubar_north(LBi:UBi) )
      END IF

to:

Code: Select all

      IF (tl_LBC(inorth,isUbar,ng)%acquire) THEN
        allocate ( BOUNDARY(ng) % tl_ubar_north(LBi:UBi) )
      END IF

fixes the problem.

Ann Kristin

[leon]

Yes, you are right. I ran into the same problem too. It's in mod_boundary.F line1100. isouth should be inorth.