Changeset 8457 for branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

Timestamp:

2017-08-23T15:27:51+02:00 (7 years ago)

Author:

dford

Message:

Option to read in a value for atmospheric xCO2.

File:

• branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

@@ -1320 +1320 @@
          f_xco2a = fq4
       endif
+#  if defined key_foam_medusa
+      IF ( xobs_xco2a > 0.0 ) THEN
+         IF(lwp) WRITE(numout,*) ' using observed atm pCO2 = ', xobs_xco2a
+         f_xco2a = xobs_xco2a
+      ELSE
+         IF(lwp) WRITE(numout,*) ' xobs_xco2a <= 0 so using default atm pCO2'
+      ENDIF
+#  endif
 #  if defined key_axy_pi_co2
       !! f_xco2a = 284.725       !! CMIP5 pre-industrial pCO2