Changeset 10049

Timestamp:

2018-08-10T17:45:44+02:00 (5 years ago)

Author:

jpalmier

Message:

correct PI CO2 flux bug

File:

• branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90 (modified) (1 diff)

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

@@ -455 +455 @@
                !! mmol/m2/s -> mmol/m3/d; correct for sea-ice; divide 
                !! through by layer thickness
-               f_co2flux(ji,jj) = (1. - fr_i(ji,jj)) * f_co2flux(ji,jj) *    &
+               f_co2flux(ji,jj) = (1. - fr_i(ji,jj)) * f_co2flux(ji,jj) *       &
                                   86400. / fse3t(ji,jj,1)
-               !!
+#   if defined key_omip_dic
+               !! PI DIC - for omip runs
+               f_pi_co2flux(ji,jj) = (1. - fr_i(ji,jj)) * f_pi_co2flux(ji,jj) * &
+                                      86400. / fse3t(ji,jj,1)
+#   endif
                !! oxygen (O2); OCMIP-2 code
                !! AXY (23/06/15): amend input list for oxygen to account 
                !!                 for common gas transfer velocity
-               CALL trc_oxy_medusa(ztmp(ji,jj),zsal(ji,jj),f_kw660(ji,jj),   &
-                                   f_pp0(ji,jj),zoxy(ji,jj),                 &
-                                   f_kwo2_dum,f_o2flux(ji,jj),               &
+               CALL trc_oxy_medusa(ztmp(ji,jj),zsal(ji,jj),f_kw660(ji,jj),     &
+                                   f_pp0(ji,jj),zoxy(ji,jj),                   &
+                                   f_kwo2_dum,f_o2flux(ji,jj),                 &
                                    f_o2sat(ji,jj))
                !!
                !! mmol/m2/s -> mol/m3/d; correct for sea-ice; divide 
                !! through by layer thickness
-               f_o2flux(ji,jj)  = (1. - fr_i(ji,jj)) * f_o2flux(ji,jj) *     &
+               f_o2flux(ji,jj)  = (1. - fr_i(ji,jj)) * f_o2flux(ji,jj) *       &
                                   86400. / fse3t(ji,jj,1)
             ENDIF