NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
Steve, Looks like the problem is solved now -- I had GPOINT=43/46/1;-96/-100/1 Changed to GPOINT=43/46/1;-100/-96/1 Thanks. Sincerely, Chirag. -----Original Message----- From: SHUKLA, CHIRAG Sent: Friday, May 13, 2005 10:18 AM To: Unidata Support Subject: 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> ><!-- 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/"${di >sp= laydate}${zvalue}${hourvalue}${height}".gif"<BR> >#set filen =3D = >"/home/gempak/web_on_lancer/climate/downloads/gptrajnew/abc.gif&qu >ot= >;<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-- > -- ************************************************************************ **** < 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.