Changeset 3982 for branches/2013/dev_r3406_CNRS_LIM3

Timestamp:

2013-07-18T12:44:32+02:00 (11 years ago)

Author:

clem

Message:

bug correction for penetrative solar flux

Location:

branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3

Files:

• limsbc.F90 (modified) (3 diffs)
• limthd.F90 (modified) (1 diff)

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

@@ -109 +109 @@
       INTEGER  ::   ifvt, i1mfr, idfr               ! some switches
       INTEGER  ::   iflt, ial, iadv, ifral, ifrdv
-      REAL(wp) ::   zinda, zindb, zfons, zpme              ! local scalars
+      REAL(wp) ::   zinda, zfons, zpme              ! local scalars
       REAL(wp) ::   zfmm             ! IOVINO freezing minus melting (F-M) 
       REAL(wp), POINTER, DIMENSION(:,:) ::   zfcm1 , zfcm2    ! solar/non solar heat fluxes
       REAL(wp), POINTER, DIMENSION(:,:,:) ::   zalb, zalbp   ! 2D/3D workspace
+      REAL(wp) ::   zzfcm1, zfscmbq ! clem: for light penetration
       !!---------------------------------------------------------------------
       
@@ -127 +128 @@
          DO ji = 1, jpi
             zinda   = 1.0 - MAX( rzero , SIGN( rone , - ( 1.0 - pfrld(ji,jj) ) ) )
-            zindb   = 1.0 - MAX( rzero , SIGN( rone , - iatte(ji,jj) ) )
             ifvt    = zinda  *  MAX( rzero , SIGN( rone, - phicif(ji,jj) ) )  !subscripts are bad here
             i1mfr   = 1.0 - MAX( rzero , SIGN( rone ,  - at_i(ji,jj) ) )
@@ -151 +151 @@
             !   computation the solar flux at ocean surface
             zfcm1(ji,jj)   = pfrld(ji,jj) * qsr(ji,jj)  + &
-                 &           zindb * ( 1. - pfrld(ji,jj) ) * fstric(ji,jj) / MAX( iatte(ji,jj), epsi20 )
+                 &           ( 1. - pfrld(ji,jj) ) * fstric(ji,jj) / ( 1.0 - zinda + zinda * iatte(ji,jj) )
             ! fstric     Solar flux transmitted trough the ice
             ! qsr        Net short wave heat flux on free ocean
             ! new line
-            fscmbq(ji,jj) = zindb * ( 1.0 - pfrld(ji,jj) ) * fstric(ji,jj) / MAX( iatte(ji,jj), epsi20 )
+            fscmbq(ji,jj) = ( 1.0 - pfrld(ji,jj) ) * fstric(ji,jj) / ( 1.0 - zinda + zinda * iatte(ji,jj) )
+
+            ! solar flux and fscmbq with light penetration (clem)
+            zzfcm1  = pfrld(ji,jj) * qsr(ji,jj) * oatte(ji,jj) + ( 1. - pfrld(ji,jj) ) * fstric(ji,jj)
+            zfscmbq = ( 1.0 - pfrld(ji,jj) ) * fstric(ji,jj)
 
             !  computation the non solar heat flux at ocean surface
-            zfcm2(ji,jj) = - zfcm1(ji,jj)                  &
-               &           + iflt    * ( fscmbq(ji,jj) )   & ! total abl -> fscmbq is given to the ocean
-               ! fscmbq and ffltbif are obsolete
-               !              &           + iflt * ffltbif(ji,jj) !!! only if one category is used
+            zfcm2(ji,jj) = - zzfcm1                  &
+               &           + iflt    * zfscmbq   & ! total abl -> fscmbq is given to the ocean
                &           + ifral   * ( ial * qcmif(ji,jj) + (1 - ial) * qldif(ji,jj) ) * r1_rdtice   &
                &           + ifrdv   * ( qfvbq(ji,jj) + qdtcn(ji,jj) )                   * r1_rdtice   &

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

@@ -185 +185 @@
             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)
+            !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"