Changeset 8062 for branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/phytoplankton.F90

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:

• branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/phytoplankton.F90 (modified) (7 diffs)

branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/phytoplankton.F90

@@ -121 +121 @@
       DO jj = 2,jpjm1
          DO ji = 2,jpim1
-            if (tmask(ji,jj,1) == 1) then
+            if (tmask(ji,jj,jk) == 1) then
                !!----------------------------------------------------------
                !! Phytoplankton light limitation
@@ -197 +197 @@
       DO jj = 2,jpjm1
          DO ji = 2,jpim1
-            if (tmask(ji,jj,1) == 1) then
+            if (tmask(ji,jj,jk) == 1) then
                !!----------------------------------------------------------
                !! Phytoplankton nutrient limitation
@@ -233 +233 @@
       DO jj = 2,jpjm1
          DO ji = 2,jpim1
-            if (tmask(ji,jj,1) == 1) then
+            if (tmask(ji,jj,jk) == 1) then
                !!----------------------------------------------------------
                !! Primary production (non-diatoms)
@@ -254 +254 @@
       DO jj = 2,jpjm1
          DO ji = 2,jpim1
-            if (tmask(ji,jj,1) == 1) then
+            if (tmask(ji,jj,jk) == 1) then
                !!----------------------------------------------------------
                !! Primary production (diatoms)
@@ -348 +348 @@
       DO jj = 2,jpjm1
          DO ji = 2,jpim1
-            if (tmask(ji,jj,1) == 1) then
+            if (tmask(ji,jj,jk) == 1) then
                !!----------------------------------------------------------
                !! Mixed layer primary production
@@ -379 +379 @@
       DO jj = 2,jpjm1
          DO ji = 2,jpim1
-            if (tmask(ji,jj,1) == 1) then
+            if (tmask(ji,jj,jk) == 1) then
                !!----------------------------------------------------------
                !! Vertical Integral --
@@ -407 +407 @@
       DO jj = 2,jpjm1
          DO ji = 2,jpim1
-            if (tmask(ji,jj,1) == 1) then
+            if (tmask(ji,jj,jk) == 1) then
                !!----------------------------------------------------------
                !! More chlorophyll calculations