Changeset 11726 for branches/NERC/dev_r5518_GO6_GeoAlk_MEDUSA/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90

Timestamp:

2019-10-18T16:49:33+02:00 (5 years ago)

Author:

jpalmier

Message:

add Changes to enable Alkalinity geo-ingeneering tests within MEDUSA

File:

• branches/NERC/dev_r5518_GO6_GeoAlk_MEDUSA/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90 (modified) (7 diffs)

branches/NERC/dev_r5518_GO6_GeoAlk_MEDUSA/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90

@@ -71 +71 @@
 # endif
                                    zpds, zphd, zphn
-      USE dom_oce,           ONLY: e3t_0, gphit, mbathy, tmask
+      USE dom_oce,           ONLY: e3t_0, gphit, mbathy, tmask, nyear          ! AXY (16/10/19)
+
 # if defined key_vvl
       USE dom_oce,           ONLY: e3t_n
@@ -232 +233 @@
                !! 
                !! riverine flux
+               !! (NOTE: this is "protected" for near-surface addition of riverine
+               !! by array friver_dep which is zero below a defined model level)
                if ( jriver_n .gt. 0 ) then
                   f_riv_loc_n(ji,jj) = f_riv_n(ji,jj) *                      &
@@ -276 +279 @@
                !! 
                !! riverine flux
+               !! (NOTE: this is "protected" for near-surface addition of riverine
+               !! by array friver_dep which is zero below a defined model level)
                if ( jriver_si .gt. 0 ) then
                   f_riv_loc_si(ji,jj) = f_riv_si(ji,jj) *                    &
@@ -362 +367 @@
                !! 
                !! riverine flux
+               !! (NOTE: this is "protected" for near-surface addition of riverine
+               !! by array friver_dep which is zero below a defined model level)
                if ( jriver_c .gt. 0 ) then
                   f_riv_loc_c(ji,jj) = f_riv_c(ji,jj) *                      &
@@ -387 +394 @@
                !! 
                !! riverine flux
-               if ( jk .eq. 1 .and. jriver_c .gt. 0 ) then
-                  fc_prod = fc_prod + f_riv_c(ji,jj)
-               endif
+               !! (NOTE: the lines immediately below are incorrect and riverine
+               !! flux is already handled above)
+               !! if ( jk .eq. 1 .and. jriver_c .gt. 0 ) then
+               !!    fc_prod = fc_prod + f_riv_c(ji,jj)
+               !! endif
                !!
                btra(ji,jj,jpdic_lc) = b0 * ( fc_prod + fc_cons )
@@ -414 +423 @@
                !! 
                !! riverine flux
+               !! (NOTE: this is "protected" for near-surface addition of riverine
+               !! by array friver_dep which is zero below a defined model level)
                if ( jriver_alk .gt. 0 ) then
                   f_riv_loc_alk(ji,jj) = f_riv_alk(ji,jj) *                  &
@@ -419 +430 @@
                                          fse3t(ji,jj,jk)
                   fa_prod = fa_prod + f_riv_loc_alk(ji,jj)
+               endif
+          !!
+               !! AXY (16/10/19)
+               !! geoengineering flux
+               !! (NOTE: this is only added to the surface ocean grid cell, only
+               !! within a particular time window, and only within a particular
+               !! sub-polar geographical range; these parameters could be passed up
+               !! to ROSE if they need different values for different experiments)
+          if ( jk .eq. 1 .and. nyear .gt. 2014 .and. nyear .lt. 2050 .and. &
+                   gphit(ji,jj) .ge .-60.0 .and. gphit(ji,jj) .le. 70.0 ) then
+                  f_geoeng_alk(ji,jj) = 1.17 / fse3t(ji,jj,jk)
+                  fa_prod = fa_prod + f_geoeng_alk(ji,jj)
                endif
                !!