Changeset 8457


Ignore:
Timestamp:
2017-08-23T15:27:51+02:00 (3 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:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90

    r8436 r8457  
    390390      & 926.4960, 935.7040 / 
    391391# endif 
     392# if defined key_foam_medusa 
     393   REAL(wp) ::   xobs_xco2a   !: Observed atmospheric xCO2, read in 
     394# endif 
    392395#endif 
    393396 
  • branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

    r8436 r8457  
    13201320         f_xco2a = fq4 
    13211321      endif 
     1322#  if defined key_foam_medusa 
     1323      IF ( xobs_xco2a > 0.0 ) THEN 
     1324         IF(lwp) WRITE(numout,*) ' using observed atm pCO2 = ', xobs_xco2a 
     1325         f_xco2a = xobs_xco2a 
     1326      ELSE 
     1327         IF(lwp) WRITE(numout,*) ' xobs_xco2a <= 0 so using default atm pCO2' 
     1328      ENDIF 
     1329#  endif 
    13221330#  if defined key_axy_pi_co2 
    13231331      !! 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

    r8131 r8457  
    8888#if defined key_roam 
    8989      NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit,        & 
     90# if defined key_foam_medusa 
     91      &    xobs_xco2a,                                        & 
     92# endif 
    9093      &    xthetarem,xo2min  
    9194#endif 
     
    569572         &   ' key_debug_medusa                                                       = INACTIVE' 
    570573# endif 
     574#if defined key_foam_medusa 
     575         WRITE(numout,*)     & 
     576         &   ' key_foam_medusa                                                        = ACTIVE' 
     577#else 
     578         WRITE(numout,*)     & 
     579         &   ' key_foam_medusa                                                        = INACTIVE' 
     580#endif 
    571581         WRITE(numout,*) ' ' 
    572582 
     
    10851095      xthetarem = 0. 
    10861096      xo2min    = 0. 
     1097# if defined key_foam_medusa 
     1098      xobs_xco2a = 0. 
     1099# endif 
    10871100 
    10881101      !READ(numnatm,natroam) 
     
    11041117!!       xthetarem :  oxygen consumption by carbon remineralisation 
    11051118!!       xo2min    :  oxygen minimum concentration 
     1119# if defined key_foam_medusa 
     1120!!       xobs_xco2a : observed atmospheric xCO2 
     1121# endif 
    11061122 
    11071123      IF(lwp) THEN 
     
    11211137          WRITE(numout,*)     & 
    11221138          &   ' oxygen minimum concentration                               xo2min      = ', xo2min 
     1139# if defined key_foam_medusa 
     1140          WRITE(numout,*)     & 
     1141          &   ' observed atmospheric xCO2                                  xobs_xco2a  = ', xobs_xco2a 
     1142# endif 
    11231143       ENDIF 
    11241144 
Note: See TracChangeset for help on using the changeset viewer.