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

20030616: Average field diffulgleties

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.

  • 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--
>