Changeset 15227 for NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_oper.F90

Timestamp:

2021-09-03T16:02:35+02:00 (3 years ago)

Author:

dford

Message:

Implement FBD/SIT (compiles but otherwise not tested).

File:

• NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_oper.F90 (modified) (5 diffs)

NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_oper.F90

@@ -18 +18 @@
    USE lib_mpp,       ONLY :   ctl_warn, ctl_stop           ! Warning and stopping routines
    USE sbcdcy,        ONLY :   sbc_dcy, nday_qsr            ! For calculation of where it is night-time
-   USE obs_grid,      ONLY :   obs_level_search     
+   USE obs_grid,      ONLY :   obs_level_search
+   USE obs_field,     ONLY : cobsname_sla, cobsname_fbd, imaxavtypes
+#if defined key_si3 || defined key_cice
+   USE phycst,        ONLY : rhos, rhoi, rhow               ! For conversion from sea ice freeboard to thickness
+#endif
    !
    USE par_kind     , ONLY :   wp   ! Precision variables
@@ -28 +32 @@
    PUBLIC   obs_prof_opt   !: Compute the model counterpart of profile obs
    PUBLIC   obs_surf_opt   !: Compute the model counterpart of surface obs
-
-   INTEGER, PARAMETER, PUBLIC ::   imaxavtypes = 20   !: Max number of daily avgd obs types
 
    !!----------------------------------------------------------------------
@@ -453 +455 @@
       &                     kit000, kdaystp, kvar, psurf, psurfmask, &
       &                     k2dint, ldnightav, plamscl, pphiscl,     &
-      &                     lindegrees, kssh, kmdt )
+      &                     lindegrees, kssh, kmdt, kfbd, ksnow )
 
       !!-----------------------------------------------------------------------
@@ -514 +516 @@
       INTEGER, OPTIONAL, INTENT(IN) :: &
          & kmdt                        ! Index of extra variable representing MDT
+      INTEGER, OPTIONAL, INTENT(IN) :: &
+         & kfbd                        ! Index of additional variable representing ice freeboard
+      INTEGER, OPTIONAL, INTENT(IN) :: &
+         & ksnow                       ! Index of extra variable representing ice snow thickness
 
       !! * Local declarations
@@ -756 +762 @@
          ENDIF
 
-         IF ( TRIM(surfdataqc%cvars(kvar)) == 'SLA' .AND. PRESENT(kssh) .AND. PRESENT(kmdt) ) THEN
+         IF ( TRIM(surfdataqc%cvars(kvar)) == cobsname_sla .AND. PRESENT(kssh) .AND. PRESENT(kmdt) ) THEN
             ! ... Remove the MDT from the SSH at the observation point to get the SLA
             surfdataqc%radd(jobs,kssh,kvar) = zext(1)
             surfdataqc%rmod(jobs,kvar) = surfdataqc%radd(jobs,kssh,kvar) - surfdataqc%rext(jobs,kmdt)
+#if defined key_si3 || defined key_cice
+         ELSE IF ( TRIM(surfdataqc%cvars(kvar)) == cobsname_fbd .AND. PRESENT(kfbd) .AND. PRESENT(ksnow) ) THEN
+            ! Convert radar freeboard to true freeboard
+            ! (add 1/4 snow depth; 1/4 based on ratio of speed of light in vacuum
+            !  compared to snow (3.0e8 vs 2.4e8 m/s))
+            surfdataqc%radd(jobs,kfbd,kvar) = surfdataqc%robs(jobs,kvar)
+            surfdataqc%robs(jobs,kvar) = surfdataqc%radd(jobs,kfbd,kvar) + 0.25*surfdataqc%rext(jobs,ksnow)
+            ! If the corrected freeboard observation is outside -0.3 to 3.0 m (CPOM) then set the QC flag to bad
+            IF ((surfdataqc%robs(jobs,kvar) < -0.3) .OR. (surfdataqc%robs(jobs,kvar) > 3.0)) THEN
+               surfdataqc%nqc(jobs) = 4
+            ENDIF           
+            ! Convert corrected freeboard to ice thickness following Tilling et al. (2016)
+            surfdataqc%robs(jobs,kvar) = (surfdataqc%robs(jobs,kvar)*rhow + surfdataqc%rext(jobs,ksnow)*rhos)/ &
+               &                         (rhow - rhoi)
+#endif
          ELSE
             surfdataqc%rmod(jobs,kvar) = zext(1)
          ENDIF
-! DO THIS FOR FBD TOO
 
          IF ( zext(1) == obfillflt ) THEN