Changeset 4768 for branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

Timestamp:

2014-09-17T11:05:03+02:00 (10 years ago)

Author:

jwhile

Message:

Adding code for bias correction to SST obs

File:

• branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90 (modified) (11 diffs)

branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

@@ -25 +25 @@
    USE obs_read_sla             ! Reading and allocation of SLA observations  
    USE obs_read_sst             ! Reading and allocation of SST observations  
+   USE obs_sstbias              ! Bias correction routine for SST 
    USE obs_readmdt              ! Reading and allocation of MDT for SLA.
    USE obs_read_seaice          ! Reading and allocation of Sea Ice observations  
@@ -87 +88 @@
    LOGICAL, PUBLIC :: ln_ignmis      !: Logical switch for ignoring missing files
    LOGICAL, PUBLIC :: ln_s_at_t      !: Logical switch to compute model S at T observations
-
+   LOGICAL, PUBLIC :: ln_sstbias     !: Logical switch for bias corection of SST 
+   
    REAL(KIND=dp), PUBLIC :: dobsini   !: Observation window start date YYYYMMDD.HHMMSS
    REAL(KIND=dp), PUBLIC :: dobsend   !: Observation window end date YYYYMMDD.HHMMSS
@@ -102 +104 @@
    !                    !: If so use endailyavtypes
       & ln_profb_enatim !: Change tim for 820 enact data set.
-   
+
+   INTEGER, DIMENSION(MaxNumFiles), PUBLIC :: sstbias_type !SST bias type    
+      
    LOGICAL, DIMENSION(MaxNumFiles) :: &
       & ln_velfb_av   !: Is the velocity feedback files daily average?
@@ -135 +139 @@
       !!        !  06-10  (A. Weaver) Cleaning and add controls
       !!        !  07-03  (K. Mogensen) General handling of profiles
+      !!        !  14-08  (J.While) Incorporated SST bias correction  
       !!----------------------------------------------------------------------
 
@@ -144 +149 @@
       CHARACTER(len=128) :: profbfiles(MaxNumFiles)
       CHARACTER(len=128) :: sstfiles(MaxNumFiles)      
-      CHARACTER(len=128) :: sstfbfiles(MaxNumFiles) 
+      CHARACTER(len=128) :: sstfbfiles(MaxNumFiles)
+      CHARACTER(len=128) :: sstbias_files(MaxNumFiles) 
       CHARACTER(len=128) :: slafilesact(MaxNumFiles)      
       CHARACTER(len=128) :: slafilespas(MaxNumFiles)      
@@ -180 +186 @@
          &            ln_velhradcp, velhradcpfiles,                   &
          &            ln_velfb, velfbfiles, ln_velfb_av,              &
-         &            ln_profb_enatim, ln_ignmis, ln_cl4
+         &            ln_profb_enatim, ln_ignmis, ln_cl4,             &
+         &            ln_sstbias, sstbias_files          
 
       INTEGER :: jprofset
@@ -193 +200 @@
       INTEGER :: jnumsst
       INTEGER :: jnumsstfb
+      INTEGER :: jnumsstbias 
       INTEGER :: jnumseaice
       INTEGER :: jnumvelavcur
@@ -208 +216 @@
       !-----------------------------------------------------------------------
 
-      velcurfiles(:) = ''
+      !Initalise all values in namelist arrays 
+      enactfiles(:) = '' 
+      coriofiles(:) = '' 
+      profbfiles(:) = '' 
+      slafilesact(:) = '' 
+      slafilespas(:) = '' 
+      slafbfiles(:) = '' 
+      sstfiles(:)   = '' 
+      sstfbfiles(:) = '' 
+      seaicefiles(:) = '' 
+      velcurfiles(:) = '' 
+      veladcpfiles(:) = '' 
+      velavcurfiles(:) = '' 
+      velhrcurfiles(:) = '' 
+      velavadcpfiles(:) = '' 
+      velhradcpfiles(:) = '' 
+      velfbfiles(:) = '' 
+      velcurfiles(:) = '' 
       veladcpfiles(:) = ''
+      sstbias_files(:) = ''  
+      endailyavtypes(:) = -1 
+      endailyavtypes(1) = 820 
+      ln_profb_ena(:) = .FALSE. 
+      ln_profb_enatim(:) = .TRUE. 
+      ln_velfb_av(:) = .FALSE. 
+      ln_ignmis = .FALSE. 
       CALL ini_date( dobsini )
       CALL fin_date( dobsend )
@@ -264 +296 @@
          lmask(:) = .FALSE.
       ENDIF
+      IF (ln_sstbias) THEN 
+         lmask(:) = .FALSE. 
+         WHERE (sstbias_files(:) /= '') lmask(:) = .TRUE. 
+         jnumsstbias = COUNT(lmask) 
+         lmask(:) = .FALSE. 
+      ENDIF      
       IF (ln_seaice) THEN
          lmask(:) = .FALSE.
@@ -315 +353 @@
          WRITE(numout,*) '             Logical switch for GHRSST observations          ln_ghrsst = ', ln_ghrsst
          WRITE(numout,*) '             Logical switch for feedback SST data             ln_sstfb = ', ln_sstfb
+         WRITE(numout,*) '             Logical switch for SST bias correction         ln_sstbias = ', ln_sstbias 
          WRITE(numout,*) '             Logical switch for night-time SST obs         ln_sstnight = ', ln_sstnight
          WRITE(numout,*) '             Logical switch for SSS observations                ln_sss = ', ln_sss
@@ -754 +793 @@
 
          ENDIF
-
+         
+         !Read in bias field and correct SST.
+         IF ( ln_sstbias ) THEN
+            IF ( jnumsstbias == 0 ) CALL ctl_stop("ln_sstbias set,"// &
+                                             "  but no bias"// &
+                                             " files to read in")   
+            CALL obs_app_sstbias( nsstsets, sstdatqc, n2dint, &
+                                  jnumsstbias, & 
+                                  sstbias_files(1:jnumsstbias) )
+         ENDIF
+         
       ENDIF