[MITgcm-devel] Beaufort experiment on mac os x

[Menemenlis, Dimitris (3248)]

Martin, does the change below get rid of segmentation fault or just help with vectorization?

Dimitris Menemenlis

On Mar 15, 2012, at 2:21 AM, Martin Losch wrote:

> Alternatively, we could replace
>        DO j=1-Oly,sNy+Oly
>         DO i=1-Olx,sNx+Olx
>          R_low(i,j,bi,bj) = rF(1)
>          DO k=Nr,1,-1
>           R_low(i,j,bi,bj) = R_low(i,j,bi,bj)
>     &                      - drF(k)*hFacC(i,j,k,bi,bj)
>          ENDDO
>         ENDDO
>        ENDDO
> 
> with
>        DO j=1-Oly,sNy+Oly
>         DO i=1-Olx,sNx+Olx
>          R_low(i,j,bi,bj) = rF(1)
>         ENDDO
>        ENDDO
>        DO k=Nr,1,-1
>         DO j=1-Oly,sNy+Oly
>          DO i=1-Olx,sNx+Olx
>           R_low(i,j,bi,bj) = R_low(i,j,bi,bj)
>     &                      - drF(k)*hFacC(i,j,k,bi,bj)
>          ENDDO
>         ENDDO
>        ENDDO
> That does work too (and makes sense, if you care about vectoriziation, although in this initialisation routine, performance is not an issue).
> 
> M.