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/trcbio_medusa.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/trcbio_medusa.F90

    r8051 r8062  
    432432         DO jj = 2,jpjm1 
    433433            DO ji = 2,jpim1 
    434                if (tmask(ji,jj,1) == 1) then 
     434               if (tmask(ji,jj,jk) == 1) then 
    435435                  !! microzooplankton 
    436436                  zzmi(ji,jj) = max(0.,trn(ji,jj,jk,jpzmi)) 
     
    452452         DO jj = 2,jpjm1 
    453453            DO ji = 2,jpim1 
    454                if (tmask(ji,jj,1) == 1) then 
     454               if (tmask(ji,jj,jk) == 1) then 
    455455                  !! detrital carbon 
    456456                  zdtc(ji,jj) = max(0.,trn(ji,jj,jk,jpdtc)) 
     
    495495         DO jj = 2,jpjm1 
    496496            DO ji = 2,jpim1 
    497                if (tmask(ji,jj,1) == 1) then 
     497               if (tmask(ji,jj,jk) == 1) then 
    498498                  !! implicit detrital carbon 
    499499                  zdtc(ji,jj) = zdet(ji,jj) * xthetad 
     
    505505         DO jj = 2,jpjm1 
    506506            DO ji = 2,jpim1 
    507                if (tmask(ji,jj,1) == 1) then 
     507               if (tmask(ji,jj,jk) == 1) then 
    508508                  if (idf.eq.1) then 
    509509                     !! AXY (15/01/10) 
     
    560560!            DO jj = 2,jpjm1 
    561561!               DO ji = 2,jpim1 
    562 !                  if (tmask(ji,jj,1) == 1) then 
     562!                  if (tmask(ji,jj,jk) == 1) then 
    563563!                     !! report state variable values 
    564564!                     IF (lwp) write (numout,*)                               & 
     
    592592!            DO jj = 2,jpjm1 
    593593!               DO ji = 2,jpim1 
    594 !                  if (tmask(ji,jj,1) == 1) then 
     594!                  if (tmask(ji,jj,jk) == 1) then 
    595595!                     IF (lwp) write (numout,*)                               & 
    596596!                       '------------------------------' 
Note: See TracChangeset for help on using the changeset viewer.