Changeset 15650

Timestamp:

2022-01-18T09:50:49+01:00 (2 years ago)

Author:

kingr

Message:

A-TSCV: reject SLA where SIC above limit.

Location:

branches/UKMO/dev_r5518_obs_oper_WISA_osse/NEMOGCM/NEMO/OPA_SRC/OBS

Files:

• diaobs.F90 (modified) (2 diffs)
• obs_oper.F90 (modified) (1 diff)
• obs_write.F90 (modified) (2 diffs)

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

@@ -959 +959 @@
                & ltype_clim = .TRUE.
 
-            IF ( (TRIM(cobstypessurf(jtype)) == 'sla') .OR. &
-               & (TRIM(cobstypessurf(jtype)) == 'sit') .OR. &
+            IF ( TRIM(cobstypessurf(jtype)) == 'sla') THEN
+               ltype_clim = .TRUE.
+               nvarssurf(jtype) = 1
+               nextrsurf(jtype) = 2
+            ELSE IF ( (TRIM(cobstypessurf(jtype)) == 'sit') .OR. &
                & (TRIM(cobstypessurf(jtype)) == 'fbd') ) THEN
                nvarssurf(jtype) = 1
@@ -1467 +1470 @@
             CASE('sla')
                zsurfvar(:,:,1) = sshn(:,:)
+               IF ( kstp == 0 ) THEN
+                  IF ( lwp .AND. surfdataqc(jtype)%nsstpmpp(1) > 0 ) THEN
+                     CALL ctl_warn( 'Sea-ice not initialised on zeroth '// &
+                        &           'time-step but some obs are valid then.' )
+                     WRITE(numout,*)surfdataqc(jtype)%nsstpmpp(1), &
+                        &           ' SLA obs will not be checked against ice coverage'
+                  ENDIF
+                  zsurfclim(:,:,1) = 0.0
+               ELSE
+                  ! Pass sea-ice concentration in through SLA climatology
+                  ! for use in rejecting SLA under ice.
+                  zsurfclim(:,:,1) = fr_i(:,:)
+               ENDIF
+
+
             CASE('sss')
                zsurfvar(:,:,1) = tsn(:,:,1,jp_sal)

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

@@ -954 +954 @@
             IF ( surfdataqc%lclim ) surfdataqc%rclm(jobs,1) = zclm(1)
 
+            IF ( surfdataqc%lclim .AND. TRIM(surfdataqc%cvars(1)) == 'SLA' ) THEN
+               IF (surfdataqc%rclm(jobs,1) > 0.01) THEN
+                   surfdataqc%nqc(jobs) = 4
+               END IF
+            ENDIF
+
             IF ( zext(1) == obfillflt ) THEN
                ! If the observation value is a fill value, set QC flag to bad

branches/UKMO/dev_r5518_obs_oper_WISA_osse/NEMOGCM/NEMO/OPA_SRC/OBS/obs_write.F90

@@ -494 +494 @@
 
          CALL alloc_obfbdata( fbdata, 1, surfdata%nsurf, 1, &
-            &                 1 + iadd_mdt + iadd_std + iadd, &
+            &                 1 + iadd_mdt + iadd_std + iadd_clm + iadd, &
             &                 1 + iext, .TRUE. )
 
@@ -516 +516 @@
          fbdata%cgrid(1)      = 'T'
          DO ja = 1, iadd
-            fbdata%caddname(1+iadd_mdt+iadd_std+ja) = padd%cdname(ja)
-            fbdata%caddlong(1+iadd_mdt+iadd_std+ja,1) = padd%cdlong(ja,1)
-            fbdata%caddunit(1+iadd_mdt+iadd_std+ja,1) = padd%cdunit(ja,1)
+            fbdata%caddname(1+iadd_mdt+iadd_std+iadd_clm+ja) = padd%cdname(ja)
+            fbdata%caddlong(1+iadd_mdt+iadd_std+iadd_clm+ja,1) = padd%cdlong(ja,1)
+            fbdata%caddunit(1+iadd_mdt+iadd_std+iadd_clm+ja,1) = padd%cdunit(ja,1)
          END DO