Changeset 5841 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcco2_medusa.F90

Timestamp:

2015-10-30T12:48:06+01:00 (8 years ago)

Author:

jpalmier

Message:

JPALM --30-10-2015-- Add MOCSY and DMS to MEDUSA-NEMO3.6

File:

• branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcco2_medusa.F90 (modified) (13 diffs)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcco2_medusa.F90

@@ -47 +47 @@
 !=======================================================================
 !
-      SUBROUTINE trc_co2_medusa( Temp, Sal, DIC, ALK, Depth, Wnd, pCO2a, &
+      SUBROUTINE trc_co2_medusa( Temp, Sal, DIC, ALK, Depth, xkw, pCO2a, &
       pH, pCO2w, h2co3, hco3, co3, om_cal, om_arg, co2flux, TDIC, TALK,  &
       dcf, henry, iters )
@@ -72 +72 @@
 ! 17/02/2010. Update calculation of K1, K2, Kb to make consistant with the OCMIP protocols.
 ! 29/07/2011. Merged into MEDUSA with a raft of changes to this subroutine; less elsewhere
+! 23/06/2015. Modified to take gas transfer velocity as an input (rather than wind speed); 
+!             alter CO2 flux to /s rather than /d for consistency with other schemes
 !
 ! Changes for MEDUSA include: 
@@ -85 +87 @@
       REAL(wp), INTENT( in )    :: ALK        ! meq  / m3
       REAL(wp), INTENT( in )    :: Depth      ! m
-      REAL(wp), INTENT( in )    :: Wnd        ! m / s
+!     REAL(wp), INTENT( in )    :: Wnd        ! m / s
+      REAL(wp), INTENT( in )    :: xkw        ! m / s
       REAL(wp), INTENT( in )    :: pCO2a      ! uatm
 !----------------------------------------------------------------------
@@ -95 +98 @@
       REAL(wp), INTENT( inout ) :: om_cal     ! normalised
       REAL(wp), INTENT( inout ) :: om_arg     ! normalised
-      REAL(wp), INTENT( inout ) :: co2flux    ! mmol / m2 / d
+      REAL(wp), INTENT( inout ) :: co2flux    ! mmol / m2 / s
       REAL(wp), INTENT( inout ) :: TDIC       ! umol / kg
       REAL(wp), INTENT( inout ) :: TALK       ! ueq  / kg
@@ -129 +132 @@
    !                 (i.e. surface calculations being performed)
    if (Depth .eq. 0.0) then
-      call Air_sea_exchange( Temp, Wnd, pCO2w, pCO2a, henry, dcf, &         ! inputs
+      call Air_sea_exchange( Temp, xkw, pCO2w, pCO2a, henry, dcf, &         ! inputs
          co2flux )                                                          ! output
    else
@@ -145 +148 @@
       IF(lwp) WRITE(numout,*) ' trc_co2_medusa: zdic    =', DIC
       IF(lwp) WRITE(numout,*) ' trc_co2_medusa: zalk    =', ALK
-      IF(lwp) WRITE(numout,*) ' trc_co2_medusa: f_wind  =', Wnd
+      IF(lwp) WRITE(numout,*) ' trc_co2_medusa: f_kw660 =', xkw
       IF(lwp) WRITE(numout,*) ' trc_co2_medusa: f_ph    =', ph
       IF(lwp) WRITE(numout,*) ' trc_co2_medusa: f_pco2w =', pCO2w
@@ -162 +165 @@
       &        ' DIC', DIC, ' ALK', ALK
       if (lwp) write (numout,'(a,a,f10.3,a,f10.3)') 'CO2FLUX-NAN', &
-      &        ' WND', Wnd, ' PH ', ph
+      &        ' XKW', xkw, ' PH ', ph
       if (lwp) write (numout,'(a,a,i6)') 'CO2FLUX-NAN', &
       &        ' ITERS', iters
@@ -196 +199 @@
 !  WRITE(*,'(A27,F10.3)') "    Omega calcite       (~) = ", om_cal
 !  WRITE(*,'(A27,F10.3)') "    Omega aragonite     (~) = ", om_arg
-!  WRITE(*,'(A27,F10.3)') "    air sea flux(mmol/m2/d) = ", flux
+!  WRITE(*,'(A27,F10.3)') "    air sea flux(mmol/m2/s) = ", flux
 !  WRITE(*,*) " "
 
@@ -287 +290 @@
 !=======================================================================
 !
-      SUBROUTINE Air_sea_exchange( T, Wnd, pco2w, pco2a, henry, dcf, &
+      SUBROUTINE Air_sea_exchange( T, xkw, pco2w, pco2a, henry, dcf, &
       flux )
 !      
@@ -302 +305 @@
 !  pCO2a    partial pressure of CO2 in the atmosphere (usually external forcing).
 !  T        temperature (C)
-!  Wnd      wind speed, metres
+!  Wnd      wind speed, metres (DELETED)
+!  xkw      gas transfer velocity
 !  Henry    henry's constant
 !  density  the density of water for conversion between mmol/m3 and umol/kg
@@ -312 +316 @@
    IMPLICIT NONE
 
-      REAL(wp), INTENT( in )    :: T, wnd, pco2w, pco2a, henry, dcf ! INPUT PARAMETERS:
+      REAL(wp), INTENT( in )    :: T, xkw, pco2w, pco2a, henry, dcf ! INPUT PARAMETERS:
 !-----------------------------------------------------------------------
       REAL(wp), INTENT( inout ) :: flux                             ! OUTPUT Variables
@@ -320 +324 @@
 ! calculate the Schmidt number and unit conversions
           sc    = 2073.1-125.62*T+3.6276*T**2.0-0.0432190*T**3.0
-          fwind = (0.222d0 * wnd**2d0 + 0.333d0 * wnd)*(sc/660.d0)**(-0.5)
+!         fwind = (0.222d0 * wnd**2d0 + 0.333d0 * wnd)*(sc/660.d0)**(-0.5)
+          fwind = xkw * (sc/660.d0)**(-0.5)
           fwind = fwind*24.d0/100.d0   ! convert to m/day
 
@@ -326 +331 @@
 ! here it is rescaled to mmol/m2/d
           flux = fwind * henry * ( pco2a - pco2w ) * dcf
+
+! AXY (23/06/15): let's get it from /d to /s
+          flux = flux / ( 86400. )
 
   RETURN