Changeset 10045 for branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90

Timestamp:

2018-08-07T14:19:25+02:00 (6 years ago)

Author:

jpalmier

Message:

Andrew's changes to add the OMIP double_DIC (activated with key_omip_dic)

File:

• branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90 (modified) (6 diffs)

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

@@ -8 +8 @@
    !!   -   ! 2017-08 (A. Yool)            Amend slow-detritus bug
    !!   -   ! 2017-08 (A. Yool)            Reformatting for clarity
+   !!   -   ! 2018-08 (A. Yool)            add OMIP preindustrial DIC
    !!----------------------------------------------------------------------
 #if defined key_medusa
@@ -38 +39 @@
                                    f_benout_si,                              &
                                    f_co2flux, f_o2flux,                      &
+# if defined key_omip_dic
+                                   f_pi_co2flux,                             & ! AXY (06/08/18)
+# endif                                   
                                    f_riv_loc_alk, f_riv_loc_c,               &
                                    f_riv_loc_n, f_riv_loc_si,                & 
@@ -395 +399 @@
                fcar_prod(ji,jj) = fcar_prod(ji,jj) + ( fse3t(ji,jj,jk) *     &
                                                        fc_prod )
+#  if defined key_omip_dic
+               !! AXY (06/08/18): OMIP PI DIC has the same BGC fluxes as
+               !!                 normal DIC, with the exception of its
+               !!                 air-sea exchange; see below
+               btra(ji,jj,jpomd_lc) = btra(ji,jj,jpdic_lc)
+#  endif               
             ENDIF
          ENDDO
@@ -573 +583 @@
                   btra(ji,jj,jpdic_lc) = btra(ji,jj,jpdic_lc) + (b0 *              &
                                                            f_co2flux(ji,jj))
+#  if defined key_omip_dic
+                  !! AXY (06/08/18): air-sea CO2 flux is the only difference
+                  !!                 between DIC and OMIP PI DIC tracers
+                  btra(ji,jj,jpomd_lc) = btra(ji,jj,jpomd_lc) + (b0 *              &
+                                                           f_pi_co2flux(ji,jj))                  
+#  endif                   
                   !!
                   !! O2 flux (mol/m3/s -> mmol/m3/d)
@@ -624 +640 @@
                   IF (lwp) write (numout,*) 'btra(ji,jj,jpoxy_lc)(',jk,')  = ', &
                                             btra(ji,jj,jpoxy_lc)
+#   if defined key_omip_dic
+                  IF (lwp) write (numout,*) 'btra(ji,jj,jpomd_lc)(',jk,')  = ', &
+                                            btra(ji,jj,jpomd_lc)
+#   endif
 #  endif
                endif
@@ -731 +751 @@
                   tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) +                &
                                         (btra(ji,jj,jpoxy_lc) / 86400.)
+#  if defined key_omip_dic
+                  tra(ji,jj,jk,jpomd) = tra(ji,jj,jk,jpomd) +                &
+                                        (btra(ji,jj,jpomd_lc) / 86400.)
+#  endif
 # endif
                endif