Changeset 8565 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90

Timestamp:

2017-09-27T12:09:10+02:00 (7 years ago)

Author:

clem

Message:

trying to respect naming convention

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90 (modified) (5 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90

@@ -378 +378 @@
       CASE( 1 )       !  constant salinity in time and space  !
          !            !---------------------------------------!
-         sz_i_1d(1:nidx,:) = rn_icesal
+         sz_i_1d(1:npti,:) = rn_icesal
          !
          !            !---------------------------------------------!
@@ -387 +387 @@
          !
          !                                      ! Slope of the linear profile 
-         WHERE( h_i_1d(1:nidx) > epsi20 )   ;   z_slope_s(1:nidx) = 2._wp * s_i_1d(1:nidx) / h_i_1d(1:nidx)
-         ELSEWHERE                           ;   z_slope_s(1:nidx) = 0._wp
+         WHERE( h_i_1d(1:npti) > epsi20 )   ;   z_slope_s(1:npti) = 2._wp * s_i_1d(1:npti) / h_i_1d(1:npti)
+         ELSEWHERE                           ;   z_slope_s(1:npti) = 0._wp
          END WHERE
          
          z1_dS = 1._wp / ( zsi1 - zsi0 )
-         DO ji = 1, nidx
+         DO ji = 1, npti
             zalpha(ji) = MAX(  0._wp , MIN(  ( zsi1 - s_i_1d(ji) ) * z1_dS , 1._wp  )  )
             !                             ! force a constant profile when SSS too low (Baltic Sea)
@@ -400 +400 @@
          ! Computation of the profile
          DO jk = 1, nlay_i
-            DO ji = 1, nidx
+            DO ji = 1, npti
                !                          ! linear profile with 0 surface value
                zs0 = z_slope_s(ji) * ( REAL(jk,wp) - 0.5_wp ) * h_i_1d(ji) * r1_nlay_i
@@ -414 +414 @@
          !            !-------------------------------------------!                                   (mean = 2.30)
          !
-         s_i_1d(1:nidx) = 2.30_wp
+         s_i_1d(1:npti) = 2.30_wp
          !
 !!gm cf remark in ice_var_salprof routine, CASE( 3 )
@@ -420 +420 @@
             zargtemp  = ( REAL(jk,wp) - 0.5_wp ) * r1_nlay_i
             zsal =  1.6_wp * ( 1._wp - COS( rpi * zargtemp**( 0.407_wp / ( 0.573_wp + zargtemp ) ) ) )
-            DO ji = 1, nidx
+            DO ji = 1, npti
                sz_i_1d(ji,jk) = zsal
             END DO