[gradsusr] Station data

Zilore Mumba zmumba at gmail.com
Tue Mar 29 01:27:51 EDT 2011


The basic problem with station data is that the length is the header
and the data are not the same. In your fortran open statement do not
specify a RECL. Instead use the option recordtype='stream'
Zilore

On 3/29/11, gradsusr-request at gradsusr.org <gradsusr-request at gradsusr.org> wrote:
> Send gradsusr mailing list submissions to
> 	gradsusr at gradsusr.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> 	http://gradsusr.org/mailman/listinfo/gradsusr
> or, via email, send a message with subject or body 'help' to
> 	gradsusr-request at gradsusr.org
>
> You can reach the person managing the list at
> 	gradsusr-owner at gradsusr.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of gradsusr digest..."
>
>
> Today's Topics:
>
>    1. the 'set ccolor color#' don't take effect for gxout	stnmark
>       (Qian Hao)
>    2. Re: station data (Qian Hao)
>    3. Re: Wind barb (yagnesh)
>    4. Plotting vector of direction Magnitude (Zablone Owiti)
>    5. Request (mani kumari)
>    6. Re: Request (sudev das M P)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 29 Mar 2011 09:33:16 +0800
> From: Qian Hao <yoru.007 at gmail.com>
> Subject: [gradsusr] the 'set ccolor color#' don't take effect for
> 	gxout	stnmark
> To: gradsusr at gradsusr.org
> Message-ID:
> 	<AANLkTik1mSAxK5CwAxo0sEyZwd7hXPeNaHBcYuvomXqr at mail.gmail.com>
> Content-Type: text/plain; charset="gb2312"
>
> Hello, everyone
>
> I found this mismatch between what the grads doc says and its actual
> behavior. The grads html doc presents that
> "For the graphics output type stnmark, the size and style of the mark are
> controled by set
> digsiz<http://www.iges.org/grads/gadoc/gradcomdsetdigsiz.html>
>  and set cmark <http://www.iges.org/grads/gadoc/gradcomdsetcmark.html>. To
> turn off rainbow colorizing and use a single color instead, use set
> ccolor<http://www.iges.org/grads/gadoc/gradcomdsetccolor.html>
>  *color#*."
> http://www.iges.org/grads/gadoc/gradcomdsetgxout.html
> <http://www.iges.org/grads/gadoc/gradcomdsetgxout.html>But after 'set ccolor
> 1', the marks drew by grads were still rainbow colorized. Although i could
> solve this problem by 'set rbcols 1', I think it's necessary to report this
> issue.
> --
> ??
>         ???
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20110329/b7b61c77/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Tue, 29 Mar 2011 09:38:39 +0800
> From: Qian Hao <yoru.007 at gmail.com>
> Subject: Re: [gradsusr] station data
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> 	<AANLkTikQ+tuZMAQcfFBb3gXkpR39RW0Xqh-WNUJek46Q at mail.gmail.com>
> Content-Type: text/plain; charset="gb2312"
>
> I see that there is a blank space in the 'DSET' section of your .ctl
> file.Maybe it's because of this.
>
> 2011/3/28 Filipa Varino <filipavarino at gmail.com>
>
>> Hi,
>>
>> I'm new grads user and i'm having problems with station data.
>>
>> i've created a .txt data with year month day station name and data
>> value(precipitation) for 16 stations during 123 days
>>  1876  10   1  PPP   41.2  -8.6     0
>>  1876  10   1  GGG   40.5  -7.2   0.4
>>  1876  10   1  CMM   39.0  -7.0     0
>>  1876  10   1  EEE   38.6  -7.9     0
>>  1876  10   1  LLL   37.1  -8.6     0
>>  1876  10   1  AHH   38.7 -27.2  17.6
>>  1876  10   1  PDD   37.8 -25.7  19.3
>>  1876  10   1  FFF   32.6 -16.9   1.6
>>  1876  10   1  LLL   38.7  -9.1     0
>>  1876  10   1  CCC   36.5  -6.3     0
>>  1876  10   1  BBB   38.9  -7.0     0
>>  ...
>> i've used the fortran code from here
>> http://www.iges.org/grads/gadoc/aboutstationdata.html  addapted for my
>> .txt file
>>
>> and created the .ctl file
>> DSET   C:\DISASTER\20CRen\dez_1876\dados\estacoes \stat.dat
>> TITLE  station data
>> DTYPE  station
>> STNMAP C:\DISASTER\20CRen\dez_1876\dados\estacoes\stat.map
>> UNDEF   -999.
>> options sequential
>> TDEF    123 LINEAR 1OCT1876 1dy
>> VARS    1
>> prec   0   99   pprecipitation
>> ENDVARS
>>
>>
>>
>> when i run stnmap -i stat.ctl in command line i have de error:
>>
>> Processing time step 1
>> Invalid sation hdr found in station binary file
>> possible causes: invalid level count in hdr
>>                          descriptor file mismatch
>>                          file not station data
>>                          invalid relative time
>>                          invalid sequential format
>> levs = 256 flag 26312 time 3.5872e-43
>>
>> Could someone help me please?!!
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> --
> ??
>         ???
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20110329/09ba3065/attachment-0001.html
>
> ------------------------------
>
> Message: 3
> Date: Tue, 29 Mar 2011 10:51:34 +0900
> From: yagnesh <yagnesh at lowtem.hokudai.ac.jp>
> Subject: Re: [gradsusr] Wind barb
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID: <4D913B26.1020701 at lowtem.hokudai.ac.jp>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi,
> is it because your resolution is very high so that vectors got packed up
> and became shaded?,
> if so, you may need to use 'skip'.
>
> eg:
> d skip(u,10);v; mag(u,v)
>
> --
> yyr
>
>
> On 03/24/2011 08:50 AM, Dr.Gamal El Afandi wrote:
>> Dear All,
>> I used ARWpost to make .ctl and .dat files from WRF output.
>> I can show all parameters but when I tried to do  graphics output
>> types |vector, stream,| and |barb through these commands:|
>> set gxout vector
>> set gxout barb
>> set gxout stream
>> |display u ; v
>> display ave(u,t=1,t=10) ; ave(v,t=1,t=5)|
>> |display u ; v ; mag(u,v)
>> display u ; v ; hcurl(u,v)|
>> To draw a 1-D time series with graphics output types vector or barb,
>> set the dimension environment so that only time is varying, then :
>> |display const(u,0); u ; v |
>> ||
>> |I did not get any barbs but only shaded graphs although I can show
>> the u and v components separately or as velocity.|
>> |Please attached is the .ctl file may be something is wrong.|
>>
>> With my kind regards
>>
>> Gamal
>>
>> _____________________________________________________
>>
>> *Dr.Gamal Salah El Afandi,Senior Researcher.
>> Data to Insight Centerof Pervasive Technology Institute*
>> *Indiana University Bloomington USA*
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20110329/86e6e1da/attachment-0001.html
>
> ------------------------------
>
> Message: 4
> Date: Mon, 28 Mar 2011 19:26:36 -0700 (PDT)
> From: Zablone Owiti <zablonowiti at yahoo.com>
> Subject: [gradsusr] Plotting vector of direction Magnitude
> To: "gradsusr at gradsusr.org" <gradsusr at gradsusr.org>
> Message-ID: <351397.4121.qm at web161310.mail.bf1.yahoo.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi Users,
>
> I have a file of points (station locations in lat, lon), bearing(the
> bearings are such that North, East, South, West
> indicate January,April,July and October respectively)and magnitude (the
> magnitude here is the amplitude of rainfall annual cycle at each station).
> I would like to make a plot where each station is depicted by an arrow and
> the direction the arrow points is based on the value of bearing
> while length of the arrow is based on the value of magnitude.
>
> What I wish to achieve is a plot showing the amplitude of the annual cycle
> at each station and its phase shown by a 12-month clock with i.e Northward
> arrow indicating maximum rainfall January.
>
> Any suggestions on how to do this with Grads.
>
> ?
> ===================================================
> ZABLONE OWITI
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20110328/6cfe2de6/attachment-0001.html
>
> ------------------------------
>
> Message: 5
> Date: Tue, 29 Mar 2011 09:41:39 +0530
> From: mani kumari <manikumarimsc at gmail.com>
> Subject: [gradsusr] Request
> To: gradsusr at gradsusr.org
> Message-ID:
> 	<AANLkTi=gx9hX53zCOGgKX1gt1mEKjUOT+NfuiFpx6N1x at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Good morning,
>
>                     I have Uwnd, Vwnd in separate files,  i want to
>  represent wind vectors, If they are in single file i can represent but it
> is difficult  with separate files, so kindly provide information about the
> command how can i represent the wind vectors with separate files.
>
> --
>  thanking you
> D.Mani Kumari
> Research Scholar
> Dept.of Meteorology&Oceanography
> Andhra University
> Visakhapatnam.
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20110329/919a6e2f/attachment-0001.html
>
> ------------------------------
>
> Message: 6
> Date: Tue, 29 Mar 2011 10:12:12 +0530
> From: sudev das M P <devdas.pnr at gmail.com>
> Subject: Re: [gradsusr] Request
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> 	<AANLkTim=Gk=RQNsdHFBmSYEjmFW=CoMnnHKiyPWK-3R_ at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hai,
>
> Try the following commands,
>
> First open the two files, if it is netcdf use the following command,
>
> sdfopen uwnd.nc
> sdfopen vwnd.nc
> ..........................................
> .........................................
> If the variables are u and v,
> .................................
> ................................
> set gxout vector
> set ccolor 1
> set arrowhead -0.2
> set arrscl 0.5 0.2
> d skip(u.1,2,2);v.2
> ...................................
> ...................................
>
> Hope this will work.
>
> sudev
>
> On Tue, Mar 29, 2011 at 9:41 AM, mani kumari <manikumarimsc at gmail.com>wrote:
>
>> Good morning,
>>
>>                     I have Uwnd, Vwnd in separate files,  i want to
>>  represent wind vectors, If they are in single file i can represent but it
>> is difficult  with separate files, so kindly provide information about the
>> command how can i represent the wind vectors with separate files.
>>
>> --
>>  thanking you
>> D.Mani Kumari
>> Research Scholar
>> Dept.of Meteorology&Oceanography
>> Andhra University
>> Visakhapatnam.
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20110329/65275744/attachment.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 13, Issue 68
> ****************************************
>



More information about the gradsusr mailing list