Changeset 5063 for branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/SAO_SRC/sao_data.F90

Timestamp:

2015-02-05T17:29:55+01:00 (9 years ago)

Author:

andrewryan

Message:

gross simplification of stand alone observation operator

File:

• branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/SAO_SRC/sao_data.F90 (modified) (6 diffs)

branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/SAO_SRC/sao_data.F90

@@ -12 +12 @@
    INTEGER, PARAMETER :: MaxNumFiles = 1000
 
-   !! Class 4 file settings
-   INTEGER :: &
-           & 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
-   REAL :: &
-      & cl4_leadtime(MaxNumFiles)      !: Lead time data
-
    !! Stand Alone Observation operator settings
    CHARACTER(len=lc) :: &
       & sao_files(MaxNumFiles)         !: model files
    INTEGER            :: &
-      & jifile, &                      !: current file list index
       & n_files, &                     !: number of files
-      & jimatch, &                     !: current match
       & nn_sao_idx(MaxNumFiles), &     !: time_counter indices
       & nn_sao_freq                    !: read frequency in time steps
-   CHARACTER(len=128) :: &
-      & alt_file                       !: altimeter file
 CONTAINS
-   SUBROUTINE sao_data_init( ld_cl4 )
+   SUBROUTINE sao_data_init()
       !!----------------------------------------------------------------------
       !!                    ***  SUBROUTINE sao_data_init ***
@@ -53 +31 @@
          & jf                           !: file dummy loop index
       LOGICAL :: lmask(MaxNumFiles)     !: Logical mask used for counting
-      LOGICAL, INTENT(IN) :: ld_cl4     !: Logical class 4 on/off
       INTEGER :: ios
 
@@ -59 +36 @@
       NAMELIST/namsao/sao_files, nn_sao_idx, nn_sao_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_fcst_idx, cl4_fcst_len, cl4_match_len
-
       ! Standard offline obs_oper initialisation
-      jimatch = 0                   !: match-up iteration variable 
-      jifile = 1                    !: input file iteration variable 
       n_files = 0                   !: number of files to cycle through
       sao_files(:) = ''             !: list of files to read in
       nn_sao_idx(:) = 0             !: list of indices inside each file
       nn_sao_freq = -1              !: input frequency in time steps
-
-      ! Class 4 initialisation
-      cl4_leadtime(:) = 0           !: Lead time axis value for each file
-      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
-      cl4_cfg = ''                  !: output file configuration
-      cl4_date = ''                 !: output file date string
-      cl4_vn = ''                   !: output file version
-      cl4_prefix = 'class4'         !: output file prefix
-      cl4_contact = ''              !: output file contact details
-      cl4_inst = ''                 !: output file institution
 
       ! Standard offline obs_oper settings
@@ -95 +51 @@
 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsao in configuration namelist', .TRUE. )
 
-      ! Read class 4 output settings
-      IF (ld_cl4) THEN
-         REWIND( numnam_ref )              ! Namelist namctl in reference namelist : Control prints & Benchmark
-         READ  ( numnam_ref, namcl4, IOSTAT = ios, ERR = 903 )
-903      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcl4 in reference namelist', .TRUE. )
-
-         REWIND( numnam_cfg )              ! Namelist namctl in confguration namelist : Control prints & Benchmark
-         READ  ( numnam_cfg, namcl4, IOSTAT = ios, ERR = 904 )
-904      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcl4 in configuration namelist', .TRUE. )
-      ENDIF
 
       ! count input files
@@ -122 +68 @@
          WRITE(numout,*) 'offline obs_oper : Initialization'
          WRITE(numout,*) '~~~~~~~~~~~~~~~~~'
-         WRITE(numout,*) '   Namelist namsao : set stand alone obs_oper parameters' 
+         WRITE(numout,*) '   Namelist namsao : set stand alone obs_oper parameters'
          DO jf = 1, n_files
             WRITE(numout,'(1X,2A)') '   Input forecast file name          forecastfile = ', &
@@ -128 +74 @@
             WRITE(numout,*) '   Input forecast file index        forecastindex = ', &
                nn_sao_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 = ', &
-               TRIM(cl4_vars(jf))
          END DO
-         WRITE(numout, '(1X,2A)') '   Input class 4 system            class 4 system = ', &
-            TRIM(cl4_sys)
-         WRITE(numout, '(1X,2A)') '   Input class 4 config            class 4 config = ', &
-            TRIM(cl4_cfg)
-         WRITE(numout, '(1X,2A)') '   Input class 4 date                class 4 date = ', &
-            TRIM(cl4_date)
-         WRITE(numout, '(1X,2A)') '   Input class 4 version          class 4 version = ', &
-            TRIM(cl4_vn)
-         WRITE(numout, '(1X,2A)') '   Input class 4 prefix            class 4 prefix = ', &
-            TRIM(cl4_prefix)
-         WRITE(numout, '(1X,2A)') '   Input class 4 contact          class 4 contact = ', &
-            TRIM(cl4_contact)
-         WRITE(numout, '(1X,2A)') '   Input class 4 institute      class 4 institute = ', &
-            TRIM(cl4_inst)
       END IF