New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
sao_data.F90 in NEMO/branches/UKMO/NEMO_4.0_GO8_package_text_diagnostics/src/SAO – NEMO

source: NEMO/branches/UKMO/NEMO_4.0_GO8_package_text_diagnostics/src/SAO/sao_data.F90 @ 10947

Last change on this file since 10947 was 10888, checked in by davestorkey, 5 years ago

branches/UKMO/NEMO_4.0_mirror : clear SVN keywords

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', .TRUE. )
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', .TRUE. )
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.