[gradsusr] Creating data files

Alexandre Santos puffy78 at gmail.com
Wed May 19 11:55:51 EDT 2010 

Anna,

send for me the data complete radar , txt, (DBZ). With the data see the
source at fortran

2010/5/19 Huddleston, John <Huddleston at cira.colostate.edu>

> Anna
>
> 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.
>
> 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.
>
> Then write out each of the records in a do loop like
> IREC=1
> DO I=1,20
> WRITE(15,REC=IREC) ((A(J,K,I),J=1,600),K=1,600)
> IREC=IREC+1
> ENDDO
>
> Your control file becomes much simpler and you have 20 layers of data to
> plot.
>
> Think about it and if you want help, I can help you later with the coding,
> but not now..am working on another project.
>
> John Huddleston, PhD
> Cooperative Institute for Research in the Atmosphere
>
> -----Original Message-----
> From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org]
> On Behalf Of Anna Cinzia Marra
> Sent: Tuesday, May 18, 2010 3:47 AM
> To: GrADS Users Forum
> Subject: Re: [gradsusr] Creating data files
>
> Hello Alexandre,
> my radar data file has this structure, as I already told you:
>
> # Fields: z y x DBZ
> # Units: km km km dBZ
> 1 -224.625 -224.625 35
> 1 -224.625 -223.875 28
> 1 -224.625 -223.125 20
> 1 -224.625 -222.375 20
> 1 -224.625 -221.625 18
>
> They are related to a single radar scan (so only one time), but they have
> 20 z levels (from 1 km to 15.25 km every 0.75 km), 600 y points and 600 x
> points (from -224.625 km to 224.625 km every 0.75 km). Both x and y are
> obviously related to longitude and latitude.
>
> I modified radar.f in order to account for my big dataset.
> I modified your control file in the following way
>
> dset ^0dbz.bin
> title  Analysis DBZ
> undef   -32
> xdef     600  linear -224.625 0.75
> ydef     600  linear -224.625 0.75
> zdef     20   levels
>    15.25
>    14.5
>    13.75
>    13
>    12.25
>    11.5
>    10.75
>    10
>     9.25
>     8.5
>     7.75
>     7
>     6.25
>     5.5
>     4.75
>     4
>     3.25
>     2.5
>     1.75
>     1
> tdef     1  linear 22oct2005 1mn
> vars 4
> z        0 999 z(km)
> y        0 999 y(km)
> x        0 999 x(km)
> dbz     20 999 DBZ
> endvars
>
> but I cannot see anything reasonable with grads. What am I wrong?
> Thank again.
> Cinzia
>
> > I Understand. Other question can consult.
> >
> > Alexandre
> >
> >
> >
> >
> >
> > 2010/5/18 Anna Cinzia Marra <a.marra at isac.cnr.it>
> >
> >> Hello John,
> >> thanks a lot, but that's not a problem of fortran. I hhad already
> >> modified
> >> radar.f according to my complete data file. Maybe I am not modifying in
> >> the right way the control file.
> >>
> >> Cinzia
> >>
> >> > Cinzia,
> >> >
> >> > Edit the radar.f file you were sent by Alexandre and change the 4x5
> >> array
> >> > to fit the size of your big dataset.
> >> >
> >> > Change the input name to the name of your input file, Alexandre used
> >> > 'dbz.txt'.
> >> >
> >> > This is not so much a GrADS issue as it is a knowledge of programming
> >> and
> >> > building software.
> >> >
> >> > After you change the array size to fit your dataset, then create the
> >> > executable. e.g. gfortran -o dbz radar.f
> >> >
> >> > Rerun the application, e.g. ./dbz
> >> >
> >> > It will create the datafile dbz.bin and you can use the control file
> >> he
> >> > sent.
> >> >
> >> >
> >> > John Huddleston, PhD
> >> >
> >> >
> >> > -----Original Message-----
> >> > From: gradsusr-bounces at gradsusr.org [mailto:
> >> gradsusr-bounces at gradsusr.org]
> >> > On Behalf Of Anna Cinzia Marra
> >> > Sent: Monday, May 17, 2010 4:09 AM
> >> > To: GrADS Users Forum
> >> > Subject: Re: [gradsusr] Creating data files
> >> >
> >> > Hello,
> >> > unfortunately I am starting now with grads. at the moment I am not
> >> able
> >> to
> >> > display data from the bin file you sent me. I don't know why. Moreover
> >> I
> >> > am not able also to re-obtain that bin file by means of radar.f and I
> >> > cannot send you the complete data file because it is too big.
> >> > I am sorry.
> >> >
> >> > Thank you in advance
> >> >
> >> > Cinzia
> >> >
> >> >> Hello,
> >> >>
> >> >> send the file radar complete with the bin+ctl, please.
> >> >>
> >> >>
> >> >>
> >> >> 2010/5/14 Anna Cinzia Marra <a.marra at isac.cnr.it>
> >> >>
> >> >>> Hello,
> >> >>> thanks a lot. I am trying to convert a complete radar data file, but
> >> I
> >> >>> have some problems with plot.
> >> >>> I will try again.
> >> >>>
> >> >>> Cinzia
> >> >>>
> >> >>> > Hello,
> >> >>> >
> >> >>> > Send the program at fortran more files binary the ctl. With the
> >> >>> binary
> >> >>> +
> >> >>> > ctl
> >> >>> > you can see in the Grads
> >> >>> >
> >> >>> >
> >> >>> >
> >> >>> > 2010/5/13 Anna Cinzia Marra <a.marra at isac.cnr.it>
> >> >>> >
> >> >>> >>
> >> >>> >> Hello,
> >> >>> >>
> >> >>> >> does anybody have a fortran program which converts measured radar
> >> >>> >> reflectivity data to binary files read by grads?
> >> >>> >>
> >> >>> >> Here is an example of my data files:
> >> >>> >>
> >> >>> >> # Fields: z y x DBZ
> >> >>> >> # Units: km km km dBZ
> >> >>> >> 1 -224.625 -224.625 35
> >> >>> >> 1 -224.625 -223.875 28
> >> >>> >> 1 -224.625 -223.125 20
> >> >>> >> 1 -224.625 -222.375 20
> >> >>> >> 1 -224.625 -221.625 18
> >> >>> >> ......................
> >> >>> >>
> >> >>> >> Thank you very much.
> >> >>> >>
> >> >>> >> Cinzia
> >> >>> >>
> >> >>> >>
> >> >>> >> _______________________________________________
> >> >>> >> gradsusr mailing list
> >> >>> >> gradsusr at gradsusr.org
> >> >>> >> http://gradsusr.org/mailman/listinfo/gradsusr
> >> >>> >>
> >> >>> > _______________________________________________
> >> >>> > gradsusr mailing list
> >> >>> > gradsusr at gradsusr.org
> >> >>> > http://gradsusr.org/mailman/listinfo/gradsusr
> >> >>> >
> >> >>>
> >> >>>
> >> >>> _______________________________________________
> >> >>> gradsusr mailing list
> >> >>> gradsusr at gradsusr.org
> >> >>> http://gradsusr.org/mailman/listinfo/gradsusr
> >> >>>
> >> >> _______________________________________________
> >> >> gradsusr mailing list
> >> >> gradsusr at gradsusr.org
> >> >> http://gradsusr.org/mailman/listinfo/gradsusr
> >> >>
> >> >
> >> >
> >> > _______________________________________________
> >> > gradsusr mailing list
> >> > gradsusr at gradsusr.org
> >> > http://gradsusr.org/mailman/listinfo/gradsusr
> >> >
> >> > _______________________________________________
> >> > gradsusr mailing list
> >> > gradsusr at gradsusr.org
> >> > http://gradsusr.org/mailman/listinfo/gradsusr
> >> >
> >> >
> >>
> >>
> >> _______________________________________________
> >> gradsusr mailing list
> >> gradsusr at gradsusr.org
> >> http://gradsusr.org/mailman/listinfo/gradsusr
> >>
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
> _______________________________________________
> 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/20100519/1eda0c56/attachment-0003.html

More information about the gradsusr mailing list