New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8062 for branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_med_diag_iomput.F90 – NEMO

Ignore:
Timestamp:
2017-05-24T10:03:18+02:00 (7 years ago)
Author:
marc
Message:

Fixes to do with tmask(ji,jj,1) and mbathy(ji,jj) outside DO loop

File:
1 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  
    367367            ENDDO 
    368368         ENDDO 
    369       ELSE IF (jk.eq.mbathy(ji,jj)) THEN 
     369      ELSE 
    370370         DO jj = 2,jpjm1 
    371371            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 
    414416                  ENDIF 
    415417               ENDIF 
Note: See TracChangeset for help on using the changeset viewer.