Changeset 4521 for branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/TOP_SRC
- Timestamp:
- 2014-03-03T19:12:24+01:00 (10 years ago)
- Location:
- branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/TOP_SRC/PISCES
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90
r4153 r4521 233 233 ! -------------------------------------------------------------------------------------- 234 234 IF( ln_fechem ) THEN 235 zfeequi = ( zFe3(ji,jj,jk) + zFe2(ji,jj,jk) ) * 1E-9235 zfeequi = ( zFe3(ji,jj,jk) + zFe2(ji,jj,jk) + zFeP(ji,jj,jk) ) * 1E-9 236 236 zfecoll = ( 0.3 * zFeL1(ji,jj,jk) + 0.5 * zFeL2(ji,jj,jk) ) * 1E-9 237 237 ELSE -
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90
r4148 r4521 203 203 204 204 zmortz2 = ztortz2 + zrespz2 205 zmortzgoc = unass2 / ( 1. - zepsherv ) * ztortz2205 zmortzgoc = unass2 / ( 1. - epsher2 ) * ztortz2 + zrespz2 206 206 tra(ji,jj,jk,jpmes) = tra(ji,jj,jk,jpmes) - zmortz2 + zepsherv * zgraztot 207 207 tra(ji,jj,jk,jpdia) = tra(ji,jj,jk,jpdia) - zgrazd -
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90
r4147 r4521 184 184 tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) + zprcaca 185 185 #if defined key_kriest 186 tra(ji,jj,jk,jpnum) = tra(ji,jj,jk,jpnum) + ( zmortz - zgrazm ) * xkr_dmicro 186 tra(ji,jj,jk,jpnum) = tra(ji,jj,jk,jpnum) + zmortz * xkr_dmicro & 187 - zgrazm * trn(ji,jj,jk,jpnum) / ( trn(ji,jj,jk,jppoc) + rtrn ) 187 188 #endif 188 189 END DO -
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90
r4230 r4521 479 479 CALL iom_close( numhydro ) 480 480 ! 481 hydrofe(:,:,:) = ( hydrofe(:,:,:) * hratio ) / ( cvol(:,:,:) * ryyss + rtrn ) / fse3t(:,:,:)481 hydrofe(:,:,:) = ( hydrofe(:,:,:) * hratio ) / ( cvol(:,:,:) * ryyss + rtrn ) / 1000._wp 482 482 ! 483 483 ENDIF -
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90
r4162 r4521 70 70 REAL(wp) :: zlim, zfact, zfactcal 71 71 REAL(wp) :: zo2, zno3, zflx, zpdenit, z1pdenit, zdenitt, zolimit 72 REAL(wp) :: zsiloss, zcaloss, zws bio3, zwsbio4, zwscal, zdep, zwstpoc72 REAL(wp) :: zsiloss, zcaloss, zws3, zws4, zwsc, zdep, zwstpoc 73 73 REAL(wp) :: ztrfer, ztrpo4, zwdust 74 74 REAL(wp) :: zrdenittot, zsdenittot, znitrpottot … … 77 77 REAL(wp), POINTER, DIMENSION(:,: ) :: zpdep, zsidep, zwork1, zwork2, zwork3, zwork4 78 78 REAL(wp), POINTER, DIMENSION(:,: ) :: zdenit2d, zironice 79 REAL(wp), POINTER, DIMENSION(:,: ) :: zwsbio3, zwsbio4, zwscal 79 80 REAL(wp), POINTER, DIMENSION(:,:,:) :: znitrpot, zirondep 80 81 !!--------------------------------------------------------------------- … … 93 94 ! Allocate temporary workspace 94 95 CALL wrk_alloc( jpi, jpj, zdenit2d, zwork1, zwork2, zwork3, zwork4 ) 96 CALL wrk_alloc( jpi, jpj, zwsbio3, zwsbio4, zwscal ) 95 97 CALL wrk_alloc( jpi, jpj, jpk, znitrpot ) 98 96 99 zdenit2d(:,:) = 0.e0 97 100 … … 196 199 ENDIF 197 200 201 202 ! OA: Warning, the following part is necessary, especially with Kriest 203 ! to avoid CFL problems above the sediments 204 ! -------------------------------------------------------------------- 205 DO jj = 1, jpj 206 DO ji = 1, jpi 207 ikt = mbkt(ji,jj) 208 zdep = fse3t(ji,jj,ikt) / xstep 209 zwsbio4(ji,jj) = MIN( 0.99 * zdep, wsbio4(ji,jj,ikt) ) 210 zwscal (ji,jj) = MIN( 0.99 * zdep, wscal (ji,jj,ikt) ) 211 zwsbio3(ji,jj) = MIN( 0.99 * zdep, wsbio3(ji,jj,ikt) ) 212 END DO 213 END DO 214 198 215 #if ! defined key_sed 199 216 ! Computation of the sediment denitrification proportion: The metamodel from midlleburg (2006) is being used … … 204 221 ikt = mbkt(ji,jj) 205 222 # if defined key_kriest 206 zflx = trn(ji,jj,ikt,jppoc) * wsbio3(ji,jj,ikt) * 1E3 * 1E6 / 1E4223 zflx = trn(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) * 1E3 * 1E6 / 1E4 207 224 # else 208 zflx = ( trn(ji,jj,ikt,jpgoc) * wsbio4(ji,jj,ikt) &209 & + trn(ji,jj,ikt,jppoc) * wsbio3(ji,jj,ikt) ) * 1E3 * 1E6 / 1E4225 zflx = ( trn(ji,jj,ikt,jpgoc) * zwsbio4(ji,jj) & 226 & + trn(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) ) * 1E3 * 1E6 / 1E4 210 227 #endif 211 228 zflx = LOG10( MAX( 1E-3, zflx ) ) … … 227 244 ikt = mbkt(ji,jj) 228 245 # if defined key_kriest 229 zwork1(ji,jj) = trn(ji,jj,ikt,jpgsi) * wscal (ji,jj,ikt)230 zwork2(ji,jj) = trn(ji,jj,ikt,jppoc) * wsbio3(ji,jj,ikt)246 zwork1(ji,jj) = trn(ji,jj,ikt,jpgsi) * zwscal (ji,jj) 247 zwork2(ji,jj) = trn(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) 231 248 # else 232 zwork1(ji,jj) = trn(ji,jj,ikt,jpgsi) * wsbio4(ji,jj,ikt)233 zwork2(ji,jj) = trn(ji,jj,ikt,jpgoc) * wsbio4(ji,jj,ikt) + trn(ji,jj,ikt,jppoc) * wsbio3(ji,jj,ikt)249 zwork1(ji,jj) = trn(ji,jj,ikt,jpgsi) * zwsbio4(ji,jj) 250 zwork2(ji,jj) = trn(ji,jj,ikt,jpgoc) * zwsbio4(ji,jj) + trn(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) 234 251 # endif 235 252 ! For calcite, burial efficiency is made a function of saturation 236 253 zfactcal = MIN( excess(ji,jj,ikt), 0.2 ) 237 254 zfactcal = MIN( 1., 1.3 * ( 0.2 - zfactcal ) / ( 0.4 - zfactcal ) ) 238 zwork3(ji,jj) = trn(ji,jj,ikt,jpcal) * wscal (ji,jj,ikt) * 2.e0 * zfactcal255 zwork3(ji,jj) = trn(ji,jj,ikt,jpcal) * zwscal(ji,jj) * 2.e0 * zfactcal 239 256 END DO 240 257 END DO … … 256 273 ikt = mbkt(ji,jj) 257 274 zdep = xstep / fse3t(ji,jj,ikt) 258 zws bio4 = wsbio4(ji,jj,ikt) * zdep259 zwsc al = wscal (ji,jj,ikt) * zdep275 zws4 = zwsbio4(ji,jj) * zdep 276 zwsc = zwscal (ji,jj) * zdep 260 277 # if defined key_kriest 261 zsiloss = trn(ji,jj,ikt,jpgsi) * zws bio4278 zsiloss = trn(ji,jj,ikt,jpgsi) * zws4 262 279 # else 263 zsiloss = trn(ji,jj,ikt,jpgsi) * zwsc al280 zsiloss = trn(ji,jj,ikt,jpgsi) * zwsc 264 281 # endif 265 zcaloss = trn(ji,jj,ikt,jpcal) * zwsc al282 zcaloss = trn(ji,jj,ikt,jpcal) * zwsc 266 283 ! 267 284 trn(ji,jj,ikt,jpgsi) = trn(ji,jj,ikt,jpgsi) - zsiloss … … 282 299 ikt = mbkt(ji,jj) 283 300 zdep = xstep / fse3t(ji,jj,ikt) 284 zws bio4 = wsbio4(ji,jj,ikt) * zdep285 zws bio3 = wsbio3(ji,jj,ikt) * zdep301 zws4 = zwsbio4(ji,jj) * zdep 302 zws3 = zwsbio3(ji,jj) * zdep 286 303 # if ! defined key_kriest 287 trn(ji,jj,ikt,jpgoc) = trn(ji,jj,ikt,jpgoc) - trn(ji,jj,ikt,jpgoc) * zws bio4288 trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc) - trn(ji,jj,ikt,jppoc) * zws bio3289 trn(ji,jj,ikt,jpbfe) = trn(ji,jj,ikt,jpbfe) - trn(ji,jj,ikt,jpbfe) * zws bio4290 trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe) - trn(ji,jj,ikt,jpsfe) * zws bio3291 zwstpoc = trn(ji,jj,ikt,jpgoc) * zwsbio4 + trn(ji,jj,ikt,jppoc) * zwsbio3304 trn(ji,jj,ikt,jpgoc) = trn(ji,jj,ikt,jpgoc) - trn(ji,jj,ikt,jpgoc) * zws4 305 trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc) - trn(ji,jj,ikt,jppoc) * zws3 306 trn(ji,jj,ikt,jpbfe) = trn(ji,jj,ikt,jpbfe) - trn(ji,jj,ikt,jpbfe) * zws4 307 trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe) - trn(ji,jj,ikt,jpsfe) * zws3 308 zwstpoc = trn(ji,jj,ikt,jpgoc) * zws4 + trn(ji,jj,ikt,jppoc) * zws3 292 309 # else 293 trn(ji,jj,ikt,jpnum) = trn(ji,jj,ikt,jpnum) - trn(ji,jj,ikt,jpnum) * zws bio4294 trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc) - trn(ji,jj,ikt,jppoc) * zws bio3295 trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe) - trn(ji,jj,ikt,jpsfe) * zws bio3296 zwstpoc = trn(ji,jj,ikt,jppoc) * zws bio3310 trn(ji,jj,ikt,jpnum) = trn(ji,jj,ikt,jpnum) - trn(ji,jj,ikt,jpnum) * zws4 311 trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc) - trn(ji,jj,ikt,jppoc) * zws3 312 trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe) - trn(ji,jj,ikt,jpsfe) * zws3 313 zwstpoc = trn(ji,jj,ikt,jppoc) * zws3 297 314 # endif 298 315 … … 383 400 ! 384 401 CALL wrk_dealloc( jpi, jpj, zdenit2d, zwork1, zwork2, zwork3, zwork4 ) 402 CALL wrk_dealloc( jpi, jpj, zwsbio3, zwsbio4, zwscal ) 385 403 CALL wrk_dealloc( jpi, jpj, jpk, znitrpot ) 386 404 ! -
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90
r4162 r4521 121 121 rdenit = 105._wp / 16._wp 122 122 rdenita = 3._wp / 5._wp 123 o2ut = 13 1._wp / 122._wp123 o2ut = 133._wp / 122._wp 124 124 125 125 ! Initialization of tracer concentration in case of no restart
Note: See TracChangeset
for help on using the changeset viewer.