Changeset 2410 for branches/nemo_v3_3_beta/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90

Timestamp:

2010-11-19T17:42:36+01:00 (13 years ago)

Author:

cetlod

Message:

Small changes to get a better behavior of NEMO on IDRIS/IBM SP6 machine with fortran compiler xlf 13.1.0.3, see ticket #761

File:

• branches/nemo_v3_3_beta/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90 (modified) (1 diff)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90

@@ -485 +485 @@
       ENDIF
 
-      ztmp(:,:) = 1. - AINT( frld(:,:), wp )                        ! return 1 as soon as there is ice
+!! ce     ztmp(:,:) = 1. - AINT( frld(:,:), wp )                        ! return 1 as soon as there is ice
+!! ce     A big warning because the model crashes on IDRIS/IBM SP6 with xlf 13.1.0.3, see ticket #761
+!! ce     We Unroll the loop and everything works fine 
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            ztmp(ji,jj) = 1. - AINT( frld(ji,jj), wp )                ! return 1 as soon as there is ice
+         END DO
+      END DO
+      !
       CALL iom_put( 'ice_pres'  , ztmp                            )   ! Ice presence                          [-]
       CALL iom_put( 'ist_ipa'   , ( sist(:,:) - rt0 ) * ztmp(:,:) )   ! Ice surface temperature               [Celius]