Changeset 3979

Timestamp:

2013-07-16T16:45:31+02:00 (11 years ago)

Author:

clem

Message:

bug correction on solar flux penetration

File:

• branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90 (modified) (1 diff)

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

@@ -184 +184 @@
             pfrld(ji,jj)   = 1.0 - at_i(ji,jj)
             zinda          = tms(ji,jj) * (1.0 - MAX( zzero , SIGN( zone , - at_i(ji,jj) ) ) )
-            !
+            ! clem: set up ocean and ice attenuation to 1 (i.e. no attenuation) when there is no ice (zinda = 0)
+            iatte(ji,jj) = (1.0 - zinda ) + zinda * iatte(ji,jj)
+            oatte(ji,jj) = (1.0 - zinda ) + zinda * oatte(ji,jj)
             !           !  solar irradiance transmission at the mixed layer bottom and used in the lead heat budget
             !           !  practically no "direct lateral ablation"