[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

20020301: GPFRONT and Satellite Overlay



Patrick,

It sounds like some memory might be getting stepped on
if the satellite file is reported "does not exist",
or the projection read from the satellite image is not 
as expected. Either are possible if the program was linked 
against a different version of the libraries.

What version of GEMPAK are you running? Type "version"
at the gpfront prompt to make sure it agrees with your
installed distribution.

Steve Chiswell
Unidata User Support




>From: "Patrick O'Reilly" <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200203012258.g21MwLx25795

>This is a multi-part message in MIME format.
>
>------=_NextPart_000_0031_01C1C142.4BAF3360
>Content-Type: text/plain;
>       charset="iso-8859-1"
>Content-Transfer-Encoding: quoted-printable
>
>Hi again,
>
>I have been fighting with gpfront to get an overlay of IR satellite and =
>fronts.  I already have a script running successfully creating an image =
>of the 6km composite radar and fronts.  When I am physically sitting at =
>the computer, entering commands and using gpfront, I can get the desired =
>image.  But when I run gpfront in a script, no dice.  I have tried two =
>approaches  (with many, many variations).
>
>I have tried running gpmap to create the satellite image, then running =
>gpfront to overlay the fronts, but no dice that way.  I get the =
>satellite image created, but get an error ( [IM -1] Image file ... does =
>not exist)   from gpfront and no fronts are drawn. =20
>
>The second approach is simply running gpfront and specifying a satellite =
>file.  This method produces no  error messages, but strange plots.  One =
>such plot is attached.  It looks like it plotted something in regards to =
>the satellite file, as the IR colorbar is there.  The clip from the =
>script:
>
># Plot the Fronts/Highs/Lows
>
>gpfront<< GPFRONT_FLAG
>satfil=3D$GEMDATA/images/sat/GOES-8/4km/IR/$satfil
>map=3D31
>garea=3D23.0;-120.0;47.0;-65.0
>proj=3Dsat//1;1;1;1
>radfil=3D
>latlon=3D0
>panel=3D0
>text=3D1.5/22/HW
>clear=3Dyes
>title=3D32/-1/Fronts at $filetime
>asfil=3D$GEMDATA/nwx/hpc/fronts/$frontfile
>device=3Dgf|$tmpgif|900;650
>ashr=3D
>asatt=3D1/N
>
>run
>
>exit
>
>GPFRONT_FLAG
>
>I can see nothing wrong with this, if you do let me know.  There aren't =
>errors finding satellite or front files, and all the variables work =
>fine.  I can't understand why it works sitting entering commands running =
>gempak interactively, but when I try the SAME thing in a script, I get a =
>hard time.  I found very little in the support email archives about =
>gpfront.
>
>Thanks once again for your help!  I am trying to get some requested =
>imagery put together for posting to the web, and have about 95% of it =
>done.=20
>
>Patrick
>
>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>Patrick O'Reilly                               Support Scientist
>The STORM Project            address@hidden
>208 Latham Hall                             ph: 319-273-3789
>University of Northern Iowa
>Cedar Falls, IA 50614            =20
>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  =20
>
>
>------=_NextPart_000_0031_01C1C142.4BAF3360
>Content-Type: text/html;
>       charset="iso-8859-1"
>Content-Transfer-Encoding: quoted-printable
>
><!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
><HTML><HEAD>
><META http-equiv=3DContent-Type content=3D"text/html; =
>charset=3Diso-8859-1">
><META content=3D"MSHTML 5.50.4522.1800" name=3DGENERATOR>
><STYLE></STYLE>
></HEAD>
><BODY bgColor=3D#efefef>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>Hi again,</FONT></DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>I have been fighting with =
>gpfront to get=20
>an overlay of IR satellite and fronts.&nbsp; I already have&nbsp;a =
>script=20
>running successfully creating an image&nbsp;of the 6km composite radar =
>and=20
>fronts.&nbsp; When I am physically sitting at the computer, entering =
>commands=20
>and using gpfront, I can get&nbsp;the desired image.&nbsp; But when I =
>run=20
>gpfront in a script, no dice.&nbsp; I have tried two approaches&nbsp; =
>(with=20
>many, many variations).</FONT></DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>I have tried running gpmap to =
>create the=20
>satellite image, then running gpfront to overlay the fronts, but no dice =
>that=20
>way.&nbsp; I get the satellite image created, but get an error ( [IM -1] =
>Image=20
>file ... does not exist)&nbsp;&nbsp;&nbsp;from gpfront and no fronts are =
>
>drawn.&nbsp; </FONT></DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>The second approach is simply =
>running=20
>gpfront and specifying a satellite file.&nbsp; This method produces =
>no&nbsp;=20
>error messages, but strange plots.&nbsp; One such plot is =
>attached.&nbsp; It=20
>looks like it plotted something in regards to the satellite file, as the =
>IR=20
>colorbar is there.&nbsp; The clip from the script:</FONT></DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2><BR># Plot the=20
>Fronts/Highs/Lows</FONT></DIV>
><DIV>&nbsp;</DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>gpfront&lt;&lt;=20
>GPFRONT_FLAG<BR>satfil=3D$GEMDATA/images/sat/GOES-8/4km/IR/$satfil<BR>map=
>=3D31<BR>garea=3D23.0;-120.0;47.0;-65.0<BR>proj=3Dsat//1;1;1;1<BR>radfil=3D=
><BR>latlon=3D0<BR>panel=3D0<BR>text=3D1.5/22/HW<BR>clear=3Dyes<BR>title=3D=
>32/-1/Fronts=20
>at=20
>$filetime<BR>asfil=3D$GEMDATA/nwx/hpc/fronts/$frontfile<BR>device=3Dgf|$t=
>mpgif|900;650<BR>ashr=3D<BR>asatt=3D1/N</FONT></DIV>
><DIV>&nbsp;</DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>run</FONT></DIV>
><DIV>&nbsp;</DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>exit</FONT></DIV>
><DIV>&nbsp;</DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>GPFRONT_FLAG<BR></FONT></DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>I can see nothing wrong with =
>this, if you=20
>do let me know.&nbsp; There aren't errors finding satellite or front =
>files, and=20
>all the variables work fine.&nbsp; I can't understand why it works =
>sitting=20
>entering commands running gempak interactively, but when I try the SAME =
>thing in=20
>a script, I get&nbsp;a hard time.&nbsp; I found very little in the =
>support email=20
>archives about gpfront.</FONT></DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>Thanks once again for your =
>help!&nbsp; I=20
>am trying to get some requested imagery put together for posting to the =
>web, and=20
>have about 95% of it done. </FONT></DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>Patrick</FONT></DIV>
><DIV><FONT face=3D"MS Sans Serif" size=3D2>&nbsp;</DIV></FONT>
><DIV><FONT face=3D"MS Sans Serif"=20
>size=3D2>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<BR>Patrick=20
>O'Reilly&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
>;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=20
>Support Scientist<BR>The STORM=20
>Project&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
> <A=20
>href=3D"mailto:address@hidden";>address@hidden</A><BR>20=
>8 Latham=20
>Hall&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
>sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
>p;&nbsp;&nbsp;&nbsp;&nbsp;=20
>ph: 319-273-3789<BR>University of Northern Iowa<BR>Cedar Falls, IA=20
>50614&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
>bsp;=20
><BR>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~&nbsp;&nbsp;=20
><BR></FONT></DIV></BODY></HTML>
>
>------=_NextPart_000_0031_01C1C142.4BAF3360--
>