Changeset 1180
- Timestamp:
- 2008-09-11T18:41:47+02:00 (16 years ago)
- Location:
- trunk/NEMO/TOP_SRC/PISCES
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/PISCES/p4zche.F90
r1152 r1180 105 105 kf4 =-0.001005 106 106 107 108 107 REAL(wp) :: & ! coeff. for 1. dissoc. of boric acid (Dickson and Goyet, 1994) 109 108 cb0 = -8966.90, & … … 129 128 cw6 = -0.01615 130 129 131 REAL(wp) :: & ! coeff. for dissoc. of phosphate (Millero (1974)132 cp10 = 115.54 , &133 cp11 = -4576.752 , &134 cp12 = -18.453 , &135 cp13 = -106.736 , &136 cp14 = 0.69171 , &137 cp15 = -0.65643 , &138 cp16 = -0.01844 , &139 cp20 = 172.1033 , &140 cp21 = -8814.715 , &141 cp22 = -27.927 , &142 cp23 = -160.340 , &143 cp24 = 1.3566 , &144 cp25 = 0.37335 , &145 cp26 = -0.05778 , &146 cp30 = -18.126 , &147 cp31 = -3070.75 , &148 cp32 = 17.27039 , &149 cp33 = 2.81197 , &150 cp34 = -44.99486 , &151 cp35 = -0.09984152 153 REAL(wp) :: & ! coeff. for dissoc. of silicates (Millero (1974)154 cs10 = 117.385 , &155 cs11 = -8904.2 , &156 cs12 = -19.334 , &157 cs13 = -458.79 , &158 cs14 = 3.5913 , &159 cs15 = 188.74 , &160 cs16 = -1.5998 , &161 cs17 = -12.1652 , &162 cs18 = 0.07871 , &163 cs19 = 0. , &164 cs20 = 1. , &165 cs21 = -0.001005166 167 130 REAL(wp) :: & ! volumetric solubility constants for o2 in ml/l (Weiss, 1974) 168 131 ox0 = -58.3877 , & -
trunk/NEMO/TOP_SRC/PISCES/p4zopt.F90
r1152 r1180 27 27 etot, enano, ediat, & !: PAR for phyto, nano and diat 28 28 emoy !: averaged PAR in the mixed layer 29 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: & !:30 heup !: Depth of the euphotic zone31 29 32 30 !! * Module variables … … 210 208 ! Computation of the euphotic depth 211 209 ! --------------------------------- 212 213 heup(:,:) = 300.e0 210 ! Euphotic layer bottom level 211 neln(:,:) = 1 ! initialisation of EL level 212 heup(:,:) = 300. 214 213 215 214 DO jk = 2, jpkm1 216 215 DO jj = 1, jpj 217 DO ji = 1, jpi 218 IF( etot(ji,jj,jk) >= 0.0043 * qsr(ji,jj) ) heup(ji,jj) = fsdepw(ji,jj,jk+1) 219 END DO 220 END DO 221 END DO 222 223 heup(:,:) = MIN( 300., heup(:,:) ) 216 DO ji = 1, jpi 217 IF( etot(ji,jj,jk) >= 0.0043 * qsr(ji,jj) ) THEN 218 neln(ji,jj) = jk+1 ! 1rst T-level strictly below EL bottom 219 ! ! nb. this is to ensure compatibility with 220 ! ! nmld_trc definition in trd_mld_trc_zint 221 heup(ji,jj) = fsdepw(ji,jj,jk+1) ! Euphotic layer depth 222 ENDIF 223 END DO 224 END DO 225 ENDDO 226 227 heup(:,:) = MIN( 300., heup(:,:) ) 224 228 225 229 ! Computation of the mean light over the mixed layer depth -
trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90
r1152 r1180 17 17 USE sms_pisces ! 18 18 USE prtctl_trc 19 USE p4zopt 19 20 USE p4zint 20 21 USE p4zlim 21 USE p4zopt22 22 23 23 IMPLICIT NONE -
trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90
r1152 r1180 63 63 !!---------------------------------------------------------------------- 64 64 !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 65 !! $ Id$65 !! $Header:$ 66 66 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 67 67 !!---------------------------------------------------------------------- … … 82 82 INTEGER :: ji, jj, jk 83 83 INTEGER :: ikt 84 #if ! defined key_sed 84 85 REAL(wp) :: zsumsedsi, zsumsedpo4, zsumsedcal 86 #endif 85 87 REAL(wp) :: zconctmp , zdenitot , znitrpottot 86 88 REAL(wp) :: zlim, zconctmp2, zstep, zfact … … 137 139 & + zirondep(:,:,jk) + ironsed(:,:,jk) * rfact2 138 140 END DO 141 142 143 #if ! defined key_sed 139 144 140 145 ! Initialisation of variables used to compute Sinking Speed … … 174 179 CALL mpp_sum( zsumsedpo4 ) ! sums over the global domain 175 180 ENDIF 181 182 #endif 176 183 177 184 ! Then this loss is scaled at each bottom grid cell for … … 192 199 193 200 trn(ji,jj,ikt,jpdsi) = trn(ji,jj,ikt,jpdsi) - zconctmp 201 202 #if ! defined key_sed 194 203 trn(ji,jj,ikt,jpsil) = trn(ji,jj,ikt,jpsil) + zconctmp & 195 204 & * ( 1.- ( sumdepsi + rivalkinput / raass / 6. ) / zsumsedsi ) 205 #endif 196 206 END DO 197 207 END DO … … 202 212 zconctmp = trn(ji,jj,ikt,jpcal) * wscal(ji,jj,ikt) * zstep / fse3t(ji,jj,ikt) 203 213 trn(ji,jj,ikt,jpcal) = trn(ji,jj,ikt,jpcal) - zconctmp 214 215 #if ! defined key_sed 204 216 trn(ji,jj,ikt,jptal) = trn(ji,jj,ikt,jptal) + zconctmp & 205 217 & * ( 1.- ( rivalkinput / raass ) / zsumsedcal ) * 2.e0 206 218 trn(ji,jj,ikt,jpdic) = trn(ji,jj,ikt,jpdic) + zconctmp & 207 219 & * ( 1.- ( rivalkinput / raass ) / zsumsedcal ) 220 #endif 208 221 END DO 209 222 END DO … … 219 232 trn(ji,jj,ikt,jpgoc) = trn(ji,jj,ikt,jpgoc) - zconctmp * wsbio4(ji,jj,ikt) * zfact 220 233 trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc) - zconctmp2 * wsbio3(ji,jj,ikt) * zfact 234 #if ! defined key_sed 221 235 trn(ji,jj,ikt,jpdoc) = trn(ji,jj,ikt,jpdoc) & 222 236 & + ( zconctmp * wsbio4(ji,jj,ikt) + zconctmp2 * wsbio3(ji,jj,ikt) ) * zfact & 223 237 & * ( 1.- rivpo4input / (raass * zsumsedpo4 ) ) 238 #endif 224 239 trn(ji,jj,ikt,jpbfe) = trn(ji,jj,ikt,jpbfe) - trn(ji,jj,ikt,jpbfe) * wsbio4(ji,jj,ikt) * zfact 225 240 trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe) - trn(ji,jj,ikt,jpsfe) * wsbio3(ji,jj,ikt) * zfact … … 233 248 trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc) & 234 249 & - zconctmp2 * wsbio3(ji,jj,ikt) * zfact 250 #if ! defined key_sed 235 251 trn(ji,jj,ikt,jpdoc) = trn(ji,jj,ikt,jpdoc) & 236 252 & + ( zconctmp2 * wsbio3(ji,jj,ikt) ) & 237 253 & * zfact * ( 1.- rivpo4input / ( raass * zsumsedpo4 ) ) 254 #endif 238 255 trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe) & 239 256 & - trn(ji,jj,ikt,jpsfe) * wsbio3(ji,jj,ikt) * zfact -
trunk/NEMO/TOP_SRC/PISCES/p4zsink.F90
r1152 r1180 27 27 28 28 !! * Module variables 29 REAL(wp), DIMENSION(jpi,jpj,jpk) :: & !:29 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: & !: 30 30 sinking, sinking2, & !: POC sinking fluxes (different meanings depending on the parameterization 31 31 sinkcal, sinksil, & !: CaCO3 and BSi sinking fluxes -
trunk/NEMO/TOP_SRC/PISCES/sms_pisces.F90
r1152 r1180 47 47 !!--------------------------------------------- 48 48 REAL(wp), DIMENSION(jpi,jpj,jpk) :: etot3 !: ??? 49 INTEGER , DIMENSION(jpi,jpj) :: neln !: number of levels in the euphotic layer 50 REAL(wp), DIMENSION(jpi,jpj) :: heup !: euphotic layer depth 49 51 50 52 !!----------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.