Changeset 12035 for branches/UKMO/dev_r5518_obs_oper_update_SlaAssim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_surf.F90

Timestamp:

2019-12-03T12:20:53+01:00 (4 years ago)

Author:

kingr

Message:

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

File:

• branches/UKMO/dev_r5518_obs_oper_update_SlaAssim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_surf.F90 (modified) (3 diffs)

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

@@ -39 +39 @@
 
    SUBROUTINE obs_rea_surf( surfdata, knumfiles, cdfilenames, &
-      &                     kvars, kextr, kstp, ddobsini, ddobsend, &
+      &                     kvars, kextr, kstp, ddobsini, ddobsend, MeanPeriodHours, &
       &                     ldignmis, ldmod, ldnightav, ldclim, cdvars )
       !!---------------------------------------------------------------------
@@ -74 +74 @@
       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 +281 @@
                inpfiles(jj)%iobsj = -1
             ENDIF
+
+            !If the 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 ( MeanPeriodHours > 0._wp ) THEN
+               DO ji = 1, inpfiles(jj)%nobs
+                  inpfiles(jj)%ptim(ji) = &
+                        & djulini(jj) + (MeanPeriodHours/24.)
+               END DO
+            ENDIF   
+
             inowin = 0
             DO ji = 1, inpfiles(jj)%nobs