Changeset 8919

Timestamp:

2017-12-06T13:52:42+01:00 (6 years ago)

Author:

mattmartin

Message:

Updated the NEMO book documentation with information on the observation averaging operator.

Location:

branches/2017/dev_r8657_UKMO_OBSoper/DOC

Files:

• Figures/Fig_OBS_avg_rad.eps (added)
• Figures/Fig_OBS_avg_rad.pdf (added)
• Figures/Fig_OBS_avg_rec.eps (added)
• Figures/Fig_OBS_avg_rec.pdf (added)
• Namelists/namobs (modified) (3 diffs)
• TexFiles/Chapters/Chap_OBS.tex (modified) (4 diffs)

branches/2017/dev_r8657_UKMO_OBSoper/DOC/Namelists/namobs

@@ -9 +9 @@
    ln_sic      = .false.             ! Logical switch for Sea Ice observations
    ln_vel3d    = .false.             ! Logical switch for velocity observations
+   ln_sss     = .false.              ! Logical swithc for SSS observations
    ln_altbias  = .false.             ! Logical switch for altimeter bias correction
    ln_nea      = .false.             ! Logical switch for rejection of observations near land
@@ -16 +17 @@
    ln_s_at_t   = .false.             ! Logical switch for computing model S at T obs if not there
    ln_sstnight = .false.             ! Logical switch for calculating night-time average for SST obs
+   ln_sla_fp_indegs = .true.         ! Logical for SLA: T=> averaging footpring is in degrees, F=> in metres
+   ln_sst_fp_indegs = .true.         ! Logical for SST: T=> averaging footpring is in degrees, F=> in metres
+   ln_sss_fp_indegs = .true.         ! Logical for SSS: T=> averaging footpring is in degrees, F=> in metres
+   ln_sic_fp_indegs = .true.         ! Logical for SIC: T=> averaging footpring is in degrees, F=> in metres
 ! All of the *files* variables below are arrays. Use namelist_cfg to add more files
    cn_profbfiles = 'profiles_01.nc'  ! Profile feedback input observation file names
    cn_slafbfiles = 'sla_01.nc'       ! SLA feedback input observation file names
    cn_sstfbfiles = 'sst_01.nc'       ! SST feedback input observation file names
+   cn_sssfbfiles = 'sss_01.nc'       ! SSS feedback input observation file names
    cn_sicfbfiles = 'sic_01.nc'       ! SIC feedback input observation file names
    cn_velfbfiles = 'vel_01.nc'       ! Velocity feedback input observation file names
@@ -27 +33 @@
    rn_dobsini  = 00010101.000000     ! Initial date in window YYYYMMDD.HHMMSS
    rn_dobsend  = 00010102.000000     ! Final date in window YYYYMMDD.HHMMSS
+   rn_sla_avglamscl = 0.             ! E/W diameter of SLA observation footprint (metres/degrees)
+   rn_sla_avgphiscl = 0.             ! N/S diameter of SLA observation footprint (metres/degrees)
+   rn_sst_avglamscl = 0.             ! E/W diameter of SST observation footprint (metres/degrees)
+   rn_sst_avgphiscl = 0.             ! N/S diameter of SST observation footprint (metres/degrees)
+   rn_sss_avglamscl = 0.             ! E/W diameter of SSS observation footprint (metres/degrees)
+   rn_sss_avgphiscl = 0.             ! N/S diameter of SSS observation footprint (metres/degrees)
+   rn_sic_avglamscl = 0.             ! E/W diameter of SIC observation footprint (metres/degrees)
+   rn_sic_avgphiscl = 0.             ! N/S diameter of SIC observation footprint (metres/degrees)
    nn_1dint    = 0                   ! Type of vertical interpolation method
    nn_2dint    = 0                   ! Type of horizontal interpolation method
    nn_msshc    = 0                   ! MSSH correction scheme
+   nn_2dint_sla = 0                    ! Horizontal interpolation method for SLA
+   nn_2dint_sst = 0                    ! Horizontal interpolation method for SST
+   nn_2dint_sss = 0                    ! Horizontal interpolation method for SSS
+   nn_2dint_sic = 0                    ! Horizontal interpolation method for SIC
    rn_mdtcorr  = 1.61                ! MDT  correction
    rn_mdtcutoff = 65.0               ! MDT cutoff for computed correction

branches/2017/dev_r8657_UKMO_OBSoper/DOC/TexFiles/Chapters/Chap_OBS.tex

@@ -27 +27 @@
 is set to true.
 
-For all data types a 2D horizontal  interpolator is needed to interpolate the model fields to
+For all data types a 2D horizontal interpolator or averager is needed to interpolate/average the model fields to
 the observation location. For {\em in situ} profiles, a 1D vertical interpolator is needed in
-addition to provide model fields at the observation depths. Currently this only works in
-z-level model configurations, but is being developed to work with a generalised vertical
+addition to provide model fields at the observation depths. This now works in a generalised vertical
 coordinate system. 
 
 Some profile observation types (e.g. tropical moored buoys) are made available as daily averaged quantities.
-The observation operator code can be set-up to calculated the equivalent daily average model temperature fields
+The observation operator code can be set-up to calculate the equivalent daily average model temperature fields
 using the \np{nn\_profdavtypes} namelist array. Some SST observations are equivalent to a night-time
 average value and the observation operator code can calculate equivalent night-time average model SST fields by
@@ -40 +39 @@
 observation time is used.
 
-The code is controlled by the namelist \textit{nam\_obs}. See the following sections for more
+The code is controlled by the namelist \textit{namobs}. See the following sections for more
 details on setting up the namelist.
 
@@ -560 +559 @@
 \label{OBS_theory}
 
-\subsection{Horizontal interpolation methods}
-
+\subsection{Horizontal interpolation and averaging methods}
+
+For most observation types, the horizontal extent of the observation is small compared to the model grid size
+and so the model equivalent of the observation is calculated by interpolating from the four surrounding grid
+points to the observation location. Some satellite observations (e.g. microwave satellite SST data, or SSS data)
+have a footprint which is similar size or larger than the model grid size (particularly when the grid size is small). 
+In those cases the model counterpart should be calculated by averaging the model grid points over the same size as the footprint.
+NEMO therefore has the capability to specify either an interpolation or an averaging (for surface observation types only). 
+
+The main namelist option associated with the interpolation/averaging is nn_2dint. This default option can be set to values from 0 to 6. 
+Values between 0 to 4 are associated with interpolation while values 5 or 6 are associated with averaging.
+\begin{itemize}
+\item n2dint=0: Distance-weighted interpolation
+\item n2dint=1: Distance-weighted interpolation (small angle)
+\item n2dint=2: Bilinear interpolation (geographical grid)
+\item n2dint=3: Bilinear remapping interpolation (general grid)
+\item n2dint=4: Polynomial interpolation
+\item n2dint=5: Radial footprint averaging with radius specified in the namelist as rn_???_avglamscl in degrees or metres (set using ln_???_fp_indegs)
+\item n2dint=6: Rectangular footprint averaging with E/W and N/S size specified in the namelist as rn_???_avglamscl and rn_???_avgphiscl in degrees or metres (set using ln_???_fp_indegs)
+\end{itemize}
+The ??? in the last two options indicate these options should be specified for each observation type for which the averaging is to be performed (see namelist example above).
+The n2dint default option can be overridden for surface observation types using namelist values nn_2dint_??? where ??? is one of sla,sst,sss,sic.
+
+Below is some more detail on the various options for interpolation and averaging available in NEMO.
+
+\subsubsection{Horizontal interpolation}
 Consider an observation point ${\rm P}$ with 
 with longitude and latitude $({\lambda_{}}_{\rm P}, \phi_{\rm P})$ and the 
@@ -660 +683 @@
   
 \end{enumerate}
+
+\subsubsection{Horizontal averaging}
+
+For each surface observation type:
+\begin{itemize}
+\item The standard grid-searching code is used to find the nearest model grid point to the observation location (see next subsection).
+\item The maximum number of grid points is calculated in the local grid domain for which the averaging is likely need to cover.
+\item The lats/longs of the grid points surrounding the nearest model grid box are extracted using existing mpi routines.
+\item The weights for each grid point associated with each observation are calculated, either for radial or rectangular footprints. For grid points completely within the footprint, the weight is one; for grid points completely outside the footprint, the weight is zero. For grid points which are partly within the footprint the ratio between the area of the footprint within the grid box and the total area of the grid box is used as the weight.
+\item The weighted average of the model grid points associated with each observation is calculated, and this is then given as the model counterpart of the observation.
+\end{itemize}
+
+%>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+\begin{figure}      \begin{center}
+\includegraphics[width=10cm,height=12cm]{Fig_OBS_avg_rec}
+\caption{      \label{fig:obsavgrec}
+Weights associated with each model grid box (blue lines and numbers) for an observation at -170.5E, 56.0N with a footprint of $1\deg \time 1\deg$} 
+\end{center}      \end{figure}
+%>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+
+
+%>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+\begin{figure}      \begin{center}
+\includegraphics[width=10cm,height=12cm]{Fig_OBS_avg_rad}
+\caption{      \label{fig:obsavgrad}
+As for figure \ref{obsavgrec} but for a radial footprint with diameter $1\deg$.
+\end{center}      \end{figure}
+%>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+
 
 \subsection{Grid search}