Ignore:
Timestamp:
2019-02-21T16:44:54+01:00 (3 years ago)
Author:
emmafiedler
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

    r10276 r10712  
    805805            surfdataqc%rext(jobs,1) = surfdataqc%robs(jobs,1)             
    806806            surfdataqc%robs(jobs,1) = surfdataqc%rext(jobs,1) + 0.25*surfdataqc%rext(jobs,2) 
     807            ! If the corrected freeboard observation is outside -0.3 to 3.0 m (CPOM) then set the QC flag to bad 
     808            IF ((surfdataqc%robs(jobs,1) < -0.3) .OR. (surfdataqc%robs(jobs,1) > 3.0)) THEN 
     809              surfdataqc%nqc(jobs) = 4 
     810            ENDIF             
    807811            ! Convert corrected freeboard to ice thickness following Tilling et al. (2016) 
    808812            surfdataqc%robs(jobs,1) = (surfdataqc%robs(jobs,1)*rhow + surfdataqc%rext(jobs,2)*rhos)/ & 
    809                                       (rhow - rhoi)             
     813                                      (rhow - rhoi) 
     814            ! Flag any negative ice thickness values as bad 
     815            IF (surfdataqc%robs(jobs,1) < 0.0) THEN 
     816              surfdataqc%nqc(jobs) = 4 
     817            ENDIF                                       
    810818         ENDIF 
    811819          
Note: See TracChangeset for help on using the changeset viewer.