Changeset 2317 for branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/PISCES

Timestamp:

2010-10-26T18:01:34+02:00 (14 years ago)

Author:

cetlod

Message:

Improve the computation of the divergence of the downward solar irradiance in traqsr.F90, see ticket #726

• change fsdepw & fse3t into fsdepw_0 & fse3t_0 in tra_qsr_init
• modify tra_qsr to recompute systematically etot3 in vvl case
• suppress the namelist parameter rn_si2 which is computed as the RGB longest depth of extinction in trc_oce.F90

File:

• branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zopt.F90 (modified) (3 diffs)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zopt.F90

@@ -28 +28 @@
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   emoy                 !: averaged PAR in the mixed layer
 
-   INTEGER  ::   nksrp   ! levels below which the light cannot penetrate ( depth larger than 391 m)
-   REAL(wp) ::   &
-      parlux = 0.43 / 3.e0
+   INTEGER  ::  nksrp   ! levels below which the light cannot penetrate ( depth larger than 391 m)
+   REAL(wp) ::  parlux = 0.43 / 3.e0
 
    REAL(wp), DIMENSION(3,61), PUBLIC ::   xkrgb  !: tabulated attenuation coefficients for RGB absorption
@@ -54 +53 @@
       !!---------------------------------------------------------------------
       INTEGER, INTENT(in) ::   kt, jnt ! ocean time step
-      INTEGER  ::   ji, jj, jk, jc
+      INTEGER  ::   ji, jj, jk
       INTEGER  ::   irgb
       REAL(wp) ::   zchl, zxsi0r
@@ -235 +234 @@
       !!----------------------------------------------------------------------
 
-      !                                ! level of light extinction
-      nksrp = trc_oce_ext_lev( rn_si2, 0.33e2 )
-      IF(lwp) THEN
-        WRITE(numout,*)
-        WRITE(numout,*) ' level max of computation of qsr = ', nksrp, ' ref depth = ', gdepw_0(nksrp+1), ' m'
-      ENDIF
-!!      CALL trc_oce_rgb( xkrgb )     ! tabulated attenuation coefficients
-      CALL trc_oce_rgb_read( xkrgb )     ! tabulated attenuation coefficients
+!!      CALL trc_oce_rgb( xkrgb )                  ! tabulated attenuation coefficients
+      CALL trc_oce_rgb_read( xkrgb )               ! tabulated attenuation coefficients
+      nksrp = trc_oce_ext_lev( r_si2, 0.33e2 )     ! max level of light extinction (Blue Chl=0.01)
+      IF(lwp) WRITE(numout,*) '        level of light extinction = ', nksrp, ' ref depth = ', gdepw_0(nksrp+1), ' m'
+      !
       etot (:,:,:) = 0.e0
       enano(:,:,:) = 0.e0
       ediat(:,:,:) = 0.e0
-      IF( ln_qsr_bio ) etot3(:,:,:) = 0.e0
+      IF( ln_qsr_bio )   etot3(:,:,:) = 0.e0
+
       ! 
    END SUBROUTINE p4z_opt_init