- Timestamp:
- 2017-01-04T17:47:47+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90
r7508 r7525 209 209 ! Exchange between organic matter compartments due to coagulation/disaggregation 210 210 ! --------------------------------------------------- 211 !$OMP PARALLEL DO schedule(static) private(jk, jj, ji, zstep, zfact, zagg, zagg1, zagg2, zagg3, zagg4, zaggfe, zaggdoc, zaggdoc2, zaggdoc3) 211 !$OMP PARALLEL DO schedule(static) private(jk, jj, ji, zstep, zfact, zagg, zagg1, zagg2) & 212 !$OMP& private(zagg3, zagg4, zaggfe, zaggdoc, zaggdoc2, zaggdoc3) 212 213 DO jk = 1, jpkm1 213 214 DO jj = 1, jpj … … 259 260 ! Total carbon export per year 260 261 IF( iom_use( "tcexp" ) .OR. ( ln_check_mass .AND. kt == nitend .AND. knt == nrdttrc ) ) & 261 & t_oce_co2_exp = glob_sum( ( sinking(:,:,ik100) + sinking2(:,:,ik100) ) * e1e2t(:,:) * tmask(:,:,1) )262 & t_oce_co2_exp = glob_sum( ( sinking(:,:,ik100) + sinking2(:,:,ik100) ) * e1e2t(:,:) * tmask(:,:,1) ) 262 263 ! 263 264 IF( lk_iomput ) THEN … … 269 270 IF( iom_use( "EPC100" ) ) THEN 270 271 !$OMP DO schedule(static) private(jj, ji) 271 DO jj = 1, jpj272 DO ji = 1, jpi273 zw2d(ji,jj) = ( sinking(ji,jj,ik100) + sinking2(ji,jj,ik100) ) * zfact * tmask(ji,jj,1) ! Export of carbon at 100m274 END DO275 END DO276 272 DO jj = 1, jpj 273 DO ji = 1, jpi 274 zw2d(ji,jj) = ( sinking(ji,jj,ik100) + sinking2(ji,jj,ik100) ) * zfact * tmask(ji,jj,1) ! Export of carbon at 100m 275 END DO 276 END DO 277 CALL iom_put( "EPC100" , zw2d ) 277 278 ENDIF 278 279 IF( iom_use( "EPFE100" ) ) THEN 279 280 !$OMP DO schedule(static) private(jj, ji) 280 DO jj = 1, jpj281 DO ji = 1, jpi282 zw2d(ji,jj) = ( sinkfer(ji,jj,ik100) + sinkfer2(ji,jj,ik100) ) * zfact * tmask(ji,jj,1) ! Export of iron at 100m283 END DO284 END DO285 281 DO jj = 1, jpj 282 DO ji = 1, jpi 283 zw2d(ji,jj) = ( sinkfer(ji,jj,ik100) + sinkfer2(ji,jj,ik100) ) * zfact * tmask(ji,jj,1) ! Export of iron at 100m 284 END DO 285 END DO 286 CALL iom_put( "EPFE100" , zw2d ) 286 287 ENDIF 287 288 IF( iom_use( "EPCAL100" ) ) THEN 288 289 !$OMP DO schedule(static) private(jj, ji) 289 DO jj = 1, jpj290 DO ji = 1, jpi291 zw2d(ji,jj) = sinkcal(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of calcite at 100m292 END DO293 END DO294 290 DO jj = 1, jpj 291 DO ji = 1, jpi 292 zw2d(ji,jj) = sinkcal(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of calcite at 100m 293 END DO 294 END DO 295 CALL iom_put( "EPCAL100" , zw2d ) 295 296 ENDIF 296 297 IF( iom_use( "EPSI100" ) ) THEN 297 298 !$OMP DO schedule(static) private(jj, ji) 298 DO jj = 1, jpj299 DO ji = 1, jpi300 zw2d(ji,jj) = sinksil(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of bigenic silica at 100m301 END DO302 END DO303 299 DO jj = 1, jpj 300 DO ji = 1, jpi 301 zw2d(ji,jj) = sinksil(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of bigenic silica at 100m 302 END DO 303 END DO 304 CALL iom_put( "EPSI100" , zw2d ) 304 305 ENDIF 305 306 IF( iom_use( "EXPC" ) ) THEN 306 307 !$OMP DO schedule(static) private(jk, jj, ji) 307 DO jk = 1, jpk308 DO jj = 1, jpj309 DO ji = 1, jpi310 zw3d(ji,jj,jk) = ( sinking(ji,jj,jk) + sinking2(ji,jj,jk) ) * zfact * tmask(ji,jj,jk) ! Export of carbon in the water column311 END DO312 END DO313 END DO314 308 DO jk = 1, jpk 309 DO jj = 1, jpj 310 DO ji = 1, jpi 311 zw3d(ji,jj,jk) = ( sinking(ji,jj,jk) + sinking2(ji,jj,jk) ) * zfact * tmask(ji,jj,jk) ! Export of carbon in the water column 312 END DO 313 END DO 314 END DO 315 CALL iom_put( "EXPC" , zw3d ) 315 316 ENDIF 316 317 IF( iom_use( "EXPFE" ) ) THEN 317 318 !$OMP DO schedule(static) private(jk, jj, ji) 318 DO jk = 1, jpk319 DO jj = 1, jpj320 DO ji = 1, jpi321 zw3d(ji,jj,jk) = ( sinkfer(ji,jj,jk) + sinkfer2(ji,jj,jk) ) * zfact * tmask(ji,jj,jk) ! Export of iron322 END DO323 END DO324 END DO325 319 DO jk = 1, jpk 320 DO jj = 1, jpj 321 DO ji = 1, jpi 322 zw3d(ji,jj,jk) = ( sinkfer(ji,jj,jk) + sinkfer2(ji,jj,jk) ) * zfact * tmask(ji,jj,jk) ! Export of iron 323 END DO 324 END DO 325 END DO 326 CALL iom_put( "EXPFE" , zw3d ) 326 327 ENDIF 327 328 IF( iom_use( "EXPCAL" ) ) THEN 328 329 !$OMP DO schedule(static) private(jk, jj, ji) 329 DO jk = 1, jpk330 DO jj = 1, jpj331 DO ji = 1, jpi332 zw3d(ji,jj,jk) = sinkcal(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of calcite333 END DO334 END DO335 END DO336 330 DO jk = 1, jpk 331 DO jj = 1, jpj 332 DO ji = 1, jpi 333 zw3d(ji,jj,jk) = sinkcal(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of calcite 334 END DO 335 END DO 336 END DO 337 CALL iom_put( "EXPCAL" , zw3d ) 337 338 ENDIF 338 339 IF( iom_use( "EXPSI" ) ) THEN 339 340 !$OMP DO schedule(static) private(jk, jj, ji) 340 DO jk = 1, jpk341 DO jj = 1, jpj342 DO ji = 1, jpi343 zw3d(ji,jj,jk) = sinksil(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of bigenic silica344 END DO345 END DO346 END DO347 341 DO jk = 1, jpk 342 DO jj = 1, jpj 343 DO ji = 1, jpi 344 zw3d(ji,jj,jk) = sinksil(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of bigenic silica 345 END DO 346 END DO 347 END DO 348 CALL iom_put( "EXPSI" , zw3d ) 348 349 ENDIF 349 350 IF( iom_use( "tcexp" ) ) CALL iom_put( "tcexp" , t_oce_co2_exp * zfact ) ! molC/s … … 526 527 zval4 = 4. + xkr_eta 527 528 528 !$OMP PARALLEL DO schedule(static) private(jk, jj, ji, znum, zeps, zfm, zsm, zdiv, zdiv1, zdiv2, zdiv3, zdiv4, zdiv5, zagg, zagg1, zagg2, zagg3, zagg4, zagg5, zfract, zaggdoc, zaggdoc1, zaggsh, zaggsi, znumdoc) 529 !$OMP PARALLEL DO schedule(static) private(jk, jj, ji, znum, zeps, zfm, zsm, zdiv, zdiv1, zdiv2, zdiv3, zdiv4) & 530 !$OMP& private(zdiv5, zagg, zagg1, zagg2, zagg3, zagg4, zagg5, zfract, zaggdoc, zaggdoc1, zaggsh, zaggsi, znumdoc) 529 531 DO jk = 1,jpkm1 530 532 DO jj = 1,jpj … … 632 634 IF( iom_use( "EPC100" ) ) THEN 633 635 !$OMP DO schedule(static) private(jj, ji) 634 DO jj = 1, jpj635 DO ji = 1, jpi636 zw2d(ji,jj) = sinking(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of carbon at 100m637 END DO638 END DO639 636 DO jj = 1, jpj 637 DO ji = 1, jpi 638 zw2d(ji,jj) = sinking(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of carbon at 100m 639 END DO 640 END DO 641 CALL iom_put( "EPC100" , zw2d ) 640 642 ENDIF 641 643 IF( iom_use( "EPN100" ) ) THEN 642 644 !$OMP DO schedule(static) private(jj, ji) 643 DO jj = 1, jpj644 DO ji = 1, jpi645 zw2d(ji,jj) = sinking2(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of number of aggregates ?646 END DO647 END DO648 645 DO jj = 1, jpj 646 DO ji = 1, jpi 647 zw2d(ji,jj) = sinking2(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of number of aggregates ? 648 END DO 649 END DO 650 CALL iom_put( "EPN100" , zw2d ) 649 651 ENDIF 650 652 IF( iom_use( "EPCAL100" ) ) THEN 651 653 !$OMP DO schedule(static) private(jj, ji) 652 DO jj = 1, jpj653 DO ji = 1, jpi654 zw2d(ji,jj) = sinkcal(ji,jj,ik100) * zfact * tmask(ji,jj,1) !Export of calcite at 100m655 END DO656 END DO657 654 DO jj = 1, jpj 655 DO ji = 1, jpi 656 zw2d(ji,jj) = sinkcal(ji,jj,ik100) * zfact * tmask(ji,jj,1) !Export of calcite at 100m 657 END DO 658 END DO 659 CALL iom_put( "EPCAL100" , zw2d ) 658 660 ENDIF 659 661 IF( iom_use( "EPSI100" ) ) THEN 660 662 !$OMP DO schedule(static) private(jj, ji) 661 DO jj = 1, jpj662 DO ji = 1, jpi663 zw2d(ji,jj) = sinksil(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of bigenic silica at 100m664 END DO665 END DO666 663 DO jj = 1, jpj 664 DO ji = 1, jpi 665 zw2d(ji,jj) = sinksil(ji,jj,ik100) * zfact * tmask(ji,jj,1) ! Export of bigenic silica at 100m 666 END DO 667 END DO 668 CALL iom_put( "EPSI100" , zw2d ) 667 669 ENDIF 668 670 IF( iom_use( "EXPC" ) ) THEN 669 671 !$OMP DO schedule(static) private(jk, jj, ji) 670 DO jk = 1, jpk671 DO jj = 1, jpj672 DO ji = 1, jpi673 zw3d(ji,jj,jk) = sinking(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of carbon in the water column674 END DO675 END DO676 END DO677 672 DO jk = 1, jpk 673 DO jj = 1, jpj 674 DO ji = 1, jpi 675 zw3d(ji,jj,jk) = sinking(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of carbon in the water column 676 END DO 677 END DO 678 END DO 679 CALL iom_put( "EXPC" , zw3d ) 678 680 ENDIF 679 681 IF( iom_use( "EXPN" ) ) THEN 680 682 !$OMP DO schedule(static) private(jk, jj, ji) 681 DO jk = 1, jpk682 DO jj = 1, jpj683 DO ji = 1, jpi684 zw3d(ji,jj,jk) = sinking(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of carbon in the water column685 END DO686 END DO687 END DO688 683 DO jk = 1, jpk 684 DO jj = 1, jpj 685 DO ji = 1, jpi 686 zw3d(ji,jj,jk) = sinking(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of carbon in the water column 687 END DO 688 END DO 689 END DO 690 CALL iom_put( "EXPN" , zw3d ) 689 691 ENDIF 690 692 IF( iom_use( "EXPCAL" ) ) THEN 691 693 !$OMP DO schedule(static) private(jk, jj, ji) 692 DO jk = 1, jpk693 DO jj = 1, jpj694 DO ji = 1, jpi695 zw3d(ji,jj,jk) = sinkcal(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of calcite696 END DO697 END DO698 END DO699 694 DO jk = 1, jpk 695 DO jj = 1, jpj 696 DO ji = 1, jpi 697 zw3d(ji,jj,jk) = sinkcal(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of calcite 698 END DO 699 END DO 700 END DO 701 CALL iom_put( "EXPCAL" , zw3d ) 700 702 ENDIF 701 703 IF( iom_use( "EXPSI" ) ) THEN 702 704 !$OMP DO schedule(static) private(jk, jj, ji) 703 DO jk = 1, jpk704 DO jj = 1, jpj705 DO ji = 1, jpi706 zw3d(ji,jj,jk) = sinksil(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of bigenic silica707 END DO708 END DO709 END DO710 705 DO jk = 1, jpk 706 DO jj = 1, jpj 707 DO ji = 1, jpi 708 zw3d(ji,jj,jk) = sinksil(ji,jj,jk) * zfact * tmask(ji,jj,jk) ! Export of bigenic silica 709 END DO 710 END DO 711 END DO 712 CALL iom_put( "EXPSI" , zw3d ) 711 713 ENDIF 712 714 IF( iom_use( "XNUM" ) ) THEN 713 715 !$OMP DO schedule(static) private(jk, jj, ji) 714 DO jk = 1, jpk715 DO jj = 1, jpj716 DO ji = 1, jpi717 zw3d(ji,jj,jk) = znum3d(ji,jj,jk) * tmask(ji,jj,jk) ! Number of particles on aggregats718 END DO719 END DO720 END DO721 716 DO jk = 1, jpk 717 DO jj = 1, jpj 718 DO ji = 1, jpi 719 zw3d(ji,jj,jk) = znum3d(ji,jj,jk) * tmask(ji,jj,jk) ! Number of particles on aggregats 720 END DO 721 END DO 722 END DO 723 CALL iom_put( "XNUM" , zw3d ) 722 724 ENDIF 723 725 IF( iom_use( "WSC" ) ) THEN 724 726 !$OMP DO schedule(static) private(jk, jj, ji) 725 DO jk = 1, jpk726 DO jj = 1, jpj727 DO ji = 1, jpi728 zw3d(ji,jj,jk) = wsbio3(ji,jj,jk) * tmask(ji,jj,jk) ! Sinking speed of carbon particles729 END DO730 END DO731 END DO732 727 DO jk = 1, jpk 728 DO jj = 1, jpj 729 DO ji = 1, jpi 730 zw3d(ji,jj,jk) = wsbio3(ji,jj,jk) * tmask(ji,jj,jk) ! Sinking speed of carbon particles 731 END DO 732 END DO 733 END DO 734 CALL iom_put( "WSC" , zw3d ) 733 735 ENDIF 734 736 IF( iom_use( "WSN" ) ) THEN 735 737 !$OMP DO schedule(static) private(jk, jj, ji) 736 DO jk = 1, jpk737 DO jj = 1, jpj738 DO ji = 1, jpi739 zw3d(ji,jj,jk) = wsbio4(ji,jj,jk) * tmask(ji,jj,jk) ! Sinking speed of particles number740 END DO741 END DO742 END DO743 738 DO jk = 1, jpk 739 DO jj = 1, jpj 740 DO ji = 1, jpi 741 zw3d(ji,jj,jk) = wsbio4(ji,jj,jk) * tmask(ji,jj,jk) ! Sinking speed of particles number 742 END DO 743 END DO 744 END DO 745 CALL iom_put( "WSN" , zw3d ) 744 746 ENDIF 745 747 !
Note: See TracChangeset
for help on using the changeset viewer.