See the link below. See the <a href="http://drawbox.gs">drawbox.gs</a>. It'll help you.<br><br><a href="http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html">http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html</a><br>
<br>Guilherme.<br><br><div class="gmail_quote">2010/5/19 <span dir="ltr"><<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Send gradsusr mailing list submissions to<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of gradsusr digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. How to use draw function (sushant puranik)<br>
2. Re: How to use draw function (Heiner K?rnich)<br>
3. Re: How to use draw function (hersala)<br>
4. Re: How to use draw function (sushant puranik)<br>
5. Re: Creating data files (Huddleston, John)<br>
6. Re: Creating data files (Alexandre Santos)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 19 May 2010 12:35:37 +0530<br>
From: sushant puranik <<a href="mailto:sushantpuranik@gmail.com">sushantpuranik@gmail.com</a>><br>
Subject: [gradsusr] How to use draw function<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<<a href="mailto:AANLkTikDALtl3-T9v96_3PkfL8F1d8rJmdWZ7UO32cyU@mail.gmail.com">AANLkTikDALtl3-T9v96_3PkfL8F1d8rJmdWZ7UO32cyU@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi<br>
I am trying to make a rectangle over the region 65 deg to 75 deg east and<br>
equator to 10 deg north.<br>
As per instructions in the grads document first of all i made set line 1 1 1<br>
then use draw rec 0 65 10 75 command. but there is no rectangle drawn on the<br>
plot.<br>
<br>
In the draw command xlo ylo xhi yhi corresponds to grid point or waht else.<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100519/e9d2858b/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100519/e9d2858b/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 19 May 2010 09:36:18 +0200<br>
From: Heiner K?rnich <<a href="mailto:heiner@misu.su.se">heiner@misu.su.se</a>><br>
Subject: Re: [gradsusr] How to use draw function<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<<a href="mailto:AANLkTinA3fpr9Kh7M5lBQdoAlSz7GgUBuffMnTfDG4bU@mail.gmail.com">AANLkTinA3fpr9Kh7M5lBQdoAlSz7GgUBuffMnTfDG4bU@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi,<br>
<br>
it is actually a little confusing. The draw command uses the coordinates of<br>
your drawing window. Those coordinates range from 0 to 8.5 and 0 to 11 inch.<br>
The coordinates that you were using are the world coordinates. Use the<br>
function "q w2xy" to convert from world to xy coordinates, eg. for 65 deg<br>
east and 10 deg north<br>
q w2xy 65 10<br>
<br>
Regards,<br>
Heiner<br>
<br>
On Wed, May 19, 2010 at 9:05 AM, sushant puranik<br>
<<a href="mailto:sushantpuranik@gmail.com">sushantpuranik@gmail.com</a>>wrote:<br>
<br>
> Hi<br>
> I am trying to make a rectangle over the region 65 deg to 75 deg east and<br>
> equator to 10 deg north.<br>
> As per instructions in the grads document first of all i made set line 1 1<br>
> 1 then use draw rec 0 65 10 75 command. but there is no rectangle drawn on<br>
> the plot.<br>
><br>
> In the draw command xlo ylo xhi yhi corresponds to grid point or waht else.<br>
><br>
><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
><br>
<br>
<br>
--<br>
Heiner K?rnich<br>
Dept. of Meteorology Tel: +46 8 164333<br>
Stockholms University, SE-106 91 Stockholm, Sweden<br>
Email: <a href="mailto:heiner@misu.su.se">heiner@misu.su.se</a> <a href="http://www.misu.su.se/%7Eheiner/" target="_blank">www.misu.su.se/~heiner/</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100519/084da8bc/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100519/084da8bc/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Wed, 19 May 2010 08:39:48 -0300<br>
From: hersala <<a href="mailto:hersala@gmail.com">hersala@gmail.com</a>><br>
Subject: Re: [gradsusr] How to use draw function<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<<a href="mailto:AANLkTikwNR_o30B06j-yxVPpsfJZR8Lz3P5yhFB7A3Yj@mail.gmail.com">AANLkTikwNR_o30B06j-yxVPpsfJZR8Lz3P5yhFB7A3Yj@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi Sushant<br>
I send you a little example that allow to draw a rectangle using world<br>
coordinates, so you can change the size or change slightly the geographic<br>
region selected for your plot and the rectangle will be "anchored" to its<br>
geographic position. This small code shoul by used in a *.gs script.<br>
Best regards.<br>
<br>
Lic. Hern?n Sala<br>
Instituto Ant?rtico Argentino<br>
<br>
<br>
*Draw rectangle<br>
'set line 7 1 6'<br>
lon1= 135<br>
lat1= -78<br>
'q w2xy 'lon1' 'lat1<br>
say result<br>
x1 = subwrd(result,3)<br>
y1 = subwrd(result,6)<br>
lon2= 239<br>
lat2= -70<br>
'q w2xy 'lon2' 'lat2<br>
say result<br>
x2 = subwrd(result,3)<br>
y2 = subwrd(result,6)<br>
'draw rec 'x1' 'y1' 'x2' 'y2<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100519/91f73a41/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100519/91f73a41/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Wed, 19 May 2010 18:16:07 +0530<br>
From: sushant puranik <<a href="mailto:sushantpuranik@gmail.com">sushantpuranik@gmail.com</a>><br>
Subject: Re: [gradsusr] How to use draw function<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<<a href="mailto:AANLkTil730cx4h1IfibVnsrD07NWdvShE7ioLf0FhIaF@mail.gmail.com">AANLkTil730cx4h1IfibVnsrD07NWdvShE7ioLf0FhIaF@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Thank you very much. Problem is solved.<br>
<br>
<br>
On Wed, May 19, 2010 at 5:09 PM, hersala <<a href="mailto:hersala@gmail.com">hersala@gmail.com</a>> wrote:<br>
<br>
> Hi Sushant<br>
> I send you a little example that allow to draw a rectangle using world<br>
> coordinates, so you can change the size or change slightly the geographic<br>
> region selected for your plot and the rectangle will be "anchored" to its<br>
> geographic position. This small code shoul by used in a *.gs script.<br>
> Best regards.<br>
><br>
> Lic. Hern?n Sala<br>
> Instituto Ant?rtico Argentino<br>
><br>
><br>
> *Draw rectangle<br>
> 'set line 7 1 6'<br>
> lon1= 135<br>
> lat1= -78<br>
> 'q w2xy 'lon1' 'lat1<br>
> say result<br>
> x1 = subwrd(result,3)<br>
> y1 = subwrd(result,6)<br>
> lon2= 239<br>
> lat2= -70<br>
> 'q w2xy 'lon2' 'lat2<br>
> say result<br>
> x2 = subwrd(result,3)<br>
> y2 = subwrd(result,6)<br>
> 'draw rec 'x1' 'y1' 'x2' 'y2<br>
><br>
><br>
><br>
><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
><br>
<br>
<br>
--<br>
Sushant Puranik<br>
Junior Research Fellow<br>
Dept. of Atmospheric & Space Sciences,<br>
University of Pune,<br>
Pune-07,<br>
India.<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100519/929469b6/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100519/929469b6/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Wed, 19 May 2010 08:20:24 -0700<br>
From: "Huddleston, John" <<a href="mailto:Huddleston@cira.colostate.edu">Huddleston@cira.colostate.edu</a>><br>
Subject: Re: [gradsusr] Creating data files<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<<a href="mailto:F367071F539DC64FB75E7D7690B8160D0105BCA7173F@EXVMBX003-2.exch003intermedia.net">F367071F539DC64FB75E7D7690B8160D0105BCA7173F@EXVMBX003-2.exch003intermedia.net</a>><br>
<br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Anna<br>
<br>
You may want to think about the architecture as 20 600x600 grids. So, instead of opening the binary file with LRECL=4, open it with a record length of 36000.<br>
<br>
Set up a real array of A(600,600,20), initialize it to -999, then read in all the data, compute the Z level 1:20 corresponding to your 20 levels, compute the (ipos,jpos) given the input (km,km) units, and set the value equal to the dBZ.<br>
<br>
Then write out each of the records in a do loop like<br>
IREC=1<br>
DO I=1,20<br>
WRITE(15,REC=IREC) ((A(J,K,I),J=1,600),K=1,600)<br>
IREC=IREC+1<br>
ENDDO<br>
<br>
Your control file becomes much simpler and you have 20 layers of data to plot.<br>
<br>
Think about it and if you want help, I can help you later with the coding, but not now..am working on another project.<br>
<br>
John Huddleston, PhD<br>
Cooperative Institute for Research in the Atmosphere<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>] On Behalf Of Anna Cinzia Marra<br>
Sent: Tuesday, May 18, 2010 3:47 AM<br>
To: GrADS Users Forum<br>
Subject: Re: [gradsusr] Creating data files<br>
<br>
Hello Alexandre,<br>
my radar data file has this structure, as I already told you:<br>
<br>
# Fields: z y x DBZ<br>
# Units: km km km dBZ<br>
1 -224.625 -224.625 35<br>
1 -224.625 -223.875 28<br>
1 -224.625 -223.125 20<br>
1 -224.625 -222.375 20<br>
1 -224.625 -221.625 18<br>
<br>
They are related to a single radar scan (so only one time), but they have<br>
20 z levels (from 1 km to 15.25 km every 0.75 km), 600 y points and 600 x<br>
points (from -224.625 km to 224.625 km every 0.75 km). Both x and y are<br>
obviously related to longitude and latitude.<br>
<br>
I modified radar.f in order to account for my big dataset.<br>
I modified your control file in the following way<br>
<br>
dset ^0dbz.bin<br>
title Analysis DBZ<br>
undef -32<br>
xdef 600 linear -224.625 0.75<br>
ydef 600 linear -224.625 0.75<br>
zdef 20 levels<br>
15.25<br>
14.5<br>
13.75<br>
13<br>
12.25<br>
11.5<br>
10.75<br>
10<br>
9.25<br>
8.5<br>
7.75<br>
7<br>
6.25<br>
5.5<br>
4.75<br>
4<br>
3.25<br>
2.5<br>
1.75<br>
1<br>
tdef 1 linear 22oct2005 1mn<br>
vars 4<br>
z 0 999 z(km)<br>
y 0 999 y(km)<br>
x 0 999 x(km)<br>
dbz 20 999 DBZ<br>
endvars<br>
<br>
but I cannot see anything reasonable with grads. What am I wrong?<br>
Thank again.<br>
Cinzia<br>
<br>
> I Understand. Other question can consult.<br>
><br>
> Alexandre<br>
><br>
><br>
><br>
><br>
><br>
> 2010/5/18 Anna Cinzia Marra <<a href="mailto:a.marra@isac.cnr.it">a.marra@isac.cnr.it</a>><br>
><br>
>> Hello John,<br>
>> thanks a lot, but that's not a problem of fortran. I hhad already<br>
>> modified<br>
>> radar.f according to my complete data file. Maybe I am not modifying in<br>
>> the right way the control file.<br>
>><br>
>> Cinzia<br>
>><br>
>> > Cinzia,<br>
>> ><br>
>> > Edit the radar.f file you were sent by Alexandre and change the 4x5<br>
>> array<br>
>> > to fit the size of your big dataset.<br>
>> ><br>
>> > Change the input name to the name of your input file, Alexandre used<br>
>> > 'dbz.txt'.<br>
>> ><br>
>> > This is not so much a GrADS issue as it is a knowledge of programming<br>
>> and<br>
>> > building software.<br>
>> ><br>
>> > After you change the array size to fit your dataset, then create the<br>
>> > executable. e.g. gfortran -o dbz radar.f<br>
>> ><br>
>> > Rerun the application, e.g. ./dbz<br>
>> ><br>
>> > It will create the datafile dbz.bin and you can use the control file<br>
>> he<br>
>> > sent.<br>
>> ><br>
>> ><br>
>> > John Huddleston, PhD<br>
>> ><br>
>> ><br>
>> > -----Original Message-----<br>
>> > From: <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto:<br>
>> <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>]<br>
>> > On Behalf Of Anna Cinzia Marra<br>
>> > Sent: Monday, May 17, 2010 4:09 AM<br>
>> > To: GrADS Users Forum<br>
>> > Subject: Re: [gradsusr] Creating data files<br>
>> ><br>
>> > Hello,<br>
>> > unfortunately I am starting now with grads. at the moment I am not<br>
>> able<br>
>> to<br>
>> > display data from the bin file you sent me. I don't know why. Moreover<br>
>> I<br>
>> > am not able also to re-obtain that bin file by means of radar.f and I<br>
>> > cannot send you the complete data file because it is too big.<br>
>> > I am sorry.<br>
>> ><br>
>> > Thank you in advance<br>
>> ><br>
>> > Cinzia<br>
>> ><br>
>> >> Hello,<br>
>> >><br>
>> >> send the file radar complete with the bin+ctl, please.<br>
>> >><br>
>> >><br>
>> >><br>
>> >> 2010/5/14 Anna Cinzia Marra <<a href="mailto:a.marra@isac.cnr.it">a.marra@isac.cnr.it</a>><br>
>> >><br>
>> >>> Hello,<br>
>> >>> thanks a lot. I am trying to convert a complete radar data file, but<br>
>> I<br>
>> >>> have some problems with plot.<br>
>> >>> I will try again.<br>
>> >>><br>
>> >>> Cinzia<br>
>> >>><br>
>> >>> > Hello,<br>
>> >>> ><br>
>> >>> > Send the program at fortran more files binary the ctl. With the<br>
>> >>> binary<br>
>> >>> +<br>
>> >>> > ctl<br>
>> >>> > you can see in the Grads<br>
>> >>> ><br>
>> >>> ><br>
>> >>> ><br>
>> >>> > 2010/5/13 Anna Cinzia Marra <<a href="mailto:a.marra@isac.cnr.it">a.marra@isac.cnr.it</a>><br>
>> >>> ><br>
>> >>> >><br>
>> >>> >> Hello,<br>
>> >>> >><br>
>> >>> >> does anybody have a fortran program which converts measured radar<br>
>> >>> >> reflectivity data to binary files read by grads?<br>
>> >>> >><br>
>> >>> >> Here is an example of my data files:<br>
>> >>> >><br>
>> >>> >> # Fields: z y x DBZ<br>
>> >>> >> # Units: km km km dBZ<br>
>> >>> >> 1 -224.625 -224.625 35<br>
>> >>> >> 1 -224.625 -223.875 28<br>
>> >>> >> 1 -224.625 -223.125 20<br>
>> >>> >> 1 -224.625 -222.375 20<br>
>> >>> >> 1 -224.625 -221.625 18<br>
>> >>> >> ......................<br>
>> >>> >><br>
>> >>> >> Thank you very much.<br>
>> >>> >><br>
>> >>> >> Cinzia<br>
>> >>> >><br>
>> >>> >><br>
>> >>> >> _______________________________________________<br>
>> >>> >> gradsusr mailing list<br>
>> >>> >> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>> >>> >> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>> >>> >><br>
>> >>> > _______________________________________________<br>
>> >>> > gradsusr mailing list<br>
>> >>> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>> >>> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>> >>> ><br>
>> >>><br>
>> >>><br>
>> >>> _______________________________________________<br>
>> >>> gradsusr mailing list<br>
>> >>> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>> >>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>> >>><br>
>> >> _______________________________________________<br>
>> >> gradsusr mailing list<br>
>> >> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>> >> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>> >><br>
>> ><br>
>> ><br>
>> > _______________________________________________<br>
>> > gradsusr mailing list<br>
>> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>> ><br>
>> > _______________________________________________<br>
>> > gradsusr mailing list<br>
>> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>> ><br>
>> ><br>
>><br>
>><br>
>> _______________________________________________<br>
>> gradsusr mailing list<br>
>> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Wed, 19 May 2010 12:55:51 -0300<br>
From: Alexandre Santos <<a href="mailto:puffy78@gmail.com">puffy78@gmail.com</a>><br>
Subject: Re: [gradsusr] Creating data files<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<<a href="mailto:AANLkTimnAms7m_XIaYEt9xyrcQeCFxO-kj8liMUlWcSf@mail.gmail.com">AANLkTimnAms7m_XIaYEt9xyrcQeCFxO-kj8liMUlWcSf@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Anna,<br>
<br>
send for me the data complete radar , txt, (DBZ). With the data see the<br>
source at fortran<br>
<br>
2010/5/19 Huddleston, John <<a href="mailto:Huddleston@cira.colostate.edu">Huddleston@cira.colostate.edu</a>><br>
<br>
> Anna<br>
><br>
> You may want to think about the architecture as 20 600x600 grids. So,<br>
> instead of opening the binary file with LRECL=4, open it with a record<br>
> length of 36000.<br>
><br>
> Set up a real array of A(600,600,20), initialize it to -999, then read in<br>
> all the data, compute the Z level 1:20 corresponding to your 20 levels,<br>
> compute the (ipos,jpos) given the input (km,km) units, and set the value<br>
> equal to the dBZ.<br>
><br>
> Then write out each of the records in a do loop like<br>
> IREC=1<br>
> DO I=1,20<br>
> WRITE(15,REC=IREC) ((A(J,K,I),J=1,600),K=1,600)<br>
> IREC=IREC+1<br>
> ENDDO<br>
><br>
> Your control file becomes much simpler and you have 20 layers of data to<br>
> plot.<br>
><br>
> Think about it and if you want help, I can help you later with the coding,<br>
> but not now..am working on another project.<br>
><br>
> John Huddleston, PhD<br>
> Cooperative Institute for Research in the Atmosphere<br>
><br>
> -----Original Message-----<br>
> From: <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>]<br>
> On Behalf Of Anna Cinzia Marra<br>
> Sent: Tuesday, May 18, 2010 3:47 AM<br>
> To: GrADS Users Forum<br>
> Subject: Re: [gradsusr] Creating data files<br>
><br>
> Hello Alexandre,<br>
> my radar data file has this structure, as I already told you:<br>
><br>
> # Fields: z y x DBZ<br>
> # Units: km km km dBZ<br>
> 1 -224.625 -224.625 35<br>
> 1 -224.625 -223.875 28<br>
> 1 -224.625 -223.125 20<br>
> 1 -224.625 -222.375 20<br>
> 1 -224.625 -221.625 18<br>
><br>
> They are related to a single radar scan (so only one time), but they have<br>
> 20 z levels (from 1 km to 15.25 km every 0.75 km), 600 y points and 600 x<br>
> points (from -224.625 km to 224.625 km every 0.75 km). Both x and y are<br>
> obviously related to longitude and latitude.<br>
><br>
> I modified radar.f in order to account for my big dataset.<br>
> I modified your control file in the following way<br>
><br>
> dset ^0dbz.bin<br>
> title Analysis DBZ<br>
> undef -32<br>
> xdef 600 linear -224.625 0.75<br>
> ydef 600 linear -224.625 0.75<br>
> zdef 20 levels<br>
> 15.25<br>
> 14.5<br>
> 13.75<br>
> 13<br>
> 12.25<br>
> 11.5<br>
> 10.75<br>
> 10<br>
> 9.25<br>
> 8.5<br>
> 7.75<br>
> 7<br>
> 6.25<br>
> 5.5<br>
> 4.75<br>
> 4<br>
> 3.25<br>
> 2.5<br>
> 1.75<br>
> 1<br>
> tdef 1 linear 22oct2005 1mn<br>
> vars 4<br>
> z 0 999 z(km)<br>
> y 0 999 y(km)<br>
> x 0 999 x(km)<br>
> dbz 20 999 DBZ<br>
> endvars<br>
><br>
> but I cannot see anything reasonable with grads. What am I wrong?<br>
> Thank again.<br>
> Cinzia<br>
><br>
> > I Understand. Other question can consult.<br>
> ><br>
> > Alexandre<br>
> ><br>
> ><br>
> ><br>
> ><br>
> ><br>
> > 2010/5/18 Anna Cinzia Marra <<a href="mailto:a.marra@isac.cnr.it">a.marra@isac.cnr.it</a>><br>
> ><br>
> >> Hello John,<br>
> >> thanks a lot, but that's not a problem of fortran. I hhad already<br>
> >> modified<br>
> >> radar.f according to my complete data file. Maybe I am not modifying in<br>
> >> the right way the control file.<br>
> >><br>
> >> Cinzia<br>
> >><br>
> >> > Cinzia,<br>
> >> ><br>
> >> > Edit the radar.f file you were sent by Alexandre and change the 4x5<br>
> >> array<br>
> >> > to fit the size of your big dataset.<br>
> >> ><br>
> >> > Change the input name to the name of your input file, Alexandre used<br>
> >> > 'dbz.txt'.<br>
> >> ><br>
> >> > This is not so much a GrADS issue as it is a knowledge of programming<br>
> >> and<br>
> >> > building software.<br>
> >> ><br>
> >> > After you change the array size to fit your dataset, then create the<br>
> >> > executable. e.g. gfortran -o dbz radar.f<br>
> >> ><br>
> >> > Rerun the application, e.g. ./dbz<br>
> >> ><br>
> >> > It will create the datafile dbz.bin and you can use the control file<br>
> >> he<br>
> >> > sent.<br>
> >> ><br>
> >> ><br>
> >> > John Huddleston, PhD<br>
> >> ><br>
> >> ><br>
> >> > -----Original Message-----<br>
> >> > From: <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto:<br>
> >> <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>]<br>
> >> > On Behalf Of Anna Cinzia Marra<br>
> >> > Sent: Monday, May 17, 2010 4:09 AM<br>
> >> > To: GrADS Users Forum<br>
> >> > Subject: Re: [gradsusr] Creating data files<br>
> >> ><br>
> >> > Hello,<br>
> >> > unfortunately I am starting now with grads. at the moment I am not<br>
> >> able<br>
> >> to<br>
> >> > display data from the bin file you sent me. I don't know why. Moreover<br>
> >> I<br>
> >> > am not able also to re-obtain that bin file by means of radar.f and I<br>
> >> > cannot send you the complete data file because it is too big.<br>
> >> > I am sorry.<br>
> >> ><br>
> >> > Thank you in advance<br>
> >> ><br>
> >> > Cinzia<br>
> >> ><br>
> >> >> Hello,<br>
> >> >><br>
> >> >> send the file radar complete with the bin+ctl, please.<br>
> >> >><br>
> >> >><br>
> >> >><br>
> >> >> 2010/5/14 Anna Cinzia Marra <<a href="mailto:a.marra@isac.cnr.it">a.marra@isac.cnr.it</a>><br>
> >> >><br>
> >> >>> Hello,<br>
> >> >>> thanks a lot. I am trying to convert a complete radar data file, but<br>
> >> I<br>
> >> >>> have some problems with plot.<br>
> >> >>> I will try again.<br>
> >> >>><br>
> >> >>> Cinzia<br>
> >> >>><br>
> >> >>> > Hello,<br>
> >> >>> ><br>
> >> >>> > Send the program at fortran more files binary the ctl. With the<br>
> >> >>> binary<br>
> >> >>> +<br>
> >> >>> > ctl<br>
> >> >>> > you can see in the Grads<br>
> >> >>> ><br>
> >> >>> ><br>
> >> >>> ><br>
> >> >>> > 2010/5/13 Anna Cinzia Marra <<a href="mailto:a.marra@isac.cnr.it">a.marra@isac.cnr.it</a>><br>
> >> >>> ><br>
> >> >>> >><br>
> >> >>> >> Hello,<br>
> >> >>> >><br>
> >> >>> >> does anybody have a fortran program which converts measured radar<br>
> >> >>> >> reflectivity data to binary files read by grads?<br>
> >> >>> >><br>
> >> >>> >> Here is an example of my data files:<br>
> >> >>> >><br>
> >> >>> >> # Fields: z y x DBZ<br>
> >> >>> >> # Units: km km km dBZ<br>
> >> >>> >> 1 -224.625 -224.625 35<br>
> >> >>> >> 1 -224.625 -223.875 28<br>
> >> >>> >> 1 -224.625 -223.125 20<br>
> >> >>> >> 1 -224.625 -222.375 20<br>
> >> >>> >> 1 -224.625 -221.625 18<br>
> >> >>> >> ......................<br>
> >> >>> >><br>
> >> >>> >> Thank you very much.<br>
> >> >>> >><br>
> >> >>> >> Cinzia<br>
> >> >>> >><br>
> >> >>> >><br>
> >> >>> >> _______________________________________________<br>
> >> >>> >> gradsusr mailing list<br>
> >> >>> >> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> >> >>> >> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> >> >>> >><br>
> >> >>> > _______________________________________________<br>
> >> >>> > gradsusr mailing list<br>
> >> >>> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> >> >>> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> >> >>> ><br>
> >> >>><br>
> >> >>><br>
> >> >>> _______________________________________________<br>
> >> >>> gradsusr mailing list<br>
> >> >>> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> >> >>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> >> >>><br>
> >> >> _______________________________________________<br>
> >> >> gradsusr mailing list<br>
> >> >> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> >> >> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> >> >><br>
> >> ><br>
> >> ><br>
> >> > _______________________________________________<br>
> >> > gradsusr mailing list<br>
> >> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> >> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> >> ><br>
> >> > _______________________________________________<br>
> >> > gradsusr mailing list<br>
> >> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> >> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> >> ><br>
> >> ><br>
> >><br>
> >><br>
> >> _______________________________________________<br>
> >> gradsusr mailing list<br>
> >> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> >> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> >><br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> ><br>
><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20100519/1eda0c56/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20100519/1eda0c56/attachment.html</a><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
End of gradsusr Digest, Vol 3, Issue 41<br>
***************************************<br>
</blockquote></div><br><br clear="all"><br>-- <br>-----------------------------------------------------------<br>Guilherme Martins<br>Instituto Nacional de Pesquisas Espaciais (INPE)<br>Centro de Previsão de Tempo e Estudos Climáticos (CPTEC)<br>
Grupo de Modelagem da Atmosfera e Interfaces (GMAI)<br>Telefone(CP): +55 12 3186-8502<br>Celular(SJC): +55 12 8152-8580 / +55 12 3208-7276<br>E-mail: <a href="mailto:guilherme.martins@inpe.br">guilherme.martins@inpe.br</a><br>
E-mail alternativo: <a href="mailto:jgmsantos@gmail.com">jgmsantos@gmail.com</a><br>Homepage: <a href="http://guilhermemartins.hd1.com.br/">http://guilhermemartins.hd1.com.br/</a><br>-----------------------------------------------------------<br>
<div style="visibility: hidden; display: inline;" id="avg_ls_inline_popup"></div><style type="text/css">#avg_ls_inline_popup { position:absolute; z-index:9999; padding: 0px 0px; margin-left: 0px; margin-top: 0px; width: 240px; overflow: hidden; word-wrap: break-word; color: black; font-size: 10px; text-align: left; line-height: 13px;}</style>