Changeset 2317 for branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

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/OPA_SRC/trc_oce.F90 (modified) (2 diffs)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

@@ -21 +21 @@
    PUBLIC   trc_oce_ext_lev    ! function called by traqsr.F90 at least
  
+   REAL(wp), PUBLIC                          ::   r_si2   !: largest depth of extinction (blue & 0.01 mg.m-3)  (RGB)
    REAL(wp), PUBLIC , DIMENSION(jpi,jpj,jpk) ::   etot3   !: light absortion coefficient
 
@@ -138 +139 @@
       prgb(:,:) = zrgb(2:4,:)
       !
+      r_si2 = 1.e0 / zrgb(2, 1)        ! blue with the smallest chlorophyll concentration)
+      IF(lwp) WRITE(numout,*) '      RGB longest depth of extinction    r_si2 = ', r_si2
+      !
       DO jc = 1, 61                         ! check
          zchl = zrgb(1,jc)