source: NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/SAO/sao_data.F90 @ 11317

Last change on this file since 11317 was 11317, checked in by smasson, 14 months ago

dev_r10984_HPC-13 : improve error handling, see #2307 and #2285

  • Property svn:keywords set to Id
File size: 3.6 KB
Line 
1MODULE sao_data
2   !!======================================================================
3   !!                       ***  MODULE sao_data  ***
4   !!======================================================================
5   !! History :  3.6  ! 2015-12  (A. Ryan)  Original code
6   !!----------------------------------------------------------------------
7   USE par_kind, ONLY: lc
8   USE lib_mpp         ! distributed memory computing
9   USE in_out_manager
10
11   IMPLICIT NONE
12
13   INTEGER, PARAMETER :: MaxNumFiles = 1000
14
15   !! Stand Alone Observation operator settings
16   CHARACTER(len=lc) ::   sao_files(MaxNumFiles)   !: model files
17   INTEGER           ::   n_files                  !: number of files
18   INTEGER           :: nn_sao_idx(MaxNumFiles)    !: time_counter indices
19   INTEGER           :: nn_sao_freq                !: read frequency in time steps
20   
21   !!----------------------------------------------------------------------
22   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
23   !! $Id$
24   !! Software governed by the CeCILL license (see ./LICENSE)
25   !!----------------------------------------------------------------------
26CONTAINS
27
28   SUBROUTINE sao_data_init()
29      !!----------------------------------------------------------------------
30      !!                    ***  SUBROUTINE sao_data_init ***
31      !!
32      !! ** Purpose : To read namelists and initialise offline_oper run.
33      !!
34      !!----------------------------------------------------------------------
35      INTEGER ::   jf                   ! file dummy loop index
36      LOGICAL ::   lmask(MaxNumFiles)   ! Logical mask used for counting
37      INTEGER ::   ios
38      !!
39      NAMELIST/namsao/sao_files, nn_sao_idx, nn_sao_freq
40      !!----------------------------------------------------------------------
41      IF(lwp) THEN
42         WRITE(numout,*)
43         WRITE(numout,*) 'sao_data_init : offline obs operator initialization'
44         WRITE(numout,*) '~~~~~~~~~~~~~'
45      ENDIF
46
47      ! Standard offline obs_oper initialisation
48      n_files       =  0         ! number of files to cycle through
49      sao_files(:)  = ''         ! list of files to read in
50      nn_sao_idx(:) =  0         ! list of indices inside each file
51      nn_sao_freq   = -1         ! input frequency in time steps
52
53      ! Standard offline obs_oper settings
54      REWIND( numnam_ref )              ! Namelist namctl in reference namelist : Control prints & Benchmark
55      READ  ( numnam_ref, namsao, IOSTAT = ios, ERR = 901 )
56901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namsao in reference namelist' )
57      REWIND( numnam_cfg )              ! Namelist namctl in confguration namelist : Control prints & Benchmark
58      READ  ( numnam_cfg, namsao, IOSTAT = ios, ERR = 902 )
59902   IF( ios >  0 )   CALL ctl_nam ( ios , 'namsao in configuration namelist' )
60     
61      lmask(:) = .FALSE.               ! count input files
62      WHERE( sao_files(:) /= '' )   lmask(:) = .TRUE.
63      n_files = COUNT(lmask)
64      !
65      IF(nn_sao_freq == -1) THEN      ! Initialise sub obs window frequency
66         nn_sao_freq = nitend - nit000 + 1      ! Run length
67      ENDIF
68      !
69      IF(lwp) THEN                     ! Print summary of settings
70         WRITE(numout,*) '   Namelist namsao : set stand alone obs_oper parameters'
71         DO jf = 1, n_files
72            WRITE(numout,'(1X,2A)') '      Input forecast file name         forecastfile  = ', TRIM(sao_files(jf))
73            WRITE(numout,*)         '      Input forecast file index        forecastindex = ', nn_sao_idx(jf)
74         END DO
75      END IF
76      !
77   END SUBROUTINE sao_data_init
78
79   !!======================================================================
80END MODULE sao_data
81
Note: See TracBrowser for help on using the repository browser.