Changeset 1445 for trunk/NEMO/OPA_SRC/trc_oce.F90

Timestamp:

2009-05-13T16:35:02+02:00 (15 years ago)

Author:

cetlod

Message:

add the use of bio-optical retroaction on dynamics when coupling with PISCES, see ticket:428

File:

• trunk/NEMO/OPA_SRC/trc_oce.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/trc_oce.F90

@@ -10 +10 @@
    !!   trc_oce_rgb : tabulated attenuation coefficients for RGB light penetration         
    !!----------------------------------------------------------------------
+   USE par_oce
    USE in_out_manager  ! I/O manager
    USE dom_oce         ! ocean space and time domain
 
-#if defined key_top && defined key_pisces
-   !!----------------------------------------------------------------------
-   !!   'key_top'   &   'key_pisces'                       PISCES bio-model          
-   !!----------------------------------------------------------------------
-   USE sms_pisces , ONLY :   etot3    =>   etot3   !:  bio-model light absorption
-
    IMPLICIT NONE
    PRIVATE
 
-   PUBLIC   trc_oce_rgb   ! routine called by p4zopt.F90
-   
-   LOGICAL, PUBLIC, PARAMETER ::   lk_qsr_bio = .TRUE.   !: bio-model light absorption flag
-   
-#else
-   !!----------------------------------------------------------------------
-   !! Default option                          No bio-model light absorption      
-   !!----------------------------------------------------------------------
-   USE par_oce
-
-   IMPLICIT NONE
-   PRIVATE
-
-   PUBLIC   trc_oce_rgb   ! routine called by traqsr.F90
-   
-   LOGICAL, PUBLIC, PARAMETER ::   lk_qsr_bio = .FALSE.   !: bio-model light absorption flag
+   PUBLIC   trc_oce_rgb        ! routine called by traqsr.F90
+   PUBLIC   trc_oce_rgb_read   ! routine called by traqsr.F90
+   PUBLIC   trc_oce_ext_lev    ! function called by traqsr.F90 at least
    
    REAL(wp), PUBLIC , DIMENSION(jpi,jpj,jpk) ::   etot3   !: light absortion coefficient
+
+#if defined key_top && defined key_pisces
+   !!----------------------------------------------------------------------
+   !!   'key_top'   &   'key_pisces'                       PISCES bio-model          
+   !!----------------------------------------------------------------------
+   LOGICAL, PUBLIC, PARAMETER ::   lk_qsr_bio = .TRUE.   !: bio-model light absorption flag
+#else
+   !!----------------------------------------------------------------------
+   !! Default option                          No bio-model light absorption      
+   !!----------------------------------------------------------------------
+   LOGICAL, PUBLIC, PARAMETER ::   lk_qsr_bio = .FALSE.   !: bio-model light absorption flag
 #endif
-
-   PUBLIC   trc_oce_ext_lev    ! function called by traqsr.F90 at least
-
-   INTEGER, PUBLIC ::   nksr   ! levels below which the light cannot penetrate ( depth larger than 391 m)
-
-   REAL(wp), DIMENSION(3,61), PUBLIC ::   rkrgb   !: tabulated attenuation coefficients for RGB absorption
 
    !! * Substitutions
@@ -52 +39 @@
    !!----------------------------------------------------------------------
    !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) 
-   !! $Id:$ 
+   !! $Id$ 
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
    !!----------------------------------------------------------------------
@@ -72 +59 @@
       !! Reference  : Lengaigne et al. 2007, Clim. Dyn., V28, 5, 503-516.
       !!----------------------------------------------------------------------
-      REAL(wp), DIMENSION(3,61), INTENT(inout) ::   prgb   ! tabulated attenuation coefficient
+      REAL(wp), DIMENSION(3,61), INTENT(out) ::   prgb   ! tabulated attenuation coefficient
       !!
       INTEGER  ::   jc     ! dummy loop indice
@@ -175 +162 @@
       !!                          attenuation coefficient (from JM Andre)
       !!----------------------------------------------------------------------
-      REAL(wp), DIMENSION(3,61), INTENT(inout) ::   prgb   ! tabulated attenuation coefficient
+      REAL(wp), DIMENSION(3,61), INTENT(out) ::   prgb   ! tabulated attenuation coefficient
       !!
       INTEGER  ::   jchl, jband   ! dummy loop indices
       INTEGER  ::   numlight
-      REAL(wp) ::   ztoto
+      REAL(wp) ::   zchl
       CHARACTER(LEN=20) :: clname
       !!----------------------------------------------------------------------
@@ -186 +173 @@
       CALL ctlopn( numlight, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', 1, numout, .TRUE., 1 )
       DO jchl = 1, 61
-         READ(numlight,*) ztoto, ( prgb(jband,jchl), jband=1,3 )
+         READ(numlight,*) zchl, ( prgb(jband,jchl), jband=1,3 )
       END DO
       CLOSE( numlight )