Changeset 8075 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

Timestamp:

2017-05-25T18:58:42+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- update CFCs - add SF6 and update gas transfert param

File:

• branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90 (modified) (10 diffs)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

@@ -7 +7 @@
    !!  NEMO      1.0  !  2004-03  (C. Ethe) free form + modularity
    !!            2.0  !  2007-12  (C. Ethe, G. Madec)  reorganisation
+   !!                 !  2016-06  (J. Palmieri)  update for UKESM1
+   !!                 !  2017-04  (A. Yool)  update to add SF6, fix coefficients
    !!----------------------------------------------------------------------
 #if defined key_cfc
@@ -43 +45 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   patm     ! atmospheric function
 
-   REAL(wp), DIMENSION(4,2) ::   soa   ! coefficient for solubility of CFC [mol/l/atm]
-   REAL(wp), DIMENSION(3,2) ::   sob   !    "               "
-   REAL(wp), DIMENSION(4,2) ::   sca   ! coefficients for schmidt number in degre Celcius
+   REAL(wp), DIMENSION(4,3) ::   soa   ! coefficient for solubility of CFC [mol/l/atm]
+   REAL(wp), DIMENSION(3,3) ::   sob   !    "               "
+   REAL(wp), DIMENSION(5,3) ::   sca   ! coefficients for schmidt number in degre Celcius
       
    !                          ! coefficients for conversion
@@ -85 +87 @@
       INTEGER  ::   im1, im2, ierr
       REAL(wp) ::   ztap, zdtap        
-      REAL(wp) ::   zt1, zt2, zt3, zv2
+      REAL(wp) ::   zt1, zt2, zt3, zt4, zv2
       REAL(wp) ::   zsol      ! solubility
       REAL(wp) ::   zsch      ! schmidt number 
@@ -106 +108 @@
       ! Temporal interpolation
       ! ----------------------
-      !! JPALM -- 15-06-2016 -- define 2 kind of CFC run.
-      !!                     we want to make cycle experiments, 
-      !!                     to periodically compare the ocean dynamic within
-      !!                     1- the SPIN-UP and 2- Hincast/Projections
+      !! JPALM -- 15-06-2016 -- define 2 kinds of CFC run:
+      !!                     1- the SPIN-UP and 2- Hindcast/Projections
       !!                     -- main difference is the way to define the year of
       !!                     simulation, that determine the atm pCFC.
       !!                     1-- Spin-up: our atm forcing is of 30y we cycle on.
       !!                     So we do 90y CFC cycles to be in good
-      !!                     correspondance with the atmosphere
+      !!                     correspondence with the atmosphere
       !!                     2-- Hindcast/proj, instead of nyear-1900 we keep
       !!                     the 2 last digit, and enable 3 cycle from 1800 to 2100.  
       !!----------------------------------------------------------------------
-      !! 1 -- SPIN-UP
       IF (simu_type==1) THEN
+         !! 1 -- SPIN-UP
          iyear_tmp = nyear - nyear_res  !! JPALM -- in our spin-up, nyear_res is 1000
          iyear_beg = MOD( iyear_tmp , 90 )
@@ -133 +133 @@
             iyear_beg = iyear_beg + 30
          ENDIF
-      !!
-      !! 2 -- Hindcast/proj
       ELSEIF (simu_type==2) THEN
+         !! 2 -- Hindcast/proj
          iyear_beg = MOD(nyear, 100)
          IF (iyear_beg < 9)  iyear_beg = iyear_beg + 100
@@ -191 +190 @@
                zt2  = zt1 * zt1 
                zt3  = zt1 * zt2
-               zsch = sca(1,jl) + sca(2,jl) * zt1 + sca(3,jl) * zt2 + sca(4,jl) * zt3
+               zt4  = zt1 * zt3
+               zsch = sca(1,jl) + sca(2,jl) * zt1 + sca(3,jl) * zt2 + sca(4,jl) * zt3 + sca(5,jl) * zt4
 
                !    speed transfert : formulae of wanninkhof 1992
                zv2     = wndm(ji,jj) * wndm(ji,jj)
                zsch    = zsch / 660.
-               zak_cfc = ( 0.39 * xconv2 * zv2 / SQRT(zsch) ) * tmask(ji,jj,1)
+               ! AXY (25/04/17): OMIP protocol specifies lower Wanninkhof (2014) value
+               ! zak_cfc = ( 0.39 * xconv2 * zv2 / SQRT(zsch) ) * tmask(ji,jj,1)
+               zak_cfc = ( 0.251 * xconv2 * zv2 / SQRT(zsch) ) * tmask(ji,jj,1)
 
                ! Input function  : speed *( conc. at equil - concen at surface )
@@ -226 +228 @@
                elseif (jl.EQ.2) then
                    WRITE(NUMOUT,*) 'Traceur = CFC12: '
+               elseif (jl.EQ.3) then
+                   WRITE(NUMOUT,*) 'Traceur = SF6: '
                endif
              WRITE(NUMOUT,*) 'nyear    = ', nyear
@@ -256 +260 @@
          CALL iom_put( "qtrCFC11"  , qtr_cfc (:,:,1) )
          CALL iom_put( "qintCFC11" , qint_cfc(:,:,1) )
+         CALL iom_put( "qtrCFC12"  , qtr_cfc (:,:,2) )
+         CALL iom_put( "qintCFC12" , qint_cfc(:,:,2) )
+         CALL iom_put( "qtrSF6"    , qtr_cfc (:,:,3) )
+         CALL iom_put( "qintSF6"   , qint_cfc(:,:,4) )
       ELSE
          IF( ln_diatrc ) THEN
             trc2d(:,:,jp_cfc0_2d    ) = qtr_cfc (:,:,1)
             trc2d(:,:,jp_cfc0_2d + 1) = qint_cfc(:,:,1)
+            trc2d(:,:,jp_cfc0_2d + 2) = qtr_cfc (:,:,2)
+            trc2d(:,:,jp_cfc0_2d + 3) = qint_cfc(:,:,2)
+            trc2d(:,:,jp_cfc0_2d + 4) = qtr_cfc (:,:,3)
+            trc2d(:,:,jp_cfc0_2d + 5) = qint_cfc(:,:,4)
          END IF
       END IF
@@ -293 +305 @@
       soa(2,1) =  319.6552
       soa(3,1) =  119.4471
-      soa(4,1) =  -1.39165
-
-      sob(1,1) =  -0.142382
-      sob(2,1) =   0.091459
-      sob(3,1) =  -0.0157274
-
-      ! Schmidt number 
-      sca(1,1) = 3501.8
-      sca(2,1) = -210.31
-      sca(3,1) =  6.1851
-      sca(4,1) = -0.07513
+      soa(4,1) =   -1.39165
+
+      sob(1,1) = -0.142382
+      sob(2,1) =  0.091459
+      sob(3,1) = -0.0157274
+
+      ! Schmidt number          AXY (25/04/17)
+      sca(1,1) = 3579.2       ! = 3501.8
+      sca(2,1) = -222.63      ! = -210.31
+      sca(3,1) =    7.5749    ! =    6.1851
+      sca(4,1) =   -0.14595   ! =   -0.07513
+      sca(5,1) =    0.0011874 ! = absent
 
       ! coefficient for CFC12 
@@ -312 +325 @@
       soa(2,2) =  298.9702
       soa(3,2) =  113.8049
-      soa(4,2) =  -1.39165
-
-      sob(1,2) =  -0.143566
-      sob(2,2) =   0.091015
-      sob(3,2) =  -0.0153924
-
-      ! schmidt number 
-      sca(1,2) =  3845.4 
-      sca(2,2) =  -228.95
-      sca(3,2) =  6.1908 
-      sca(4,2) =  -0.067430
+      soa(4,2) =   -1.39165
+
+      sob(1,2) = -0.143566
+      sob(2,2) =  0.091015
+      sob(3,2) = -0.0153924
+
+      ! schmidt number         AXY (25/04/17)
+      sca(1,2) = 3828.1      ! = 3845.4 
+      sca(2,2) = -249.86     ! = -228.95
+      sca(3,2) =    8.7603   ! =    6.1908 
+      sca(4,2) =   -0.1716   ! =   -0.067430
+      sca(5,2) =    0.001408 ! = absent
+
+      ! coefficients for SF6   AXY (25/04/17)
+      !---------------------
+      
+      ! Solubility
+      soa(1,3) =  -80.0343
+      soa(2,3) =  117.232
+      soa(3,3) =   29.5817
+      soa(4,3) =    0.0
+
+      sob(1,3) =  0.0335183
+      sob(2,3) = -0.0373942
+      sob(3,3) =  0.00774862
+
+      ! Schmidt number
+      sca(1,3) = 3177.5
+      sca(2,3) = -200.57
+      sca(3,3) =    6.8865
+      sca(4,3) =   -0.13335
+      sca(5,3) =    0.0010877
 
       !!---------------------------------------------