- Timestamp:
- 2016-11-28T17:04:10+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_INGV_UKMO_2016/DOC/TexFiles/Chapters/Chap_OBS.tex
r4245 r7351 1 \documentclass[NEMO_book]{subfiles} 2 \begin{document} 1 3 % ================================================================ 2 4 % Chapter observation operator (OBS) … … 18 20 location and nearest model timestep. The resulting data are saved in a ``feedback'' file (or 19 21 files). The code was originally developed for use with the NEMOVAR data assimilation code, but 20 can be used for validation or verification of model orany other data assimilation system.21 22 The OBS code is called from \mdl{nemogcm .F90} for model initialisation and to calculate the model22 can be used for validation or verification of the model or with any other data assimilation system. 23 24 The OBS code is called from \mdl{nemogcm} for model initialisation and to calculate the model 23 25 equivalent values for observations on the 0th timestep. The code is then called again after 24 each timestep from \mdl{step .F90}. To build with the OBS code active \key{diaobs} must be25 set.26 each timestep from \mdl{step}. The code is only activated if the namelist logical \np{ln\_diaobs} 27 is set to true. 26 28 27 29 For all data types a 2D horizontal interpolator is needed to interpolate the model fields to … … 29 31 addition to provide model fields at the observation depths. Currently this only works in 30 32 z-level model configurations, but is being developed to work with a generalised vertical 31 coordinate system. Temperature data from moored buoys (TAO, TRITON, PIRATA) in the 32 ENACT/ENSEMBLES data-base are available as daily averaged quantities. For this type of 33 observation the observation operator will compare such observations to the model temperature 34 fields averaged over one day. The relevant observation type may be specified in the namelist 35 using \np{endailyavtypes}. Otherwise the model value from the nearest timestep to the 33 coordinate system. 34 35 Some profile observation types (e.g. tropical moored buoys) are made available as daily averaged quantities. 36 The observation operator code can be set-up to calculated the equivalent daily average model temperature fields 37 using the \np{nn\_profdavtypes} namelist array. Some SST observations are equivalent to a night-time 38 average value and the observation operator code can calculate equivalent night-time average model SST fields by 39 setting the namelist value \np{ln\_sstnight} to true. Otherwise the model value from the nearest timestep to the 36 40 observation time is used. 37 41 38 42 The code is controlled by the namelist \textit{nam\_obs}. See the following sections for more 39 details on setting up the namelist. 43 details on setting up the namelist. 40 44 41 45 Section~\ref{OBS_example} introduces a test example of the observation operator code including … … 59 63 60 64 \begin{enumerate} 61 \item Compile NEMO with \key{diaobs} set.62 63 \item Download some EN SEMBLES EN3data from64 \href{http://www. hadobs.org}{http://www.hadobs.org}. Choose observations which are65 \item Compile NEMO. 66 67 \item Download some EN4 data from 68 \href{http://www.metoffice.gov.uk/hadobs}{http://www.metoffice.gov.uk/hadobs}. Choose observations which are 65 69 valid for the period of your test run because the observation operator compares 66 70 the model and observations for a matching date and time. 67 71 68 \item Add the following to the NEMO namelist to run the observation 69 operator on this data. Set the \np{enactfiles} namelist variable to the 70 observation file name: 72 \item Compile the OBSTOOLS code using: 73 \begin{verbatim} 74 ./maketools -n OBSTOOLS -m [ARCH]. 75 \end{verbatim} 76 77 \item Convert the EN4 data into feedback format: 78 \begin{verbatim} 79 enact2fb.exe profiles_01.nc EN.4.1.1.f.profiles.g10.YYYYMM.nc 80 \end{verbatim} 81 82 \item Include the following in the NEMO namelist to run the observation 83 operator on this data: 71 84 \end{enumerate} 72 85 … … 77 90 Options are defined through the \ngn{namobs} namelist variables. 78 91 The options \np{ln\_t3d} and \np{ln\_s3d} switch on the temperature and salinity 79 profile observation operator code. The \np{ln\_ena} switch turns on the reading 80 of ENACT/ENSEMBLES type profile data. The filename or array of filenames are 81 specified using the \np{enactfiles} variable. The model grid points for a 92 profile observation operator code. The filename or array of filenames are 93 specified using the \np{cn\_profbfiles} variable. The model grid points for a 82 94 particular observation latitude and longitude are found using the grid 83 95 searching part of the code. This can be expensive, particularly for large … … 92 104 A number of utilities are now provided to plot the feedback files, convert and 93 105 recombine the files. These are explained in more detail in section~\ref{OBS_obsutils}. 106 Utilites to convert other input data formats into the feedback format are also 107 described in section~\ref{OBS_obsutils}. 94 108 95 109 \section{Technical details} … … 104 118 %------------------------------------------------------------------------------------------------------------- 105 119 106 Th is name list uses the "feedback" type observation file inputformat for107 profile, sea level anomaly and sea surface temperature data. All the120 The observation operator code uses the "feedback" observation file format for 121 all data types. All the 108 122 observation files must be in NetCDF format. Some example headers (produced using 109 123 \mbox{\textit{ncdump~-h}}) for profile … … 732 746 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 733 747 \begin{figure} \begin{center} 734 \includegraphics[width=10cm,height=12cm,angle=-90.]{ ./TexFiles/Figures/Fig_ASM_obsdist_local}748 \includegraphics[width=10cm,height=12cm,angle=-90.]{Fig_ASM_obsdist_local} 735 749 \caption{ \label{fig:obslocal} 736 750 Example of the distribution of observations with the geographical distribution of observational data.} … … 759 773 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 760 774 \begin{figure} \begin{center} 761 \includegraphics[width=10cm,height=12cm,angle=-90.]{ ./TexFiles/Figures/Fig_ASM_obsdist_global}775 \includegraphics[width=10cm,height=12cm,angle=-90.]{Fig_ASM_obsdist_global} 762 776 \caption{ \label{fig:obsglobal} 763 777 Example of the distribution of observations with the round-robin distribution of observational data.} … … 1376 1390 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 1377 1391 \begin{figure} \begin{center} 1378 %\includegraphics[width=10cm,height=12cm,angle=-90.]{ ./TexFiles/Figures/Fig_OBS_dataplot_main}1379 \includegraphics[width=9cm,angle=-90.]{ ./TexFiles/Figures/Fig_OBS_dataplot_main}1392 %\includegraphics[width=10cm,height=12cm,angle=-90.]{Fig_OBS_dataplot_main} 1393 \includegraphics[width=9cm,angle=-90.]{Fig_OBS_dataplot_main} 1380 1394 \caption{ \label{fig:obsdataplotmain} 1381 1395 Main window of dataplot.} … … 1388 1402 %>>>>>>>>>>>>>>>>>>>>>>>>>>>> 1389 1403 \begin{figure} \begin{center} 1390 %\includegraphics[width=10cm,height=12cm,angle=-90.]{ ./TexFiles/Figures/Fig_OBS_dataplot_prof}1391 \includegraphics[width=7cm,angle=-90.]{ ./TexFiles/Figures/Fig_OBS_dataplot_prof}1404 %\includegraphics[width=10cm,height=12cm,angle=-90.]{Fig_OBS_dataplot_prof} 1405 \includegraphics[width=7cm,angle=-90.]{Fig_OBS_dataplot_prof} 1392 1406 \caption{ \label{fig:obsdataplotprofile} 1393 1407 Profile plot from dataplot produced by right clicking on a point in the main window.} … … 1398 1412 1399 1413 1414 \end{document}
Note: See TracChangeset
for help on using the changeset viewer.