[4849] | 1 | MODULE sao_data |
---|
[7646] | 2 | !!====================================================================== |
---|
| 3 | !! *** MODULE sao_data *** |
---|
| 4 | !!====================================================================== |
---|
| 5 | !! History : 3.6 ! 2015-12 (A. Ryan) Original code |
---|
| 6 | !!---------------------------------------------------------------------- |
---|
[4030] | 7 | USE par_kind, ONLY: lc |
---|
[4844] | 8 | USE lib_mpp ! distributed memory computing |
---|
[7646] | 9 | USE in_out_manager |
---|
[4110] | 10 | |
---|
[4030] | 11 | IMPLICIT NONE |
---|
[4110] | 12 | |
---|
[4030] | 13 | INTEGER, PARAMETER :: MaxNumFiles = 1000 |
---|
[4048] | 14 | |
---|
[4846] | 15 | !! Stand Alone Observation operator settings |
---|
[7646] | 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 | !!---------------------------------------------------------------------- |
---|
[9598] | 22 | !! NEMO/OCE 4.0 , NEMO Consortium (2018) |
---|
[10069] | 23 | !! $Id$ |
---|
[10068] | 24 | !! Software governed by the CeCILL license (see ./LICENSE) |
---|
[7646] | 25 | !!---------------------------------------------------------------------- |
---|
[4030] | 26 | CONTAINS |
---|
[7646] | 27 | |
---|
[5063] | 28 | SUBROUTINE sao_data_init() |
---|
[4030] | 29 | !!---------------------------------------------------------------------- |
---|
[4849] | 30 | !! *** SUBROUTINE sao_data_init *** |
---|
[4030] | 31 | !! |
---|
| 32 | !! ** Purpose : To read namelists and initialise offline_oper run. |
---|
| 33 | !! |
---|
| 34 | !!---------------------------------------------------------------------- |
---|
[7646] | 35 | INTEGER :: jf ! file dummy loop index |
---|
| 36 | LOGICAL :: lmask(MaxNumFiles) ! Logical mask used for counting |
---|
| 37 | INTEGER :: ios |
---|
| 38 | !! |
---|
[4846] | 39 | NAMELIST/namsao/sao_files, nn_sao_idx, nn_sao_freq |
---|
[7646] | 40 | !!---------------------------------------------------------------------- |
---|
[9169] | 41 | IF(lwp) THEN |
---|
| 42 | WRITE(numout,*) |
---|
| 43 | WRITE(numout,*) 'sao_data_init : offline obs operator initialization' |
---|
| 44 | WRITE(numout,*) '~~~~~~~~~~~~~' |
---|
| 45 | ENDIF |
---|
[4030] | 46 | |
---|
| 47 | ! Standard offline obs_oper initialisation |
---|
[9169] | 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 |
---|
[4030] | 52 | |
---|
| 53 | ! Standard offline obs_oper settings |
---|
[4844] | 54 | REWIND( numnam_ref ) ! Namelist namctl in reference namelist : Control prints & Benchmark |
---|
[4846] | 55 | READ ( numnam_ref, namsao, IOSTAT = ios, ERR = 901 ) |
---|
[11536] | 56 | 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsao in reference namelist' ) |
---|
[4844] | 57 | REWIND( numnam_cfg ) ! Namelist namctl in confguration namelist : Control prints & Benchmark |
---|
[4846] | 58 | READ ( numnam_cfg, namsao, IOSTAT = ios, ERR = 902 ) |
---|
[11536] | 59 | 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namsao in configuration namelist' ) |
---|
[7646] | 60 | |
---|
| 61 | lmask(:) = .FALSE. ! count input files |
---|
[9169] | 62 | WHERE( sao_files(:) /= '' ) lmask(:) = .TRUE. |
---|
[4048] | 63 | n_files = COUNT(lmask) |
---|
[7646] | 64 | ! |
---|
| 65 | IF(nn_sao_freq == -1) THEN ! Initialise sub obs window frequency |
---|
| 66 | nn_sao_freq = nitend - nit000 + 1 ! Run length |
---|
[4048] | 67 | ENDIF |
---|
[7646] | 68 | ! |
---|
| 69 | IF(lwp) THEN ! Print summary of settings |
---|
[5063] | 70 | WRITE(numout,*) ' Namelist namsao : set stand alone obs_oper parameters' |
---|
[4048] | 71 | DO jf = 1, n_files |
---|
[9169] | 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) |
---|
[4030] | 74 | END DO |
---|
| 75 | END IF |
---|
[7646] | 76 | ! |
---|
[4849] | 77 | END SUBROUTINE sao_data_init |
---|
[4030] | 78 | |
---|
[7646] | 79 | !!====================================================================== |
---|
[4849] | 80 | END MODULE sao_data |
---|
[4030] | 81 | |
---|