Changeset 8457

Timestamp:

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

Author:

dford

Message:

Option to read in a value for atmospheric xCO2.

Location:

branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA

Files:

• sms_medusa.F90 (modified) (1 diff)
• trcbio_medusa.F90 (modified) (1 diff)
• trcnam_medusa.F90 (modified) (5 diffs)

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

@@ -390 +390 @@
       & 926.4960, 935.7040 /
 # endif
+# if defined key_foam_medusa
+   REAL(wp) ::   xobs_xco2a   !: Observed atmospheric xCO2, read in
+# endif
 #endif
 

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

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

@@ -88 +88 @@
 #if defined key_roam
       NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit,        &
+# if defined key_foam_medusa
+      &    xobs_xco2a,                                        &
+# endif
       &    xthetarem,xo2min 
 #endif
@@ -569 +572 @@
          &   ' key_debug_medusa                                                       = INACTIVE'
 # endif
+#if defined key_foam_medusa
+         WRITE(numout,*)     &
+         &   ' key_foam_medusa                                                        = ACTIVE'
+#else
+         WRITE(numout,*)     &
+         &   ' key_foam_medusa                                                        = INACTIVE'
+#endif
          WRITE(numout,*) ' '
 
@@ -1085 +1095 @@
       xthetarem = 0.
       xo2min    = 0.
+# if defined key_foam_medusa
+      xobs_xco2a = 0.
+# endif
 
       !READ(numnatm,natroam)
@@ -1104 +1117 @@
 !!       xthetarem :  oxygen consumption by carbon remineralisation
 !!       xo2min    :  oxygen minimum concentration
+# if defined key_foam_medusa
+!!       xobs_xco2a : observed atmospheric xCO2
+# endif
 
       IF(lwp) THEN
@@ -1121 +1137 @@
           WRITE(numout,*)     &
           &   ' oxygen minimum concentration                               xo2min      = ', xo2min
+# if defined key_foam_medusa
+          WRITE(numout,*)     &
+          &   ' observed atmospheric xCO2                                  xobs_xco2a  = ', xobs_xco2a
+# endif
        ENDIF