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 5845 for branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90 – NEMO

Ignore:
Timestamp:
2015-10-31T08:40:45+01:00 (8 years ago)
Author:
gm
Message:

#1613: vvl by default: suppression of domzgr_substitute.h90

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

    r5385 r5845  
    6565#endif 
    6666 
    67    !! * Substitutions 
    68 #  include "domzgr_substitute.h90" 
    6967   !!---------------------------------------------------------------------- 
    7068   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    … …  
    104102      !!---------------------------------------------------------------------- 
    105103      REAL(wp), DIMENSION(3,61), INTENT(out) ::   prgb   ! tabulated attenuation coefficient 
    106       !! 
     104      ! 
    107105      INTEGER  ::   jc     ! dummy loop indice 
    108106      INTEGER  ::   irgb   ! temporary integer 
    … …  
    188186         zchl = zrgb(1,jc) 
    189187         irgb = NINT( 41 + 20.* LOG10( zchl ) + 1.e-15 ) 
    190          IF(lwp) WRITE(numout,*) '    jc =', jc, '  Chl = ', zchl, '  irgb = ', irgb 
     188         IF(lwp .AND. nn_print >= 1 ) WRITE(numout,*) '    jc =', jc, '  Chl = ', zchl, '  irgb = ', irgb 
    191189         IF( irgb /= jc ) THEN 
    192190            IF(lwp) WRITE(numout,*) '    jc =', jc, '  Chl = ', zchl, '  Chl class = ', irgb 
    … …  
    210208      !!---------------------------------------------------------------------- 
    211209      REAL(wp), DIMENSION(3,61), INTENT(out) ::   prgb   ! tabulated attenuation coefficient 
    212       !! 
     210      ! 
    213211      INTEGER  ::   jc, jb ! dummy loop indice 
    214212      INTEGER  ::   irgb   ! temporary integer 
    … …  
    262260      REAL(wp), INTENT(in) ::   prldex    ! longest depth of extinction 
    263261      REAL(wp), INTENT(in) ::   pqsr_frc  ! frac. solar radiation which penetrates  
    264       !! 
     262      ! 
    265263      INTEGER  ::   jk, pjl            ! levels 
    266264      REAL(wp) ::   zhext              ! deepest level till which light penetrates 
    … …  
    276274      DO jk = jpkm1, 1, -1 
    277275         IF(SUM(tmask(:,:,jk)) > 0 ) THEN 
    278             zem = MAXVAL( fsdepw(:,:,jk+1) * tmask(:,:,jk) ) 
     276            zem = MAXVAL( gdepw_0(:,:,jk+1) * tmask(:,:,jk) ) 
    279277            IF( zem >= zhext )   pjl = jk                       ! last T-level reached by Qsr 
    280278         ELSE 
Note: See TracChangeset for help on using the changeset viewer.