Changeset 1445 for trunk/NEMO/OFF_SRC/TRA/traqsr.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/OFF_SRC/TRA/traqsr.F90 (modified) (3 diffs)

trunk/NEMO/OFF_SRC/TRA/traqsr.F90

@@ -23 +23 @@
 
    !! * Shared module variables
-   LOGICAL, PUBLIC ::   ln_traqsr = .TRUE.   !: qsr flag (Default=T)
-
-   !! * Module variables
-   REAL(wp), PUBLIC ::  & !!! * penetrative solar radiation namelist *
-      xsi1 = 0.35_wp    ! first depth of extinction 
-   LOGICAL , PUBLIC ::   ln_qsr_sms = .false. ! flag to use or not the biological fluxes for light
+   LOGICAL , PUBLIC ::   ln_traqsr  = .TRUE.    !: light absorption (qsr) flag
+   LOGICAL , PUBLIC ::   ln_qsr_bio = .FALSE.   !: bio-optical retroaction
+   REAL(wp), PUBLIC ::   rn_abs     = 0.58_wp   !: fraction absorbed in the very near surface (RGB & 2 bands)
+   REAL(wp), PUBLIC ::   rn_si0     = 0.35_wp   !: very near surface depth of extinction      (RGB & 2 bands)
+   REAL(wp), PUBLIC ::   rn_si2     = 61.8_wp   !: deepest depth of extinction (blue & 0.01 mg.m-3)     (RGB)
 
 
@@ -55 +54 @@
       !!----------------------------------------------------------------------
       !! * Local declarations
-
-      NAMELIST/namqsr/ ln_traqsr, xsi1
+      NAMELIST/namqsr/ ln_qsr_bio, rn_abs, rn_si0, rn_si2
       !!----------------------------------------------------------------------
 
@@ -64 +62 @@
       READ   ( numnam, namqsr )
 
+      ln_qsr_bio = .FALSE. ! Offline mode : No retroaction on dynamics 
       ! Parameter control and print
       ! ---------------------------
-      IF( ln_traqsr  ) THEN
-        IF ( lwp ) THEN
+      IF(lwp) THEN                ! control print
          WRITE(numout,*)
          WRITE(numout,*) 'tra_qsr_init : penetration of the surface solar radiation'
          WRITE(numout,*) '~~~~~~~~~~~~'
          WRITE(numout,*) '    Namelist namqsr : set the parameter of penetration'
-         WRITE(numout,*) '        first depth of extinction        xsi1        = ',xsi1
-         WRITE(numout,*) '     Biological fluxes for light(Y/N) ln_qsr_sms  = ',ln_qsr_sms
-         WRITE(numout,*) ' '
-        END IF
-      ELSE
-        IF ( lwp ) THEN
-         WRITE(numout,*)
-         WRITE(numout,*) 'tra_qsr_init : NO solar flux penetration'
-         WRITE(numout,*) '~~~~~~~~~~~~'
-        END IF
+         WRITE(numout,*) '        bio-model            light penetration   ln_qsr_bio = ', ln_qsr_bio
+         WRITE(numout,*) '        RGB & 2 bands: fraction of light (rn_si1)    rn_abs = ', rn_abs
+         WRITE(numout,*) '        RGB & 2 bands: shortess depth of extinction  rn_si0 = ', rn_si0
+         WRITE(numout,*) '        3 bands: longest depth of extinction         rn_si2 = ', rn_si2
       ENDIF
-
-      IF( xsi1 < 0.e0  ) &
-         CALL ctl_stop( '              0<xsi1 not satisfied' )
-
+      
 
    END SUBROUTINE tra_qsr_init