1 | MODULE 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/OPA 3.7 , NEMO Consortium (2015) |
---|
23 | !! $Id: trazdf_imp.F90 6140 2015-12-21 11:35:23Z timgraham $ |
---|
24 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
---|
25 | !!---------------------------------------------------------------------- |
---|
26 | CONTAINS |
---|
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 | |
---|
42 | ! Standard offline obs_oper initialisation |
---|
43 | n_files = 0 ! number of files to cycle through |
---|
44 | sao_files(:) = '' ! list of files to read in |
---|
45 | nn_sao_idx(:) = 0 ! list of indices inside each file |
---|
46 | nn_sao_freq = -1 ! input frequency in time steps |
---|
47 | |
---|
48 | ! Standard offline obs_oper settings |
---|
49 | REWIND( numnam_ref ) ! Namelist namctl in reference namelist : Control prints & Benchmark |
---|
50 | READ ( numnam_ref, namsao, IOSTAT = ios, ERR = 901 ) |
---|
51 | 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsao in reference namelist', .TRUE. ) |
---|
52 | ! |
---|
53 | REWIND( numnam_cfg ) ! Namelist namctl in confguration namelist : Control prints & Benchmark |
---|
54 | READ ( numnam_cfg, namsao, IOSTAT = ios, ERR = 902 ) |
---|
55 | 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsao in configuration namelist', .TRUE. ) |
---|
56 | |
---|
57 | lmask(:) = .FALSE. ! count input files |
---|
58 | WHERE (sao_files(:) /= '') lmask(:) = .TRUE. |
---|
59 | n_files = COUNT(lmask) |
---|
60 | ! |
---|
61 | IF(nn_sao_freq == -1) THEN ! Initialise sub obs window frequency |
---|
62 | nn_sao_freq = nitend - nit000 + 1 ! Run length |
---|
63 | ENDIF |
---|
64 | ! |
---|
65 | IF(lwp) THEN ! Print summary of settings |
---|
66 | WRITE(numout,*) |
---|
67 | WRITE(numout,*) 'offline obs_oper : Initialization' |
---|
68 | WRITE(numout,*) '~~~~~~~~~~~~~~~~~' |
---|
69 | WRITE(numout,*) ' Namelist namsao : set stand alone obs_oper parameters' |
---|
70 | DO jf = 1, n_files |
---|
71 | WRITE(numout,'(1X,2A)') ' Input forecast file name forecastfile = ', TRIM(sao_files(jf)) |
---|
72 | WRITE(numout,*) ' Input forecast file index forecastindex = ', nn_sao_idx(jf) |
---|
73 | END DO |
---|
74 | END IF |
---|
75 | ! |
---|
76 | END SUBROUTINE sao_data_init |
---|
77 | |
---|
78 | !!====================================================================== |
---|
79 | END MODULE sao_data |
---|
80 | |
---|