Changeset 12140 for branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

Timestamp:

2019-12-10T12:38:10+01:00 (4 years ago)

Author:

kingr

Message:

Added option to remove tides from SLA bkg by taking average over 24h50m.

File:

• branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90 (modified) (7 diffs)

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

@@ -55 +55 @@
    LOGICAL :: ln_sic_fp_indegs     !: T=> sea-ice obs footprint size specified in degrees, F=> in metres
    LOGICAL :: ln_output_clim       !: Logical switch for interpolating and writing T/S climatology
+   LOGICAL :: ln_time_mean_sla_bkg !: Logical switch for applying time mean of SLA background to remove tidal signal
 
    REAL(wp) :: rn_default_avglamscl !: Default E/W diameter of observation footprint
@@ -66 +67 @@
    REAL(wp) :: rn_sic_avglamscl     !: E/W diameter of sea-ice observation footprint
    REAL(wp) :: rn_sic_avgphiscl     !: N/S diameter of sea-ice observation footprint
+   REAL(wp), PUBLIC :: &
+      &        MeanPeriodHours = 24. + (5./6.) !: Meaning period for surface data.
+
 
    INTEGER :: nn_1dint         !: Vertical interpolation method
@@ -262 +266 @@
          &            ln_ignmis, ln_s_at_t, ln_bound_reject,          &
          &            ln_sstnight,  ln_output_clim,                   &
-         &            ln_default_fp_indegs,                           &
+         &            ln_time_mean_sla_bkg, ln_default_fp_indegs,     &
          &            ln_sla_fp_indegs, ln_sst_fp_indegs,             &
          &            ln_sss_fp_indegs, ln_sic_fp_indegs,             &
@@ -427 +431 @@
          WRITE(numout,*) '             Logical switch for night-time SST obs               ln_sstnight = ', ln_sstnight
          WRITE(numout,*) '             Logical switch for writing climat. at obs points ln_output_clim = ', ln_output_clim
+         WRITE(numout,*) '             Logical switch for time-mean of SLA        ln_time_mean_sla_bkg = ', ln_time_mean_sla_bkg
       ENDIF
       !-----------------------------------------------------------------------
@@ -912 +917 @@
                &               clsurffiles(jtype,1:ifilessurf(jtype)), &
                &               nvarssurf(jtype), nextrsurf(jtype), nitend-nit000+2, &
-               &               rn_dobsini, rn_dobsend, ln_ignmis, .FALSE., &
-               &               llnightav(jtype), ltype_clim, clvars )
+               &               rn_dobsini, rn_dobsend, MeanPeriodHours, ln_ignmis, .FALSE., &
+               &               llnightav(jtype), ltype_clim, ln_time_mean_sla_bkg, clvars )
 
             CALL obs_pre_surf( surfdata(jtype), surfdataqc(jtype), ln_nea, ln_bound_reject )
@@ -1031 +1036 @@
       !! * Local declarations
       INTEGER :: idaystp           ! Number of timesteps per day
+      INTEGER :: imeanstp          ! Number of timesteps for sla averaging
       INTEGER :: jtype             ! Data loop variable
       INTEGER :: jvar              ! Variable number
@@ -1624 +1630 @@
             ENDIF
 
-            CALL obs_surf_opt( surfdataqc(jtype), kstp, jpi, jpj,       &
-               &               nit000, idaystp, zsurfvar,               &
-               &               zsurfclim, zsurfmask,                    &
-               &               n2dintsurf(jtype), llnightav(jtype),     &
-               &               ravglamscl(jtype), ravgphiscl(jtype),    &
-               &               lfpindegs(jtype) )
+            IF ( TRIM(cobstypessurf(jtype)) == 'sla' .AND.                 &
+                  &  ln_time_mean_sla_bkg ) THEN
+               !Number of time-steps in meaning period
+               imeanstp = NINT( ( MeanPeriodHours * 60. * 60. ) / rdt )
+               CALL obs_surf_opt( surfdataqc(jtype), kstp, jpi, jpj,       &
+                  &               nit000, idaystp, zsurfvar,               &
+                  &               zsurfclim, zsurfmask,                    &
+                  &               n2dintsurf(jtype), llnightav(jtype),     &
+                  &               ravglamscl(jtype), ravgphiscl(jtype),    &
+                  &               lfpindegs(jtype), kmeanstp = imeanstp )
+
+            ELSE
+               CALL obs_surf_opt( surfdataqc(jtype), kstp, jpi, jpj,       &
+                  &               nit000, idaystp, zsurfvar,               &
+                  &               zsurfclim, zsurfmask,                    &
+                  &               n2dintsurf(jtype), llnightav(jtype),     &
+                  &               ravglamscl(jtype), ravgphiscl(jtype),    &
+                  &               lfpindegs(jtype) )
+            ENDIF
 
          END DO