New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7351 for branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90 – NEMO

Ignore:
Timestamp:
2016-11-28T17:04:10+01:00 (7 years ago)
Author:
emanuelaclementi
Message:

ticket #1805 step 3: /2016/dev_INGV_UKMO_2016 aligned to the trunk at revision 7161

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

    r5385 r7351  
    88 
    99   !!---------------------------------------------------------------------- 
    10    !!   trc_oce_rgb : tabulated attenuation coefficients for RGB light penetration          
    11    !!---------------------------------------------------------------------- 
    12    USE par_oce 
    13    USE in_out_manager  ! I/O manager 
    14    USE dom_oce         ! ocean space and time domain 
    15    USE lib_mpp         ! MPP library 
     10   !!   trc_oce_rgb   : tabulated attenuation coefficients for RGB light penetration          
     11   !!---------------------------------------------------------------------- 
     12   USE par_oce        ! ocean parameters 
     13   USE dom_oce        ! ocean space and time domain 
     14   ! 
     15   USE in_out_manager ! I/O manager 
     16   USE lib_mpp        ! MPP library 
    1617 
    1718   IMPLICIT NONE 
     
    4142   LOGICAL, PUBLIC, PARAMETER ::   lk_qsr_bio = .FALSE.   !: bio-model light absorption flag 
    4243#endif 
    43  
    4444#if defined key_offline 
    4545   !!---------------------------------------------------------------------- 
     
    6464   LOGICAL, PUBLIC, PARAMETER ::   lk_degrad = .FALSE.   !: degradation flag 
    6565#endif 
    66  
    67    !! * Substitutions 
    68 #  include "domzgr_substitute.h90" 
    6966   !!---------------------------------------------------------------------- 
    7067   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
     
    104101      !!---------------------------------------------------------------------- 
    105102      REAL(wp), DIMENSION(3,61), INTENT(out) ::   prgb   ! tabulated attenuation coefficient 
    106       !! 
     103      ! 
    107104      INTEGER  ::   jc     ! dummy loop indice 
    108105      INTEGER  ::   irgb   ! temporary integer 
     
    188185         zchl = zrgb(1,jc) 
    189186         irgb = NINT( 41 + 20.* LOG10( zchl ) + 1.e-15 ) 
    190          IF(lwp) WRITE(numout,*) '    jc =', jc, '  Chl = ', zchl, '  irgb = ', irgb 
     187         IF(lwp .AND. nn_print >= 1 ) WRITE(numout,*) '    jc =', jc, '  Chl = ', zchl, '  irgb = ', irgb 
    191188         IF( irgb /= jc ) THEN 
    192189            IF(lwp) WRITE(numout,*) '    jc =', jc, '  Chl = ', zchl, '  Chl class = ', irgb 
     
    210207      !!---------------------------------------------------------------------- 
    211208      REAL(wp), DIMENSION(3,61), INTENT(out) ::   prgb   ! tabulated attenuation coefficient 
    212       !! 
     209      ! 
    213210      INTEGER  ::   jc, jb ! dummy loop indice 
    214211      INTEGER  ::   irgb   ! temporary integer 
     
    262259      REAL(wp), INTENT(in) ::   prldex    ! longest depth of extinction 
    263260      REAL(wp), INTENT(in) ::   pqsr_frc  ! frac. solar radiation which penetrates  
    264       !! 
     261      ! 
    265262      INTEGER  ::   jk, pjl            ! levels 
    266263      REAL(wp) ::   zhext              ! deepest level till which light penetrates 
     
    276273      DO jk = jpkm1, 1, -1 
    277274         IF(SUM(tmask(:,:,jk)) > 0 ) THEN 
    278             zem = MAXVAL( fsdepw(:,:,jk+1) * tmask(:,:,jk) ) 
     275            zem = MAXVAL( gdepw_0(:,:,jk+1) * tmask(:,:,jk) ) 
    279276            IF( zem >= zhext )   pjl = jk                       ! last T-level reached by Qsr 
    280277         ELSE 
Note: See TracChangeset for help on using the changeset viewer.