Changeset 11468 for branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/obs_write.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_write.F90 (modified) (22 diffs)

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

@@ -95 +95 @@
       INTEGER :: je
       INTEGER :: iadd
+      INTEGER :: iadd_clm ! 1 if climatology present
       INTEGER :: iext
       REAL(wp) :: zpres
 
+
+      iadd_clm = 0 
+      IF ( profdata%lclim ) iadd_clm = 1
+      
       IF ( PRESENT( padd ) ) THEN
          iadd = padd%inum
@@ -123 +128 @@
          clfiletype='profb'
          CALL alloc_obfbdata( fbdata, 2, profdata%nprof, ilevel, &
-            &                 1 + iadd, 1 + iext, .TRUE. )
+            &                 1 + iadd_clm + iadd, 1 + iext, .TRUE. )
          fbdata%cname(1)      = profdata%cvars(1)
          fbdata%cname(2)      = profdata%cvars(2)
@@ -137 +142 @@
          fbdata%caddunit(1,1) = 'Degrees centigrade'
          fbdata%caddunit(1,2) = 'PSU'
+         IF ( profdata%lclim ) THEN
+            fbdata%caddlong(2,1) = 'Climatology interpolated potential temperature'
+            fbdata%caddlong(2,2) = 'Climatology interpolated practical salinity'
+            fbdata%caddunit(2,1) = 'Degrees centigrade'
+            fbdata%caddunit(2,2) = 'PSU'
+         ENDIF
          fbdata%cgrid(:)      = 'T'
          DO je = 1, iext
@@ -144 +155 @@
          END DO
          DO ja = 1, iadd
-            fbdata%caddname(1+ja) = padd%cdname(ja)
+            fbdata%caddname(1+iadd_clm+ja) = padd%cdname(ja)
             DO jvar = 1, 2
-               fbdata%caddlong(1+ja,jvar) = padd%cdlong(ja,jvar)
-               fbdata%caddunit(1+ja,jvar) = padd%cdunit(ja,jvar)
+               fbdata%caddlong(1+iadd_clm+ja,jvar) = padd%cdlong(ja,jvar)
+               fbdata%caddunit(1+iadd_clm+ja,jvar) = padd%cdunit(ja,jvar)
             END DO
          END DO
@@ -154 +165 @@
 
          clfiletype='velfb'
-         CALL alloc_obfbdata( fbdata, 2, profdata%nprof, ilevel, 1, 0, .TRUE. )
+         CALL alloc_obfbdata( fbdata, 2, profdata%nprof, ilevel, &
+            &                 1 + iadd_clm + iadd, 0, .TRUE. )
          fbdata%cname(1)      = profdata%cvars(1)
          fbdata%cname(2)      = profdata%cvars(2)
@@ -170 +182 @@
          fbdata%caddunit(1,1) = 'm/s'
          fbdata%caddunit(1,2) = 'm/s'
+         IF ( profdata%lclim ) THEN
+            fbdata%caddlong(2,1) = 'Climatology interpolated zonal velocity'
+            fbdata%caddlong(2,2) = 'Climatology interpolated meridional velocity'
+            fbdata%caddunit(2,1) = 'm/s'
+            fbdata%caddunit(2,2) = 'm/s'
+         ENDIF         
          fbdata%cgrid(1)      = 'U' 
          fbdata%cgrid(2)      = 'V'
          DO ja = 1, iadd
-            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)
+            fbdata%caddname(1+iadd_clm+ja) = padd%cdname(ja)
+            fbdata%caddlong(1+iadd_clm+ja,1) = padd%cdlong(ja,1)
+            fbdata%caddunit(1+iadd_clm+ja,1) = padd%cdunit(ja,1)
          END DO
 
@@ -246 +264 @@
          & ( TRIM(profdata%cvars(1)) /= 'UVEL' ) ) THEN
          CALL alloc_obfbdata( fbdata, 1, profdata%nprof, ilevel, &
-            &                 1 + iadd, iext, .TRUE. )
+            &                 1 + iadd_clm + iadd, iext, .TRUE. )
          fbdata%cname(1)      = profdata%cvars(1)
          fbdata%coblong(1)    = cllongname
@@ -252 +270 @@
          fbdata%caddlong(1,1) = 'Model interpolated ' // TRIM(cllongname)
          fbdata%caddunit(1,1) = clunits
+         IF ( profdata%lclim ) THEN
+            fbdata%caddlong(2,1) = 'Climatological interpolated ' // TRIM(cllongname)
+            fbdata%caddunit(2,1) = clunits
+         ENDIF         
          fbdata%cgrid(:)      = clgrid
          DO je = 1, iext
@@ -259 +281 @@
          END DO
          DO ja = 1, iadd
-            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)
+            fbdata%caddname(1+iadd_clm+ja) = padd%cdname(ja)
+            fbdata%caddlong(1+iadd_clm+ja,1) = padd%cdlong(ja,1)
+            fbdata%caddunit(1+iadd_clm+ja,1) = padd%cdunit(ja,1)
          END DO
       ENDIF
 
       fbdata%caddname(1)   = 'Hx'
-
+      IF ( profdata%lclim ) fbdata%caddname(1+iadd_clm)   = 'CLM'
+      
       WRITE(clfname, FMT="(A,'_fdbk_',I4.4,'.nc')") TRIM(clfiletype), nproc
 
@@ -319 +342 @@
             DO jk = profdata%npvsta(jo,jvar), profdata%npvend(jo,jvar)
                ik = profdata%var(jvar)%nvlidx(jk)
-               fbdata%padd(ik,jo,1,jvar) = profdata%var(jvar)%vmod(jk)
                fbdata%pob(ik,jo,jvar)    = profdata%var(jvar)%vobs(jk)
                fbdata%pdep(ik,jo)        = profdata%var(jvar)%vdep(jk)
@@ -333 +355 @@
                ENDIF
                fbdata%iobsk(ik,jo,jvar)  = profdata%var(jvar)%mvk(jk)
+               
+               fbdata%padd(ik,jo,1,jvar) = profdata%var(jvar)%vmod(jk)
+               IF ( profdata%lclim ) THEN           
+                  fbdata%padd(ik,jo,1+iadd_clm,jvar) = profdata%var(jvar)%vclm(jk)     
+               ENDIF              
                DO ja = 1, iadd
-                  fbdata%padd(ik,jo,1+ja,jvar) = &
+                  fbdata%padd(ik,jo,1+iadd_clm+ja,jvar) = &
                      & profdata%var(jvar)%vext(jk,padd%ipoint(ja))
                END DO
@@ -420 +447 @@
       INTEGER :: indx_std
       INTEGER :: iadd_std
-
+      INTEGER :: iadd_clm     
+      INTEGER :: iadd_mdt 
+
+      IF ( PRESENT( pext ) ) THEN
+         iext = pext%inum
+      ELSE
+         iext = 0
+      ENDIF
+
+
+      ! Set up number of additional variables to be ouput:
+      ! Hx, CLM, STD, MDT...
+ 
       IF ( PRESENT( padd ) ) THEN
          iadd = padd%inum
@@ -426 +465 @@
          iadd = 0
       ENDIF
-
-      IF ( PRESENT( pext ) ) THEN
-         iext = pext%inum
-      ELSE
-         iext = 0
-      ENDIF
-
+      
       iadd_std = 0
       indx_std = -1
@@ -444 +477 @@
       ENDIF
       
+      iadd_clm = 0
+      IF ( surfdata%lclim ) iadd_clm = 1
+      
+      iadd_mdt = 0
+      IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) iadd_mdt = 1
+      
       CALL init_obfbdata( fbdata )
 
       SELECT CASE ( TRIM(surfdata%cvars(1)) )
       CASE('SLA')
-
+         
          ! SLA needs special treatment because of MDT, so is all done here
          ! Other variables are done more generically
+         ! No climatology for SLA, MDT is our best estimate of that and is already output.
 
          CALL alloc_obfbdata( fbdata, 1, surfdata%nsurf, 1, &
-            &                 2 + iadd + iadd_std, 1 + iext, .TRUE. )
+            &                 1 + iadd_mdt + iadd_std + iadd, &
+            &                 1 + iext, .TRUE. )
 
          clfiletype = 'slafb'
@@ -474 +515 @@
          fbdata%cgrid(1)      = 'T'
          DO ja = 1, iadd
-            fbdata%caddname(2+iadd_std+ja) = padd%cdname(ja)
-            fbdata%caddlong(2+iadd_std+ja,1) = padd%cdlong(ja,1)
-            fbdata%caddunit(2+iadd_std+ja,1) = padd%cdunit(ja,1)
+            fbdata%caddname(1+iadd_mdt+iadd_std+ja) = padd%cdname(ja)
+            fbdata%caddlong(1+iadd_mdt+iadd_std+ja,1) = padd%cdlong(ja,1)
+            fbdata%caddunit(1+iadd_mdt+iadd_std+ja,1) = padd%cdunit(ja,1)
          END DO
 
@@ -485 +526 @@
          clunits    = 'Degree centigrade'
          clgrid     = 'T'
-
+         
       CASE('ICECONC')
 
@@ -499 +540 @@
          clunits    = 'psu'
          clgrid     = 'T'
-
+         
       CASE('SLCHLTOT','LOGCHL','LogChl','logchl')
 
@@ -610 +651 @@
       
          CALL alloc_obfbdata( fbdata, 1, surfdata%nsurf, 1, &
-            &                 1 + iadd + iadd_std, iext, .TRUE. )
+            &                 1 + iadd_std + iadd_clm + iadd, iext, .TRUE. )
 
          fbdata%cname(1)      = surfdata%cvars(1)
@@ -628 +669 @@
          fbdata%cgrid(1)      = clgrid
          DO ja = 1, iadd
-            fbdata%caddname(1+iadd_std+ja) = padd%cdname(ja)
-            fbdata%caddlong(1+iadd_std+ja,1) = padd%cdlong(ja,1)
-            fbdata%caddunit(1+iadd_std+ja,1) = padd%cdunit(ja,1)
+            fbdata%caddname(1+iadd_mdt+iadd_std+iadd_clm+ja) = padd%cdname(ja)
+            fbdata%caddlong(1+iadd_mdt+iadd_std+iadd_clm+ja,1) = padd%cdlong(ja,1)
+            fbdata%caddunit(1+iadd_mdt+iadd_std+iadd_clm+ja,1) = padd%cdunit(ja,1)
          END DO
 
@@ -637 +678 @@
       fbdata%caddname(1)   = 'Hx'
       IF ( indx_std /= -1 ) THEN
-         IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) iadd_std = iadd_std + 1
-         fbdata%caddname(1+iadd_std)   = surfdata%cext(indx_std)
-         fbdata%caddlong(1+iadd_std,1) = 'Obs error standard deviation'
-         fbdata%caddunit(1+iadd_std,1) = fbdata%cobunit(1)
-      ENDIF
+         fbdata%caddname(1+iadd_mdt+iadd_std)   = surfdata%cext(indx_std)
+         fbdata%caddlong(1+iadd_mdt+iadd_std,1) = 'Obs error standard deviation'
+         fbdata%caddunit(1+iadd_mdt+iadd_std,1) = fbdata%cobunit(1)
+      ENDIF
+      
+      IF ( surfdata%lclim ) THEN
+         fbdata%caddname(1+iadd_mdt+iadd_std+iadd_clm)   = 'CLM'
+         fbdata%caddlong(1+iadd_mdt+iadd_std+iadd_clm,1) = 'Climatology'
+         fbdata%caddunit(1+iadd_mdt+iadd_std+iadd_clm,1) = fbdata%cobunit(1)
+      ENDIF
+      
       
       WRITE(clfname, FMT="(A,'_fdbk_',I4.4,'.nc')") TRIM(clfiletype), nproc
@@ -689 +736 @@
             &           fbdata%ptim(jo),   &
             &           krefdate = 19500101 )
-         fbdata%padd(1,jo,1,1) = surfdata%rmod(jo,1)
-         IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) fbdata%padd(1,jo,2,1) = surfdata%rext(jo,1)
+                    
          fbdata%pob(1,jo,1)    = surfdata%robs(jo,1) 
          fbdata%pdep(1,jo)     = 0.0
@@ -706 +752 @@
          ENDIF
          fbdata%iobsk(1,jo,1)  = 0
-         IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) fbdata%pext(1,jo,1) = surfdata%rext(jo,2)
+ 
+         ! Additional variables.
+         ! Hx is always the first additional variable
+         fbdata%padd(1,jo,1,1) = surfdata%rmod(jo,1)
+         ! MDT is output as an additional variable if SLA obs type
+         IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) THEN
+            fbdata%padd(1,jo,1+iadd_mdt,1) = surfdata%rext(jo,1)
+         ENDIF    
+         ! STD is output as an additional variable if available
          IF ( indx_std /= -1 ) THEN
-            fbdata%padd(1,jo,1+iadd_std,1) = surfdata%rext(jo,indx_std)
+            fbdata%padd(1,jo,1+iadd_mdt+iadd_std,1) = surfdata%rext(jo,indx_std)
          ENDIF
+         ! CLM is output as an additional variable if available
+         IF ( surfdata%lclim ) THEN
+            fbdata%padd(1,jo,1+iadd_mdt+iadd_std+iadd_clm,1) = surfdata%rclm(jo,1)
+         ENDIF
+         ! Then other additional variables are output
+         DO ja = 1, iadd
+            fbdata%padd(1,jo,1+iadd_mdt+iadd_std+iadd_clm+ja,1) = &
+               & surfdata%rext(jo,padd%ipoint(ja))
+         END DO
          
-         DO ja = 1, iadd
-            fbdata%padd(1,jo,2+iadd_std+ja,1) = &
-               & surfdata%rext(jo,padd%ipoint(ja))
-         END DO
+         ! Extra variables
+         IF ( TRIM(surfdata%cvars(1)) == 'SLA' ) fbdata%pext(1,jo,1) = surfdata%rext(jo,2)         
          DO je = 1, iext
             fbdata%pext(1,jo,1+je) = &