Changeset 5992 for branches

Timestamp:

2015-12-03T12:16:20+01:00 (8 years ago)

Author:

timgraham

Message:

Merged branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION into branch

Location:

branches/2015/dev_MetOffice_merge_2015/NEMOGCM

Files:

• CONFIG/SHARED/namelist_ref (modified) (1 diff)
• NEMO/OPA_SRC/OBS/diaobs.F90 (modified) (11 diffs)
• NEMO/OPA_SRC/OBS/obs_sstbias.F90 (copied) (copied from branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/NEMO/OPA_SRC/OBS/obs_sstbias.F90)

branches/2015/dev_MetOffice_merge_2015/NEMOGCM/CONFIG/SHARED/namelist_ref

@@ -1214 +1214 @@
    ln_grid_global = .true.
    ln_grid_search_lookup = .false.
+   ln_sstbias = .false. 
+   sstbias_files = 'sstbias.nc'
 /
 !-----------------------------------------------------------------------

branches/2015/dev_MetOffice_merge_2015/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?
@@ -137 +141 @@
       !!        !  06-10  (A. Weaver) Cleaning and add controls
       !!        !  07-03  (K. Mogensen) General handling of profiles
+      !!        !  14-08  (J.While) Incorporated SST bias correction  
       !!----------------------------------------------------------------------
 
@@ -146 +151 @@
       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)      
@@ -182 +188 @@
          &            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
@@ -195 +202 @@
       INTEGER :: jnumsst
       INTEGER :: jnumsstfb
+      INTEGER :: jnumsstbias 
       INTEGER :: jnumseaice
       INTEGER :: jnumvelavcur
@@ -210 +218 @@
       !-----------------------------------------------------------------------
       
-      !Initalise all values in namelist arrays
-      enactfiles(:) = ''
-      coriofiles(:) = ''
-      profbfiles(:) = ''
-      slafilesact(:) = ''
-      slafilespas(:) = ''
-      slafbfiles(:) = ''
-      sstfiles(:)   = ''
-      sstfbfiles(:) = ''
-      seaicefiles(:) = ''
-      velcurfiles(:) = ''
+      !Initalise all values in namelist arrays 
+      enactfiles(:) = '' 
+      coriofiles(:) = '' 
+      profbfiles(:) = '' 
+      slafilesact(:) = '' 
+      slafilespas(:) = '' 
+      slafbfiles(:) = '' 
+      sstfiles(:)   = '' 
+      sstfbfiles(:) = '' 
+      seaicefiles(:) = '' 
+      velcurfiles(:) = '' 
+      velavcurfiles(:) = '' 
+      velhrcurfiles(:) = '' 
+      velavadcpfiles(:) = '' 
+      velhradcpfiles(:) = '' 
+      velfbfiles(:) = '' 
       veladcpfiles(:) = ''
-      velavcurfiles(:) = ''
-      velhrcurfiles(:) = ''
-      velavadcpfiles(:) = ''
-      velhradcpfiles(:) = ''
-      velfbfiles(:) = ''
-      velcurfiles(:) = ''
-      veladcpfiles(:) = ''
-      endailyavtypes(:) = -1
-      endailyavtypes(1) = 820
-      ln_profb_ena(:) = .FALSE.
-      ln_profb_enatim(:) = .TRUE.
-      ln_velfb_av(:) = .FALSE.
-      ln_ignmis = .FALSE.
+      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 )
@@ -290 +297 @@
          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.
@@ -341 +354 @@
          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
@@ -776 +790 @@
 
          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