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

Timestamp:

2019-08-23T10:37:22+02:00 (5 years ago)

Author:

mattmartin

Message:

Merged changes to allow writing of climatological information to feedback files.

File:

• branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_surf.F90 (modified) (4 diffs)

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

@@ -40 +40 @@
    SUBROUTINE obs_rea_surf( surfdata, knumfiles, cdfilenames, &
       &                     kvars, kextr, kstp, ddobsini, ddobsend, &
-      &                     ldignmis, ldmod, ldnightav )
+      &                     ldignmis, ldmod, ldnightav, ldclim )
       !!---------------------------------------------------------------------
       !!
@@ -71 +71 @@
       LOGICAL, INTENT(IN) :: ldmod      ! Initialize model from input data
       LOGICAL, INTENT(IN) :: ldnightav  ! Observations represent a night-time average
+      LOGICAL, INTENT(IN) :: ldclim     ! Will include climatology at obs points.
       REAL(dp), INTENT(IN) :: ddobsini   ! Obs. ini time in YYYYMMDD.HHMMSS
       REAL(dp), INTENT(IN) :: ddobsend   ! Obs. end time in YYYYMMDD.HHMMSS
@@ -359 +360 @@
          &               iindx   )
 
-      CALL obs_surf_alloc( surfdata, iobs, kvars, iextr, kstp, jpi, jpj )
+      CALL obs_surf_alloc( surfdata, iobs, kvars, iextr, kstp, jpi, jpj, ldclim )
 
       ! Read obs/positions, QC, all variable and assign to surfdata
@@ -462 +463 @@
                   IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) surfdata%rext(iobs,:) = fbrmdi
                ENDIF
-
+               
+               ! Initialise climatology if set
+               IF ( surfdata%lclim ) surfdata%rclm(iobs,1) = fbrmdi
+               
                ! STD (obs error standard deviation) read from file and passed through obs operator
                IF ( iadd_std(jj) /= -1 ) THEN