Changeset 834 for trunk/NEMO/LIM_SRC_3/limthd_dif.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/limthd_dif.F90 (modified) (16 diffs)

trunk/NEMO/LIM_SRC_3/limthd_dif.F90

@@ -1 +1 @@
 MODULE limthd_dif
 #if defined key_lim3
+   !!----------------------------------------------------------------------
+   !!   'key_lim3'                                      LIM3 sea-ice model
+   !!----------------------------------------------------------------------
    !!======================================================================
    !!                       ***  MODULE limthd_dif ***
@@ -14 +17 @@
    USE thd_ice
    USE iceini
-   USE limicepoints
    USE limistate
    USE in_out_manager
@@ -34 +36 @@
 
    !!----------------------------------------------------------------------
-   !!   LIM 3.0,  UCL-ASTR-LOCEAN-IPSL (2005)
+   !!   LIM 3.0,  UCL-ASTR-LOCEAN-IPSL (2008)
    !!   (c) UCL-ASTR and Martin Vancoppenolle
    !!----------------------------------------------------------------------
@@ -88 +90 @@
         !! ** History :
         !!           (02-2003) Martin Vancoppenolle, Louvain-la-Neuve, Belgium
-        !!           (06-2005) Martin Vancoppenolle still!!! for the 3d version...
+        !!           (06-2005) Martin Vancoppenolle, 3d version
         !!           (11-2006) Vectorized by Xavier Fettweis (UCL-ASTR)
         !!           (04-2007) Energy conservation tested by M. Vancoppenolle
@@ -106 +108 @@
                     layer,    &   ! vertical dummy loop index 
                     nconv,    &   ! number of iterations in iterative procedure
-!                   nconvmax = 50, &! maximum number of iterations for temperature computations (50)
-                    index,    &   !
                     minnumeqmin, & !
                     maxnumeqmax
@@ -158 +158 @@
           zfsw     ,   & !solar radiation absorbed at the surface
           zf       ,   & ! surface flux function
-          dzf      ,   & ! derivative of the surface flux function
-          zksn           !effective snow conductivity
+          dzf            ! derivative of the surface flux function
 
        REAL(wp)  ::           &  ! constant values
@@ -168 +167 @@
           zbeta     =  0.117,     & !: for thermal conductivity (could be 0.13)
           zraext_s  =  1.0e08,    & !: extinction coefficient of radiation in the snow
-!         zraext_i  =  1.5,       & !extinction coefficient of radiation in the ice
-!         zraext_i  =  1.0,       & !extinction coefficient of radiation in the ice
-                                    ! MY 0.87, FY 1.2 (Grenfell, Perovich and Light
-                                    ! 07)
-!         zerrmax   =  1.0e-4 ,   & !: error criterion for convergence of the iterative process
-          zfc_int           ,     & !: conductive heat flux at the interface
           zkimin    =  0.10 ,     & !: minimum ice thermal conductivity
           zht_smin  =  1.0e-4       !: minimum snow depth
 
         REAL(wp)  ::          &  ! local variables 
-          zheshth,            &  ! = zhe(ji) / thth
           ztmelt_i,           &  ! ice melting temperature
-          zerritmax,          &  ! current maximal error on temperature 
-          zexp
+          zerritmax              ! current maximal error on temperature 
 
         REAL(wp), DIMENSION(jpij)  :: &
           zerrit,             &  ! current error on temperature 
           zdifcase,           &  ! case of the equation resolution (1->4)
-          zghe,               &  ! correction factor of the thermal conductivity
-          zhe,                &  ! effective thickness for compu. of equ. thermal conductivity
           zftrice,            &  ! solar radiation transmitted through the ice
-          zihic, zihe, zhsu
-
-        CHARACTER (len=50) ::   charout
+          zihic, zhsu
+
 !!-- End of declarations
 !!----------------------------------------------------------------------------------------------
 
-       ! new namelist parameters
-!      maxer_i_thd  = maxer_i_thd
-!      thcon_i_swi = thcon_i_swi
-!      kappa_i = kappa_i
-!      nconv_i_thd = nconv_i_thd
-
        IF(lwp) WRITE(numout,*)'lim_thd_dif : Heat diffusion in sea ice for cat :', jl
 
-       thcon_i_swi = 1 ! Pringle et al. (2007)
 !
 !------------------------------------------------------------------------------!
@@ -339 +320 @@
                         zradtr_i(ji,nlay_i)
        END DO
+       ! +++++
+
        DO layer = 1, nlay_i
           DO ji = kideb , kiut
@@ -345 +328 @@
        END DO
 
-       ! +++++
        
 !
@@ -379 +361 @@
        DO WHILE ((zerritmax > maxer_i_thd).AND.(nconv < nconv_i_thd))
 
-!      ! +++++ 
-!         DO ji = kideb , kiut
-!         zji                 = MOD( npb(ji) - 1, jpi ) + 1
-!         zjj                 = ( npb(ji) - 1 ) / jpi + 1
-!         IF ( (zji.EQ.jiindex) .AND. (zjj.EQ.jjindex) ) THEN
-!            WRITE(numout,*) ' Iteration : ', nconv
-!            WRITE(numout,*) ' t_su     : ', t_su_b(ji)
-!            WRITE(numout,*) ' t_s      : ', t_s_b(ji,1)
-!            WRITE(numout,*) ' t_i      : ', t_i_b(ji,:)
-!            WRITE(numout,*) ' difcase  : ', zdifcase(ji)
-!         ENDIF
-!         END DO
-!      ! +++++ 
-
        nconv   =  nconv+1
 
@@ -410 +378 @@
 
        IF ( thcon_i_swi .EQ. 1 ) THEN
-       ! Pringle et al formula included, should be tested
+       ! Pringle et al formula included,
        ! 2.11 + 0.09 S/T - 0.011.T
        DO ji = kideb , kiut
@@ -489 +457 @@
           zkappa_i(ji,0)        = ztcond_i(ji,0)/MAX(zeps,zh_i(ji))
           zkappa_i(ji,nlay_i)   = ztcond_i(ji,nlay_i) / MAX(zeps,zh_i(ji))
-
        !-- Interface
           zkappa_s(ji,nlay_s)   = 2.0*rcdsn*ztcond_i(ji,0)/MAX(zeps, &
@@ -616 +583 @@
                    numeqmin(ji)    =  1
                    numeqmax(ji)    =  nlay_i + nlay_s + 1
+
                    !!surface equation
                    ztrid(ji,1,1) = 0.0
                    ztrid(ji,1,2) = dzf(ji) - zg1s*zkappa_s(ji,0)
                    ztrid(ji,1,3) = zg1s*zkappa_s(ji,0)
-                   ! bug again ???
-!                  zindterm(ji,1) = dzf(ji)*ztsuold(ji)   - zf(ji)
                    zindterm(ji,1) = dzf(ji)*t_su_b(ji)   - zf(ji)
 
@@ -853 +819 @@
        END DO
 
+       !-------------------------!
+       ! Heat conservation       !
+       !-------------------------!
        IF ( con_i ) THEN
-       !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-       ! ++++ heat conservation ++++ '
-       ! Internal fluxes ... goes to trash!!!
 
        DO ji = kideb, kiut
@@ -884 +850 @@
           zji                 = MOD( npb(ji) - 1, jpi ) + 1
           zjj                 = ( npb(ji) - 1 ) / jpi + 1
-!         IF ( ( zji.EQ.jiindex ) .AND. ( zjj.EQ.jjindex) ) THEN
-!            WRITE(numout,*) ' layer    : ', layer
-!            WRITE(numout,*) ' fc_i     : ', fc_i(ji,layer)
-!            WRITE(numout,*) ' zkappa_i : ', zkappa_i(ji,layer)
-!            WRITE(numout,*) ' ztcond_i : ', ztcond_i(ji,layer)
-!            WRITE(numout,*) ' zspeche_i: ', zspeche_i(ji,layer)
-!            WRITE(numout,*) ' ti2      : ', t_i_b(ji,layer+1)
-!            WRITE(numout,*) ' ti1      : ', t_i_b(ji,layer)
-!         ENDIF
        END DO
        END DO
@@ -903 +860 @@
 
        ENDIF
-
-       !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
     END SUBROUTINE lim_thd_dif