Changeset 8954

Timestamp:

2017-12-08T12:32:59+01:00 (7 years ago)

Author:

dford

Message:

Import changes from reanalysis branch.

Location:

branches/UKMO/dev_r4650_general_vert_coord_obsoper_logchl_std/NEMOGCM/NEMO/OPA_SRC/OBS

Files:

• diaobs.F90 (modified) (1 diff)
• obs_read_logchl.F90 (modified) (5 diffs)
• obs_write.F90 (modified) (3 diffs)

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

@@ -1166 +1166 @@
          nlogchlvars = 1
 
-         ! Set the number of extra variables for logchl to 0
-         nlogchlextr = 0
+         ! Set the number of extra variables for logchl to 1 (STD)
+         nlogchlextr = 1
          
          IF ( ln_logchlfb ) THEN

branches/UKMO/dev_r4650_general_vert_coord_obsoper_logchl_std/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_logchl.F90

@@ -82 +82 @@
       INTEGER :: jj
       INTEGER :: jk
+      INTEGER :: jadd
       INTEGER :: iflag
       INTEGER :: inobf
@@ -115 +116 @@
       TYPE(obfbdata), POINTER, DIMENSION(:) :: &
          & inpfiles
+      LOGICAL :: ll_get_std   ! Logical for getting STD
       real(wp), DIMENSION(knumfiles) :: &
          & djulini, &
@@ -120 +122 @@
       INTEGER :: iobs
       INTEGER :: iobstot
+      INTEGER :: iaddref_std
       INTEGER :: ios
       INTEGER :: ioserrcount
@@ -205 +208 @@
                CALL ctl_stop( 'File format unknown' )
             ENDIF
+
+            !------------------------------------------------------------------ 
+            ! Find the references to additional entries, if any 
+            !------------------------------------------------------------------ 
+            ! STD 
+            iaddref_std = 0 
+            ll_get_std = .FALSE.
+            DO jadd = 1,inpfiles(jj) % nadd  
+               IF (  TRIM(inpfiles(jj) % caddname(jadd)) == "STD" ) THEN 
+                  iaddref_std = jadd
+                  ll_get_std = .TRUE. 
+                  EXIT 
+               ENDIF 
+            END DO
 
             !------------------------------------------------------------------
@@ -417 +434 @@
                   logchldata%rmod(iobs,1) = fbrmdi
                ENDIF
+                  
+               ! Copy in STD 
+               IF ( ll_get_std ) THEN 
+                  logchldata%rext(iobs,1) = inpfiles(jj)%padd(1,ji,iaddref_std,1) 
+               ELSE
+                  IF ( jk == 1 )   &
+                  &   CALL ctl_warn(   &
+                  &   "No STD values for logchl observations, setting to zero.") 
+                  logchldata%rext(iobs,1) = 0. 
+               ENDIF 
             ENDIF
          ENDIF

branches/UKMO/dev_r4650_general_vert_coord_obsoper_logchl_std/NEMOGCM/NEMO/OPA_SRC/OBS/obs_write.F90

@@ -974 +974 @@
          nadd = padd%inum
       ELSE
-         nadd = 0
+         nadd = 1
       ENDIF
 
@@ -1000 +1000 @@
       fbdata%caddunit(1,1) = 'mg/m3'
       fbdata%cgrid(1)      = 'T'
-      DO ja = 1, nadd
-         fbdata%caddname(1+ja) = padd%cdname(ja)
-         fbdata%caddlong(1+ja,1) = padd%cdlong(ja,1)
-         fbdata%caddunit(1+ja,1) = padd%cdunit(ja,1)
-      END DO
+      IF ( PRESENT( padd ) ) THEN
+         DO ja = 1, nadd
+            fbdata%caddname(1+ja) = padd%cdname(ja)
+            fbdata%caddlong(1+ja,1) = padd%cdlong(ja,1)
+            fbdata%caddunit(1+ja,1) = padd%cdunit(ja,1)
+         END DO
+      ELSE
+         fbdata%caddname(2) = 'STD'
+         fbdata%caddlong(2) = 'Standard deviation of LOGCHL'
+         fbdata%caddunit(2) = 'mg/m3'
+      ENDIF
 
       WRITE(cfname, FMT="(A,'_fdbk_',I4.4,'.nc')") TRIM(cprefix), nproc
@@ -1067 +1073 @@
          fbdata%iobsk(1,jo,1)  = 0
          DO ja = 1, nadd
-            fbdata%padd(1,jo,1+ja,1) = &
-               & logchldata%rext(jo,padd%ipoint(ja))
+            IF ( PRESENT( padd ) ) THEN
+               fbdata%padd(1,jo,1+ja,1) = &
+                  & logchldata%rext(jo,padd%ipoint(ja))
+            ELSE
+               fbdata%padd(1,jo,1+ja,1) = &
+                  & logchldata%rext(jo,1)
+            ENDIF
          END DO
          DO je = 1, next