Changeset 7479 for branches/UKMO

Timestamp:

2016-12-08T14:25:38+01:00 (7 years ago)

Author:

dford

Message:

Changes to get chlorophyll from FABM.

File:

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

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

@@ -1446 +1446 @@
          & MEDUSA_FILL_FLT
 #elif defined key_fabm
-      !USE ???                           ! ERSEM chlorophyll, fCO2 and pCO2
+      USE fabm
+      USE trcsms_fabm, ONLY: model      ! ERSEM chlorophyll, fCO2 and pCO2
 #endif
 #if defined key_spm
@@ -1487 +1488 @@
          maskpco2                       ! array for special pco2 mask
       INTEGER :: jn                     ! loop index
+#if defined key_fabm
+      INTEGER :: chl_index
+      REAL(wp), DIMENSION(jpi,jpj,jpk) :: logchl_3d
+#endif
       CHARACTER(LEN=20) :: datestr=" ",timestr=" "
  
@@ -1599 +1604 @@
 
 #if defined key_hadocc
-         logchl(:,:)  = HADOCC_CHL(:,:,1)    ! (not log) chlorophyll from HadOCC
+         logchl(:,:) = HADOCC_CHL(:,:,1)    ! (not log) chlorophyll from HadOCC
 #elif defined key_medusa && defined key_foam_medusa
-         logchl(:,:)  = MEDUSA_CHL(:,:,1)    ! (not log) chlorophyll from HadOCC
+         logchl(:,:) = MEDUSA_CHL(:,:,1)    ! (not log) chlorophyll from HadOCC
 #elif defined key_fabm
-         !logchl(:,:)  =  ???                 ! (not log) chlorophyll from ERSEM
-         CALL ctl_stop( ' Trying to run logchl observation operator', &
-            &           ' but not properly implemented for FABM-ERSEM yet' )
+         DO jn = 1, SIZE(model%diagnostic_variables)
+            IF (TRIM(model%diagnostic_variables(jn)%standard_variable%name) == 'total_chlorophyll') THEN
+               chl_index = jn
+               EXIT
+            ENDIF
+         END DO
+         logchl_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, chl_index)
+         logchl(:,:) = logchl_3d(:,:,1)
 #else
          CALL ctl_stop( ' Trying to run logchl observation operator', &
@@ -1758 +1768 @@
       INTEGER :: jfbini
       CHARACTER(LEN=20) :: datestr=" ",timestr=" "
-      CHARACTER(LEN=10) :: cdtmp
+      CHARACTER(LEN=20) :: cdtmp
       !-----------------------------------------------------------------------
       ! Depending on switches call various observation output routines