Changeset 8856 for branches/UKMO/dev_r5518_obs_oper_update_obserr/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_surf.F90

Timestamp:

2017-11-30T14:06:43+01:00 (7 years ago)

Author:

mattmartin

Message:

Add first version of code to include obs error standard deviations in output files and read in from input files.

File:

• branches/UKMO/dev_r5518_obs_oper_update_obserr/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_surf.F90 (modified) (10 diffs)

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

@@ -97 +97 @@
       INTEGER :: ios
       INTEGER :: ioserrcount
+      INTEGER :: iextr
       INTEGER, PARAMETER :: jpsurfmaxtype = 1024
       INTEGER, DIMENSION(knumfiles) :: irefdate
@@ -108 +109 @@
          & iindx,    &
          & ifileidx, &
-         & isurfidx
+         & isurfidx, &
+         & iadd_std
       REAL(wp), DIMENSION(:), ALLOCATABLE :: &
          & zphi, &
@@ -122 +124 @@
 
       ! Local initialization
-      iobs = 0
-
+      iobs  = 0
+      iextr = kextr
       !-----------------------------------------------------------------------
       ! Count the number of files needed and allocate the obfbdata type
@@ -131 +133 @@
 
       ALLOCATE( inpfiles(inobf) )
+      ALLOCATE( iadd_std(inobf) )
 
       surf_files : DO jj = 1, inobf
@@ -195 +198 @@
                   ENDIF
                END DO
+            ENDIF
+
+            iadd_std(jj) = -1
+            IF ( inpfiles(jj)%nadd > 0 ) THEN
+               DO ji = 1, inpfiles(jj)%nadd
+                  IF ( TRIM( inpfiles(jj)%caddname(ji) ) == 'STD' ) THEN
+                     iextr = kextr + 1
+                     iadd_std(jj) = ji
+                  ENDIF
+               END DO
+            ENDIF
+
+            IF(lwp) THEN
+               IF ( iadd_std(jj) /= -1 ) THEN
+                  WRITE(numout,*) ' STD variable available in input file so passing it through the obs oper'
+                  WRITE(numout,*)
+               ENDIF
             ENDIF
 
@@ -339 +359 @@
          &               iindx   )
 
-      CALL obs_surf_alloc( surfdata, iobs, kvars, kextr, kstp, jpi, jpj )
+      CALL obs_surf_alloc( surfdata, iobs, kvars, iextr, kstp, jpi, jpj )
 
       ! Read obs/positions, QC, all variable and assign to surfdata
@@ -345 +365 @@
       iobs = 0
       surfdata%cvars(:)  = clvars(:)
+      IF ( ldmod .AND. ( TRIM( surfdata%cvars(1) ) == 'SLA' ) ) THEN
+         surfdata%cext(1) = 'SSH'
+         surfdata%cext(2) = 'MDT'
+      ENDIF
+      IF ( iextra > kextr ) surfdata%cext(iextr) = 'STD'
 
       ityp   (:) = 0
@@ -426 +451 @@
                surfdata%robs(iobs,1) = inpfiles(jj)%pob(1,ji,1)
 
-
                ! Model and MDT is set to fbrmdi unless read from file
                IF ( ldmod ) THEN
@@ -437 +461 @@
                   surfdata%rmod(iobs,1) = fbrmdi
                   IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) surfdata%rext(iobs,:) = fbrmdi
+               ENDIF
+
+               ! STD (obs error standard deviation) read from file and passed through obs operator
+               IF ( iadd_std(jj) /= -1 ) THEN
+                  surfdata%rext(iobs,iextr) = inpfiles(jj)%padd(1,ji,iadd_std(jj),1)
                ENDIF
             ENDIF
@@ -488 +517 @@
       END DO
       DEALLOCATE( inpfiles )
+      DEALLOCATE( iadd_std )
 
    END SUBROUTINE obs_rea_surf