Changeset 7625 for branches

Timestamp:

2017-01-31T13:06:14+01:00 (7 years ago)

Author:

dford

Message:

Update the retrieval of pCO2 from FABM.

File:

• branches/UKMO/dev_r4650_general_vert_coord_obsoper_surf_bgc/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90 (modified) (5 diffs)

branches/UKMO/dev_r4650_general_vert_coord_obsoper_surf_bgc/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

@@ -1423 +1423 @@
          & tmask, umask, vmask                            
       USE phycst, ONLY : &              ! Physical constants
-         & rday                         
+         & rday, &
+         & rt0
       USE oce, ONLY : &                 ! Ocean dynamics and tracers variables
          & tsn,  &             
@@ -1455 +1456 @@
       USE par_spm, ONLY: &              ! ERSEM/SPM sediments
          & jp_spm
+      USE trc, ONLY :  &
+         & trn
 #endif
       IMPLICIT NONE
@@ -1493 +1496 @@
 #if defined key_fabm
       REAL(wp), DIMENSION(jpi,jpj,jpk) :: logchl_3d
+      REAL(wp), DIMENSION(jpi,jpj,jpk) :: pco2_3d
 #endif
       CHARACTER(LEN=20) :: datestr=" ",timestr=" "
@@ -1674 +1678 @@
          ENDIF
 #elif defined key_fabm
-         !fco2(:,:)  =  ???                 ! fCO2 from ERSEM
-         CALL ctl_stop( ' Trying to run fco2 observation operator', &
-            &           ' but not properly implemented for FABM-ERSEM yet' )
+         ! First, get pCO2 from FABM
+         pco2_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabm_o3pc)
+         pco2(:,:) = pco2_3d(:,:,1)
+         ! Now, convert pCO2 to fCO2, based on SST in K
+         ! This follows the method used by SOCAT to convert pCO2 to fCO2, see:
+         ! http://www.socat.info/Uniform_format_surface_fCO2_database.pdf
+         fco2(:,:) = pco2(:,:) * EXP((-1636.75                                                                               + &
+            &                         12.0408      * (tsn(:,:,1,jp_tem)+rt0)                                                 - &
+            &                         0.0327957    * (tsn(:,:,1,jp_tem)+rt0)*(tsn(:,:,1,jp_tem)+rt0)                         + &
+            &                         0.0000316528 * (tsn(:,:,1,jp_tem)+rt0)*(tsn(:,:,1,jp_tem)+rt0)*(tsn(:,:,1,jp_tem)+rt0) + &
+            &                         2.0 * (57.7 - 0.118 * (tsn(:,:,1,jp_tem)+rt0)))                                        / &
+            &                        (82.0578 * (tsn(:,:,1,jp_tem)+rt0)))
 #else
          CALL ctl_stop( ' Trying to run fco2 observation operator', &
@@ -1710 +1723 @@
          ENDIF
 #elif defined key_fabm
-         !pco2(:,:)  =  ???                 ! pCO2 from ERSEM
-         CALL ctl_stop( ' Trying to run pCO2 observation operator', &
-            &           ' but not properly implemented for FABM-ERSEM yet' )
+         pco2_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabm_o3pc)
+         pco2(:,:) = pco2_3d(:,:,1)
 #else
          CALL ctl_stop( ' Trying to run pCO2 observation operator', &