Commit b1a9775d authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

Outputgmt: docu update

parent 35fce878
......@@ -71,6 +71,8 @@ doc/tex/eca_int.tex -text
doc/tex/eca_ref.tex -text
doc/tex/figures/collgrid.pdf -text
doc/tex/figures/distgrid.pdf -text
doc/tex/figures/gmtgrfill.pdf -text svneol=unset#application/pdf
doc/tex/figures/gmtshaded.pdf -text svneol=unset#application/pdf
doc/tex/figures/maggrfill.pdf -text svneol=unset#application/pdf
doc/tex/figures/magshaded.pdf -text svneol=unset#application/pdf
doc/tex/figures/plotcollgrid -text
......
......@@ -501,6 +501,8 @@ Operator catalog:
Output outputsrv SERVICE ASCII output
Output outputext EXTRA ASCII output
Outputtab outputtab Table output
Outputgmt outputcenter Output grid center
Outputgmt outputbounds Output grid bounds
-------------------------------------------------------------
Miscellaneous
-------------------------------------------------------------
......
......@@ -116,6 +116,7 @@ Importamsr Import/Export
Input Import/Export
Output Import/Export
Outputtab Import/Export
Outputgmt Import/Export
Gradsdes Miscellaneous
Afterburner Miscellaneous
Filter Miscellaneous
......
......@@ -26,3 +26,38 @@
@EndDescription
@EndOperator
@BeginExample
1) GMT shaded contour plot of a global temperature field with a resolution of 4 degree.
The contour interval is 3 with a rainbow color table.
@BeginVerbatim
cdo outputcenter -temp,global_4 > data.gmt
makecpt -T213/318/3 -Crainbow > gmt.cpt
pscontour -K -JQ0/10i -Rd -I -Cgmt.cpt data.gmt > gmtplot.ps
pscoast -O -J -R -Dc -W -B40g20 >> gmtplot.ps
@EndVerbatim
@IfDoc
\vspace{1mm}
\includegraphics[width=\textwidth]{gmtshaded.pdf}
\vspace{1mm}
@EndifDoc
2) GMT shaded gridfill plot of a global temperature field with a resolution of 4 degree.
The contour interval is 3 with a rainbow color table.
@BeginVerbatim
cdo outputbounds -temp,global_4 > data.gmt
makecpt -T213/318/3 -Crainbow > gmt.cpt
psxy -K -JQ0/10i -Rd -L -Cgmt.cpt -m data.gmt > gmtplot.ps
pscoast -O -J -R -Dc -W -B40g20 >> gmtplot.ps
@EndVerbatim
@IfDoc
\vspace{1mm}
\includegraphics[width=\textwidth]{gmtgrfill.pdf}
\vspace{1mm}
@EndifDoc
@EndExample
......@@ -509,10 +509,9 @@ void *Outputgmt(void *argument)
double minval, maxval, meanval;
array_stat(gridsize, array, missval, &minval, &maxval, &meanval);
double range = maxval - minval;
fprintf(stderr, "minval, maxval, meanval, range %g %g %g %g\n", minval, maxval, meanval, range);
fprintf(stderr, "makecpt -T%g/%g/%g -Crainbow > gmt.cpt\n", minval, maxval, range/20);
fprintf(stderr, "pscontour -K -Jx0.02id -Rd -I -Cgmt.cpt data.gmt > gmtplot.ps\n");
fprintf(stderr, "pscoast -O -J -R -Dc -W -B30g30 >> gmtplot.ps\n");
fprintf(stderr, "pscontour -K -JQ0/10i -Rd -I -Cgmt.cpt data.gmt > gmtplot.ps\n");
fprintf(stderr, "pscoast -O -J -R -Dc -W -B40g20 >> gmtplot.ps\n");
}
for ( i = 0; i < nvals; i++ )
......@@ -613,12 +612,11 @@ void *Outputgmt(void *argument)
double minval, maxval, meanval;
array_stat(gridsize, array, missval, &minval, &maxval, &meanval);
double range = maxval - minval;
fprintf(stderr, "minval, maxval, meanval, range %g %g %g %g\n", minval, maxval, meanval, range);
fprintf(stderr, "makecpt -T%g/%g/%g -Crainbow > gmt.cpt\n", minval, maxval, range/20);
fprintf(stderr, "psxy -K -JQ0/10i -Rd -L -Cgmt.cpt -m data.gmt > gmtplot.ps\n");
// fprintf(stderr, "psxy -K -Jx0.028id -Rd -L -Cgmt.cpt -m data.gmt > gmtplot.ps\n");
// fprintf(stderr, "psxy -K -JN0/10i -Rd -L -Cgmt.cpt -m data.gmt > gmtplot.ps\n");
fprintf(stderr, "pscoast -O -J -R -Dc -W -B30g30 >> gmtplot.ps\n");
fprintf(stderr, "pscoast -O -J -R -Dc -W -B40g20 >> gmtplot.ps\n");
fprintf(stderr, "ps2pdf gmtplot.ps\n");
}
......
......@@ -651,7 +651,7 @@ static modules_t Modules[] =
{ Nmltest, NULL, NmltestOperators, 0, CDI_REAL, 0, 0 },
{ Output, OutputHelp, OutputOperators, 1, CDI_REAL, -1, 0 },
{ Output, OutputtabHelp, OutputtabOperators, 1, CDI_REAL, -1, 0 },
{ Outputgmt, NULL, OutputgmtOperators, 1, CDI_REAL, 1, 0 },
{ Outputgmt, OutputgmtHelp, OutputgmtOperators, 1, CDI_REAL, 1, 0 },
{ Pack, NULL, PackOperators, 1, CDI_REAL, 1, 1 },
{ Pardup, NULL, PardupOperators, 1, CDI_REAL, 1, 1 },
{ Pinfo, NULL, PinfoOperators, 1, CDI_REAL, 1, 1 },
......
......@@ -4244,6 +4244,21 @@ static const char *OutputtabHelp[] = {
NULL
};
static const char *OutputgmtHelp[] = {
"NAME",
" outputcenter, outputbounds - Formatted GMT output",
"",
"SYNOPSIS",
" <operator> ifiles",
"",
"DESCRIPTION",
"",
"OPERATORS",
" outputcenter Output grid center",
" outputbounds Output grid bounds",
NULL
};
static const char *GradsdesHelp[] = {
"NAME",
" gradsdes - GrADS data descriptor file",
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment