- Timestamp:
- 2016-11-28T17:04:10+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90
r5385 r7351 8 8 9 9 !!---------------------------------------------------------------------- 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 16 17 17 18 IMPLICIT NONE … … 41 42 LOGICAL, PUBLIC, PARAMETER :: lk_qsr_bio = .FALSE. !: bio-model light absorption flag 42 43 #endif 43 44 44 #if defined key_offline 45 45 !!---------------------------------------------------------------------- … … 64 64 LOGICAL, PUBLIC, PARAMETER :: lk_degrad = .FALSE. !: degradation flag 65 65 #endif 66 67 !! * Substitutions68 # include "domzgr_substitute.h90"69 66 !!---------------------------------------------------------------------- 70 67 !! NEMO/OPA 3.3 , NEMO Consortium (2010) … … 104 101 !!---------------------------------------------------------------------- 105 102 REAL(wp), DIMENSION(3,61), INTENT(out) :: prgb ! tabulated attenuation coefficient 106 ! !103 ! 107 104 INTEGER :: jc ! dummy loop indice 108 105 INTEGER :: irgb ! temporary integer … … 188 185 zchl = zrgb(1,jc) 189 186 irgb = NINT( 41 + 20.* LOG10( zchl ) + 1.e-15 ) 190 IF(lwp ) WRITE(numout,*) ' jc =', jc, ' Chl = ', zchl, ' irgb = ', irgb187 IF(lwp .AND. nn_print >= 1 ) WRITE(numout,*) ' jc =', jc, ' Chl = ', zchl, ' irgb = ', irgb 191 188 IF( irgb /= jc ) THEN 192 189 IF(lwp) WRITE(numout,*) ' jc =', jc, ' Chl = ', zchl, ' Chl class = ', irgb … … 210 207 !!---------------------------------------------------------------------- 211 208 REAL(wp), DIMENSION(3,61), INTENT(out) :: prgb ! tabulated attenuation coefficient 212 ! !209 ! 213 210 INTEGER :: jc, jb ! dummy loop indice 214 211 INTEGER :: irgb ! temporary integer … … 262 259 REAL(wp), INTENT(in) :: prldex ! longest depth of extinction 263 260 REAL(wp), INTENT(in) :: pqsr_frc ! frac. solar radiation which penetrates 264 ! !261 ! 265 262 INTEGER :: jk, pjl ! levels 266 263 REAL(wp) :: zhext ! deepest level till which light penetrates … … 276 273 DO jk = jpkm1, 1, -1 277 274 IF(SUM(tmask(:,:,jk)) > 0 ) THEN 278 zem = MAXVAL( fsdepw(:,:,jk+1) * tmask(:,:,jk) )275 zem = MAXVAL( gdepw_0(:,:,jk+1) * tmask(:,:,jk) ) 279 276 IF( zem >= zhext ) pjl = jk ! last T-level reached by Qsr 280 277 ELSE
Note: See TracChangeset
for help on using the changeset viewer.