[MITgcm-devel] Beaufort experiment on mac os x

[Martin Losch]

It get's rid of the segmentation fault (probably because it resolves a situation that gfortran -O3 cannot deal with)

Martin

On Mar 15, 2012, at 3:19 PM, Menemenlis, Dimitris (3248) wrote:

> 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.
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel