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

RE: 20050505: gptraj, gpanot, gf driver - is there a better way to overlay images?

Subject: RE: 20050505: gptraj, gpanot, gf driver - is there a better way to overlay images?
Date: Fri, 13 May 2005 10:18:07 -0500

Title: RE: 20050505: gptraj, gpanot, gf driver - is there a better way to overlay images?

Steve,

All of a sudden, I am getting an error in gptraj:
[TG 0] F0-F12
Error specifying lat/lon/increment:
43/46/1;-96/-100/1

What could be wrong? When I see the processes, I dont see gplt process and everything else seems alright. A few days back everything was running okay and nothing changed on the system. What would you suggest?

Thanks.

Sincerely,
Chirag Shukla

Here is the output of the cshell script:

MAP Map color/dash/width/filter flag 10
GAREA Graphics area uslcc
PROJ Map projection/angles/margins|dr lcc
SATFIL Satellite image filename(s)
RADFIL Radar image filename(s)
LUTFIL Enhancement lookup table filenam
IMCBAR Color/ornt/anch/x;y/ln;wd/freq|t
LATLON Line color/dash/width/freq/inc/l 11
PANEL Panel loc/color/dash/width/regn 0
TITLE Title color/line/title 12/0/File: 20050513 00Z 12hr forecast at 850mb
TEXT Size/fnt/wdth/brdr/N-rot/just/hw 1/21/1/hw
CLEAR Clear screen flag YES
DEVICE Device|name|x size;y size|color GIF|/home/gempak/scripts/gptraj/eastriver/05130012850.gif|750;550
GDFILE Grid file /home/gempak/unidataldm/ldm/gempak/model/2005051300_eta212.gem
GVECT Vector grid WND
GPOINT Grid point 43/46/1;-96/-100/1
GLEVEL Grid level 850
GVCORD Grid vertical coordinate PRES
GDATTIM Grid date/time F0-F12
TSTEP Time step (minutes)
MARKER Marker color/type/size/width/hw 0
LINE Color/type/width/label/smth/fltr 4
RTRAJ Back/Reverse trajectory flag YES
Parameters requested: MAP,GAREA,PROJ,SATFIL,RADFIL,LUTFIL,IMCBAR,LATLON,
PANEL,TITLE,TEXT,CLEAR,DEVICE,GDFILE,GVECT,GPOINT,GLEVEL,GVCORD,GDATTIM,
TSTEP,MARKER,LINE,RTRAJ.
GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GEMPAK-GPTRAJ> GPMAP PARAMETERS:

Device: GIF|/home/gempak/scripts/gptraj/eastriver/05130012850.gif|750;550
Projection: lcc
Graphics area name: uslcc
Map: 10
Title: 12/0/File: 20050513 00Z 12hr forecast at 850mb
Panel: 0
Latlon: 11
Grid time range: F0-F12
Grid file: /home/gempak/unidataldm/ldm/gempak/model/2005051300_eta212.gem
Grid Vector: WND
Grid level: 850
Vertical Coordinate: PRES
Back Trajectory: YES
Time Step:
Clear: YES
[TG 0] F0-F12
Error specifying lat/lon/increment:
43/46/1;-96/-100/1

-----Original Message-----
From: Unidata Support [mailto:address@hidden]
Sent: Fri 5/6/2005 3:44 PM
To: SHUKLA, CHIRAG
Cc: address@hidden
Subject: 20050505: gptraj, gpanot, gf driver - is there a better way to overlay images?

Chirag,

The "southward" facing arrow is a new bug introduced by the
gemlib routine pd_drct() calling sequence being changed
and wasn't updated in the ggtraj() routine of gptraj.
I have fixed this now that you called my attention to the problem.

The arrow should be pointing in the direction of the last vector
in the trajectory.

If you built the distribution locally, I can send you the changes.
Otherwise, your need a new binary.

Steve Chiswell
Unidata User Support

>From: "SHUKLA, CHIRAG" <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200505060434.j464Y9P3010877

>This is a multi-part message in MIME format.
>
>------_=_NextPart_001_01C551F4.D463BBD9
>Content-Type: text/plain;
> charset="iso-8859-1"
>Content-Transfer-Encoding: quoted-printable
>
>Hello,
>
>I have been using gptraj to create reverse trajectory profiles. Here is =
>an example:
>http://climate.sdstate.edu/downloads/gptrajnew/043000361000.gif
>
>I was wondering if we had an arrow of some kind that would indicate the =
>path of the trajectory. I tried using markers. It looked too crowded =
>sometimes and an arrow facing south appeared at each gridpoint. We didnt =
>want the arrow facing south. Other option was pursued - the use of =
>gpanot recursively to put a little star at each gridpoint location so as =
>to indicate where the end-points were. Here is an example of the image =
>we wanted:=20
>http://climate.sdstate.edu/downloads/gptrajnew/abc.gif
>
>The script I wrote follows soon. The question is -- is there a better =
>way of doing this? Another question is -- when I tried using GIF driver, =
>trajectories were created successfully but the stars (using gpanot) were =
>not created on the map. All I did was replace GF with GIF. Was I missing =
>something --or-- was GIF driver not appropriate for overlays? When I =
>used GF driver, everything seemed to fall in place.
>
>Here goes the script:
>
>#!/bin/csh
>
># trajectory
># USAGE: ./gptrajnew 925 12 00 for 925mb 12hour forecast and 00Z file
># USAGE: ./gptrajnew 0 48 12 for 0mb 48hour forecast and 12Z file
>
>source $GEMPAKSTART/Gemenviron
>cd /home/gempak/unidataldm/ldm/gempak/model/
>
>setenv GEMCOLTBL coltbl.xwp
>
>set utcdate =3D `date +%Y%m%d`
>#set utcdate =3D 20050503
>set displaydate =3D `echo $utcdate | awk '{print substr($0,5)}'`
>set zvalue =3D $3
>set hourvalue =3D $2
>set height =3D $1
>set title =3D "File: ${utcdate} ${zvalue}Z ${hourvalue}hr forecast at =
>${height}mb"
>
>set filen =3D =
>"/home/gempak/web_on_lancer/climate/downloads/gptrajnew/"${displaydate}${=
>zvalue}${hourvalue}${height}".gif"
>#set filen =3D =
>"/home/gempak/web_on_lancer/climate/downloads/gptrajnew/abc.gif"
>
>set DEVICE =3D "GF|"${filen}"|750;550"
>
>if ( -e ${filen} ) then
> rm ${filen}
>endif
>
>gptraj << EOF > /home/gempak/scripts/gptraj/gptrajnew_out
> MAP =3D 10
> GAREA =3D uslcc
> PROJ =3D lcc
> SATFILE =3D
> RADFIL =3D
> LUTFIL =3D
> IMCBAR =3D
> LATLON =3D 11
> PANEL =3D 0
> TITLE =3D 12/0/${title}
> TEXT =3D 1/21/1/hw
> CLEAR =3D YES
> !DEVICE =3D ${DEVICE}
> GDFILE =3D =
>/home/gempak/unidataldm/ldm/gempak/model/${utcdate}${3}_eta212.gem
> GVECT =3D WND
> GPOINT =3D 43/46/1;-98/-96/1
> GLEVEL =3D $1
> GVCORD =3D PRES
> GDATTIM =3D F0-F$2
> TSTEP =3D
> MARKER =3D 0
> LINE =3D 4
> RTRAJ =3D YES
> run
>
> exit
>EOF
>
>foreach i (43 44 45 46)
> foreach j (-98 -97 -96)
> set loci1 =3D "#"${i}";"${j}
> gpanot << EOF1 >> /home/gempak/scripts/gptraj/gptrajnew
> _out
> GDFILE =3D=09
> SATFIL =3D
> RADFIL =3D
> PROJ =3D lcc
> GAREA =3D uslcc
> PANEL =3D 0
> DEVICE =3D ${DEVICE}
> CLEAR =3D no
> SHAPE =3D marker
> INFO =3D 14/1
> LOCI =3D ${loci1}
> LINE =3D 4
> CTYPE =3D=09
> r
> =09
> e
> EOF1
>
> end
>end
>
>gpend
>
>
>
>Thank you very much.
>
>Sincerely,
>Chirag Shukla
>Climate Data Specialist
>Box 2120, Ag Engineering
>South Dakota State University
>Brookings, SD 57007
>(605) 688-5682
>http://climate.sdstate.edu
>
>
>
>
>
>
>
>------_=_NextPart_001_01C551F4.D463BBD9
>Content-Type: text/html;
> charset="iso-8859-1"
>Content-Transfer-Encoding: quoted-printable
>
><!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
><HTML>
><HEAD>
><META HTTP-EQUIV=3D"Content-Type" CONTENT=3D"text/html; =
>charset=3Diso-8859-1">
><META NAME=3D"Generator" CONTENT=3D"MS Exchange Server version =
>6.5.7232.36">
><TITLE>gptraj, gpanot, gf driver - is there a better way to overlay =
>images?</TITLE>
></HEAD>
><BODY>
>
>
>Hello, 
> 
>I have been using gptraj to create reverse trajectory profiles. Here is =
>an example: 
><A =
>HREF="" href="http://climate.sdstate.edu/downloads/gptrajnew/043000361000.gif">http://climate.sdstate.edu/downloads/gptrajnew/043000361000.gif">=
>http://climate.sdstate.edu/downloads/gptrajnew/043000361000.gif</A> 
> 
>I was wondering if we had an arrow of some kind that would indicate the =
>path of the trajectory. I tried using markers. It looked too crowded =
>sometimes and an arrow facing south appeared at each gridpoint. We didnt =
>want the arrow facing south. Other option was pursued - the use of =
>gpanot recursively to put a little star at each gridpoint location so as =
>to indicate where the end-points were. Here is an example of the image =
>we wanted: 
><A =
>HREF="" href="http://climate.sdstate.edu/downloads/gptrajnew/abc.gif">http://climate.sdstate.edu/downloads/gptrajnew/abc.gif">http://cl=
>imate.sdstate.edu/downloads/gptrajnew/abc.gif</A> 
> 
>The script I wrote follows soon. The question is -- is there a better =
>way of doing this? Another question is -- when I tried using GIF driver, =
>trajectories were created successfully but the stars (using gpanot) were =
>not created on the map. All I did was replace GF with GIF. Was I missing =
>something --or-- was GIF driver not appropriate for overlays? When I =
>used GF driver, everything seemed to fall in place. 
> 
>Here goes the script: 
> 
>#!/bin/csh 
> 
># trajectory 
># USAGE: ./gptrajnew 925 12 00 for 925mb 12hour forecast and 00Z =
>file 
># USAGE: ./gptrajnew 0 48 12 for 0mb 48hour forecast and 12Z file 
> 
>source $GEMPAKSTART/Gemenviron 
>cd /home/gempak/unidataldm/ldm/gempak/model/ 
> 
>setenv GEMCOLTBL coltbl.xwp 
> 
>set utcdate =3D `date +%Y%m%d` 
>#set utcdate =3D 20050503 
>set displaydate =3D `echo $utcdate | awk '{print substr($0,5)}'` 
>set zvalue =3D $3 
>set hourvalue =3D $2 
>set height =3D $1 
>set title =3D "File: ${utcdate} ${zvalue}Z ${hourvalue}hr forecast =
>at ${height}mb" 
> 
>set filen =3D =
>"/home/gempak/web_on_lancer/climate/downloads/gptrajnew/"${disp=
>laydate}${zvalue}${hourvalue}${height}".gif" 
>#set filen =3D =
>"/home/gempak/web_on_lancer/climate/downloads/gptrajnew/abc.gif&quot=
>; 
> 
>set DEVICE =3D "GF|"${filen}"|750;550" 
> 
>if ( -e ${filen} ) then 
>        rm ${filen} 
>endif 
> 
>gptraj << EOF > /home/gempak/scripts/gptraj/gptrajnew_out 
>        MAP     =
>        =
>=3D       10 
>        GAREA   =
>=3D       uslcc 
>        PROJ    =
>=3D       lcc 
>        SATFILE =3D 
>        RADFIL  =3D 
>        LUTFIL  =3D 
>        IMCBAR  =3D 
>        LATLON  =
>=3D       11 
>        PANEL   =
>=3D       0 
>        TITLE   =
>=3D       12/0/${title} 
>        TEXT    =
>=3D       1/21/1/hw 
>        CLEAR   =
>=3D       YES 
>        !DEVICE =
>=3D       ${DEVICE} 
>        GDFILE  =
>=3D       =
>/home/gempak/unidataldm/ldm/gempak/model/${utcdate}${3}_eta212.gem 
>        GVECT   =
>=3D       WND 
>        GPOINT  =
>=3D       43/46/1;-98/-96/1 
>        GLEVEL  =
>=3D       $1 
>        GVCORD  =
>=3D       PRES 
>        GDATTIM =
>=3D       F0-F$2 
>        TSTEP   =3D 
>        MARKER  =
>=3D       0 
>        LINE    =
>=3D       4 
>        RTRAJ   =
>=3D       YES 
>        run 
> 
>        exit 
>EOF 
> 
>foreach i (43 44 45 46) 
>        foreach j (-98 -97 -96) 
>        =
>        set loci1 =3D =
>"#"${i}";"${j} 
>        =
>        =
>        gpanot << EOF1 >> =
>/home/gempak/scripts/gptraj/gptrajnew_out 
>        =
>        =
>        =
>        GDFILE  =
>=3D       
>        =
>        =
>        =
>        SATFIL  =3D 
>        =
>        =
>        =
>        RADFIL  =3D 
>        =
>        =
>        =
>        PROJ    =
>=3D       lcc 
>        =
>        =
>        =
>        GAREA   =
>=3D       uslcc 
>        =
>        =
>        =
>        PANEL   =
>=3D       0 
>        =
>        =
>        =
>        DEVICE  =
>=3D       ${DEVICE} 
>        =
>        =
>        =
>        CLEAR   =
>=3D       no 
>        =
>        =
>        =
>        SHAPE   =
>=3D       marker 
>        =
>        =
>        =
>        INFO    =
>=3D       14/1 
>        =
>        =
>        =
>        LOCI    =
>=3D       ${loci1} 
>        =
>        =
>        =
>        LINE    =
>=3D       4 
>        =
>        =
>        =
>        CTYPE   =
>=3D       
>        =
>        =
>        =
>        r 
>        =
>        =
>        =
>        
>        =
>        =
>        =
>        e 
>        =
>        =
>        EOF1 
> 
>        end 
>end 
> 
>gpend 
> 
> 
> 
>Thank you very much. 
> 
>Sincerely, 
>Chirag Shukla 
>Climate Data Specialist 
>Box 2120, Ag Engineering 
>South Dakota State University 
>Brookings, SD 57007 
>(605) 688-5682 
><A =
>HREF="" href="http://climate.sdstate.edu">http://climate.sdstate.edu">http://climate.sdstate.edu</A> 
> 
> 
> 
> 
> 
> 
>
>
>
></BODY>
></HTML>
>------_=_NextPart_001_01C551F4.D463BBD9--
>
--
**************************************************************************** <
Unidata User Support UCAR Unidata Program <
(303)497-8643 P.O. Box 3000 <
address@hidden Boulder, CO 80307 <
---------------------------------------------------------------------------- <
Unidata WWW Service http://my.unidata.ucar.edu/content/support <
---------------------------------------------------------------------------- <
NOTE: All email exchanges with Unidata User Support are recorded in the
Unidata inquiry tracking system and then made publicly available
through the web. If you do not want to have your interactions made
available in this way, you must let us know in each email you send to us.

Prev by Date: 20050512: Unidata U.S. Radar Reflectivity Composite Data
Next by Date: RE: 20050505: gptraj, gpanot, gf driver - is there a better way to overlay images?
Previous by thread: 20050505: gptraj, gpanot, gf driver - is there a better way to overlay images?
Next by thread: RE: 20050505: gptraj, gpanot, gf driver - is there a better way to overlay images?
Index(es):
- Date
- Thread