- 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/P4Z/p4zsed.F90
r7617 r8003 243 243 IF( tmask(ji,jj,1) == 1 ) THEN 244 244 ikt = mbkt(ji,jj) 245 # if defined key_kriest246 zflx = trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) * 1E3 * 1E6 / 1E4247 # else248 245 zflx = ( trb(ji,jj,ikt,jpgoc) * zwsbio4(ji,jj) & 249 246 & + trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) ) * 1E3 * 1E6 / 1E4 250 #endif251 247 zflx = LOG10( MAX( 1E-3, zflx ) ) 252 248 zo2 = LOG10( MAX( 10. , trb(ji,jj,ikt,jpoxy) * 1E6 ) ) … … 272 268 IF( tmask(ji,jj,1) == 1 ) THEN 273 269 ikt = mbkt(ji,jj) 274 # if defined key_kriest275 zwork1(ji,jj) = trb(ji,jj,ikt,jpgsi) * zwscal (ji,jj)276 zwork2(ji,jj) = trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj)277 # else278 270 zwork1(ji,jj) = trb(ji,jj,ikt,jpgsi) * zwsbio4(ji,jj) 279 271 zwork2(ji,jj) = trb(ji,jj,ikt,jpgoc) * zwsbio4(ji,jj) + trb(ji,jj,ikt,jppoc) * zwsbio3(ji,jj) 280 # endif281 272 ! For calcite, burial efficiency is made a function of saturation 282 273 zfactcal = MIN( excess(ji,jj,ikt), 0.2 ) … … 304 295 zws4 = zwsbio4(ji,jj) * zdep 305 296 zwsc = zwscal (ji,jj) * zdep 306 # if defined key_kriest307 zsiloss = trb(ji,jj,ikt,jpgsi) * zws4308 # else309 297 zsiloss = trb(ji,jj,ikt,jpgsi) * zwsc 310 # endif311 298 zcaloss = trb(ji,jj,ikt,jpcal) * zwsc 312 299 ! … … 336 323 #endif 337 324 zrivno3 = 1. - zbureff(ji,jj) 338 # if ! defined key_kriest339 325 tra(ji,jj,ikt,jpgoc) = tra(ji,jj,ikt,jpgoc) - trb(ji,jj,ikt,jpgoc) * zws4 340 326 tra(ji,jj,ikt,jppoc) = tra(ji,jj,ikt,jppoc) - trb(ji,jj,ikt,jppoc) * zws3 … … 345 331 tra(ji,jj,ikt,jpfep) = tra(ji,jj,ikt,jpfep) - trn(ji,jj,ikt,jpfep) * zwssfep 346 332 # endif 347 # else348 tra(ji,jj,ikt,jpnum) = tra(ji,jj,ikt,jpnum) - trb(ji,jj,ikt,jpnum) * zws4349 tra(ji,jj,ikt,jppoc) = tra(ji,jj,ikt,jppoc) - trb(ji,jj,ikt,jppoc) * zws3350 tra(ji,jj,ikt,jpsfe) = tra(ji,jj,ikt,jpsfe) - trb(ji,jj,ikt,jpsfe) * zws3351 zwstpoc = trb(ji,jj,ikt,jppoc) * zws3352 # if defined key_ligand353 tra(ji,jj,ikt,jpfep) = tra(ji,jj,ikt,jpfep) - trn(ji,jj,ikt,jpfep) * zwssfep354 # endif355 # endif356 333 357 334 #if ! defined key_sed … … 384 361 zlim = ( 1.- xnanono3(ji,jj,jk) - xnanonh4(ji,jj,jk) ) 385 362 IF( zlim <= 0.2 ) zlim = 0.01 386 #if defined key_degrad387 zfact = zlim * rfact2 * facvol(ji,jj,jk)388 #else389 363 zfact = zlim * rfact2 390 #endif391 364 ztrfer = biron(ji,jj,jk) / ( concfediaz + biron(ji,jj,jk) ) 392 365 ztrpo4 = trb (ji,jj,jk,jppo4) / ( concnnh4 + trb (ji,jj,jk,jppo4) ) … … 431 404 IF( iom_use("Sdenit" ) ) CALL iom_put( "Sdenit", sdenit (:,:) * 1.e+3 * rno3 ) 432 405 ENDIF 433 ELSE434 IF( ln_diatrc ) &435 & trc2d(:,:,jp_pcs0_2d + 12) = nitrpot(:,:,1) * nitrfix * rno3 * 1.e+3 * rfact2r * fse3t(:,:,1) * tmask(:,:,1)436 406 ENDIF 437 407 !
Note: See TracChangeset
for help on using the changeset viewer.