Changeset 4072 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limcat_1D.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_1D.F90 (modified) (4 diffs)

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

@@ -40 +40 @@
       !! Local variables
       INTEGER  :: ji, jk, jl             ! dummy loop indices
-      INTEGER  :: zjpi, i_fill, jl0, ztest_1, ztest_2, ztest_3, ztest_4, ztests  
+      INTEGER  :: ijpij, 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
@@ -50 +50 @@
       ! initialisation of variables
       !--------------------------------------------------------------------
-      !WRITE(numout,*) 'size(zht_i)=',SIZE(zht_i)
-      !WRITE(numout,*) 'size(zht_s)=',SIZE(zht_s)
-      !WRITE(numout,*) 'size(za_i)=',SIZE(za_i)
-      !WRITE(numout,*) 'size(zhti)=',SIZE(zhti)
-      !WRITE(numout,*) 'size(zhts)=',SIZE(zhts)
-      !WRITE(numout,*) 'size(zai)=',SIZE(zai)
-
-      zjpi = SIZE(zhti,1)
-      zht_i(1:zjpi,1:jpl) = 0.d0
-      zht_s(1:zjpi,1:jpl) = 0.d0
-      za_i (1:zjpi,1:jpl) = 0.d0
+      ijpij = SIZE(zhti,1)
+      zht_i(1:ijpij,1:jpl) = 0.d0
+      zht_s(1:ijpij,1:jpl) = 0.d0
+      za_i (1:ijpij,1:jpl) = 0.d0
 
       !------------------------------------------------------------------------------------
@@ -73 +66 @@
       !         fulfills ice volume concervation between input and output (ztests=4) 
       !--------------------------------------------------------------------------------------
-      IF(lwp) THEN
-         !WRITE(numout,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ '
-         !WRITE(numout,*) 'lim_cat_1D : Distribute ice categories '
-         !WRITE(numout,*) '  of thick. = ', ( hi_max(jl), jl = 0, jpl )
-         !WRITE(numout,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ '
-         !WRITE(numout,*) 'zjpi=',zjpi
-      END IF
  
       ! ----------------------------------------
       ! distribution over the jpl ice categories
       ! ----------------------------------------
-      DO ji = 1, zjpi
+      DO ji = 1, ijpij
         ! snow thickness in each category 
          zht_s(ji,1:jpl) = zhts(ji)
@@ -181 +167 @@
       END DO ! i loop
 
-      !WRITE(numout,*) 'clem zai=',MAXVAL(zai)
-      !WRITE(numout,*) 'clem za_i(:,jpl)=',(MAXVAL(za_i(:,jl)),jl=1,jpl)
-      !WRITE(numout,*) 'clem zhti=',MAXVAL(zhti)
-      !WRITE(numout,*) 'clem zht_i(:,jpl)=',(MAXVAL(zht_i(:,jl)),jl=1,jpl)
       IF( nn_timing == 1 )  CALL timing_stop('limcat_1D')