Changeset 4072 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limcat_2D.F90

Timestamp:

2013-10-17T15:10:59+02:00 (11 years ago)

Author:

clem

Message:

few rewritings

File:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limcat_2D.F90 (modified) (3 diffs)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limcat_2D.F90

@@ -39 +39 @@
       !! Local variables
       INTEGER  :: ji, jj, jk, jl             ! dummy loop indices
-      INTEGER  :: zjpi, zjpj, i_fill, jl0, ztest_1, ztest_2, ztest_3, ztest_4, ztests  
+      INTEGER  :: ijpi, ijpj, i_fill, jl0, ztest_1, ztest_2, ztest_3, ztest_4, ztests  
       REAL(wp) :: zarg, zV, zconv
       REAL(wp), DIMENSION(:,:),   INTENT(in)  ::   zhti, zhts, zai    ! input ice/snow variables
@@ -48 +48 @@
       ! initialisation of variables
       !--------------------------------------------------------------------
-      zjpi = SIZE(zhti,1)
-      zjpj = SIZE(zhti,2)
-      zht_i(1:zjpi,1:zjpj,1:jpl) = 0.d0
-      zht_s(1:zjpi,1:zjpj,1:jpl) = 0.d0
-      za_i (1:zjpi,1:zjpj,1:jpl) = 0.d0
+      ijpi = SIZE(zhti,1)
+      ijpj = SIZE(zhti,2)
+      zht_i(1:ijpi,1:ijpj,1:jpl) = 0.d0
+      zht_s(1:ijpi,1:ijpj,1:jpl) = 0.d0
+      za_i (1:ijpi,1:ijpj,1:jpl) = 0.d0
 
       !------------------------------------------------------------------------------------
@@ -65 +65 @@
       !         fulfills ice volume concervation between input and output (ztests=4) 
       !--------------------------------------------------------------------------------------
-      IF(lwp) THEN
-         WRITE(numout,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ '
-         WRITE(numout,*) 'lim_cat_2D : Distribute ice categories '
-         WRITE(numout,*) '  of thick. = ', ( hi_max(jl), jl = 0, jpl )
-         WRITE(numout,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ '
-      END IF
-      !WRITE(numout,*) ' * zht_i_ini : ', ( zht_i_ini(1,1,jl), jl = 1, jpl )
-      !WRITE(numout,*) ' * za_i_ini :  ', ( za_i_ini(1,1,jl), jl = 1, jpl )
  
       ! ----------------------------------------
       ! distribution over the jpl ice categories
       ! ----------------------------------------
-      DO jj = 1, zjpj
-         DO ji = 1, zjpi
+      DO jj = 1, ijpj
+         DO ji = 1, ijpi
 
             ! snow thickness in each category