Changeset 8457 for branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_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/trcnam_medusa.F90 (modified) (5 diffs)

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