- Timestamp:
- 2017-05-09T12:14:45+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zsed.F90
r7617 r8003 244 244 IF( tmask(ji,jj,1) == 1 ) THEN 245 245 ikt = mbkt(ji,jj) 246 # if defined key_kriest247 zflx = trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) * 1E3 * 1E6 / 1E4248 # else249 246 zflx = ( trb(ji,jj,ikt,jpgoc) * zwsbio4(ji,jj) & 250 247 & + trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) ) * 1E3 * 1E6 / 1E4 251 #endif252 248 zflx = LOG10( MAX( 1E-3, zflx ) ) 253 249 zo2 = LOG10( MAX( 10. , trb(ji,jj,ikt,jpoxy) * 1E6 ) ) … … 273 269 IF( tmask(ji,jj,1) == 1 ) THEN 274 270 ikt = mbkt(ji,jj) 275 # if defined key_kriest276 zwork1(ji,jj) = trb(ji,jj,ikt,jpgsi) * zwscal (ji,jj)277 zwork2(ji,jj) = trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj)278 # else279 271 zwork1(ji,jj) = trb(ji,jj,ikt,jpgsi) * zwsbio4(ji,jj) 280 272 zwork2(ji,jj) = trb(ji,jj,ikt,jpgoc) * zwsbio4(ji,jj) + trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) 281 # endif282 273 ! For calcite, burial efficiency is made a function of saturation 283 274 zfactcal = MIN( excess(ji,jj,ikt), 0.2 ) … … 305 296 zws4 = zwsbio4(ji,jj) * zdep 306 297 zwsc = zwscal (ji,jj) * zdep 307 # if defined key_kriest308 zsiloss = trb(ji,jj,ikt,jpgsi) * zws4309 # else310 298 zsiloss = trb(ji,jj,ikt,jpgsi) * zwsc 311 # endif312 299 zcaloss = trb(ji,jj,ikt,jpcal) * zwsc 313 300 ! … … 336 323 #endif 337 324 zrivno3 = 1. - zbureff(ji,jj) 338 # if ! defined key_kriest339 325 trb(ji,jj,ikt,jpgoc) = trb(ji,jj,ikt,jpgoc) - trb(ji,jj,ikt,jpgoc) * zws4 340 326 trb(ji,jj,ikt,jppoc) = trb(ji,jj,ikt,jppoc) - trb(ji,jj,ikt,jppoc) * zws3 … … 351 337 trn(ji,jj,ikt,jpfep) = trn(ji,jj,ikt,jpfep) - trn(ji,jj,ikt,jpfep) * zwssfep 352 338 # endif 353 # else354 trb(ji,jj,ikt,jpnum) = trb(ji,jj,ikt,jpnum) - trb(ji,jj,ikt,jpnum) * zws4355 trb(ji,jj,ikt,jppoc) = trb(ji,jj,ikt,jppoc) - trb(ji,jj,ikt,jppoc) * zws3356 trb(ji,jj,ikt,jppon) = trb(ji,jj,ikt,jppon) - trb(ji,jj,ikt,jppon) * zws3357 trb(ji,jj,ikt,jppop) = trb(ji,jj,ikt,jppop) - trb(ji,jj,ikt,jppop) * zws3358 trb(ji,jj,ikt,jpsfe) = trb(ji,jj,ikt,jpsfe) - trb(ji,jj,ikt,jpsfe) * zws3359 zwstpop = trb(ji,jj,ikt,jppop) * zws3360 zwstpon = trb(ji,jj,ikt,jppon) * zws3361 zwstpoc = trb(ji,jj,ikt,jppoc) * zws3362 # if defined key_ligand363 trn(ji,jj,ikt,jpfep) = trn(ji,jj,ikt,jpfep) - trn(ji,jj,ikt,jpfep) * zwssfep364 # endif365 # endif366 339 367 340 #if ! defined key_sed … … 405 378 zlim = ( 1.- xdiano3 - xdianh4 ) 406 379 IF( zlim <= 0.1 ) zlim = 0.01 407 #if defined key_degrad408 zfact = zlim * rfact2 * facvol(ji,jj,jk)409 #else410 380 zfact = zlim * rfact2 411 #endif412 381 ztrfer = biron(ji,jj,jk) / ( concfediaz + biron(ji,jj,jk) ) 413 382 ztrpo4(ji,jj,jk) = trb(ji,jj,jk,jppo4) / ( 1E-6 + trb(ji,jj,jk,jppo4) )
Note: See TracChangeset
for help on using the changeset viewer.