Changeset 5206 for branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

Timestamp:

2015-04-13T12:13:10+02:00 (9 years ago)

Author:

cetlod

Message:

First improvment of diurnal cycle for passive tracers, see ticket #1484

File:

• branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90 (modified) (4 diffs)

branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

@@ -126 +126 @@
       REAL(wp), POINTER, DIMENSION(:)     ::   zdifcase    ! case of the equation resolution (1->4)
       REAL(wp), POINTER, DIMENSION(:)     ::   zftrice     ! solar radiation transmitted through the ice
+      REAL(wp), POINTER, DIMENSION(:)     ::   zftrice_mean ! daily mean solar radiation transmitted through the ice
       REAL(wp), POINTER, DIMENSION(:)     ::   zihic
       
@@ -174 +175 @@
       CALL wrk_alloc( jpij,nlay_i+3, zindterm, zindtbis, zdiagbis  )
       CALL wrk_alloc( jpij,nlay_i+3,3, ztrid )
+      IF( l_trcdm2dc )  CALL wrk_alloc( jpij, zftrice_mean )
 
       CALL wrk_alloc( jpij, zdq, zq_ini, zhfx_err )
@@ -249 +251 @@
       END DO
 
+      IF( l_trcdm2dc ) THEN
+         DO ji = kideb , kiut
+            zftrice_mean(ji) =  qsr_ice_mean_1d(ji) * i0(ji) ! Solar radiation transmitted below the surface layer
+         END DO
+      ENDIF
+
       !---------------------------------------------------------
       ! Transmission - absorption of solar radiation in the ice
@@ -282 +290 @@
          ftr_ice_1d(ji) = zradtr_i(ji,nlay_i) 
       END DO
+
+
+      IF( l_trcdm2dc ) THEN
+         DO ji = kideb , kiut
+             ftr_ice_mean_1d(ji) =  ftr_ice_mean_1d(ji) & 
+                 &                + a_i_1d(ji) * zftrice_mean(ji) &
+                 &                             * EXP( - rn_kappa_i * ( MAX ( 0._wp , h_i_1d(ji) ) ) ) &
+                                               * EXP( - zraext_s   * ( MAX ( 0._wp , h_s_1d(ji) ) ) )
+         END DO
+      ENDIF
 
       !