Changeset 12380 for branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

Timestamp:

2020-02-13T15:10:57+01:00 (4 years ago)

Author:

dcarneir

Message:

Updating obs_oper with obs_oper version in the apps trunk

File:

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

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

@@ -57 +57 @@
    LOGICAL :: ln_sit_fp_indegs     !: T=> SIT 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
@@ -70 +71 @@
    REAL(wp) :: rn_sit_avglamscl     !: E/W diameter of SIT observation footprint
    REAL(wp) :: rn_sit_avgphiscl     !: N/S diameter of SIT observation footprint
+   REAL(wp), PUBLIC :: &
+      &        MeanPeriodHours = 24. + (5./6.) !: Meaning period for surface data.
+
 
    INTEGER :: nn_1dint         !: Vertical interpolation method
@@ -274 +278 @@
          &            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,             &
@@ -448 +452 @@
          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
       !-----------------------------------------------------------------------
@@ -952 +957 @@
                &               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, ln_seaicetypes )
@@ -1077 +1082 @@
       !! * 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
@@ -1686 +1692 @@
             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) )
-
-
-           ! Change label of data from FBD ("freeboard") to SIT ("Sea Ice
-           ! Thickness")
+            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
+
+            ! Change label of data from FBD ("freeboard") to SIT ("Sea Ice
+            ! Thickness")
             IF ( TRIM(surfdataqc(jtype)%cvars(1)) == 'FBD' ) THEN
                  surfdata(jtype)%cvars(1) = 'SIT'