Changeset 1445 for trunk/NEMO/OPA_SRC
- Timestamp:
- 2009-05-13T16:35:02+02:00 (15 years ago)
- Location:
- trunk/NEMO/OPA_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRA/traqsr.F90
r1425 r1445 43 43 REAL(wp), PUBLIC :: rn_si2 = 61.8_wp !: deepest depth of extinction (blue & 0.01 mg.m-3) (RGB) 44 44 45 ! Module variables 45 46 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_chl ! structure of input Chl (file informations, fields read) 47 INTEGER :: nksr ! levels below which the light cannot penetrate ( depth larger than 391 m) 48 REAL(wp), DIMENSION(3,61) :: rkrgb !: tabulated attenuation coefficients for RGB absorption 46 49 47 50 !! * Substitutions … … 111 114 112 115 ! ! ============================================== ! 113 IF( lk_qsr_bio ) THEN! bio-model fluxes : all vertical coordinates !116 IF( lk_qsr_bio .AND. ln_qsr_bio ) THEN ! bio-model fluxes : all vertical coordinates ! 114 117 ! ! ============================================== ! 115 118 DO jk = 1, jpkm1 … … 266 269 WRITE(numout,*) ' 3 bands: longest depth of extinction rn_si2 = ', rn_si2 267 270 ENDIF 268 ! ! control consistency269 IF( lk_qsr_bio .AND. .NOT.ln_qsr_bio ) THEN270 ln_qsr_bio = .true.271 CALL ctl_warn( 'Force bio-model light penetraton ln_qsr_bio = TRUE ' )272 ENDIF273 271 274 272 ! ! ===================================== ! -
trunk/NEMO/OPA_SRC/trc_oce.F90
r1423 r1445 10 10 !! trc_oce_rgb : tabulated attenuation coefficients for RGB light penetration 11 11 !!---------------------------------------------------------------------- 12 USE par_oce 12 13 USE in_out_manager ! I/O manager 13 14 USE dom_oce ! ocean space and time domain 14 15 15 #if defined key_top && defined key_pisces16 !!----------------------------------------------------------------------17 !! 'key_top' & 'key_pisces' PISCES bio-model18 !!----------------------------------------------------------------------19 USE sms_pisces , ONLY : etot3 => etot3 !: bio-model light absorption20 21 16 IMPLICIT NONE 22 17 PRIVATE 23 18 24 PUBLIC trc_oce_rgb ! routine called by p4zopt.F90 25 26 LOGICAL, PUBLIC, PARAMETER :: lk_qsr_bio = .TRUE. !: bio-model light absorption flag 27 28 #else 29 !!---------------------------------------------------------------------- 30 !! Default option No bio-model light absorption 31 !!---------------------------------------------------------------------- 32 USE par_oce 33 34 IMPLICIT NONE 35 PRIVATE 36 37 PUBLIC trc_oce_rgb ! routine called by traqsr.F90 38 39 LOGICAL, PUBLIC, PARAMETER :: lk_qsr_bio = .FALSE. !: bio-model light absorption flag 19 PUBLIC trc_oce_rgb ! routine called by traqsr.F90 20 PUBLIC trc_oce_rgb_read ! routine called by traqsr.F90 21 PUBLIC trc_oce_ext_lev ! function called by traqsr.F90 at least 40 22 41 23 REAL(wp), PUBLIC , DIMENSION(jpi,jpj,jpk) :: etot3 !: light absortion coefficient 24 25 #if defined key_top && defined key_pisces 26 !!---------------------------------------------------------------------- 27 !! 'key_top' & 'key_pisces' PISCES bio-model 28 !!---------------------------------------------------------------------- 29 LOGICAL, PUBLIC, PARAMETER :: lk_qsr_bio = .TRUE. !: bio-model light absorption flag 30 #else 31 !!---------------------------------------------------------------------- 32 !! Default option No bio-model light absorption 33 !!---------------------------------------------------------------------- 34 LOGICAL, PUBLIC, PARAMETER :: lk_qsr_bio = .FALSE. !: bio-model light absorption flag 42 35 #endif 43 44 PUBLIC trc_oce_ext_lev ! function called by traqsr.F90 at least45 46 INTEGER, PUBLIC :: nksr ! levels below which the light cannot penetrate ( depth larger than 391 m)47 48 REAL(wp), DIMENSION(3,61), PUBLIC :: rkrgb !: tabulated attenuation coefficients for RGB absorption49 36 50 37 !! * Substitutions … … 52 39 !!---------------------------------------------------------------------- 53 40 !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) 54 !! $Id :$41 !! $Id$ 55 42 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 56 43 !!---------------------------------------------------------------------- … … 72 59 !! Reference : Lengaigne et al. 2007, Clim. Dyn., V28, 5, 503-516. 73 60 !!---------------------------------------------------------------------- 74 REAL(wp), DIMENSION(3,61), INTENT( inout) :: prgb ! tabulated attenuation coefficient61 REAL(wp), DIMENSION(3,61), INTENT(out) :: prgb ! tabulated attenuation coefficient 75 62 !! 76 63 INTEGER :: jc ! dummy loop indice … … 175 162 !! attenuation coefficient (from JM Andre) 176 163 !!---------------------------------------------------------------------- 177 REAL(wp), DIMENSION(3,61), INTENT( inout) :: prgb ! tabulated attenuation coefficient164 REAL(wp), DIMENSION(3,61), INTENT(out) :: prgb ! tabulated attenuation coefficient 178 165 !! 179 166 INTEGER :: jchl, jband ! dummy loop indices 180 167 INTEGER :: numlight 181 REAL(wp) :: z toto168 REAL(wp) :: zchl 182 169 CHARACTER(LEN=20) :: clname 183 170 !!---------------------------------------------------------------------- … … 186 173 CALL ctlopn( numlight, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', 1, numout, .TRUE., 1 ) 187 174 DO jchl = 1, 61 188 READ(numlight,*) z toto, ( prgb(jband,jchl), jband=1,3 )175 READ(numlight,*) zchl, ( prgb(jband,jchl), jband=1,3 ) 189 176 END DO 190 177 CLOSE( numlight )
Note: See TracChangeset
for help on using the changeset viewer.