Changeset 8062
- Timestamp:
- 2017-05-24T10:03:18+02:00 (7 years ago)
- Location:
- branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_med_diag_iomput.F90
r8025 r8062 367 367 ENDDO 368 368 ENDDO 369 ELSE IF (jk.eq.mbathy(ji,jj)) THEN369 ELSE 370 370 DO jj = 2,jpjm1 371 371 DO ji = 2,jpim1 372 IF (tmask(ji,jj,jk) == 1) THEN 373 IF( med_diag%IBEN_N%dgsave ) THEN 374 iben_n2d(ji,jj) = f_sbenin_n(ji,jj) + & 375 f_fbenin_n(ji,jj) 376 ENDIF 377 IF( med_diag%IBEN_FE%dgsave ) THEN 378 iben_fe2d(ji,jj) = f_sbenin_fe(ji,jj) + & 379 f_fbenin_fe(ji,jj) 380 ENDIF 381 IF( med_diag%IBEN_C%dgsave ) THEN 382 iben_c2d(ji,jj) = f_sbenin_c(ji,jj) + & 383 f_fbenin_c(ji,jj) 384 ENDIF 385 IF( med_diag%IBEN_SI%dgsave ) THEN 386 iben_si2d(ji,jj) = f_fbenin_si(ji,jj) 387 ENDIF 388 IF( med_diag%IBEN_CA%dgsave ) THEN 389 iben_ca2d(ji,jj) = f_fbenin_ca(ji,jj) 390 ENDIF 391 IF( med_diag%OBEN_N%dgsave ) THEN 392 oben_n2d(ji,jj) = f_benout_n(ji,jj) 393 ENDIF 394 IF( med_diag%OBEN_FE%dgsave ) THEN 395 oben_fe2d(ji,jj) = f_benout_fe(ji,jj) 396 ENDIF 397 IF( med_diag%OBEN_C%dgsave ) THEN 398 oben_c2d(ji,jj) = f_benout_c(ji,jj) 399 ENDIF 400 IF( med_diag%OBEN_SI%dgsave ) THEN 401 oben_si2d(ji,jj) = f_benout_si(ji,jj) 402 ENDIF 403 IF( med_diag%OBEN_CA%dgsave ) THEN 404 oben_ca2d(ji,jj) = f_benout_ca(ji,jj) 405 ENDIF 406 IF( med_diag%SFR_OCAL%dgsave ) THEN 407 sfr_ocal2d(ji,jj) = f3_omcal(ji,jj,jk) 408 ENDIF 409 IF( med_diag%SFR_OARG%dgsave ) THEN 410 sfr_oarg2d(ji,jj) = f3_omarg(ji,jj,jk) 411 ENDIF 412 IF( med_diag%LYSO_CA%dgsave ) THEN 413 lyso_ca2d(ji,jj) = f_benout_lyso_ca(ji,jj) 372 IF (jk.eq.mbathy(ji,jj)) THEN 373 IF (tmask(ji,jj,jk) == 1) THEN 374 IF( med_diag%IBEN_N%dgsave ) THEN 375 iben_n2d(ji,jj) = f_sbenin_n(ji,jj) + & 376 f_fbenin_n(ji,jj) 377 ENDIF 378 IF( med_diag%IBEN_FE%dgsave ) THEN 379 iben_fe2d(ji,jj) = f_sbenin_fe(ji,jj) + & 380 f_fbenin_fe(ji,jj) 381 ENDIF 382 IF( med_diag%IBEN_C%dgsave ) THEN 383 iben_c2d(ji,jj) = f_sbenin_c(ji,jj) + & 384 f_fbenin_c(ji,jj) 385 ENDIF 386 IF( med_diag%IBEN_SI%dgsave ) THEN 387 iben_si2d(ji,jj) = f_fbenin_si(ji,jj) 388 ENDIF 389 IF( med_diag%IBEN_CA%dgsave ) THEN 390 iben_ca2d(ji,jj) = f_fbenin_ca(ji,jj) 391 ENDIF 392 IF( med_diag%OBEN_N%dgsave ) THEN 393 oben_n2d(ji,jj) = f_benout_n(ji,jj) 394 ENDIF 395 IF( med_diag%OBEN_FE%dgsave ) THEN 396 oben_fe2d(ji,jj) = f_benout_fe(ji,jj) 397 ENDIF 398 IF( med_diag%OBEN_C%dgsave ) THEN 399 oben_c2d(ji,jj) = f_benout_c(ji,jj) 400 ENDIF 401 IF( med_diag%OBEN_SI%dgsave ) THEN 402 oben_si2d(ji,jj) = f_benout_si(ji,jj) 403 ENDIF 404 IF( med_diag%OBEN_CA%dgsave ) THEN 405 oben_ca2d(ji,jj) = f_benout_ca(ji,jj) 406 ENDIF 407 IF( med_diag%SFR_OCAL%dgsave ) THEN 408 sfr_ocal2d(ji,jj) = f3_omcal(ji,jj,jk) 409 ENDIF 410 IF( med_diag%SFR_OARG%dgsave ) THEN 411 sfr_oarg2d(ji,jj) = f3_omarg(ji,jj,jk) 412 ENDIF 413 IF( med_diag%LYSO_CA%dgsave ) THEN 414 lyso_ca2d(ji,jj) = f_benout_lyso_ca(ji,jj) 415 ENDIF 414 416 ENDIF 415 417 ENDIF -
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_med_diag_trc.F90
r8051 r8062 529 529 530 530 531 if (jk.eq.mbathy(ji,jj)) then532 DO j j = 2,jpjm1533 DO ji = 2,jpim1534 IF (tmask(ji,jj,jk) == 1) THEN531 DO jj = 2,jpjm1 532 DO ji = 2,jpim1 533 IF (tmask(ji,jj,jk) == 1) THEN 534 if (jk.eq.mbathy(ji,jj)) then 535 535 !! N sedimentation flux 536 536 trc2d(ji,jj,69) = fsedn(ji,jj) … … 543 543 !! Ca sedimentation flux 544 544 trc2d(ji,jj,73) = fsedca(ji,jj) 545 ENDIF546 END DO547 ENDDO 548 endif545 endif 546 ENDIF 547 ENDDO 548 ENDDO 549 549 550 550 if (jk.eq.1) then … … 729 729 endif 730 730 731 if (jk .eq. mbathy(ji,jj)) then732 DO j j = 2,jpjm1733 DO ji = 2,jpim1734 IF (tmask(ji,jj,jk) == 1) THEN731 DO jj = 2,jpjm1 732 DO ji = 2,jpim1 733 IF (tmask(ji,jj,jk) == 1) THEN 734 if (jk .eq. mbathy(ji,jj)) then 735 735 !! seafloor omega calcite 736 736 trc2d(ji,jj,109) = f3_omcal(ji,jj,jk) 737 737 !! seafloor omega aragonite 738 738 trc2d(ji,jj,110) = f3_omarg(ji,jj,jk) 739 ENDIF740 END DO741 ENDDO 742 endif739 endif 740 ENDIF 741 ENDDO 742 ENDDO 743 743 744 744 if (jk.eq.i0100) then … … 779 779 endif 780 780 781 if (jk.eq.mbathy(ji,jj)) then782 DO j j = 2,jpjm1783 DO ji = 2,jpim1784 IF (tmask(ji,jj,jk) == 1) THEN781 DO jj = 2,jpjm1 782 DO ji = 2,jpim1 783 IF (tmask(ji,jj,jk) == 1) THEN 784 if (jk.eq.mbathy(ji,jj)) then 785 785 !! AXY (18/01/12): benthic flux diagnostics 786 786 trc2d(ji,jj,121) = f_sbenin_n(ji,jj) + f_fbenin_n(ji,jj) … … 794 794 trc2d(ji,jj,129) = f_benout_si(ji,jj) 795 795 trc2d(ji,jj,130) = f_benout_ca(ji,jj) 796 ENDIF797 END DO798 ENDDO 799 endif796 endif 797 ENDIF 798 ENDDO 799 ENDDO 800 800 801 801 DO jj = 2,jpjm1 … … 876 876 !! oxygen inventory 877 877 trc2d(ji,jj,148) = trc2d(ji,jj,148) + ftot_o2(ji,jj) 878 ENDIF 879 ENDDO 880 ENDDO 881 882 if (jk.eq.mbathy(ji,jj)) then 883 DO jj = 2,jpjm1 884 DO ji = 2,jpim1 885 IF (tmask(ji,jj,jk) == 1) THEN 878 879 if (jk.eq.mbathy(ji,jj)) then 886 880 trc2d(ji,jj,149) = f_benout_lyso_ca(ji,jj) 887 ENDIF 888 ENDDO 889 ENDDO 890 endif 891 892 DO jj = 2,jpjm1 893 DO ji = 2,jpim1 894 IF (tmask(ji,jj,jk) == 1) THEN 881 endif 882 895 883 !! community respiration 896 884 trc2d(ji,jj,150) = fcomm_resp(ji,jj) * fse3t(ji,jj,jk) -
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_diag.F90
r8025 r8062 78 78 DO jj = 2,jpjm1 79 79 DO ji = 2,jpim1 80 if (tmask(ji,jj, 1) == 1) then80 if (tmask(ji,jj,jk) == 1) then 81 81 !! sum tracers for inventory checks 82 82 IF ( med_diag%INVTN%dgsave ) THEN … … 105 105 DO jj = 2,jpjm1 106 106 DO ji = 2,jpim1 107 if (tmask(ji,jj, 1) == 1) then107 if (tmask(ji,jj,jk) == 1) then 108 108 IF ( med_diag%INVTC%dgsave ) THEN 109 109 ftot_c(ji,jj) = ftot_c(ji,jj) + & … … 128 128 DO jj = 2,jpjm1 129 129 DO ji = 2,jpim1 130 if (tmask(ji,jj, 1) == 1) then130 if (tmask(ji,jj,jk) == 1) then 131 131 IF ( med_diag%INVTC%dgsave ) THEN 132 132 !! … … 155 155 DO jj = 2,jpjm1 156 156 DO ji = 2,jpim1 157 if (tmask(ji,jj, 1) == 1) then157 if (tmask(ji,jj,jk) == 1) then 158 158 IF ( med_diag%O2min%dgsave ) THEN 159 159 if ( zoxy(ji,jj) < o2min(ji,jj) ) then -
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/detritus_fast_sink.F90
r8039 r8062 116 116 DO jj = 2,jpjm1 117 117 DO ji = 2,jpim1 118 IF (tmask(ji,jj, 1) == 1) THEN118 IF (tmask(ji,jj,jk) == 1) THEN 119 119 !! 120 120 !! this variable integrates the creation of slow sinking … … 167 167 DO jj = 2,jpjm1 168 168 DO ji = 2,jpim1 169 IF (tmask(ji,jj, 1) == 1) THEN169 IF (tmask(ji,jj,jk) == 1) THEN 170 170 !! 171 171 !! nitrogen -
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/iron_chem_scav.F90
r8045 r8062 87 87 DO ji = 2,jpim1 88 88 !! OPEN wet point IF..THEN loop 89 if (tmask(ji,jj, 1) == 1) then89 if (tmask(ji,jj,jk) == 1) then 90 90 !! 91 91 !! total iron concentration (mmol Fe / m3 -> umol Fe / m3) -
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/phytoplankton.F90
r8039 r8062 121 121 DO jj = 2,jpjm1 122 122 DO ji = 2,jpim1 123 if (tmask(ji,jj, 1) == 1) then123 if (tmask(ji,jj,jk) == 1) then 124 124 !!---------------------------------------------------------- 125 125 !! Phytoplankton light limitation … … 197 197 DO jj = 2,jpjm1 198 198 DO ji = 2,jpim1 199 if (tmask(ji,jj, 1) == 1) then199 if (tmask(ji,jj,jk) == 1) then 200 200 !!---------------------------------------------------------- 201 201 !! Phytoplankton nutrient limitation … … 233 233 DO jj = 2,jpjm1 234 234 DO ji = 2,jpim1 235 if (tmask(ji,jj, 1) == 1) then235 if (tmask(ji,jj,jk) == 1) then 236 236 !!---------------------------------------------------------- 237 237 !! Primary production (non-diatoms) … … 254 254 DO jj = 2,jpjm1 255 255 DO ji = 2,jpim1 256 if (tmask(ji,jj, 1) == 1) then256 if (tmask(ji,jj,jk) == 1) then 257 257 !!---------------------------------------------------------- 258 258 !! Primary production (diatoms) … … 348 348 DO jj = 2,jpjm1 349 349 DO ji = 2,jpim1 350 if (tmask(ji,jj, 1) == 1) then350 if (tmask(ji,jj,jk) == 1) then 351 351 !!---------------------------------------------------------- 352 352 !! Mixed layer primary production … … 379 379 DO jj = 2,jpjm1 380 380 DO ji = 2,jpim1 381 if (tmask(ji,jj, 1) == 1) then381 if (tmask(ji,jj,jk) == 1) then 382 382 !!---------------------------------------------------------- 383 383 !! Vertical Integral -- … … 407 407 DO jj = 2,jpjm1 408 408 DO ji = 2,jpim1 409 if (tmask(ji,jj, 1) == 1) then409 if (tmask(ji,jj,jk) == 1) then 410 410 !!---------------------------------------------------------- 411 411 !! More chlorophyll calculations -
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/plankton.F90
r8012 r8062 84 84 DO jj = 2,jpjm1 85 85 DO ji = 2,jpim1 86 if (tmask(ji,jj, 1) == 1) then86 if (tmask(ji,jj,jk) == 1) then 87 87 !!---------------------------------------------------------- 88 88 !! Plankton mortality losses … … 112 112 DO jj = 2,jpjm1 113 113 DO ji = 2,jpim1 114 if (tmask(ji,jj, 1) == 1) then114 if (tmask(ji,jj,jk) == 1) then 115 115 !! 116 116 !! diatom phytoplankton … … 135 135 DO jj = 2,jpjm1 136 136 DO ji = 2,jpim1 137 if (tmask(ji,jj, 1) == 1) then137 if (tmask(ji,jj,jk) == 1) then 138 138 !! 139 139 !! microzooplankton … … 157 157 DO jj = 2,jpjm1 158 158 DO ji = 2,jpim1 159 if (tmask(ji,jj, 1) == 1) then159 if (tmask(ji,jj,jk) == 1) then 160 160 !! 161 161 !! mesozooplankton -
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r8051 r8062 432 432 DO jj = 2,jpjm1 433 433 DO ji = 2,jpim1 434 if (tmask(ji,jj, 1) == 1) then434 if (tmask(ji,jj,jk) == 1) then 435 435 !! microzooplankton 436 436 zzmi(ji,jj) = max(0.,trn(ji,jj,jk,jpzmi)) … … 452 452 DO jj = 2,jpjm1 453 453 DO ji = 2,jpim1 454 if (tmask(ji,jj, 1) == 1) then454 if (tmask(ji,jj,jk) == 1) then 455 455 !! detrital carbon 456 456 zdtc(ji,jj) = max(0.,trn(ji,jj,jk,jpdtc)) … … 495 495 DO jj = 2,jpjm1 496 496 DO ji = 2,jpim1 497 if (tmask(ji,jj, 1) == 1) then497 if (tmask(ji,jj,jk) == 1) then 498 498 !! implicit detrital carbon 499 499 zdtc(ji,jj) = zdet(ji,jj) * xthetad … … 505 505 DO jj = 2,jpjm1 506 506 DO ji = 2,jpim1 507 if (tmask(ji,jj, 1) == 1) then507 if (tmask(ji,jj,jk) == 1) then 508 508 if (idf.eq.1) then 509 509 !! AXY (15/01/10) … … 560 560 ! DO jj = 2,jpjm1 561 561 ! DO ji = 2,jpim1 562 ! if (tmask(ji,jj, 1) == 1) then562 ! if (tmask(ji,jj,jk) == 1) then 563 563 ! !! report state variable values 564 564 ! IF (lwp) write (numout,*) & … … 592 592 ! DO jj = 2,jpjm1 593 593 ! DO ji = 2,jpim1 594 ! if (tmask(ji,jj, 1) == 1) then594 ! if (tmask(ji,jj,jk) == 1) then 595 595 ! IF (lwp) write (numout,*) & 596 596 ! '------------------------------' -
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/zooplankton.F90
r8039 r8062 122 122 DO jj = 2,jpjm1 123 123 DO ji = 2,jpim1 124 if (tmask(ji,jj, 1) == 1) then124 if (tmask(ji,jj,jk) == 1) then 125 125 !! 126 126 !! which translates to these incoming N and C fluxes … … 175 175 DO jj = 2,jpjm1 176 176 DO ji = 2,jpim1 177 if (tmask(ji,jj, 1) == 1) then177 if (tmask(ji,jj,jk) == 1) then 178 178 !!---------------------------------------------------------- 179 179 !! Mesozooplankton second … … 230 230 DO jj = 2,jpjm1 231 231 DO ji = 2,jpim1 232 if (tmask(ji,jj, 1) == 1) then232 if (tmask(ji,jj,jk) == 1) then 233 233 !! 234 234 !! the ideal food C:N ratio for mesozooplankton … … 281 281 DO jj = 2,jpjm1 282 282 DO ji = 2,jpim1 283 if (tmask(ji,jj, 1) == 1) then283 if (tmask(ji,jj,jk) == 1) then 284 284 fzmi_i(ji,jj) = fzmi_i(ji,jj) + fse3t(ji,jj,jk) * & 285 285 ( fgmipn(ji,jj) + fgmid(ji,jj) )
Note: See TracChangeset
for help on using the changeset viewer.