New operator: sealevelpressure - sea level pressure

.gitattributes

@@ -98,6 +98,7 @@ doc/tex/mod/Consecstat -text
 doc/tex/mod/Copy -text
 doc/tex/mod/Daypctl -text
 doc/tex/mod/Daystat -text
+doc/tex/mod/Derivepar -text
 doc/tex/mod/Detrend -text
 doc/tex/mod/Diff -text
 doc/tex/mod/Duplicate -text
@@ -195,7 +196,6 @@ doc/tex/mod/Rhopot -text
 doc/tex/mod/Rotuv -text
 doc/tex/mod/Runpctl -text
 doc/tex/mod/Runstat -text
-doc/tex/mod/Sealevelpressure -text
 doc/tex/mod/Seaspctl -text
 doc/tex/mod/Seasstat -text
 doc/tex/mod/Selbox -text

ChangeLog

@@ -3,6 +3,10 @@
 	* using CDI library version 1.6.4
 	* Version 1.6.4 released
 
+2014-04-14  Uwe Schulzweida
+
+	* New operator: sealevelpressure - sea level pressure
+
 2014-03-14  Uwe Schulzweida
 
 	* remapcon: check for missing values in normalisation

NEWS

@@ -5,6 +5,10 @@ Version 1.6.4 (20 May 2014):
 
    New features:
      * Option -H: Do not append to netCDF "history" global attribute
+   New operators:
+     * sealevelpressure: sea level pressure
+   Changed operators:
+     * diff: print number of different values
 
 Version 1.6.3 (18 February 2014):
 

OPERATORS

@@ -488,6 +488,7 @@ Operator catalog:
    Vargen        stdatm          Create values for pressure and temperature for hydrostatic atmosphere
    Rotuv         rotuvb          Backward rotation
    Mastrfu       mastrfu         Mass stream function
+   Derivepar     sealevelpressure Sea level pressure
    Adisit        adisit          Potential temperature to in-situ temperature
    Adisit        adipot          In-situ temperature to potential temperature
    Rhopot        rhopot          Calculates potential density

doc/tex/Modules

@@ -113,6 +113,7 @@ Timsort       Miscellaneous
 Vargen        Miscellaneous
 Rotuv         Miscellaneous
 Mastrfu       Miscellaneous
+Derivepar     Miscellaneous
 Adisit        Miscellaneous
 Rhopot        Miscellaneous
 Histogram     Miscellaneous

doc/tex/appendix.tex

@@ -88,14 +88,30 @@ CDO\_PCTL\_NBINS            &     101    & Number of histogram bins. \\ \hline
 CDO\_RESET\_HISTORY       &         0  & Set to 1 to reset the netCDF {\it history} global attribute.\\ \hline
 CDO\_REMAP\_SEARCH\_RADIUS &  180 & Remap search radius in degree. Used by the operators \\
                                                      &        & remapdis and remapnn. \\ \hline
-
-
 \end{tabular}
 
 
 \input{cdo_openmp}
 
 
+\chapter{\label{stdnametable}Standard name table}
+
+The following CF standard names are supported by {\CDO}.
+
+\begin{tabular}[t]{|>{\columncolor{pcolor1}}l|r|r|l|}
+\hline
+\rowcolor{pcolor2}
+%\cellcolor{pcolor2}
+{\bf CF standard name}       &  Units & {\bf GRIB 1 code} & {\bf variable name} \\ \hline
+surface\_geopotential           &  m2 s-2 &  129   &  geosp \\ \hline
+air\_temperature                    & K           & 130  &  ta \\ \hline
+specific\_humidity                 & 1           & 133  & hus  \\ \hline
+surface\_air\_pressure           & Pa          & 134  & aps  \\ \hline
+air\_pressure\_at\_sea\_level & Pa         & 151  & psl  \\ \hline
+geopotential\_height              & m         & 156  &  zg \\ \hline
+\end{tabular}
+
+
 \chapter{\label{appendixgrid}Grid description examples}
 
 \section{Example of a curvilinear grid description}

doc/tex/mod/Sealevelpressure → doc/tex/mod/Derivepar

 @BeginModule
 @NewPage
-@Name      = Sealevelpressure
+@Name      = Derivepar
 @Title     = Sea level pressure
 @Section   = Miscellaneous
 @Arguments = ifile ofile

src/Derivepar.c

@@ -380,7 +380,7 @@ void *Derivepar(void *argument)
 	cdoAbort("Spectral data unsupported!");
     }
 
-  if ( tempID == -1 ) cdoAbort("Air temperature not found!");
+  if ( tempID == -1 ) cdoAbort("%s not found!", var_stdname(air_temperature));
 
   array  = (double*) malloc(ngp*sizeof(double));
   geop   = (double*) malloc(ngp*sizeof(double));

src/modules.c

@@ -562,7 +562,7 @@ static modules_t Modules[] =
   { Consecstat,     ConsecstatHelp,    ConsecstatOperators,    CDI_REAL,  1,  1 },
   { Copy,           CopyHelp,          CopyOperators,          CDI_REAL, -1,  1 },
   { Deltime,        NULL,              DeltimeOperators,       CDI_REAL,  1,  1 },
-  { Derivepar,      NULL,              DeriveparOperators,     CDI_REAL,  1,  1 },
+  { Derivepar,      DeriveparHelp,     DeriveparOperators,     CDI_REAL,  1,  1 },
   { Detrend,        DetrendHelp,       DetrendOperators,       CDI_REAL,  1,  1 },
   { Diff,           DiffHelp,          DiffOperators,          CDI_REAL,  2,  0 },
   { Duplicate,      DuplicateHelp,     DuplicateOperators,     CDI_REAL,  1,  1 },

src/operator_help.h

@@ -4070,6 +4070,35 @@ static char *MastrfuHelp[] = {
     NULL
 };
 
+static char *DeriveparHelp[] = {
+    "NAME",
+    "    sealevelpressure - Sea level pressure",
+    "",
+    "SYNOPSIS",
+    "    sealevelpressure  ifile ofile",
+    "",
+    "DESCRIPTION",
+    "    This operator computes the sea level pressure (air_pressure_at_sea_level). Required input fields",
+    "    are surface_air_pressure and air_temperature on hybrid sigma pressure levels.",
+    "",
+    "NOTE",
+    "    The netCDF CF convention for vertical hybrid coordinates is not supported, yet!",
+    "    The vertical coordinate definition of a netCDF file has to follow the ECHAM convention.",
+    "    This means:",
+    "    - the dimension of the full level coordinate and the corresponding variable is called mlev,",
+    "    - the dimension of the half level coordinate and the corresponding variable is called ilev (ilev must have one element more than mlev)",
+    "    - the hybrid vertical coefficient a is given in units of Pa and called hyai (hyam for level midpoints)",
+    "    - the hybrid vertical coefficient b is given in units of 1 and called hybi (hybm for level midpoints)",
+    "    - the mlev variable has a borders attribute containing the character string 'ilev'",
+    "    ",
+    "    Use the sinfo command to test if your vertical coordinate system is recognized as hybrid system.",
+    "    ",
+    "    In case this operator complains about not finding any data on hybrid model levels you may wish",
+    "    to use the setzaxis command to generate a zaxis description which conforms to the ECHAM convention.",
+    "    See section \"1.4 Z-axis description\" for an example how to define a hybrid Z-axis.",
+    NULL
+};
+
 static char *AdisitHelp[] = {
     "NAME",
     "    adisit, adipot - Potential temperature to in-situ temperature and vice versa",