Changeset 4117
- Timestamp:
- 2013-10-23T22:32:11+02:00 (10 years ago)
- Location:
- branches/2013/dev_r3987_UKMO4_OBS
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3987_UKMO4_OBS/DOC/TexFiles/Chapters/Chap_OBS.tex
r4115 r4117 869 869 ! namoff Offline obs_oper namelist 870 870 !---------------------------------------------------------------------- 871 ! o ff_files specifies the files containing the model counterpart872 ! nn_o ff_idx specifies the time_counter index within the model file871 ! ooo_files specifies the files containing the model counterpart 872 ! nn_ooo_idx specifies the time_counter index within the model file 873 873 &namoff 874 o ff_files = "foo.nc"875 nn_o ff_idx = 2874 ooo_files = "foo.nc" 875 nn_ooo_idx = 2 876 876 / 877 877 \end{verbatim} … … 880 880 \subsubsection{Multiple fields per run} 881 881 882 Model field iteration is controlled via \textbf{nn\_o ff\_freq} which specifies882 Model field iteration is controlled via \textbf{nn\_ooo\_freq} which specifies 883 883 the number of model steps at which the next field gets read. For example, if 884 884 12 hourly fields are to be interpolated in a setup where 288 steps equals 24 hours. … … 890 890 ! namoff Offline obs_oper namelist 891 891 !---------------------------------------------------------------------- 892 ! o ff_files specifies the files containing the model counterpart893 ! nn_o ff_idx specifies the time_counter index within the model file894 ! nn_o ff_freq specifies number of time steps between read operations892 ! ooo_files specifies the files containing the model counterpart 893 ! nn_ooo_idx specifies the time_counter index within the model file 894 ! nn_ooo_freq specifies number of time steps between read operations 895 895 &namoff 896 o ff_files = "foo.nc" "foo.nc"897 nn_o ff_idx = 1 2898 nn_o ff_freq = 144896 ooo_files = "foo.nc" "foo.nc" 897 nn_ooo_idx = 1 2 898 nn_ooo_freq = 144 899 899 / 900 900 \end{verbatim} … … 1046 1046 \tiny 1047 1047 \begin{verbatim} 1048 o ff_files = "F.1.nc" "F.2.nc" "F.3.nc" "P.1.nc" "P.2.nc" "P.3.nc" "A.1.nc"1049 nn_o ff_idx = 1 1 1 1 1 1 11048 ooo_files = "F.1.nc" "F.2.nc" "F.3.nc" "P.1.nc" "P.2.nc" "P.3.nc" "A.1.nc" 1049 nn_ooo_idx = 1 1 1 1 1 1 1 1050 1050 \end{verbatim} 1051 1051 \end{alltt} … … 1060 1060 ! namoff Offline obs_oper namelist 1061 1061 !---------------------------------------------------------------------- 1062 ! o ff_files specifies the files containing the model counterpart1063 ! nn_o ff_idx specifies the time_counter index within the model file1064 ! nn_o ff_freq specifies number of time steps between read operations1062 ! ooo_files specifies the files containing the model counterpart 1063 ! nn_ooo_idx specifies the time_counter index within the model file 1064 ! nn_ooo_freq specifies number of time steps between read operations 1065 1065 &namoff 1066 o ff_files = "F.1.nc" "F.2.nc" "F.3.nc" "P.1.nc" "P.2.nc" "P.3.nc" "A.1.nc"1067 nn_o ff_idx = 1 1 1 1 1 1 11066 ooo_files = "F.1.nc" "F.2.nc" "F.3.nc" "P.1.nc" "P.2.nc" "P.3.nc" "A.1.nc" 1067 nn_ooo_idx = 1 1 1 1 1 1 1 1068 1068 / 1069 1069 !---------------------------------------------------------------------- … … 1122 1122 In certain cases it may be desirable to include both multiple model fields per 1123 1123 observation window with multiple match ups per observation. This can be achieved 1124 by specifying \textbf{nn\_o ff\_freq} as well as the class 4 settings. Care must1125 be taken in generating the o ff\_files list such that the files are arranged into1124 by specifying \textbf{nn\_ooo\_freq} as well as the class 4 settings. Care must 1125 be taken in generating the ooo\_files list such that the files are arranged into 1126 1126 consecutive blocks of single match ups. For example, 2 forecast fields 1127 1127 of 12 hourly data would result in 4 separate read operations but only 2 write … … 1131 1131 \tiny 1132 1132 \begin{verbatim} 1133 o ff_files = "F1.nc" "F1.nc" "F2.nc" "F2.nc"1133 ooo_files = "F1.nc" "F1.nc" "F2.nc" "F2.nc" 1134 1134 ... 1135 1135 cl4_fcst_idx = 1 2 -
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/namelist
r4095 r4117 3 3 ! namoff Offline obs_oper namelist 4 4 !---------------------------------------------------------------------- 5 ! o ff_files specifies the files containing the model counterpart6 ! nn_o ff_idx specifies the index within the model file7 ! nn_o ff_freq specifies the number of timesteps between file reads5 ! ooo_files specifies the files containing the model counterpart 6 ! nn_ooo_idx specifies the index within the model file 7 ! nn_ooo_freq specifies the number of timesteps between file reads 8 8 &namoff 9 o ff_files = 'fcst.0.nc'10 nn_o ff_idx = 19 ooo_files = 'fcst.0.nc' 10 nn_ooo_idx = 1 11 11 / 12 12 !---------------------------------------------------------------------- -
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/obs_fbm.F90
r4110 r4117 1564 1564 !! ** Purpose : Write an obfbdata structure into a class 4 file. 1565 1565 !! 1566 !! ** Method : 1. Allocate memory needed by o ff_wri1566 !! ** Method : 1. Allocate memory needed by ooo_write 1567 1567 !! 2. Map obfbdata into allocated memory 1568 !! 3. Pass mapped data to o ff_wri1568 !! 3. Pass mapped data to ooo_write 1569 1569 !! 4. Deallocate memory 1570 1570 !!---------------------------------------------------------------------- -
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/ooo_data.F90
r4110 r4117 31 31 !! Offline obs_oper settings 32 32 CHARACTER(len=lc) :: & 33 & o ff_files(MaxNumFiles) !: model files33 & ooo_files(MaxNumFiles) !: model files 34 34 INTEGER :: & 35 35 & jifile, & !: current file list index 36 36 & n_files, & !: number of files 37 37 & jimatch, & !: current match 38 & nn_o ff_idx(MaxNumFiles), & !: time_counter indices39 & nn_o ff_freq !: read frequency in time steps38 & nn_ooo_idx(MaxNumFiles), & !: time_counter indices 39 & nn_ooo_freq !: read frequency in time steps 40 40 CHARACTER(len=128) :: & 41 41 & alt_file !: altimeter file … … 55 55 56 56 ! Standard offline obs_oper information 57 NAMELIST/namoff/o ff_files, nn_off_idx, nn_off_freq57 NAMELIST/namoff/ooo_files, nn_ooo_idx, nn_ooo_freq 58 58 59 59 ! Class 4 file specifiers … … 66 66 jifile = 1 !: input file iteration variable 67 67 n_files = 0 !: number of files to cycle through 68 o ff_files(:) = '' !: list of files to read in69 nn_o ff_idx(:) = 0 !: list of indices inside each file70 nn_o ff_freq = -1 !: input frequency in time steps68 ooo_files(:) = '' !: list of files to read in 69 nn_ooo_idx(:) = 0 !: list of indices inside each file 70 nn_ooo_freq = -1 !: input frequency in time steps 71 71 72 72 ! Class 4 initialisation … … 94 94 ! count input files 95 95 lmask(:) = .FALSE. 96 WHERE (o ff_files(:) /= '') lmask(:) = .TRUE.96 WHERE (ooo_files(:) /= '') lmask(:) = .TRUE. 97 97 n_files = COUNT(lmask) 98 98 99 99 !! Initialise sub obs window frequency 100 IF (nn_o ff_freq == -1) THEN100 IF (nn_ooo_freq == -1) THEN 101 101 !! Run length 102 nn_o ff_freq = nitend - nit000 + 1102 nn_ooo_freq = nitend - nit000 + 1 103 103 ENDIF 104 104 … … 111 111 DO jf = 1, n_files 112 112 WRITE(numout,'(1X,2A)') ' Input forecast file name forecastfile = ', & 113 TRIM(o ff_files(jf))113 TRIM(ooo_files(jf)) 114 114 WRITE(numout,*) ' Input forecast file index forecastindex = ', & 115 nn_o ff_idx(jf)115 nn_ooo_idx(jf) 116 116 WRITE(numout,*) ' Output forecast leadtime index leadtimeindex = ', & 117 117 cl4_fcst_idx(jf) -
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/ooo_read.F90
r4110 r4117 42 42 43 43 !! Filename, index and match-up kind 44 cdfilename = TRIM(o ff_files(kfile))44 cdfilename = TRIM(ooo_files(kfile)) 45 45 cmatchname = TRIM(cl4_vars(kfile)) 46 kindex = nn_o ff_idx(kfile)46 kindex = nn_ooo_idx(kfile) 47 47 48 48 !! Update model fields -
branches/2013/dev_r3987_UKMO4_OBS/NEMOGCM/NEMO/OOO_SRC/ooomod.F90
r4110 r4117 75 75 DO WHILE ( istp <= nitend .AND. nstop == 0 ) 76 76 IF (jifile <= n_files + 1) THEN 77 IF ( MOD(istp, nn_o ff_freq) == nit000 ) THEN77 IF ( MOD(istp, nn_ooo_freq) == nit000 ) THEN 78 78 !! Read next model counterpart 79 79 CALL ooo_rea_dri(jifile)
Note: See TracChangeset
for help on using the changeset viewer.