Changeset 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zopt.F90

Timestamp:

2010-12-27T18:33:53+01:00 (13 years ago)

Author:

rblod

Message:

Update NEMOGCM from branch nemo_v3_3_beta

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zopt.F90 (modified) (7 diffs, 1 prop)

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zopt.F90

@@ -16 +16 @@
    USE trc            ! tracer variables
    USE oce_trc        ! tracer-ocean share variables
-   USE trc_oce        ! ocean-tracer share variables
    USE sms_pisces     ! Source Minus Sink of PISCES
    USE iom
@@ -23 +22 @@
    PRIVATE
 
-   PUBLIC   p4z_opt   ! called in p4zbio.F90 module
+   PUBLIC   p4z_opt        ! called in p4zbio.F90 module
+   PUBLIC   p4z_opt_init   ! called in trcsms_pisces.F90 module
 
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   etot, enano, ediat   !: PAR for phyto, nano and diat 
    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
@@ -37 +36 @@
 #  include "top_substitute.h90"
    !!----------------------------------------------------------------------
-   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 
+   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
    !! $Id$ 
-   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 
 CONTAINS
 
-   SUBROUTINE p4z_opt(kt, jnt)
+   SUBROUTINE p4z_opt( kt, jnt )
       !!---------------------------------------------------------------------
       !!                     ***  ROUTINE p4z_opt  ***
@@ -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
@@ -64 +63 @@
 
 
-      !                                        !* tabulated attenuation coef. 
-      IF( kt * jnt == nittrc000 ) THEN
-         !                                ! 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
-         etot (:,:,:) = 0.e0
-         enano(:,:,:) = 0.e0
-         ediat(:,:,:) = 0.e0
-         IF( ln_qsr_bio ) etot3(:,:,:) = 0.e0
-      ENDIF
-
-
-!     Initialisation of variables used to compute PAR
-!     -----------------------------------------------
+      !     Initialisation of variables used to compute PAR
+      !     -----------------------------------------------
       ze1 (:,:,jpk) = 0.e0
       ze2 (:,:,jpk) = 0.e0
@@ -227 +209 @@
       END DO
 
-#if defined key_trc_diaadd
+#if defined key_diatrc
 # if ! defined key_iomput
       ! save for outputs
@@ -243 +225 @@
    END SUBROUTINE p4z_opt
 
+   SUBROUTINE p4z_opt_init
+      !!----------------------------------------------------------------------
+      !!                  ***  ROUTINE p4z_opt_init  ***
+      !!
+      !! ** Purpose :   Initialization of tabulated attenuation coef
+      !!
+      !!
+      !!----------------------------------------------------------------------
+
+      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
+      ! 
+   END SUBROUTINE p4z_opt_init
 #else
    !!----------------------------------------------------------------------