Changeset 4672 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90

Timestamp:

2014-06-17T17:06:59+02:00 (10 years ago)

Author:

clem

Message:

 

File:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90 (modified) (12 diffs)

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

@@ -43 +43 @@
 CONTAINS
 
-   SUBROUTINE lim_thd_dh( kideb, kiut, jl )
+   SUBROUTINE lim_thd_dh( kideb, kiut )
       !!------------------------------------------------------------------
       !!                ***  ROUTINE lim_thd_dh  ***
@@ -68 +68 @@
       !!------------------------------------------------------------------
       INTEGER , INTENT(in) ::   kideb, kiut   ! Start/End point on which the  the computation is applied
-      INTEGER , INTENT(in) ::   jl            ! Thickness cateogry number
       !! 
       INTEGER  ::   ji , jk        ! dummy loop indices
@@ -151 +150 @@
       zintermelt(:)   = 0._wp
       icount    (:)   = 0
-
-      ! debug
-      dq_i(:) = 0._wp
-      dq_s(:) = 0._wp
 
       ! initialize layer thicknesses and enthalpies
@@ -273 +268 @@
          zh_s  (ji) = ht_s_b(ji) / REAL( nlay_s )
 
-         ! clem debug: variation of enthalpy (J.m-2)
-         dq_s(ji) = dq_s(ji) + ( zdh_s_pre(ji) + zdh_s_mel(ji) ) * zqprec(ji)  
          ENDIF
       END DO
@@ -297 +290 @@
             ht_s_b(ji) = MAX( 0._wp , ht_s_b(ji) + zdeltah(ji,jk) )
 
-            ! clem debug: variation of enthalpy (J.m-2)
-            dq_s(ji) = dq_s(ji) + zdeltah(ji,jk) * q_s_b(ji,jk)  
          END DO
       END DO
@@ -325 +316 @@
             ! new snow thickness
             ht_s_b(ji)     =  MAX( 0._wp , ht_s_b(ji) + zdh_s_sub(ji) )
-            ! clem debug: variation of enthalpy (J.m-2)
-            dq_s(ji) = dq_s(ji) + zdh_s_sub(ji) * q_s_b(ji,1)  
          END DO
       ENDIF
@@ -397 +386 @@
             icount(ji)  = icount(ji) + zindh
             zh_i(ji,jk) = MAX( 0._wp , zh_i(ji,jk) + zdeltah(ji,jk) )
-
-            ! clem debug: variation of enthalpy (J.m-2)
-            dq_i(ji) = dq_i(ji) + zdeltah(ji,jk) * q_i_b(ji,jk)  
 
             ! update heat content (J.m-2) and layer thickness
@@ -513 +499 @@
             wfx_bog_1d(ji) =  wfx_bog_1d(ji) - rhoic * a_i_b(ji) * dh_i_bott(ji) * r1_rdtice
 
-            ! clem debug: variation of enthalpy (J.m-2)
-            dq_i(ji) = dq_i(ji) + dh_i_bott(ji) * q_i_b(ji,nlay_i+1)  
-
             ! update heat content (J.m-2) and layer thickness
             qh_i_old(ji,nlay_i+1) = qh_i_old(ji,nlay_i+1) + dh_i_bott(ji) * q_i_b(ji,nlay_i+1)
@@ -557 +540 @@
                   ! Contribution to mass flux
                   wfx_res_1d(ji) =  wfx_res_1d(ji) - rhoic * a_i_b(ji) * zdeltah(ji,jk) * r1_rdtice
-
-                  ! clem debug: variation of enthalpy (J.m-2)
-                  dq_i(ji) = dq_i(ji) + zdeltah(ji,jk) * q_i_b(ji,jk)  
 
                   ! update heat content (J.m-2) and layer thickness
@@ -598 +578 @@
                   ! Contribution to mass flux
                   wfx_bom_1d(ji) =  wfx_bom_1d(ji) - rhoic * a_i_b(ji) * zdeltah(ji,jk) * r1_rdtice
-
-                  ! clem debug: variation of enthalpy (J.m-2)
-                  dq_i(ji) = dq_i(ji) + zdeltah(ji,jk) * q_i_b(ji,jk)  
 
                   ! update heat content (J.m-2) and layer thickness
@@ -664 +641 @@
 !         ! Contribution to mass flux
 !         wfx_snw_1d(ji)  =  wfx_snw_1d(ji) - rhosn * a_i_b(ji) * zdeltah(ji,1) * r1_rdtice
-!         ! clem debug: variation of enthalpy (J.m-2)
-!         dq_s(ji) = dq_s(ji) + zdeltah(ji,1) * q_s_b(ji,1)  
 !    
          ii = MOD( npb(ji) - 1, jpi ) + 1 ; ij = ( npb(ji) - 1 ) / jpi + 1
@@ -721 +696 @@
          wfx_snw_1d(ji) = wfx_snw_1d(ji) + a_i_b(ji) * dh_snowice(ji) * rhosn * r1_rdtice
 
-         ! clem debug: variation of enthalpy (J.m-2)
-         dq_s(ji) = dq_s(ji) - dh_snowice(ji) * q_s_b(ji,1)  
-         dq_i(ji) = dq_i(ji) + dh_snowice(ji) * q_s_b(ji,1) + zfmdt * zEw  
-
          ! update heat content (J.m-2) and layer thickness
          qh_i_old(ji,0) = qh_i_old(ji,0) + dh_snowice(ji) * q_s_b(ji,1) + zfmdt * zEw