[
Date Prev][
Date Next][
Thread Prev][
Thread Next][
Date Index][
Thread Index]
20050505: gptraj, gpanot, gf driver - is there a better way to overlay images?
- Subject: 20050505: gptraj, gpanot, gf driver - is there a better way to overlay images?
- Date: Fri, 06 May 2005 14:44:54 -0600
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>
><!-- Converted from text/plain format -->
>
><P><FONT SIZE=3D2>Hello,<BR>
><BR>
>I have been using gptraj to create reverse trajectory profiles. Here is =
>an example:<BR>
><A =
>HREF=3D"http://climate.sdstate.edu/downloads/gptrajnew/043000361000.gif";>=
>http://climate.sdstate.edu/downloads/gptrajnew/043000361000.gif</A><BR>
><BR>
>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:<BR>
><A =
>HREF=3D"http://climate.sdstate.edu/downloads/gptrajnew/abc.gif";>http://cl=
>imate.sdstate.edu/downloads/gptrajnew/abc.gif</A><BR>
><BR>
>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.<BR>
><BR>
>Here goes the script:<BR>
><BR>
>#!/bin/csh<BR>
><BR>
># trajectory<BR>
># USAGE: ./gptrajnew 925 12 00 for 925mb 12hour forecast and 00Z =
>file<BR>
># USAGE: ./gptrajnew 0 48 12 for 0mb 48hour forecast and 12Z file<BR>
><BR>
>source $GEMPAKSTART/Gemenviron<BR>
>cd /home/gempak/unidataldm/ldm/gempak/model/<BR>
><BR>
>setenv GEMCOLTBL coltbl.xwp<BR>
><BR>
>set utcdate =3D `date +%Y%m%d`<BR>
>#set utcdate =3D 20050503<BR>
>set displaydate =3D `echo $utcdate | awk '{print substr($0,5)}'`<BR>
>set zvalue =3D $3<BR>
>set hourvalue =3D $2<BR>
>set height =3D $1<BR>
>set title =3D "File: ${utcdate} ${zvalue}Z ${hourvalue}hr forecast =
>at ${height}mb"<BR>
><BR>
>set filen =3D =
>"/home/gempak/web_on_lancer/climate/downloads/gptrajnew/"${disp=
>laydate}${zvalue}${hourvalue}${height}".gif"<BR>
>#set filen =3D =
>"/home/gempak/web_on_lancer/climate/downloads/gptrajnew/abc.gif"=
>;<BR>
><BR>
>set DEVICE =3D "GF|"${filen}"|750;550"<BR>
><BR>
>if ( -e ${filen} ) then<BR>
> rm ${filen}<BR>
>endif<BR>
><BR>
>gptraj << EOF > /home/gempak/scripts/gptraj/gptrajnew_out<BR>
> MAP =
> =
>=3D 10<BR>
> GAREA =
>=3D uslcc<BR>
> PROJ =
>=3D lcc<BR>
> SATFILE =3D<BR>
> RADFIL =3D<BR>
> LUTFIL =3D<BR>
> IMCBAR =3D<BR>
> LATLON =
>=3D 11<BR>
> PANEL =
>=3D 0<BR>
> TITLE =
>=3D 12/0/${title}<BR>
> TEXT =
>=3D 1/21/1/hw<BR>
> CLEAR =
>=3D YES<BR>
> !DEVICE =
>=3D ${DEVICE}<BR>
> GDFILE =
>=3D =
>/home/gempak/unidataldm/ldm/gempak/model/${utcdate}${3}_eta212.gem<BR>
> GVECT =
>=3D WND<BR>
> GPOINT =
>=3D 43/46/1;-98/-96/1<BR>
> GLEVEL =
>=3D $1<BR>
> GVCORD =
>=3D PRES<BR>
> GDATTIM =
>=3D F0-F$2<BR>
> TSTEP =3D<BR>
> MARKER =
>=3D 0<BR>
> LINE =
>=3D 4<BR>
> RTRAJ =
>=3D YES<BR>
> run<BR>
><BR>
> exit<BR>
>EOF<BR>
><BR>
>foreach i (43 44 45 46)<BR>
> foreach j (-98 -97 -96)<BR>
> =
> set loci1 =3D =
>"#"${i}";"${j}<BR>
> =
> =
> gpanot << EOF1 >> =
>/home/gempak/scripts/gptraj/gptrajnew_out<BR>
> =
> =
> =
> GDFILE =
>=3D <BR>
> =
> =
> =
> SATFIL =3D<BR>
> =
> =
> =
> RADFIL =3D<BR>
> =
> =
> =
> PROJ =
>=3D lcc<BR>
> =
> =
> =
> GAREA =
>=3D uslcc<BR>
> =
> =
> =
> PANEL =
>=3D 0<BR>
> =
> =
> =
> DEVICE =
>=3D ${DEVICE}<BR>
> =
> =
> =
> CLEAR =
>=3D no<BR>
> =
> =
> =
> SHAPE =
>=3D marker<BR>
> =
> =
> =
> INFO =
>=3D 14/1<BR>
> =
> =
> =
> LOCI =
>=3D ${loci1}<BR>
> =
> =
> =
> LINE =
>=3D 4<BR>
> =
> =
> =
> CTYPE =
>=3D <BR>
> =
> =
> =
> r<BR>
> =
> =
> =
> <BR>
> =
> =
> =
> e<BR>
> =
> =
> EOF1<BR>
><BR>
> end<BR>
>end<BR>
><BR>
>gpend<BR>
><BR>
><BR>
><BR>
>Thank you very much.<BR>
><BR>
>Sincerely,<BR>
>Chirag Shukla<BR>
>Climate Data Specialist<BR>
>Box 2120, Ag Engineering<BR>
>South Dakota State University<BR>
>Brookings, SD 57007<BR>
>(605) 688-5682<BR>
><A =
>HREF=3D"http://climate.sdstate.edu";>http://climate.sdstate.edu</A><BR>
><BR>
><BR>
><BR>
><BR>
><BR>
><BR>
></FONT>
></P>
>
></BODY>
></HTML>
>------_=_NextPart_001_01C551F4.D463BBD9--
>
--
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.
![http://climate.sdstate.edu/downloads/gptrajnew/043000361000.gif"](http://climate.sdstate.edu/downloads/gptrajnew/043000361000.gif"){kind=link}
![http://climate.sdstate.edu/downloads/gptrajnew/abc.gif"](http://climate.sdstate.edu/downloads/gptrajnew/abc.gif"){kind=link}