Changeset 868 for trunk/NEMO/LIM_SRC_3/limvar.F90

Timestamp:

2008-03-14T19:53:00+01:00 (16 years ago)

Author:

rblod

Message:

First optimisation of LIM3, limrhg optimisation induces computation change

File:

• trunk/NEMO/LIM_SRC_3/limvar.F90 (modified) (11 diffs)

trunk/NEMO/LIM_SRC_3/limvar.F90

@@ -236 +236 @@
       ! Ice thickness, snow thickness, ice salinity, ice age
       !-------------------------------------------------------
+!CDIR NOVERRCHK
       DO jl = 1, jpl
+!CDIR NOVERRCHK
          DO jj = 1, jpj
+!CDIR NOVERRCHK
             DO ji = 1, jpi
                zindb          = 1.0-MAX(0.0,SIGN(1.0,- a_i(ji,jj,jl))) !0 if no ice and 1 if yes
@@ -249 +252 @@
       IF ( ( num_sal .EQ. 2 ) .OR. ( num_sal .EQ. 4 ) )THEN
 
+!CDIR NOVERRCHK
       DO jl = 1, jpl
+!CDIR NOVERRCHK
          DO jj = 1, jpj
+!CDIR NOVERRCHK
             DO ji = 1, jpi
                zindb          = 1.0-MAX(0.0,SIGN(1.0,-a_i(ji,jj,jl))) !0 if no ice and 1 if yes
@@ -266 +272 @@
       ! Ice temperatures
       !-------------------
+!CDIR NOVERRCHK
       DO jl = 1, jpl
+!CDIR NOVERRCHK
         DO jk = 1, nlay_i
+!CDIR NOVERRCHK
           DO jj = 1, jpj
+!CDIR NOVERRCHK
             DO ji = 1, jpi
               !Energy of melting q(S,T) [J.m-3]
@@ -298 +308 @@
       zfac1 = 1. / ( rhosn * cpic )
       zfac2 = lfus / cpic  
+!CDIR NOVERRCHK
       DO jl = 1, jpl
+!CDIR NOVERRCHK
         DO jk = 1, nlay_s
+!CDIR NOVERRCHK
           DO jj = 1, jpj
+!CDIR NOVERRCHK
             DO ji = 1, jpi
               !Energy of melting q(S,T) [J.m-3]
@@ -321 +335 @@
       !-------------------
       tm_i(:,:) = 0.0
+!CDIR NOVERRCHK
       DO jl = 1, jpl
+!CDIR NOVERRCHK
          DO jk = 1, nlay_i
+!CDIR NOVERRCHK
             DO jj = 1, jpj
+!CDIR NOVERRCHK
                DO ji = 1, jpi
                   zindb          = 1.0-MAX(0.0,SIGN(1.0,-a_i(ji,jj,jl)))
@@ -462 +480 @@
          zalpha(:,:,:) = 0.0
 
+!CDIR NOVERRCHK
          DO jl = 1, jpl
+!CDIR NOVERRCHK
             DO jj = 1, jpj
+!CDIR NOVERRCHK
                DO ji = 1, jpi
                   ! zind0 = 1 if sm_i le s_i_0 and 0 otherwise
@@ -511 +532 @@
          sm_i(:,:,:) = 2.30
 
+!CDIR NOVERRCHK
          DO jl = 1, jpl
+!CDIR NOVERRCHK
             DO jk = 1, nlay_i
+!CDIR NOVERRCHK
                DO jj = 1, jpj
+!CDIR NOVERRCHK
                   DO ji = 1, jpi
                      zargtemp  = ( jk - 0.5 ) / nlay_i
@@ -567 +592 @@
        zeps = 1.0e-13
        bv_i(:,:) = 0.0
+!CDIR NOVERRCHK
        DO jl = 1, jpl
+!CDIR NOVERRCHK
           DO jk = 1, nlay_i
+!CDIR NOVERRCHK
              DO jj = 1, jpj
+!CDIR NOVERRCHK
                 DO ji = 1, jpi
                    zindb          = 1.0-MAX(0.0,SIGN(1.0,-a_i(ji,jj,jl))) !0 if no ice and 1 if yes
@@ -639 +668 @@
          ! Slope of the linear profile zs_zero
          !-------------------------------------
+!CDIR NOVERRCHK
          DO ji = kideb, kiut 
                z_slope_s(ji) = 2.0 * sm_i_b(ji) / MAX( 0.01      &
@@ -650 +680 @@
          dummy_fac2 = 1. / nlay_i
 
+!CDIR NOVERRCHK
          DO jk = 1, nlay_i
+!CDIR NOVERRCHK
             DO ji = kideb, kiut
                zji    =  MOD( npb(ji) - 1, jpi ) + 1
@@ -688 +720 @@
          sm_i_b(:) = 2.30
 
+!CDIR NOVERRCHK
          DO ji = kideb, kiut
+!CDIR NOVERRCHK
             DO jk = 1, nlay_i
                zargtemp  = ( jk - 0.5 ) / nlay_i