Anna,<br><br>send for me the data complete radar , txt, (DBZ). With the data see the source at fortran<br><br><div class="gmail_quote">2010/5/19 Huddleston, John <span dir="ltr"><<a href="mailto:Huddleston@cira.colostate.edu">Huddleston@cira.colostate.edu</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;">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>
<font color="#888888"><br>
John Huddleston, PhD<br>
Cooperative Institute for Research in the Atmosphere<br>
</font><div class="im"><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>
</div><div><div></div><div class="h5">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>
</div></div></blockquote></div><br>