Changeset 4048 for branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/off_data.F90

Timestamp:

2013-10-02T18:32:18+02:00 (11 years ago)

Author:

djlea

Message:

Cleaning and debugging of the observation operator. Turn off the night time averaging of SST data by default, but add a namelist option to switch it on.

File:

• branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/off_data.F90 (modified) (4 diffs)

branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/off_data.F90

@@ -3 +3 @@
    IMPLICIT NONE
    INTEGER, PARAMETER :: MaxNumFiles = 1000
+
    !! Class 4 file settings
    INTEGER :: &
-           &  cl4_fclen, &            !: number of forecast days
-           &  jimatch, &              !: current match
-           &  n_match                 !: number of matches
+           & cl4_fcst_idx(MaxNumFiles), & !: forecast indices
+           & cl4_match_len, &             !: number of match types
+           & cl4_fcst_len                 !: number of forecast days
    CHARACTER(len=lc) :: &
-           & cl4_vars(MaxNumFiles), & !: class 4 variables
-           & cl4_sys, &               !: class 4 system
-           & cl4_cfg, &               !: class 4 configuration
-           & cl4_date, &              !: class 4 date
-           & cl4_vn,  &               !: class 4 version
-           & cl4_prefix, &            !: class 4 prefix
-           & cl4_contact, &           !: class 4 contact
-           & cl4_inst                 !: class 4 institute
-   REAL ::   cl4_modjuld              !: model Julian day
+           & cl4_vars(MaxNumFiles), &  !: class 4 variables
+           & cl4_sys, &                !: class 4 system
+           & cl4_cfg, &                !: class 4 configuration
+           & cl4_date, &               !: class 4 date
+           & cl4_vn,  &                !: class 4 version
+           & cl4_prefix, &             !: class 4 prefix
+           & cl4_contact, &            !: class 4 contact
+           & cl4_inst                  !: class 4 institute
+   REAL ::   cl4_modjuld               !: model Julian day
    REAL :: &
-      & cl4_leadtime(MaxNumFiles)     !: Lead time data
+      & cl4_leadtime(MaxNumFiles)      !: Lead time data
+
    !! Offline obs_oper settings
    CHARACTER(len=lc) :: &
-      & model_files(MaxNumFiles)      !: model files
+      & off_files(MaxNumFiles)         !: model files
    INTEGER            :: &
-      & nn_modindex(MaxNumFiles), &   !: model file indices
-      & nn_forecast(MaxNumFiles)      !: forecast indices
+      & jifile, &                      !: current file list index
+      & n_files, &                     !: number of files
+      & jimatch, &                     !: current match
+      & nn_off_idx(MaxNumFiles), &     !: time_counter indices
+      & nn_off_freq                    !: read frequency in time steps
    CHARACTER(len=128) :: &
-      & alt_file                      !: altimeter file
+      & alt_file                       !: altimeter file
 CONTAINS
    SUBROUTINE off_data_init( ld_cl4 )
@@ -44 +49 @@
 
       ! Standard offline obs_oper information
-      NAMELIST/namoff/model_files, nn_modindex, nn_forecast
+      NAMELIST/namoff/off_files, nn_off_idx, nn_off_freq
 
       ! Class 4 file specifiers
       NAMELIST/namcl4/cl4_vars, cl4_sys, cl4_cfg, cl4_date, cl4_vn, &
          &            cl4_prefix, cl4_contact, cl4_inst, cl4_leadtime, &
-         &            cl4_fclen
+         &            cl4_fcst_idx, cl4_fcst_len, cl4_match_len
 
       ! Standard offline obs_oper initialisation
-      jimatch = 0                   !: match iteration variable 
-      n_match = 0                   !: number of matches to perform
-      model_files(:) = ''           !: list of files to read in
-      nn_modindex(:) = 0            !: list of indices inside each file
+      jimatch = 0                   !: match-up iteration variable 
+      jifile = 1                    !: input file iteration variable 
+      n_files = 0                   !: number of files to cycle through
+      off_files(:) = ''             !: list of files to read in
+      nn_off_idx(:) = 0             !: list of indices inside each file
+      nn_off_freq = -1              !: input frequency in time steps
 
       ! Class 4 initialisation
       cl4_leadtime(:) = 0           !: Lead time axis value for each file
-      cl4_fclen = 0                 !: Length of the forecast dimension
+      cl4_fcst_len = 0              !: Length of the forecast dimension
+      cl4_match_len = 1             !: Number of match types
+      cl4_fcst_idx(:) = 0           !: output file forecast index
       cl4_vars(:) = ''              !: output file variable names
       cl4_sys = ''                  !: output file system
@@ -77 +86 @@
       ENDIF
 
-      ! count forecast/persistence files
+      ! count input files
       lmask(:) = .FALSE.
-      WHERE (model_files(:) /= '') lmask(:) = .TRUE.
-      n_match = COUNT(lmask)
+      WHERE (off_files(:) /= '') lmask(:) = .TRUE.
+      n_files = COUNT(lmask)
 
+      !! Initialise sub obs window frequency
+      IF (nn_off_freq == -1) THEN
+         !! Run length
+         nn_off_freq = nitend - nit000 + 1
+      ENDIF
+
+      !! Print summary of settings
       IF(lwp) THEN
          WRITE(numout,*)
@@ -87 +103 @@
          WRITE(numout,*) '~~~~~~~~~~~~~~~~~'
          WRITE(numout,*) '   Namelist namoff : set offline obs_oper parameters' 
-         DO jf = 1, n_match
+         DO jf = 1, n_files
             WRITE(numout,'(1X,2A)') '   Input forecast file name          forecastfile = ', &
-               TRIM(model_files(jf))
+               TRIM(off_files(jf))
             WRITE(numout,*) '   Input forecast file index        forecastindex = ', &
-               nn_modindex(jf)
-            WRITE(numout,*) '   Input forecast leadtime index    leadtimeindex = ', &
-               nn_forecast(jf)
-            WRITE(numout,*) '   Input forecast leadtime value    leadtimevalue = ', &
+               nn_off_idx(jf)
+            WRITE(numout,*) '   Output forecast leadtime index   leadtimeindex = ', &
+               cl4_fcst_idx(jf)
+            WRITE(numout,*) '   Output forecast leadtime value   leadtimevalue = ', &
                cl4_leadtime(jf)
             WRITE(numout,'(1X,2A)') '   Input class 4 variable       class 4 parameter = ', &