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

Timestamp:

2008-03-07T18:11:35+01:00 (16 years ago)

Author:

ctlod

Message:

Clean comments and useless lines, see ticket:#72

File:

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

trunk/NEMO/LIM_SRC_3/limvar.F90

@@ -1 +1 @@
 MODULE limvar
 #if defined key_lim3
+   !!----------------------------------------------------------------------
+   !!   'key_lim3'                                      LIM3 sea-ice model
+   !!----------------------------------------------------------------------
    !!======================================================================
    !!                       ***  MODULE limvar ***
@@ -38 +41 @@
    USE ice_oce          ! ice variables
    USE thd_ice
-   USE limicepoints
    USE in_out_manager
    USE ice
@@ -62 +64 @@
 
    !!----------------------------------------------------------------------
-   !!   LIM-@ 4.0,  UCL-ASTR (2006)
+   !!   LIM3.0,  UCL-ASTR-LOCEAN-IPSL (2008)
    !!   (c) UCL-ASTR and Martin Vancoppenolle
    !!----------------------------------------------------------------------
@@ -91 +93 @@
         !! ** History :
         !!           (01-2006) Martin Vancoppenolle, UCL-ASTR
-        !!                     Take it easy man
-        !!                     Life is just a simple game, between ups / and downs \ :@)
         !!
         !! note : you could add an argument when you need only at_i, vt_i
@@ -103 +103 @@
                      jj,       &   ! spatial dummy loop index
                      jk,       &   ! vertical layering dummy loop index
-                     jl,       &   ! ice category dummy loop index
-                     jm,       &   ! ice types    dummy loop index
-                     index
+                     jl            ! ice category dummy loop index
  
         REAL ::      zeps, epsi16, zinda, epsi06
@@ -116 +114 @@
        epsi16 = 1.0e-16
        epsi06 = 1.0e-6
-
-      IF(lwp) THEN
-         WRITE(numout,*)
-         WRITE(numout,*) 'lim_var_agg : Aggregate ice variables '
-         WRITE(numout,*) '~~~~~~~~~~~'
-      ENDIF
 
        !------------------
@@ -228 +220 @@
                      jj,       &   ! spatial dummy loop index
                      jk,       &   ! vertical layering dummy loop index
-                     jl,       &   ! ice category dummy loop index
-                     jm,       &   ! ice types    dummy loop index
-                     index
+                     jl            ! ice category dummy loop index
 
         REAL :: zq_i, zaaa, zbbb, zccc, zdiscrim, &
@@ -242 +232 @@
 !!-- End of declarations
 !!------------------------------------------------------------------------------
-!     IF (lwp) THEN
-!     WRITE(numout,*)
-!     WRITE(numout,*) 'lim_var_glo2eqv : transform global variables into equivalent ice variables '
-!     WRITE(numout,*) '~~~~~~~~~~~~~~~'
-!     ENDIF
 
       !-------------------------------------------------------
@@ -377 +362 @@
         !!
         !!------------------------------------------------------------------
-        !! * Arguments
-
-        !! * Local variables
-        INTEGER ::   ji,       &   ! spatial dummy loop index
-                     jj,       &   ! spatial dummy loop index
-                     jk,       &   ! vertical layering dummy loop index
-                     jl,       &   ! ice category dummy loop index
-                     jm,       &   ! ice types    dummy loop index
-                     index
-
-        REAL :: zq_i, zaaa, zbbb, zccc, zdiscrim, &
-                ztmelts
-
-        REAL :: zeps
-
-        !Ice specific heat, massive latent heat
-        zeps    = 1.0e-10
-
-!!-- End of declarations
-!!------------------------------------------------------------------------------
-
-!      IF (lwp) THEN
-!      WRITE(numout,*)
-!      WRITE(numout,*) 'lim_var_eqv2glo : transform equivalent variables into global ice variables '
-!      WRITE(numout,*) '~~~~~~~~~~~~~~~'
-!      ENDIF
 
        v_i(:,:,:)   = ht_i(:,:,:) * a_i(:,:,:)
@@ -408 +367 @@
        smv_i(:,:,:) = sm_i(:,:,:) * v_i(:,:,:)
        oa_i (:,:,:) = o_i (:,:,:) * a_i(:,:,:)
-
-!      DO jl = 1, jpl
-!        DO jk = 1, nlay_i
-!          DO jj = 1, jpj
-!            DO ji = 1, jpi
-!               ! heat content per unit volume
-!               e_i(ji,jj,jk,jl) = zidto * rhoic * &
-!               (   cpic    * ( ztmelts - t_i(ji,jj,jk,jl) ) &
-!               + lfus  * ( 1.0 - (ztmelts-rtt) / MIN((t_i(ji,jj,jk,jl)-rtt),-zeps) ) &
-!               - rcp      * ( ztmelts - rtt ) &
-!               )
-
-!               ! Correct dimensions to avoid big values
-!                    e_i(ji,jj,jk,jl) = e_i(ji,jj,jk,jl) / unit_fac
-
-!               ! Mutliply by ice volume, and divide by number of layers to get heat content in 10^9 J
-!                    e_i(ji,jj,jk,jl) = e_i(ji,jj,jk,jl) * &
-!                                    area(ji,jj) * a_i(ji,jj,jl) * ht_i(ji,jj,jl) / &
-!                                    FLOAT(nlay_i)
-
-!           END DO
-!         END DO
-!       END DO
-!     END DO
 
     END SUBROUTINE lim_var_eqv2glo
@@ -787 +722 @@
 
 #endif
- END MODULE limvar
+END MODULE limvar