Changeset 4979 for branches/2014/dev_4728_CNRS04_coupled_interface

Timestamp:

2014-12-08T17:39:09+01:00 (9 years ago)

Author:

vancop

Message:

bugfixes found during the merge party

Location:

branches/2014/dev_4728_CNRS04_coupled_interface/NEMOGCM/NEMO/LIM_SRC_3

Files:

• limthd_dif.F90 (modified) (1 diff)
• limtrp.F90 (modified) (4 diffs)

branches/2014/dev_4728_CNRS04_coupled_interface/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

@@ -108 +108 @@
       REAL(wp) ::   zgamma    =  18009._wp    ! for specific heat
       REAL(wp) ::   zbeta     =  0.117_wp     ! for thermal conductivity (could be 0.13)
-      REAL(wp) ::   zraext_s  =  1.e+8_wp     ! extinction coefficient of radiation in the snow
+      REAL(wp) ::   zraext_s  =  10.0_wp      ! extinction coefficient of radiation in the snow
       REAL(wp) ::   zkimin    =  0.10_wp      ! minimum ice thermal conductivity
       REAL(wp) ::   ztsu_err  =  1.e-5_wp     ! range around which t_su is considered as 0°C 

branches/2014/dev_4728_CNRS04_coupled_interface/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90

@@ -144 +144 @@
             zs0oi (:,:,jl)   = oa_i (:,:,jl) * area(:,:)    ! Age content
             zs0c0 (:,:,jl)   = e_s  (:,:,1,jl)              ! Snow heat content
-            zs0e  (:,:,:,jl) = e_i  (:,:,:,jl)              ! Ice  heat content
+            zs0e  (:,:,:,jl) = e_i  (:,:,1:nlay_i,jl)       ! Ice  heat content
          END DO
 
@@ -393 +393 @@
                      zsmv = smv_i(ji,jj,jl)
                      zes  = e_s  (ji,jj,1,jl)
-                     zei  = SUM( e_i(ji,jj,:,jl) )
+                     zei  = SUM( e_i(ji,jj,1:nlay_i,jl) )
                      zdv  = v_i(ji,jj,jl) - zviold(ji,jj,jl)   
                      !zda = a_i(ji,jj,jl) - zaiold(ji,jj,jl)   
@@ -414 +414 @@
                      smv_i(ji,jj,jl) = zindh * smv_i(ji,jj,jl)
                      e_s(ji,jj,1,jl) = zindh * e_s(ji,jj,1,jl)
-                     e_i(ji,jj,:,jl) = zindh * e_i(ji,jj,:,jl)
+                     e_i(ji,jj,1:nlay_i,jl) = zindh * e_i(ji,jj,1:nlay_i,jl)
 
                      ! Update mass fluxes
@@ -421 +421 @@
                      sfx_res(ji,jj) = sfx_res(ji,jj) - ( smv_i(ji,jj,jl) - zsmv ) * rhoic * r1_rdtice 
                      hfx_res(ji,jj) = hfx_res(ji,jj) + ( e_s(ji,jj,1,jl) - zes ) * unit_fac / area(ji,jj) * r1_rdtice ! W.m-2 <0
-                     hfx_res(ji,jj) = hfx_res(ji,jj) + ( SUM( e_i(ji,jj,:,jl) ) - zei ) * unit_fac / area(ji,jj) * r1_rdtice ! W.m-2 <0
+                     hfx_res(ji,jj) = hfx_res(ji,jj) + ( SUM( e_i(ji,jj,1:nlay_i,jl) ) - zei ) * unit_fac / area(ji,jj) * r1_rdtice ! W.m-2 <0
 
                   ENDIF