[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
20030616: Average field diffulgleties
- Subject: 20030616: Average field diffulgleties
- Date: Mon, 16 Jun 2003 10:38:17 -0600
John,
Your assignment parameter lines below are lacking the "=" character,
>gdpfun mag(vecr((avn_urel,avn_vrel)))//s !s !vecr((avn_urel,avn_vrel))
>type c !f !s
>cint 10 !10 !0
>fint 0 !10 !0
I'm currently out of the office, so haven't looked further past this problem.
Let me know if you encouter other problems.
Steve
>From: John Merrill <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200306161456.h5GEuULd015098
>This is a multi-part message in MIME format.
>--------------080901050305040905020004
>Content-Type: text/plain; charset=us-ascii; format=flowed
>Content-Transfer-Encoding: 7bit
>
>Hello, Steve:
> I am embarrassed to say that I've been frustrated to my limit trying
>to plot maps using a .gem file created from monthly
>average data from the Reanalysis set. You helped me with elementary
>syntax on defining a vector grid recently, and that
>got me as far as isotach contours. (I'd gotten mean height fields on my
>own.) Well, trying to use gdplot2 to create a streamline-
>isotach chart I get innumerable error messages and nothing more than an
>isotach contour plot. Herewith are the script and the resulting
>terminal log. I've used gdplot2 before and don't see what I'm doing
>wrong. I assume that syntax errors are behind this, too. Thanks. John
>
>--------------080901050305040905020004
>Content-Type: text/plain; x-mac-type="0"; x-mac-creator="0";
> name="plot2meanSI200"
>Content-Transfer-Encoding: 7bit
>Content-Disposition: inline;
> filename="plot2meanSI200"
>
>#!/bin/csh
>#
># Create 200 hPa Streamline-Isotach map for April mean, and save.
>#
>source /ljuka/gemadm/gempak/Gemenviron
>#
>set gfile=/afer/rlp/monthlyrean/2001.all.gem
>#
>set gtim=010401/0000f00
>set gout=SI200.ps
>echo "File " $gfile "Output " $gout "Time " $gtim
>gdplot2 << EOF
>gdfile=$gfile
>gdattim=$gtim
>device=ps|$gout
>garea=20;85;60;150
>latlon=32/1/1//15;15
>map=32/1/1
>gvcord=pres
>glevel=200
>skip=0
>wnd=vecr((avn_urel,avn_vrel))
>TITLE=32/-1/~ 200 hPa Height and Isotachs, April 2001 Mean
>TEXT=2/1/3/sw
>border=//3
>CLEAR=yes
>gdpfun mag(vecr((avn_urel,avn_vrel)))//s !s !vecr((avn_urel,avn_vrel))
>type c !f !s
>cint 10 !10 !0
>fint 0 !10 !0
>run
>
>exit
>
>EOF
>#
>gpend
>#
>exit(0)
>
>
>--------------080901050305040905020004
>Content-Type: text/plain; x-mac-type="0"; x-mac-creator="0";
> name="ng_SI.log"
>Content-Transfer-Encoding: 7bit
>Content-Disposition: inline;
filename="ng_SI.log"
>
>aspre being added to access control list
>File /afer/rlp/monthlyrean/2001.all.gem Output SI200.ps Time 010401/0000f00
> GDFILE Grid file /afer/rlp/monthlyrean/2001.all.ge
> m
> GDATTIM Grid date/time 010401/0000f00
> GLEVEL Grid level 200
> GVCORD Grid vertical coordinate pres
> PANEL Panel loc/color/dash/width/regn 0
> SKIP Skip_cntr/skip_plt_x;skip_plt_y 0
> SCALE Scalar scale / vector scale 0
> GDPFUN Scalar grid or vector grid funct mag(vecr(avn_urel,avn_vrel))
> TYPE GDPLOT2 function processing type c
> CONTUR Subbox/smooth 2/2
> CINT Contour interval/min/max 10
> LINE Color/type/width/label/smth/fltr 32/2/4/20;40
> FINT Fill interval/min/max 0
> FLINE Fill colors/fill types 10-20
> HILO Color/symbol/rng/rad/cnt/intp
> HLSYM HILO txt size/posn/font/wdth/hw
> CLRBAR Color/ornt/anch/x;y/ln;wd/freq|t
> WIND Wind symbol/siz/wdth/typ/hdsz BM1
> REFVEC Mag;x;y;txtsiz/font/wdth/HW;labl
> TITLE Title color/line/title 32/-1/~ 200 hPa Height and Isotac
> hs, April 2001 Mean
> TEXT Size/fnt/wdth/brdr/N-rot/just/hw 2/1/3/sw
> CLEAR Clear screen flag yes
> GAREA Graphics area 20;85;60;150
> PROJ Map projection/angles/margins|dr ced
> MAP Map color/dash/width/filter flag 32/1/1
> LATLON Line color/dash/width/freq/inc/l 32/1/1//15;15
> DEVICE Device|name|x size;y size|color ps|SI200.ps
> STNPLT Txtc/txt attr|marker attr|stnfil
> SATFIL Satellite image filename(s)
> RADFIL Radar image filename(s)
> LUTFIL Enhancement lookup table filenam
> STREAM lines/arrows/stop/slow/scale
> POSN Position 0
> COLORS Color list 1
> MARKER Marker color/type/size/width/hw 0
> GRDLBL Grid point label color 0
> FILTER Filter data factor YES
> IMBAR Color/ornt/anch/x;y/ln;wd/freq|t
> Parameters requested: GDFILE,GDATTIM,GLEVEL,GVCORD,PANEL,SKIP,SCALE,GDPFUN,
> TYPE,CONTUR,CINT,LINE,FINT,FLINE,HILO,HLSYM,CLRBAR,WIND,REFVEC,TITLE,TEXT,
> CLEAR,GAREA,PROJ,MAP,LATLON,DEVICE,STNPLT,SATFIL,RADFIL,LUTFIL,STREAM,POSN,
> COLORS,MARKER,GRDLBL,FILTER,IMBAR.
> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2> GEMPAK-GDPLOT
> 2> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2> GEMPAK-GDP
> LOT2> [IP -7] WND is an unrecognized parameter.
> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2> GEMPAK-GDPLOT
> 2> [IP -5] gdpfun mag(vecr((avn_urel,avn_vrel)))//s !s !vecr((avn_ur
> el,avn_vrel)) is an unrecognized command.
> GEMPAK-GDPLOT2> [IP -5] type c !f
> !s is an unrecognized command.
> GEMPAK-GDPLOT2> [IP -5] cint 10 !10
> !0 is an unrecognized command.
> GEMPAK-GDPLOT2> [IP -5] fint 0 !10
> !0 is an unrecognized command.
> GEMPAK-GDPLOT2> Number of times = 1
> 1 010401/0000F000
>
> PARAMETERS FOR GDPLOT : FRAME NUMBER 1 INITIAL PLOT SETTINGS -
> GDFILE=/afer/rlp/monthlyrean/2001.all.gem, GDATTIM=010401/0000F000
> GAREA=20;85;60;150, PROJ=CED, MAP=32/1/1, LATLON=32/1/1//15;15
> PANEL=0, TITLE=32/-1/~ 200 hPa Height and Isotachs, April 2001 Mean
> OVERLAY TYPE LEVEL VCORD GRID FUNCTION
> 1 C 200 PRES MAG(VECR(AVN_UREL,AVN_VREL
> ))
> Parameters requested: GDFILE,GDATTIM,GLEVEL,GVCORD,PANEL,SKIP,SCALE,GDPFUN,
> TYPE,CONTUR,CINT,LINE,FINT,FLINE,HILO,HLSYM,CLRBAR,WIND,REFVEC,TITLE,TEXT,
> CLEAR,GAREA,PROJ,MAP,LATLON,DEVICE,STNPLT,SATFIL,RADFIL,LUTFIL,STREAM,POSN,
> COLORS,MARKER,GRDLBL,FILTER,IMBAR.
> GEMPAK-GDPLOT2> GEMPAK-GDPLOT2>
>--------------080901050305040905020004--
>