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

Timestamp:

2019-12-10T12:38:10+01:00 (3 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/obs_read_surf.F90 (modified) (3 diffs)

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

@@ -39 +39 @@
 
    SUBROUTINE obs_rea_surf( surfdata, knumfiles, cdfilenames, &
-      &                     kvars, kextr, kstp, ddobsini, ddobsend, &
-      &                     ldignmis, ldmod, ldnightav, ldclim, cdvars )
+      &                     kvars, kextr, kstp, ddobsini, ddobsend, MeanPeriodHours, &
+      &                     ldignmis, ldmod, ldnightav, ldclim, ln_time_mean_sla_bkg, cdvars )
       !!---------------------------------------------------------------------
       !!
@@ -72 +72 @@
       LOGICAL, INTENT(IN) :: ldnightav  ! Observations represent a night-time average
       LOGICAL, INTENT(IN) :: ldclim     ! Will include climatology at obs points.
+      LOGICAL, INTENT(IN) :: ln_time_mean_sla_bkg     ! Will reset times to end of averaging period.
       REAL(dp), INTENT(IN) :: ddobsini   ! Obs. ini time in YYYYMMDD.HHMMSS
       REAL(dp), INTENT(IN) :: ddobsend   ! Obs. end time in YYYYMMDD.HHMMSS
+      REAL(wp), INTENT(IN) :: MeanPeriodHours ! Averaging period in hours
       CHARACTER(len=8), DIMENSION(kvars), INTENT(IN) :: cdvars
 
@@ -280 +282 @@
                inpfiles(jj)%iobsj = -1
             ENDIF
+
+            !If SLA observations are representing a time mean then set the time
+            !of the obs to the end of that meaning period relative to the start of the run
+            IF ( ln_time_mean_sla_bkg .AND. ( TRIM( clvarsin(1) ) == 'SLA' ) ) THEN
+               DO ji = 1, inpfiles(jj)%nobs
+                  ! Only do this for obs within time window
+                  IF ( ( inpfiles(jj)%ptim(ji) > djulini(jj) ) .AND. &
+                     & ( inpfiles(jj)%ptim(ji) <= djulend(jj) ) ) THEN
+                     inpfiles(jj)%ptim(ji) = &
+                           & djulini(jj) + (MeanPeriodHours/24.)
+                  ENDIF      
+               END DO
+            ENDIF   
+
             inowin = 0
             DO ji = 1, inpfiles(jj)%nobs