Changeset 10712 for branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

Timestamp:

2019-02-21T16:44:54+01:00 (5 years ago)

Author:

emmafiedler

Message:

Correct flagging of ice data at zeroth timestep, QC for sea ice thickness

File:

• branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

@@ -805 +805 @@
             surfdataqc%rext(jobs,1) = surfdataqc%robs(jobs,1)            
             surfdataqc%robs(jobs,1) = surfdataqc%rext(jobs,1) + 0.25*surfdataqc%rext(jobs,2)
+            ! 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,1) < -0.3) .OR. (surfdataqc%robs(jobs,1) > 3.0)) THEN
+              surfdataqc%nqc(jobs) = 4
+            ENDIF            
             ! Convert corrected freeboard to ice thickness following Tilling et al. (2016)
             surfdataqc%robs(jobs,1) = (surfdataqc%robs(jobs,1)*rhow + surfdataqc%rext(jobs,2)*rhos)/ &
-                                      (rhow - rhoi)            
+                                      (rhow - rhoi)
+            ! Flag any negative ice thickness values as bad
+            IF (surfdataqc%robs(jobs,1) < 0.0) THEN
+              surfdataqc%nqc(jobs) = 4
+            ENDIF                                      
          ENDIF