[Met_help] [rt.rap.ucar.edu #49784] History for Error in converting TRMM 3-hr rainfall into MET netcdf format
John Halley Gotway via RT
met_help at ucar.edu
Wed Oct 19 09:40:16 MDT 2011
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear Methelp,
I tried to convert the TRMM 3-hrly accumulated data (from
http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml) using
trmm2nc.R file.
I copied the ASCII rainfall data and attached for your reference.
When i am running, it is giving error that
*Error: unexpected symbol in "Rscript trmm2nc.R"*
Can you help me why it is like that.
Is this with NETCDF library?
For your information, I have netcdf library and already export that path in
my .bashrc file.
Will the Rscript automatically takes or should I mentioned somewhere when
using this Rscript.
Can you suggest what went wrong?
Thanks in advance,
Krishna.
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #49784] Error in converting TRMM 3-hr rainfall into MET netcdf format
From: John Halley Gotway
Time: Mon Sep 26 12:20:26 2011
Krishna,
I used the trmm2nc.R script to convert the sample data file you sent,
and everything worked fine.
So there is likely a problem in the environment somewhere. I have
never seen that sort of error message from R before. Let's try
looking how R is installed on your system:
Check what version of R you're running:
R --version
Check that you have the Rscript command in your path:
which Rscript
Try running the Rscript command:
Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt TRMM_rain_20070420-
9_12UTC.nc
Please note that this script expects that the R "ncdf" package is
installed on your system. To check if it's there, just try executing
the following commands:
R
library(ncdf)
If it loads without error, everything is fine. If you get an error,
you'll need to install that "ncdf" package. To do so, execute the
following commands:
R
install.packages("ncdf", configure.args="--with-
netcdf_include=/usr/local/netcdf/include
-with-netcdf_lib=/usr/local/netcdf/lib")
Please let me know what you learn from these steps.
Thanks,
John Halley Gotway
met_help at ucar.edu
On 09/23/2011 08:01 PM, RAL HelpDesk {for krishna kishore} wrote:
>
> Fri Sep 23 20:01:29 2011: Request 49784 was acted upon.
> Transaction: Ticket created by osurikishore at gmail.com
> Queue: met_help
> Subject: Error in converting TRMM 3-hr rainfall into MET netcdf
format
> Owner: Nobody
> Requestors: osurikishore at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>
>
> Dear Methelp,
> I tried to convert the TRMM 3-hrly accumulated data (from
> http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml)
using
> trmm2nc.R file.
> I copied the ASCII rainfall data and attached for your reference.
> When i am running, it is giving error that
> *Error: unexpected symbol in "Rscript trmm2nc.R"*
>
> Can you help me why it is like that.
> Is this with NETCDF library?
> For your information, I have netcdf library and already export that
path in
> my .bashrc file.
> Will the Rscript automatically takes or should I mentioned somewhere
when
> using this Rscript.
>
> Can you suggest what went wrong?
> Thanks in advance,
> Krishna.
>
>
------------------------------------------------
Subject: Error in converting TRMM 3-hr rainfall into MET netcdf format
From: krishna kishore
Time: Tue Oct 04 16:23:37 2011
Hi John,
I Installed the ncdf library in R.
For your more information....I Installed R software in my ubuntu
linux
machine (through sudo apt-get install r-base) ....and installed the
ncdf
library also.
Then I tried with the commands
R for R prompt and then try to run the trmm2nc.R script
> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
TRMM_rain_20070420-9_12UTC.nc
Error: unexpected symbol in "Rscript trmm2nc.R"
> R --version
Error: object 'R' not found
>which Rscript
Error: unexpected symbol in "which Rscript"
>library(ncdf)
nothing is displaying....means installed the library correctly..
But why the above errors are giving when i am executin Rscript and R
--version commands.
Thats means, R software is not installed in my account properly?
What I should do to pass through this problem..
Thanks,
Krishna.
On Mon, Sep 26, 2011 at 11:50 PM, RAL HelpDesk {for John Halley
Gotway} <
met_help at ucar.edu> wrote:
> Krishna,
>
> I used the trmm2nc.R script to convert the sample data file you
sent, and
> everything worked fine.
>
> So there is likely a problem in the environment somewhere. I have
never
> seen that sort of error message from R before. Let's try looking
how R is
> installed on your system:
>
> Check what version of R you're running:
> R --version
> Check that you have the Rscript command in your path:
> which Rscript
> Try running the Rscript command:
> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
> TRMM_rain_20070420-9_12UTC.nc
>
> Please note that this script expects that the R "ncdf" package is
installed
> on your system. To check if it's there, just try executing the
following
> commands:
> R
> library(ncdf)
>
> If it loads without error, everything is fine. If you get an error,
you'll
> need to install that "ncdf" package. To do so, execute the
following
> commands:
> R
> install.packages("ncdf",
> configure.args="--with-netcdf_include=/usr/local/netcdf/include
> -with-netcdf_lib=/usr/local/netcdf/lib")
>
> Please let me know what you learn from these steps.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
> On 09/23/2011 08:01 PM, RAL HelpDesk {for krishna kishore} wrote:
> >
> > Fri Sep 23 20:01:29 2011: Request 49784 was acted upon.
> > Transaction: Ticket created by osurikishore at gmail.com
> > Queue: met_help
> > Subject: Error in converting TRMM 3-hr rainfall into MET
netcdf
> format
> > Owner: Nobody
> > Requestors: osurikishore at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
> >
> >
> > Dear Methelp,
> > I tried to convert the TRMM 3-hrly accumulated data (from
> > http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml)
using
> > trmm2nc.R file.
> > I copied the ASCII rainfall data and attached for your reference.
> > When i am running, it is giving error that
> > *Error: unexpected symbol in "Rscript trmm2nc.R"*
> >
> > Can you help me why it is like that.
> > Is this with NETCDF library?
> > For your information, I have netcdf library and already export
that path
> in
> > my .bashrc file.
> > Will the Rscript automatically takes or should I mentioned
somewhere when
> > using this Rscript.
> >
> > Can you suggest what went wrong?
> > Thanks in advance,
> > Krishna.
> >
> >
>
>
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #49784] Error in converting TRMM 3-hr rainfall into MET netcdf format
From: John Halley Gotway
Time: Wed Oct 05 09:20:48 2011
Krishna,
Ah, I see the confusion. Those commands I sent should be run from the
command line - NOT FROM WITHIN R. Please try running the commands
again, but do NOT call R to start up an R session first.
Hopefully, that'll do it.
Thanks,
John
On 10/04/2011 04:23 PM, krishna kishore via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>
> Hi John,
> I Installed the ncdf library in R.
> For your more information....I Installed R software in my ubuntu
linux
> machine (through sudo apt-get install r-base) ....and installed the
ncdf
> library also.
> Then I tried with the commands
> R for R prompt and then try to run the trmm2nc.R script
>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
> TRMM_rain_20070420-9_12UTC.nc
> Error: unexpected symbol in "Rscript trmm2nc.R"
>
>
>> R --version
> Error: object 'R' not found
>
>> which Rscript
> Error: unexpected symbol in "which Rscript"
>
>> library(ncdf)
> nothing is displaying....means installed the library correctly..
>
> But why the above errors are giving when i am executin Rscript and R
> --version commands.
>
> Thats means, R software is not installed in my account properly?
> What I should do to pass through this problem..
>
> Thanks,
> Krishna.
>
>
> On Mon, Sep 26, 2011 at 11:50 PM, RAL HelpDesk {for John Halley
Gotway} <
> met_help at ucar.edu> wrote:
>
>> Krishna,
>>
>> I used the trmm2nc.R script to convert the sample data file you
sent, and
>> everything worked fine.
>>
>> So there is likely a problem in the environment somewhere. I have
never
>> seen that sort of error message from R before. Let's try looking
how R is
>> installed on your system:
>>
>> Check what version of R you're running:
>> R --version
>> Check that you have the Rscript command in your path:
>> which Rscript
>> Try running the Rscript command:
>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
>> TRMM_rain_20070420-9_12UTC.nc
>>
>> Please note that this script expects that the R "ncdf" package is
installed
>> on your system. To check if it's there, just try executing the
following
>> commands:
>> R
>> library(ncdf)
>>
>> If it loads without error, everything is fine. If you get an
error, you'll
>> need to install that "ncdf" package. To do so, execute the
following
>> commands:
>> R
>> install.packages("ncdf",
>> configure.args="--with-netcdf_include=/usr/local/netcdf/include
>> -with-netcdf_lib=/usr/local/netcdf/lib")
>>
>> Please let me know what you learn from these steps.
>>
>> Thanks,
>> John Halley Gotway
>> met_help at ucar.edu
>>
>> On 09/23/2011 08:01 PM, RAL HelpDesk {for krishna kishore} wrote:
>>>
>>> Fri Sep 23 20:01:29 2011: Request 49784 was acted upon.
>>> Transaction: Ticket created by osurikishore at gmail.com
>>> Queue: met_help
>>> Subject: Error in converting TRMM 3-hr rainfall into MET
netcdf
>> format
>>> Owner: Nobody
>>> Requestors: osurikishore at gmail.com
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>>>
>>>
>>> Dear Methelp,
>>> I tried to convert the TRMM 3-hrly accumulated data (from
>>> http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml)
using
>>> trmm2nc.R file.
>>> I copied the ASCII rainfall data and attached for your reference.
>>> When i am running, it is giving error that
>>> *Error: unexpected symbol in "Rscript trmm2nc.R"*
>>>
>>> Can you help me why it is like that.
>>> Is this with NETCDF library?
>>> For your information, I have netcdf library and already export
that path
>> in
>>> my .bashrc file.
>>> Will the Rscript automatically takes or should I mentioned
somewhere when
>>> using this Rscript.
>>>
>>> Can you suggest what went wrong?
>>> Thanks in advance,
>>> Krishna.
>>>
>>>
>>
>>
>
>
------------------------------------------------
Subject: Error in converting TRMM 3-hr rainfall into MET netcdf format
From: krishna kishore
Time: Wed Oct 05 09:48:48 2011
Hi John,
Thank you very much..
Now I tried to convert the TRMM data....Still giving error..
**************************************************************************
/home/osuri/Desktop/TRMM2NC_for_MET>Rscript trmm2nc.R
TRMM_rain_20070420-12_15UTC.txt TRMM_rain_20070420-12_15UTC.nc
Reading: TRMM_rain_20070420-12_15UTC.txt
Error in dim.def.ncdf("lat", "degrees_north", lat, create_dimvar =
TRUE, :
unused argument(s) (longname = "latitude")
Execution halted
**************************************************************************
Can you check what is the problem?
Thanks,
Krishna.
On Wed, Oct 5, 2011 at 8:50 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Krishna,
>
> Ah, I see the confusion. Those commands I sent should be run from
the
> command line - NOT FROM WITHIN R. Please try running the commands
again,
> but do NOT call R to start up an R session first.
>
> Hopefully, that'll do it.
>
> Thanks,
> John
>
> On 10/04/2011 04:23 PM, krishna kishore via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
> >
> > Hi John,
> > I Installed the ncdf library in R.
> > For your more information....I Installed R software in my ubuntu
linux
> > machine (through sudo apt-get install r-base) ....and installed
the ncdf
> > library also.
> > Then I tried with the commands
> > R for R prompt and then try to run the trmm2nc.R script
> >> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
> > TRMM_rain_20070420-9_12UTC.nc
> > Error: unexpected symbol in "Rscript trmm2nc.R"
> >
> >
> >> R --version
> > Error: object 'R' not found
> >
> >> which Rscript
> > Error: unexpected symbol in "which Rscript"
> >
> >> library(ncdf)
> > nothing is displaying....means installed the library correctly..
> >
> > But why the above errors are giving when i am executin Rscript and
R
> > --version commands.
> >
> > Thats means, R software is not installed in my account properly?
> > What I should do to pass through this problem..
> >
> > Thanks,
> > Krishna.
> >
> >
> > On Mon, Sep 26, 2011 at 11:50 PM, RAL HelpDesk {for John Halley
Gotway} <
> > met_help at ucar.edu> wrote:
> >
> >> Krishna,
> >>
> >> I used the trmm2nc.R script to convert the sample data file you
sent,
> and
> >> everything worked fine.
> >>
> >> So there is likely a problem in the environment somewhere. I
have never
> >> seen that sort of error message from R before. Let's try looking
how R
> is
> >> installed on your system:
> >>
> >> Check what version of R you're running:
> >> R --version
> >> Check that you have the Rscript command in your path:
> >> which Rscript
> >> Try running the Rscript command:
> >> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
> >> TRMM_rain_20070420-9_12UTC.nc
> >>
> >> Please note that this script expects that the R "ncdf" package is
> installed
> >> on your system. To check if it's there, just try executing the
> following
> >> commands:
> >> R
> >> library(ncdf)
> >>
> >> If it loads without error, everything is fine. If you get an
error,
> you'll
> >> need to install that "ncdf" package. To do so, execute the
following
> >> commands:
> >> R
> >> install.packages("ncdf",
> >> configure.args="--with-netcdf_include=/usr/local/netcdf/include
> >> -with-netcdf_lib=/usr/local/netcdf/lib")
> >>
> >> Please let me know what you learn from these steps.
> >>
> >> Thanks,
> >> John Halley Gotway
> >> met_help at ucar.edu
> >>
> >> On 09/23/2011 08:01 PM, RAL HelpDesk {for krishna kishore} wrote:
> >>>
> >>> Fri Sep 23 20:01:29 2011: Request 49784 was acted upon.
> >>> Transaction: Ticket created by osurikishore at gmail.com
> >>> Queue: met_help
> >>> Subject: Error in converting TRMM 3-hr rainfall into MET
netcdf
> >> format
> >>> Owner: Nobody
> >>> Requestors: osurikishore at gmail.com
> >>> Status: new
> >>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784>
> >>>
> >>>
> >>> Dear Methelp,
> >>> I tried to convert the TRMM 3-hrly accumulated data (from
> >>> http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml)
using
> >>> trmm2nc.R file.
> >>> I copied the ASCII rainfall data and attached for your
reference.
> >>> When i am running, it is giving error that
> >>> *Error: unexpected symbol in "Rscript trmm2nc.R"*
> >>>
> >>> Can you help me why it is like that.
> >>> Is this with NETCDF library?
> >>> For your information, I have netcdf library and already export
that
> path
> >> in
> >>> my .bashrc file.
> >>> Will the Rscript automatically takes or should I mentioned
somewhere
> when
> >>> using this Rscript.
> >>>
> >>> Can you suggest what went wrong?
> >>> Thanks in advance,
> >>> Krishna.
> >>>
> >>>
> >>
> >>
> >
> >
>
>
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #49784] Error in converting TRMM 3-hr rainfall into MET netcdf format
From: John Halley Gotway
Time: Wed Oct 05 12:10:02 2011
Krishna,
It would probably be easiest for you just to send me the file that's
giving the error. Can you please send me this file:
TRMM_rain_20070420-12_15UTC.txt
I'll run it here and try to diagnose what's going on. You can either
send it via email or post it to our anonymous ftp site, as described
here:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
If you do post it to our ftp site, please write me back to tell me
when it's there.
Hopefully we can get this working for you.
Thanks,
John
On 10/05/2011 09:48 AM, krishna kishore via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>
> Hi John,
> Thank you very much..
> Now I tried to convert the TRMM data....Still giving error..
>
**************************************************************************
> /home/osuri/Desktop/TRMM2NC_for_MET>Rscript trmm2nc.R
> TRMM_rain_20070420-12_15UTC.txt TRMM_rain_20070420-12_15UTC.nc
> Reading: TRMM_rain_20070420-12_15UTC.txt
> Error in dim.def.ncdf("lat", "degrees_north", lat, create_dimvar =
TRUE, :
> unused argument(s) (longname = "latitude")
> Execution halted
>
**************************************************************************
> Can you check what is the problem?
> Thanks,
> Krishna.
>
>
> On Wed, Oct 5, 2011 at 8:50 PM, John Halley Gotway via RT
<met_help at ucar.edu
>> wrote:
>
>> Krishna,
>>
>> Ah, I see the confusion. Those commands I sent should be run from
the
>> command line - NOT FROM WITHIN R. Please try running the commands
again,
>> but do NOT call R to start up an R session first.
>>
>> Hopefully, that'll do it.
>>
>> Thanks,
>> John
>>
>> On 10/04/2011 04:23 PM, krishna kishore via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>>>
>>> Hi John,
>>> I Installed the ncdf library in R.
>>> For your more information....I Installed R software in my ubuntu
linux
>>> machine (through sudo apt-get install r-base) ....and installed
the ncdf
>>> library also.
>>> Then I tried with the commands
>>> R for R prompt and then try to run the trmm2nc.R script
>>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
>>> TRMM_rain_20070420-9_12UTC.nc
>>> Error: unexpected symbol in "Rscript trmm2nc.R"
>>>
>>>
>>>> R --version
>>> Error: object 'R' not found
>>>
>>>> which Rscript
>>> Error: unexpected symbol in "which Rscript"
>>>
>>>> library(ncdf)
>>> nothing is displaying....means installed the library correctly..
>>>
>>> But why the above errors are giving when i am executin Rscript and
R
>>> --version commands.
>>>
>>> Thats means, R software is not installed in my account properly?
>>> What I should do to pass through this problem..
>>>
>>> Thanks,
>>> Krishna.
>>>
>>>
>>> On Mon, Sep 26, 2011 at 11:50 PM, RAL HelpDesk {for John Halley
Gotway} <
>>> met_help at ucar.edu> wrote:
>>>
>>>> Krishna,
>>>>
>>>> I used the trmm2nc.R script to convert the sample data file you
sent,
>> and
>>>> everything worked fine.
>>>>
>>>> So there is likely a problem in the environment somewhere. I
have never
>>>> seen that sort of error message from R before. Let's try looking
how R
>> is
>>>> installed on your system:
>>>>
>>>> Check what version of R you're running:
>>>> R --version
>>>> Check that you have the Rscript command in your path:
>>>> which Rscript
>>>> Try running the Rscript command:
>>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
>>>> TRMM_rain_20070420-9_12UTC.nc
>>>>
>>>> Please note that this script expects that the R "ncdf" package is
>> installed
>>>> on your system. To check if it's there, just try executing the
>> following
>>>> commands:
>>>> R
>>>> library(ncdf)
>>>>
>>>> If it loads without error, everything is fine. If you get an
error,
>> you'll
>>>> need to install that "ncdf" package. To do so, execute the
following
>>>> commands:
>>>> R
>>>> install.packages("ncdf",
>>>> configure.args="--with-netcdf_include=/usr/local/netcdf/include
>>>> -with-netcdf_lib=/usr/local/netcdf/lib")
>>>>
>>>> Please let me know what you learn from these steps.
>>>>
>>>> Thanks,
>>>> John Halley Gotway
>>>> met_help at ucar.edu
>>>>
>>>> On 09/23/2011 08:01 PM, RAL HelpDesk {for krishna kishore} wrote:
>>>>>
>>>>> Fri Sep 23 20:01:29 2011: Request 49784 was acted upon.
>>>>> Transaction: Ticket created by osurikishore at gmail.com
>>>>> Queue: met_help
>>>>> Subject: Error in converting TRMM 3-hr rainfall into MET
netcdf
>>>> format
>>>>> Owner: Nobody
>>>>> Requestors: osurikishore at gmail.com
>>>>> Status: new
>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784>
>>>>>
>>>>>
>>>>> Dear Methelp,
>>>>> I tried to convert the TRMM 3-hrly accumulated data (from
>>>>> http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml)
using
>>>>> trmm2nc.R file.
>>>>> I copied the ASCII rainfall data and attached for your
reference.
>>>>> When i am running, it is giving error that
>>>>> *Error: unexpected symbol in "Rscript trmm2nc.R"*
>>>>>
>>>>> Can you help me why it is like that.
>>>>> Is this with NETCDF library?
>>>>> For your information, I have netcdf library and already export
that
>> path
>>>> in
>>>>> my .bashrc file.
>>>>> Will the Rscript automatically takes or should I mentioned
somewhere
>> when
>>>>> using this Rscript.
>>>>>
>>>>> Can you suggest what went wrong?
>>>>> Thanks in advance,
>>>>> Krishna.
>>>>>
>>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>
>
------------------------------------------------
Subject: Error in converting TRMM 3-hr rainfall into MET netcdf format
From: krishna kishore
Time: Wed Oct 05 12:13:20 2011
Thank you,
I am attaching the same file with this mail.
Please explain what is the error...
Krishna.
On Wed, Oct 5, 2011 at 11:40 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Krishna,
>
> It would probably be easiest for you just to send me the file that's
giving
> the error. Can you please send me this file:
> TRMM_rain_20070420-12_15UTC.txt
>
> I'll run it here and try to diagnose what's going on. You can
either send
> it via email or post it to our anonymous ftp site, as described
here:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
> If you do post it to our ftp site, please write me back to tell me
when
> it's there.
>
> Hopefully we can get this working for you.
>
> Thanks,
> John
>
>
>
> On 10/05/2011 09:48 AM, krishna kishore via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
> >
> > Hi John,
> > Thank you very much..
> > Now I tried to convert the TRMM data....Still giving error..
> >
>
**************************************************************************
> > /home/osuri/Desktop/TRMM2NC_for_MET>Rscript trmm2nc.R
> > TRMM_rain_20070420-12_15UTC.txt TRMM_rain_20070420-12_15UTC.nc
> > Reading: TRMM_rain_20070420-12_15UTC.txt
> > Error in dim.def.ncdf("lat", "degrees_north", lat, create_dimvar =
TRUE,
> :
> > unused argument(s) (longname = "latitude")
> > Execution halted
> >
>
**************************************************************************
> > Can you check what is the problem?
> > Thanks,
> > Krishna.
> >
> >
> > On Wed, Oct 5, 2011 at 8:50 PM, John Halley Gotway via RT <
> met_help at ucar.edu
> >> wrote:
> >
> >> Krishna,
> >>
> >> Ah, I see the confusion. Those commands I sent should be run
from the
> >> command line - NOT FROM WITHIN R. Please try running the
commands
> again,
> >> but do NOT call R to start up an R session first.
> >>
> >> Hopefully, that'll do it.
> >>
> >> Thanks,
> >> John
> >>
> >> On 10/04/2011 04:23 PM, krishna kishore via RT wrote:
> >>>
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
> >>>
> >>> Hi John,
> >>> I Installed the ncdf library in R.
> >>> For your more information....I Installed R software in my
ubuntu linux
> >>> machine (through sudo apt-get install r-base) ....and installed
the
> ncdf
> >>> library also.
> >>> Then I tried with the commands
> >>> R for R prompt and then try to run the trmm2nc.R script
> >>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
> >>> TRMM_rain_20070420-9_12UTC.nc
> >>> Error: unexpected symbol in "Rscript trmm2nc.R"
> >>>
> >>>
> >>>> R --version
> >>> Error: object 'R' not found
> >>>
> >>>> which Rscript
> >>> Error: unexpected symbol in "which Rscript"
> >>>
> >>>> library(ncdf)
> >>> nothing is displaying....means installed the library correctly..
> >>>
> >>> But why the above errors are giving when i am executin Rscript
and R
> >>> --version commands.
> >>>
> >>> Thats means, R software is not installed in my account properly?
> >>> What I should do to pass through this problem..
> >>>
> >>> Thanks,
> >>> Krishna.
> >>>
> >>>
> >>> On Mon, Sep 26, 2011 at 11:50 PM, RAL HelpDesk {for John Halley
Gotway}
> <
> >>> met_help at ucar.edu> wrote:
> >>>
> >>>> Krishna,
> >>>>
> >>>> I used the trmm2nc.R script to convert the sample data file you
sent,
> >> and
> >>>> everything worked fine.
> >>>>
> >>>> So there is likely a problem in the environment somewhere. I
have
> never
> >>>> seen that sort of error message from R before. Let's try
looking how
> R
> >> is
> >>>> installed on your system:
> >>>>
> >>>> Check what version of R you're running:
> >>>> R --version
> >>>> Check that you have the Rscript command in your path:
> >>>> which Rscript
> >>>> Try running the Rscript command:
> >>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
> >>>> TRMM_rain_20070420-9_12UTC.nc
> >>>>
> >>>> Please note that this script expects that the R "ncdf" package
is
> >> installed
> >>>> on your system. To check if it's there, just try executing the
> >> following
> >>>> commands:
> >>>> R
> >>>> library(ncdf)
> >>>>
> >>>> If it loads without error, everything is fine. If you get an
error,
> >> you'll
> >>>> need to install that "ncdf" package. To do so, execute the
following
> >>>> commands:
> >>>> R
> >>>> install.packages("ncdf",
> >>>> configure.args="--with-netcdf_include=/usr/local/netcdf/include
> >>>> -with-netcdf_lib=/usr/local/netcdf/lib")
> >>>>
> >>>> Please let me know what you learn from these steps.
> >>>>
> >>>> Thanks,
> >>>> John Halley Gotway
> >>>> met_help at ucar.edu
> >>>>
> >>>> On 09/23/2011 08:01 PM, RAL HelpDesk {for krishna kishore}
wrote:
> >>>>>
> >>>>> Fri Sep 23 20:01:29 2011: Request 49784 was acted upon.
> >>>>> Transaction: Ticket created by osurikishore at gmail.com
> >>>>> Queue: met_help
> >>>>> Subject: Error in converting TRMM 3-hr rainfall into MET
netcdf
> >>>> format
> >>>>> Owner: Nobody
> >>>>> Requestors: osurikishore at gmail.com
> >>>>> Status: new
> >>>>> Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784>
> >>>>>
> >>>>>
> >>>>> Dear Methelp,
> >>>>> I tried to convert the TRMM 3-hrly accumulated data (from
> >>>>>
http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml)
> using
> >>>>> trmm2nc.R file.
> >>>>> I copied the ASCII rainfall data and attached for your
reference.
> >>>>> When i am running, it is giving error that
> >>>>> *Error: unexpected symbol in "Rscript trmm2nc.R"*
> >>>>>
> >>>>> Can you help me why it is like that.
> >>>>> Is this with NETCDF library?
> >>>>> For your information, I have netcdf library and already export
that
> >> path
> >>>> in
> >>>>> my .bashrc file.
> >>>>> Will the Rscript automatically takes or should I mentioned
somewhere
> >> when
> >>>>> using this Rscript.
> >>>>>
> >>>>> Can you suggest what went wrong?
> >>>>> Thanks in advance,
> >>>>> Krishna.
> >>>>>
> >>>>>
> >>>>
> >>>>
> >>>
> >>>
> >>
> >>
> >
> >
>
>
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
------------------------------------------------
Subject: Error in converting TRMM 3-hr rainfall into MET netcdf format
From: krishna kishore
Time: Wed Oct 05 12:13:20 2011
Selected parameter: 3-hourly TRMM 3B42(V6) Accumulated Rainfall
Selected area: lat=[18N,28N], lon=[81E,93E]
Selected time period: (12Z20Apr2007-15Z20Apr2007)
Undefined/Missing Value: -9999
Unit: (mm)
Latitude Longitude AccRain
18.000 81.000 0.0000
18.000 81.250 0.0000
18.000 81.500 0.4632
18.000 81.750 2.5835
18.000 82.000 5.1967
18.000 82.250 4.7258
18.000 82.500 4.5513
18.000 82.750 9.3033
18.000 83.000 9.0213
18.000 83.250 0.0000
18.000 83.500 0.0000
18.000 83.750 0.0000
18.000 84.000 0.0000
18.000 84.250 0.0000
18.000 84.500 0.0000
18.000 84.750 0.0000
18.000 85.000 0.0000
18.000 85.250 0.0000
18.000 85.500 0.0000
18.000 85.750 0.0000
18.000 86.000 0.0000
18.000 86.250 0.0000
18.000 86.500 0.0000
18.000 86.750 0.0000
18.000 87.000 0.0000
18.000 87.250 0.0000
18.000 87.500 0.0000
18.000 87.750 0.0000
18.000 88.000 0.0000
18.000 88.250 0.0000
18.000 88.500 0.0000
18.000 88.750 0.0000
18.000 89.000 0.0000
18.000 89.250 0.0000
18.000 89.500 0.0000
18.000 89.750 0.0000
18.000 90.000 0.0000
18.000 90.250 0.0000
18.000 90.500 0.0000
18.000 90.750 0.0000
18.000 91.000 0.0000
18.000 91.250 0.0000
18.000 91.500 0.0000
18.000 91.750 0.0000
18.000 92.000 0.0000
18.000 92.250 0.0000
18.000 92.500 0.0000
18.000 92.750 0.0000
18.000 93.000 0.0000
18.250 81.000 0.0000
18.250 81.250 0.0000
18.250 81.500 0.0000
18.250 81.750 0.0000
18.250 82.000 0.0555
18.250 82.250 0.0914
18.250 82.500 0.2913
18.250 82.750 0.8561
18.250 83.000 0.0000
18.250 83.250 0.0000
18.250 83.500 0.0000
18.250 83.750 0.0000
18.250 84.000 0.0000
18.250 84.250 0.0000
18.250 84.500 0.0000
18.250 84.750 0.0000
18.250 85.000 0.0000
18.250 85.250 0.0000
18.250 85.500 0.0000
18.250 85.750 0.0000
18.250 86.000 0.0000
18.250 86.250 0.0000
18.250 86.500 0.0000
18.250 86.750 0.0000
18.250 87.000 0.0000
18.250 87.250 0.0000
18.250 87.500 0.0000
18.250 87.750 0.0000
18.250 88.000 0.0000
18.250 88.250 0.0000
18.250 88.500 0.0000
18.250 88.750 0.0000
18.250 89.000 0.0000
18.250 89.250 0.0000
18.250 89.500 0.0000
18.250 89.750 0.0000
18.250 90.000 0.0000
18.250 90.250 0.0000
18.250 90.500 0.0000
18.250 90.750 0.0000
18.250 91.000 0.0000
18.250 91.250 0.0000
18.250 91.500 0.0000
18.250 91.750 0.0000
18.250 92.000 0.0000
18.250 92.250 0.0000
18.250 92.500 0.0000
18.250 92.750 0.0000
18.250 93.000 0.0000
18.500 81.000 0.0000
18.500 81.250 0.0000
18.500 81.500 0.0000
18.500 81.750 0.0000
18.500 82.000 0.0000
18.500 82.250 0.0000
18.500 82.500 0.0000
18.500 82.750 0.0000
18.500 83.000 0.0000
18.500 83.250 0.0000
18.500 83.500 0.0000
18.500 83.750 0.0000
18.500 84.000 0.0000
18.500 84.250 0.0000
18.500 84.500 0.0000
18.500 84.750 0.0000
18.500 85.000 0.0000
18.500 85.250 0.0000
18.500 85.500 0.0000
18.500 85.750 0.0000
18.500 86.000 0.0000
18.500 86.250 0.0000
18.500 86.500 0.0000
18.500 86.750 0.0000
18.500 87.000 0.0000
18.500 87.250 0.0000
18.500 87.500 0.0000
18.500 87.750 0.0000
18.500 88.000 0.0000
18.500 88.250 0.0000
18.500 88.500 0.0000
18.500 88.750 0.0000
18.500 89.000 0.0000
18.500 89.250 0.0000
18.500 89.500 0.0000
18.500 89.750 0.0000
18.500 90.000 0.0000
18.500 90.250 0.0000
18.500 90.500 0.0000
18.500 90.750 0.0000
18.500 91.000 0.0000
18.500 91.250 0.0000
18.500 91.500 0.0000
18.500 91.750 0.0000
18.500 92.000 0.0000
18.500 92.250 0.0000
18.500 92.500 0.0000
18.500 92.750 0.0000
18.500 93.000 0.0000
18.750 81.000 0.0000
18.750 81.250 0.0000
18.750 81.500 0.0000
18.750 81.750 0.0000
18.750 82.000 0.0000
18.750 82.250 0.0000
18.750 82.500 0.0000
18.750 82.750 0.0000
18.750 83.000 0.0000
18.750 83.250 0.0000
18.750 83.500 0.0000
18.750 83.750 0.0000
18.750 84.000 0.0000
18.750 84.250 0.0000
18.750 84.500 0.0000
18.750 84.750 0.0000
18.750 85.000 0.0000
18.750 85.250 0.0000
18.750 85.500 0.0000
18.750 85.750 0.0000
18.750 86.000 0.0000
18.750 86.250 0.0000
18.750 86.500 0.0000
18.750 86.750 0.0000
18.750 87.000 0.0000
18.750 87.250 0.0000
18.750 87.500 0.0000
18.750 87.750 0.0000
18.750 88.000 0.0000
18.750 88.250 0.0000
18.750 88.500 0.0000
18.750 88.750 0.0000
18.750 89.000 0.0000
18.750 89.250 0.0000
18.750 89.500 0.0000
18.750 89.750 0.0000
18.750 90.000 0.0000
18.750 90.250 0.0000
18.750 90.500 0.0000
18.750 90.750 0.0000
18.750 91.000 0.0000
18.750 91.250 0.0000
18.750 91.500 0.0000
18.750 91.750 0.0000
18.750 92.000 0.0000
18.750 92.250 0.0000
18.750 92.500 0.0000
18.750 92.750 0.0000
18.750 93.000 0.0000
19.000 81.000 0.0000
19.000 81.250 0.0000
19.000 81.500 0.0000
19.000 81.750 0.0000
19.000 82.000 0.0000
19.000 82.250 0.0000
19.000 82.500 0.0000
19.000 82.750 0.0000
19.000 83.000 0.0000
19.000 83.250 0.0000
19.000 83.500 0.0000
19.000 83.750 0.0000
19.000 84.000 0.0000
19.000 84.250 0.0000
19.000 84.500 0.0000
19.000 84.750 0.0000
19.000 85.000 0.0000
19.000 85.250 0.0000
19.000 85.500 0.0000
19.000 85.750 0.0000
19.000 86.000 0.0000
19.000 86.250 0.0000
19.000 86.500 0.0000
19.000 86.750 0.0000
19.000 87.000 0.0000
19.000 87.250 0.0000
19.000 87.500 0.0000
19.000 87.750 0.0000
19.000 88.000 0.0000
19.000 88.250 0.0000
19.000 88.500 0.0000
19.000 88.750 0.0000
19.000 89.000 0.0000
19.000 89.250 0.0000
19.000 89.500 0.0000
19.000 89.750 0.0000
19.000 90.000 0.0000
19.000 90.250 0.0000
19.000 90.500 0.0000
19.000 90.750 0.0000
19.000 91.000 0.0000
19.000 91.250 0.0000
19.000 91.500 0.0000
19.000 91.750 0.0000
19.000 92.000 0.0000
19.000 92.250 0.0000
19.000 92.500 0.0000
19.000 92.750 0.0000
19.000 93.000 0.0000
19.250 81.000 0.0000
19.250 81.250 0.0000
19.250 81.500 0.0000
19.250 81.750 0.0000
19.250 82.000 0.0000
19.250 82.250 0.0000
19.250 82.500 0.0000
19.250 82.750 0.0000
19.250 83.000 0.0000
19.250 83.250 0.2917
19.250 83.500 0.9825
19.250 83.750 0.9498
19.250 84.000 0.7938
19.250 84.250 0.0000
19.250 84.500 0.0000
19.250 84.750 0.0000
19.250 85.000 0.0000
19.250 85.250 0.0000
19.250 85.500 0.0000
19.250 85.750 0.0000
19.250 86.000 0.0000
19.250 86.250 0.0000
19.250 86.500 0.0000
19.250 86.750 0.0000
19.250 87.000 0.0000
19.250 87.250 0.0000
19.250 87.500 0.0000
19.250 87.750 0.0000
19.250 88.000 0.0000
19.250 88.250 0.0000
19.250 88.500 0.0000
19.250 88.750 0.0000
19.250 89.000 0.0000
19.250 89.250 0.0000
19.250 89.500 0.0000
19.250 89.750 0.0000
19.250 90.000 0.0000
19.250 90.250 0.0000
19.250 90.500 0.0000
19.250 90.750 0.0000
19.250 91.000 0.0000
19.250 91.250 0.0000
19.250 91.500 0.0000
19.250 91.750 0.0000
19.250 92.000 0.0000
19.250 92.250 0.0000
19.250 92.500 0.0000
19.250 92.750 0.0000
19.250 93.000 0.0000
19.500 81.000 0.0000
19.500 81.250 0.0000
19.500 81.500 0.0000
19.500 81.750 0.0000
19.500 82.000 0.0000
19.500 82.250 0.0000
19.500 82.500 0.0000
19.500 82.750 0.0000
19.500 83.000 0.3256
19.500 83.250 5.2225
19.500 83.500 5.3465
19.500 83.750 6.1165
19.500 84.000 0.0000
19.500 84.250 0.0000
19.500 84.500 0.0437
19.500 84.750 0.0000
19.500 85.000 0.0000
19.500 85.250 0.0000
19.500 85.500 0.0000
19.500 85.750 0.0000
19.500 86.000 0.0000
19.500 86.250 0.0000
19.500 86.500 0.0000
19.500 86.750 0.0000
19.500 87.000 0.0000
19.500 87.250 0.0000
19.500 87.500 0.0000
19.500 87.750 0.0000
19.500 88.000 0.0000
19.500 88.250 0.0000
19.500 88.500 0.0000
19.500 88.750 0.0000
19.500 89.000 0.0000
19.500 89.250 0.0000
19.500 89.500 0.0000
19.500 89.750 0.0000
19.500 90.000 0.0000
19.500 90.250 0.0000
19.500 90.500 0.0000
19.500 90.750 0.0000
19.500 91.000 0.0000
19.500 91.250 0.0000
19.500 91.500 0.0000
19.500 91.750 0.0000
19.500 92.000 0.0000
19.500 92.250 0.0000
19.500 92.500 0.0000
19.500 92.750 0.0000
19.500 93.000 0.0000
19.750 81.000 0.0000
19.750 81.250 0.0000
19.750 81.500 0.0000
19.750 81.750 0.0000
19.750 82.000 0.0000
19.750 82.250 0.0000
19.750 82.500 0.0000
19.750 82.750 0.4939
19.750 83.000 1.8164
19.750 83.250 0.6567
19.750 83.500 0.0000
19.750 83.750 0.0000
19.750 84.000 0.0000
19.750 84.250 0.0000
19.750 84.500 0.0000
19.750 84.750 0.0000
19.750 85.000 0.0000
19.750 85.250 0.0000
19.750 85.500 0.0000
19.750 85.750 0.0000
19.750 86.000 0.0000
19.750 86.250 0.0000
19.750 86.500 0.0000
19.750 86.750 0.0000
19.750 87.000 0.0000
19.750 87.250 0.0000
19.750 87.500 0.0000
19.750 87.750 0.0000
19.750 88.000 0.0000
19.750 88.250 0.0000
19.750 88.500 0.0000
19.750 88.750 0.0000
19.750 89.000 0.0000
19.750 89.250 0.0000
19.750 89.500 0.0000
19.750 89.750 0.0000
19.750 90.000 0.0000
19.750 90.250 0.0000
19.750 90.500 0.0000
19.750 90.750 0.0000
19.750 91.000 0.0000
19.750 91.250 0.0000
19.750 91.500 0.0000
19.750 91.750 0.0000
19.750 92.000 0.0000
19.750 92.250 0.0000
19.750 92.500 0.0000
19.750 92.750 0.0000
19.750 93.000 0.0000
20.000 81.000 0.0000
20.000 81.250 0.0000
20.000 81.500 0.0000
20.000 81.750 0.0000
20.000 82.000 0.0000
20.000 82.250 0.0000
20.000 82.500 4.0861
20.000 82.750 5.2817
20.000 83.000 5.1671
20.000 83.250 3.1528
20.000 83.500 0.0000
20.000 83.750 0.0000
20.000 84.000 0.0000
20.000 84.250 0.0000
20.000 84.500 0.1621
20.000 84.750 0.0000
20.000 85.000 0.0000
20.000 85.250 0.0000
20.000 85.500 0.0000
20.000 85.750 0.0000
20.000 86.000 0.0000
20.000 86.250 0.0000
20.000 86.500 0.0000
20.000 86.750 0.0000
20.000 87.000 0.0000
20.000 87.250 0.0000
20.000 87.500 0.0000
20.000 87.750 0.0000
20.000 88.000 0.0000
20.000 88.250 0.0000
20.000 88.500 0.0000
20.000 88.750 0.0000
20.000 89.000 0.0000
20.000 89.250 0.0000
20.000 89.500 0.0000
20.000 89.750 0.0000
20.000 90.000 0.0000
20.000 90.250 0.0000
20.000 90.500 0.0000
20.000 90.750 0.0000
20.000 91.000 0.0000
20.000 91.250 0.0000
20.000 91.500 0.0000
20.000 91.750 0.0000
20.000 92.000 0.0000
20.000 92.250 0.0000
20.000 92.500 0.0000
20.000 92.750 0.0000
20.000 93.000 0.0000
20.250 81.000 0.0000
20.250 81.250 0.0000
20.250 81.500 0.0000
20.250 81.750 0.0000
20.250 82.000 0.0000
20.250 82.250 1.0190
20.250 82.500 5.2732
20.250 82.750 4.9871
20.250 83.000 5.4140
20.250 83.250 5.7986
20.250 83.500 2.2281
20.250 83.750 2.4492
20.250 84.000 1.8091
20.250 84.250 1.2222
20.250 84.500 1.0645
20.250 84.750 0.0000
20.250 85.000 1.3844
20.250 85.250 0.0000
20.250 85.500 0.0000
20.250 85.750 0.0000
20.250 86.000 0.0000
20.250 86.250 0.0000
20.250 86.500 0.0000
20.250 86.750 0.0000
20.250 87.000 0.0000
20.250 87.250 0.0000
20.250 87.500 0.0000
20.250 87.750 0.0000
20.250 88.000 0.0000
20.250 88.250 0.0000
20.250 88.500 0.0000
20.250 88.750 0.0000
20.250 89.000 0.0000
20.250 89.250 0.0000
20.250 89.500 0.0000
20.250 89.750 0.0000
20.250 90.000 0.0000
20.250 90.250 0.0000
20.250 90.500 0.0000
20.250 90.750 0.0000
20.250 91.000 0.0000
20.250 91.250 0.0000
20.250 91.500 0.0000
20.250 91.750 0.0000
20.250 92.000 0.0000
20.250 92.250 0.0000
20.250 92.500 0.0000
20.250 92.750 0.0000
20.250 93.000 0.0000
20.500 81.000 0.0000
20.500 81.250 0.0000
20.500 81.500 0.0000
20.500 81.750 0.0000
20.500 82.000 0.0000
20.500 82.250 3.7129
20.500 82.500 4.4829
20.500 82.750 4.1209
20.500 83.000 1.1185
20.500 83.250 0.4919
20.500 83.500 3.4049
20.500 83.750 8.8447
20.500 84.000 4.5988
20.500 84.250 1.6684
20.500 84.500 0.9492
20.500 84.750 0.0000
20.500 85.000 0.5759
20.500 85.250 1.0255
20.500 85.500 0.0000
20.500 85.750 0.0000
20.500 86.000 0.0000
20.500 86.250 0.0000
20.500 86.500 0.0000
20.500 86.750 0.0000
20.500 87.000 0.0000
20.500 87.250 0.0000
20.500 87.500 0.0000
20.500 87.750 0.0000
20.500 88.000 0.0000
20.500 88.250 0.0000
20.500 88.500 0.0000
20.500 88.750 0.0000
20.500 89.000 0.0000
20.500 89.250 0.0000
20.500 89.500 0.0000
20.500 89.750 0.0000
20.500 90.000 0.0000
20.500 90.250 0.0000
20.500 90.500 0.0000
20.500 90.750 0.0000
20.500 91.000 0.0000
20.500 91.250 0.0000
20.500 91.500 0.0000
20.500 91.750 0.0000
20.500 92.000 0.0000
20.500 92.250 0.0000
20.500 92.500 0.0000
20.500 92.750 0.0000
20.500 93.000 0.0000
20.750 81.000 0.0000
20.750 81.250 0.0000
20.750 81.500 0.0000
20.750 81.750 0.0000
20.750 82.000 0.4198
20.750 82.250 0.0000
20.750 82.500 0.0000
20.750 82.750 0.0000
20.750 83.000 0.0000
20.750 83.250 0.1621
20.750 83.500 0.9673
20.750 83.750 9.6513
20.750 84.000 5.1845
20.750 84.250 1.4524
20.750 84.500 0.2574
20.750 84.750 0.0000
20.750 85.000 0.0000
20.750 85.250 0.0000
20.750 85.500 0.0000
20.750 85.750 0.0000
20.750 86.000 0.0000
20.750 86.250 0.0000
20.750 86.500 0.0000
20.750 86.750 0.0000
20.750 87.000 0.0000
20.750 87.250 0.0000
20.750 87.500 0.0000
20.750 87.750 0.0000
20.750 88.000 0.0000
20.750 88.250 0.0000
20.750 88.500 0.0000
20.750 88.750 0.0000
20.750 89.000 0.0000
20.750 89.250 0.0000
20.750 89.500 0.0000
20.750 89.750 0.0000
20.750 90.000 0.0000
20.750 90.250 0.0000
20.750 90.500 0.0000
20.750 90.750 0.0000
20.750 91.000 0.0000
20.750 91.250 0.0000
20.750 91.500 0.0000
20.750 91.750 0.0000
20.750 92.000 0.0000
20.750 92.250 0.0000
20.750 92.500 0.0000
20.750 92.750 0.0000
20.750 93.000 0.0000
21.000 81.000 0.0000
21.000 81.250 0.0000
21.000 81.500 0.0000
21.000 81.750 0.0000
21.000 82.000 1.5804
21.000 82.250 0.9089
21.000 82.500 0.0000
21.000 82.750 1.1782
21.000 83.000 1.4028
21.000 83.250 0.3216
21.000 83.500 0.7039
21.000 83.750 0.0000
21.000 84.000 0.0000
21.000 84.250 0.0520
21.000 84.500 0.0000
21.000 84.750 0.0000
21.000 85.000 0.0000
21.000 85.250 0.0000
21.000 85.500 0.0000
21.000 85.750 1.1791
21.000 86.000 1.0469
21.000 86.250 0.4821
21.000 86.500 0.0000
21.000 86.750 0.0000
21.000 87.000 0.0000
21.000 87.250 0.0000
21.000 87.500 0.0000
21.000 87.750 0.0000
21.000 88.000 0.0000
21.000 88.250 0.0000
21.000 88.500 0.0000
21.000 88.750 0.0000
21.000 89.000 0.0000
21.000 89.250 0.0000
21.000 89.500 0.0000
21.000 89.750 0.0000
21.000 90.000 0.0000
21.000 90.250 0.0000
21.000 90.500 0.0000
21.000 90.750 0.0000
21.000 91.000 0.0000
21.000 91.250 0.0000
21.000 91.500 0.0000
21.000 91.750 0.0000
21.000 92.000 0.0000
21.000 92.250 0.0000
21.000 92.500 0.0000
21.000 92.750 0.0000
21.000 93.000 0.0000
21.250 81.000 0.6446
21.250 81.250 0.0000
21.250 81.500 0.6624
21.250 81.750 3.7872
21.250 82.000 6.9334
21.250 82.250 4.6835
21.250 82.500 2.4336
21.250 82.750 2.6843
21.250 83.000 2.2329
21.250 83.250 0.6182
21.250 83.500 0.7137
21.250 83.750 0.4114
21.250 84.000 0.7473
21.250 84.250 0.2646
21.250 84.500 0.8019
21.250 84.750 0.4661
21.250 85.000 0.4362
21.250 85.250 0.0000
21.250 85.500 0.0000
21.250 85.750 1.3635
21.250 86.000 4.2196
21.250 86.250 3.5045
21.250 86.500 0.4692
21.250 86.750 0.0000
21.250 87.000 0.0000
21.250 87.250 0.0000
21.250 87.500 0.0000
21.250 87.750 0.0000
21.250 88.000 0.0000
21.250 88.250 0.0000
21.250 88.500 0.0000
21.250 88.750 0.0000
21.250 89.000 0.0000
21.250 89.250 0.0000
21.250 89.500 0.0000
21.250 89.750 0.0000
21.250 90.000 0.0000
21.250 90.250 0.0000
21.250 90.500 0.0000
21.250 90.750 0.0000
21.250 91.000 0.0000
21.250 91.250 0.0000
21.250 91.500 0.0000
21.250 91.750 0.0000
21.250 92.000 0.0000
21.250 92.250 0.0000
21.250 92.500 0.0000
21.250 92.750 0.0000
21.250 93.000 0.0000
21.500 81.000 2.3870
21.500 81.250 0.4604
21.500 81.500 3.1637
21.500 81.750 6.6653
21.500 82.000 7.3179
21.500 82.250 7.4931
21.500 82.500 6.7549
21.500 82.750 3.2660
21.500 83.000 3.1620
21.500 83.250 1.1821
21.500 83.500 0.5848
21.500 83.750 0.5517
21.500 84.000 1.3218
21.500 84.250 0.5931
21.500 84.500 1.9777
21.500 84.750 0.5139
21.500 85.000 0.2273
21.500 85.250 0.0000
21.500 85.500 0.0000
21.500 85.750 0.4809
21.500 86.000 1.3468
21.500 86.250 1.6364
21.500 86.500 2.2935
21.500 86.750 0.0000
21.500 87.000 0.0000
21.500 87.250 0.0000
21.500 87.500 0.0000
21.500 87.750 0.0000
21.500 88.000 0.0000
21.500 88.250 0.0000
21.500 88.500 0.0000
21.500 88.750 0.0000
21.500 89.000 0.0000
21.500 89.250 0.0000
21.500 89.500 0.0000
21.500 89.750 0.0000
21.500 90.000 0.0000
21.500 90.250 0.0000
21.500 90.500 0.0000
21.500 90.750 0.0000
21.500 91.000 0.0000
21.500 91.250 0.0000
21.500 91.500 0.0000
21.500 91.750 0.0000
21.500 92.000 0.0000
21.500 92.250 0.0000
21.500 92.500 0.0000
21.500 92.750 0.0000
21.500 93.000 0.0000
21.750 81.000 0.9251
21.750 81.250 0.9915
21.750 81.500 1.4599
21.750 81.750 3.7088
21.750 82.000 6.2689
21.750 82.250 6.3614
21.750 82.500 9.4681
21.750 82.750 3.9792
21.750 83.000 2.0967
21.750 83.250 0.9263
21.750 83.500 1.1022
21.750 83.750 0.8112
21.750 84.000 2.4276
21.750 84.250 3.8598
21.750 84.500 2.3712
21.750 84.750 1.7959
21.750 85.000 0.2839
21.750 85.250 0.0000
21.750 85.500 0.0000
21.750 85.750 0.0000
21.750 86.000 0.0000
21.750 86.250 0.5579
21.750 86.500 0.7785
21.750 86.750 0.2041
21.750 87.000 0.0000
21.750 87.250 0.0000
21.750 87.500 0.0000
21.750 87.750 0.0000
21.750 88.000 0.0000
21.750 88.250 0.0000
21.750 88.500 0.0000
21.750 88.750 0.0000
21.750 89.000 0.0000
21.750 89.250 0.0000
21.750 89.500 0.0000
21.750 89.750 0.0000
21.750 90.000 0.0000
21.750 90.250 0.0000
21.750 90.500 0.0000
21.750 90.750 0.0000
21.750 91.000 0.0000
21.750 91.250 0.0000
21.750 91.500 0.0000
21.750 91.750 0.0000
21.750 92.000 0.0000
21.750 92.250 0.0000
21.750 92.500 0.0000
21.750 92.750 0.0000
21.750 93.000 0.0000
22.000 81.000 0.7521
22.000 81.250 0.7876
22.000 81.500 0.9449
22.000 81.750 2.2642
22.000 82.000 3.8134
22.000 82.250 3.3622
22.000 82.500 3.4227
22.000 82.750 3.3661
22.000 83.000 0.7290
22.000 83.250 0.6111
22.000 83.500 1.0837
22.000 83.750 0.7547
22.000 84.000 0.8203
22.000 84.250 0.1461
22.000 84.500 0.6666
22.000 84.750 5.1010
22.000 85.000 2.2856
22.000 85.250 1.2195
22.000 85.500 0.6985
22.000 85.750 0.0000
22.000 86.000 0.0000
22.000 86.250 0.0000
22.000 86.500 0.0000
22.000 86.750 0.0000
22.000 87.000 0.0000
22.000 87.250 0.0000
22.000 87.500 0.0000
22.000 87.750 0.0000
22.000 88.000 0.0000
22.000 88.250 0.0000
22.000 88.500 0.0000
22.000 88.750 0.0000
22.000 89.000 0.0000
22.000 89.250 0.0000
22.000 89.500 0.0000
22.000 89.750 0.0000
22.000 90.000 0.0000
22.000 90.250 0.0000
22.000 90.500 0.0000
22.000 90.750 0.0000
22.000 91.000 0.0000
22.000 91.250 0.0000
22.000 91.500 0.0000
22.000 91.750 0.0000
22.000 92.000 0.0000
22.000 92.250 0.0000
22.000 92.500 0.0000
22.000 92.750 0.0000
22.000 93.000 0.0000
22.250 81.000 0.4883
22.250 81.250 0.3946
22.250 81.500 0.8132
22.250 81.750 0.7983
22.250 82.000 0.6530
22.250 82.250 0.0833
22.250 82.500 0.0906
22.250 82.750 1.0352
22.250 83.000 2.7045
22.250 83.250 6.7197
22.250 83.500 7.8344
22.250 83.750 2.8986
22.250 84.000 0.4641
22.250 84.250 1.7283
22.250 84.500 1.2827
22.250 84.750 0.9366
22.250 85.000 0.2860
22.250 85.250 0.2081
22.250 85.500 0.0000
22.250 85.750 0.0000
22.250 86.000 0.0000
22.250 86.250 0.0000
22.250 86.500 0.0000
22.250 86.750 0.0000
22.250 87.000 0.0000
22.250 87.250 0.0000
22.250 87.500 0.0000
22.250 87.750 0.0000
22.250 88.000 0.0000
22.250 88.250 0.0000
22.250 88.500 0.0000
22.250 88.750 0.0000
22.250 89.000 0.0000
22.250 89.250 0.0000
22.250 89.500 0.0000
22.250 89.750 0.0000
22.250 90.000 0.0000
22.250 90.250 0.0000
22.250 90.500 0.0000
22.250 90.750 0.0000
22.250 91.000 0.0000
22.250 91.250 0.0000
22.250 91.500 0.0000
22.250 91.750 0.0000
22.250 92.000 0.0000
22.250 92.250 0.0000
22.250 92.500 0.0000
22.250 92.750 0.0000
22.250 93.000 0.0000
22.500 81.000 1.2359
22.500 81.250 2.5590
22.500 81.500 1.5217
22.500 81.750 1.6180
22.500 82.000 0.1197
22.500 82.250 0.0797
22.500 82.500 0.0000
22.500 82.750 0.8486
22.500 83.000 2.6191
22.500 83.250 2.6061
22.500 83.500 11.2457
22.500 83.750 7.3395
22.500 84.000 0.3818
22.500 84.250 0.3108
22.500 84.500 0.4108
22.500 84.750 0.0000
22.500 85.000 0.0000
22.500 85.250 0.0000
22.500 85.500 0.0000
22.500 85.750 0.3574
22.500 86.000 1.2454
22.500 86.250 0.6831
22.500 86.500 0.0000
22.500 86.750 0.0000
22.500 87.000 0.0000
22.500 87.250 0.0000
22.500 87.500 0.0000
22.500 87.750 0.0000
22.500 88.000 0.0000
22.500 88.250 0.0000
22.500 88.500 0.0000
22.500 88.750 0.0000
22.500 89.000 0.0000
22.500 89.250 0.0000
22.500 89.500 0.0000
22.500 89.750 0.0000
22.500 90.000 0.0000
22.500 90.250 0.0000
22.500 90.500 0.0000
22.500 90.750 0.0000
22.500 91.000 0.0000
22.500 91.250 0.0000
22.500 91.500 0.0000
22.500 91.750 0.0000
22.500 92.000 0.0000
22.500 92.250 0.0000
22.500 92.500 0.0000
22.500 92.750 0.0000
22.500 93.000 0.0000
22.750 81.000 0.0000
22.750 81.250 0.0000
22.750 81.500 1.5469
22.750 81.750 0.3762
22.750 82.000 0.0000
22.750 82.250 0.0000
22.750 82.500 0.0000
22.750 82.750 0.0000
22.750 83.000 1.4566
22.750 83.250 0.3730
22.750 83.500 0.1746
22.750 83.750 0.0000
22.750 84.000 0.0000
22.750 84.250 0.0000
22.750 84.500 0.0000
22.750 84.750 0.0000
22.750 85.000 0.0000
22.750 85.250 0.9206
22.750 85.500 1.6599
22.750 85.750 1.9209
22.750 86.000 2.2685
22.750 86.250 4.5318
22.750 86.500 3.2441
22.750 86.750 0.3929
22.750 87.000 0.0000
22.750 87.250 3.3776
22.750 87.500 3.5554
22.750 87.750 1.5482
22.750 88.000 0.0000
22.750 88.250 0.0000
22.750 88.500 0.0000
22.750 88.750 0.0000
22.750 89.000 0.0000
22.750 89.250 0.0000
22.750 89.500 0.0000
22.750 89.750 0.0000
22.750 90.000 0.0000
22.750 90.250 0.0000
22.750 90.500 0.0000
22.750 90.750 0.0000
22.750 91.000 0.0000
22.750 91.250 0.0000
22.750 91.500 0.0000
22.750 91.750 0.0000
22.750 92.000 0.0000
22.750 92.250 0.0000
22.750 92.500 0.0000
22.750 92.750 0.0000
22.750 93.000 0.0000
23.000 81.000 0.0000
23.000 81.250 0.0000
23.000 81.500 0.0000
23.000 81.750 0.0799
23.000 82.000 0.0000
23.000 82.250 0.0000
23.000 82.500 0.0000
23.000 82.750 0.0000
23.000 83.000 0.0000
23.000 83.250 0.0000
23.000 83.500 0.0000
23.000 83.750 0.0000
23.000 84.000 0.0000
23.000 84.250 0.0000
23.000 84.500 0.0000
23.000 84.750 0.0000
23.000 85.000 0.0000
23.000 85.250 1.0908
23.000 85.500 5.4043
23.000 85.750 5.1232
23.000 86.000 5.2945
23.000 86.250 6.0117
23.000 86.500 1.4181
23.000 86.750 0.3871
23.000 87.000 0.6372
23.000 87.250 3.9512
23.000 87.500 4.6588
23.000 87.750 4.0816
23.000 88.000 2.1208
23.000 88.250 1.4213
23.000 88.500 0.3562
23.000 88.750 0.0000
23.000 89.000 0.0000
23.000 89.250 0.0000
23.000 89.500 0.0000
23.000 89.750 0.0000
23.000 90.000 0.0000
23.000 90.250 0.0000
23.000 90.500 0.0000
23.000 90.750 0.0000
23.000 91.000 0.0000
23.000 91.250 0.0000
23.000 91.500 0.0000
23.000 91.750 0.0000
23.000 92.000 0.0000
23.000 92.250 0.0000
23.000 92.500 0.0000
23.000 92.750 0.0000
23.000 93.000 0.0000
23.250 81.000 0.0845
23.250 81.250 0.0000
23.250 81.500 0.0000
23.250 81.750 0.0000
23.250 82.000 0.0000
23.250 82.250 0.0000
23.250 82.500 0.0000
23.250 82.750 0.0000
23.250 83.000 0.0000
23.250 83.250 0.0000
23.250 83.500 0.0000
23.250 83.750 0.0000
23.250 84.000 0.0000
23.250 84.250 0.0000
23.250 84.500 0.0000
23.250 84.750 0.0000
23.250 85.000 0.0000
23.250 85.250 3.6075
23.250 85.500 1.6726
23.250 85.750 1.1323
23.250 86.000 0.3922
23.250 86.250 2.2399
23.250 86.500 0.4699
23.250 86.750 3.0116
23.250 87.000 3.0673
23.250 87.250 4.1799
23.250 87.500 2.5479
23.250 87.750 1.4695
23.250 88.000 0.7510
23.250 88.250 0.4423
23.250 88.500 0.2105
23.250 88.750 0.3485
23.250 89.000 1.3237
23.250 89.250 1.1748
23.250 89.500 0.0000
23.250 89.750 0.0000
23.250 90.000 0.0000
23.250 90.250 0.0000
23.250 90.500 0.0000
23.250 90.750 0.0000
23.250 91.000 0.0000
23.250 91.250 0.0000
23.250 91.500 0.0000
23.250 91.750 0.0000
23.250 92.000 0.0000
23.250 92.250 0.0000
23.250 92.500 0.0000
23.250 92.750 0.0000
23.250 93.000 0.0000
23.500 81.000 0.0000
23.500 81.250 0.0000
23.500 81.500 0.0000
23.500 81.750 0.0000
23.500 82.000 0.0000
23.500 82.250 0.0000
23.500 82.500 0.0000
23.500 82.750 0.0000
23.500 83.000 0.0000
23.500 83.250 0.0000
23.500 83.500 0.0000
23.500 83.750 0.0000
23.500 84.000 0.7662
23.500 84.250 0.7470
23.500 84.500 0.0000
23.500 84.750 0.0000
23.500 85.000 0.0000
23.500 85.250 0.0000
23.500 85.500 0.0000
23.500 85.750 0.0000
23.500 86.000 0.0000
23.500 86.250 0.0000
23.500 86.500 3.5077
23.500 86.750 3.4698
23.500 87.000 4.5317
23.500 87.250 13.8581
23.500 87.500 4.2388
23.500 87.750 2.1122
23.500 88.000 0.0000
23.500 88.250 0.5257
23.500 88.500 0.7983
23.500 88.750 1.0454
23.500 89.000 1.7977
23.500 89.250 1.5289
23.500 89.500 0.8189
23.500 89.750 0.0000
23.500 90.000 0.1250
23.500 90.250 0.0000
23.500 90.500 0.0000
23.500 90.750 0.0000
23.500 91.000 0.0000
23.500 91.250 0.0000
23.500 91.500 0.0000
23.500 91.750 0.0000
23.500 92.000 0.0000
23.500 92.250 0.0000
23.500 92.500 0.0000
23.500 92.750 0.0000
23.500 93.000 0.0000
23.750 81.000 0.0000
23.750 81.250 0.0000
23.750 81.500 0.0000
23.750 81.750 0.0000
23.750 82.000 0.0000
23.750 82.250 0.0000
23.750 82.500 0.0000
23.750 82.750 0.0000
23.750 83.000 0.0000
23.750 83.250 0.0000
23.750 83.500 0.0000
23.750 83.750 0.0000
23.750 84.000 0.6893
23.750 84.250 5.1776
23.750 84.500 0.0000
23.750 84.750 0.0000
23.750 85.000 0.0000
23.750 85.250 0.0000
23.750 85.500 0.0000
23.750 85.750 0.0000
23.750 86.000 0.0000
23.750 86.250 0.0000
23.750 86.500 0.4825
23.750 86.750 0.8503
23.750 87.000 3.6547
23.750 87.250 1.8267
23.750 87.500 3.0571
23.750 87.750 4.5780
23.750 88.000 3.8249
23.750 88.250 0.0000
23.750 88.500 0.1696
23.750 88.750 0.5211
23.750 89.000 3.4791
23.750 89.250 4.5958
23.750 89.500 1.5661
23.750 89.750 0.0000
23.750 90.000 0.0000
23.750 90.250 0.0000
23.750 90.500 0.4679
23.750 90.750 0.0000
23.750 91.000 0.0000
23.750 91.250 0.0000
23.750 91.500 0.0000
23.750 91.750 0.0000
23.750 92.000 0.0000
23.750 92.250 0.0000
23.750 92.500 0.0000
23.750 92.750 0.0000
23.750 93.000 0.0000
24.000 81.000 0.0000
24.000 81.250 0.0000
24.000 81.500 0.0000
24.000 81.750 0.0000
24.000 82.000 0.0000
24.000 82.250 0.0000
24.000 82.500 0.0000
24.000 82.750 0.0000
24.000 83.000 0.0000
24.000 83.250 0.0000
24.000 83.500 0.0000
24.000 83.750 0.0000
24.000 84.000 0.7546
24.000 84.250 0.6742
24.000 84.500 0.7790
24.000 84.750 0.6740
24.000 85.000 1.0129
24.000 85.250 0.9121
24.000 85.500 1.1117
24.000 85.750 0.4746
24.000 86.000 0.0000
24.000 86.250 0.0000
24.000 86.500 0.0000
24.000 86.750 2.0795
24.000 87.000 4.0874
24.000 87.250 6.8610
24.000 87.500 1.6119
24.000 87.750 16.7975
24.000 88.000 19.9028
24.000 88.250 3.0389
24.000 88.500 0.0000
24.000 88.750 1.7491
24.000 89.000 7.7203
24.000 89.250 10.9378
24.000 89.500 7.2037
24.000 89.750 3.3923
24.000 90.000 2.4071
24.000 90.250 0.9209
24.000 90.500 0.7150
24.000 90.750 0.0000
24.000 91.000 0.0000
24.000 91.250 0.0000
24.000 91.500 0.0000
24.000 91.750 0.0000
24.000 92.000 11.0796
24.000 92.250 7.1117
24.000 92.500 5.1856
24.000 92.750 4.3372
24.000 93.000 0.0000
24.250 81.000 0.0000
24.250 81.250 0.0000
24.250 81.500 0.0000
24.250 81.750 0.0000
24.250 82.000 0.0000
24.250 82.250 0.0000
24.250 82.500 0.0000
24.250 82.750 0.0000
24.250 83.000 0.0000
24.250 83.250 0.0000
24.250 83.500 0.0000
24.250 83.750 0.0000
24.250 84.000 0.0000
24.250 84.250 0.0000
24.250 84.500 0.8655
24.250 84.750 0.0000
24.250 85.000 0.3634
24.250 85.250 0.0000
24.250 85.500 0.0000
24.250 85.750 2.2395
24.250 86.000 0.0000
24.250 86.250 0.0000
24.250 86.500 0.1107
24.250 86.750 1.7211
24.250 87.000 5.1319
24.250 87.250 5.1515
24.250 87.500 2.0358
24.250 87.750 13.2939
24.250 88.000 26.0125
24.250 88.250 8.7435
24.250 88.500 0.0000
24.250 88.750 0.6084
24.250 89.000 7.0076
24.250 89.250 12.8800
24.250 89.500 24.1701
24.250 89.750 15.9608
24.250 90.000 13.5364
24.250 90.250 5.7320
24.250 90.500 2.0490
24.250 90.750 0.0000
24.250 91.000 0.0000
24.250 91.250 0.0000
24.250 91.500 0.0000
24.250 91.750 1.9865
24.250 92.000 2.1183
24.250 92.250 2.8419
24.250 92.500 0.0000
24.250 92.750 0.0000
24.250 93.000 0.0000
24.500 81.000 0.0000
24.500 81.250 0.0000
24.500 81.500 0.0000
24.500 81.750 0.0000
24.500 82.000 0.0000
24.500 82.250 0.0000
24.500 82.500 0.0000
24.500 82.750 0.0000
24.500 83.000 0.0000
24.500 83.250 0.0000
24.500 83.500 0.0000
24.500 83.750 0.0000
24.500 84.000 0.0000
24.500 84.250 0.0000
24.500 84.500 0.0000
24.500 84.750 0.0000
24.500 85.000 0.0000
24.500 85.250 0.0000
24.500 85.500 0.0000
24.500 85.750 0.0000
24.500 86.000 0.0000
24.500 86.250 0.0000
24.500 86.500 0.0000
24.500 86.750 1.0342
24.500 87.000 1.4762
24.500 87.250 9.5031
24.500 87.500 18.9740
24.500 87.750 15.6528
24.500 88.000 15.3174
24.500 88.250 14.7985
24.500 88.500 2.8200
24.500 88.750 1.3109
24.500 89.000 4.7647
24.500 89.250 9.6441
24.500 89.500 26.1013
24.500 89.750 13.6605
24.500 90.000 9.6030
24.500 90.250 7.3814
24.500 90.500 5.3098
24.500 90.750 2.7171
24.500 91.000 1.6197
24.500 91.250 0.1139
24.500 91.500 0.0000
24.500 91.750 0.0000
24.500 92.000 0.0000
24.500 92.250 0.0000
24.500 92.500 0.0000
24.500 92.750 0.0000
24.500 93.000 0.0000
24.750 81.000 0.0000
24.750 81.250 0.0000
24.750 81.500 0.0000
24.750 81.750 0.0000
24.750 82.000 0.0000
24.750 82.250 0.0000
24.750 82.500 0.0000
24.750 82.750 0.0000
24.750 83.000 0.0000
24.750 83.250 0.0000
24.750 83.500 0.0000
24.750 83.750 0.0000
24.750 84.000 0.0000
24.750 84.250 0.0000
24.750 84.500 0.0000
24.750 84.750 0.0000
24.750 85.000 0.0000
24.750 85.250 0.0000
24.750 85.500 0.0000
24.750 85.750 2.0508
24.750 86.000 0.0000
24.750 86.250 0.0000
24.750 86.500 0.0000
24.750 86.750 0.0000
24.750 87.000 0.0000
24.750 87.250 1.7051
24.750 87.500 2.3009
24.750 87.750 3.9556
24.750 88.000 7.6299
24.750 88.250 6.1145
24.750 88.500 17.9022
24.750 88.750 30.5539
24.750 89.000 26.8550
24.750 89.250 28.6321
24.750 89.500 15.2418
24.750 89.750 4.5243
24.750 90.000 5.7112
24.750 90.250 7.6150
24.750 90.500 5.5750
24.750 90.750 3.7909
24.750 91.000 2.0234
24.750 91.250 0.6670
24.750 91.500 0.0000
24.750 91.750 0.0000
24.750 92.000 0.0000
24.750 92.250 0.0000
24.750 92.500 0.0000
24.750 92.750 0.0000
24.750 93.000 0.0000
25.000 81.000 0.0000
25.000 81.250 0.0000
25.000 81.500 0.0000
25.000 81.750 0.0000
25.000 82.000 0.0000
25.000 82.250 0.0000
25.000 82.500 0.0000
25.000 82.750 0.0000
25.000 83.000 0.0000
25.000 83.250 0.0000
25.000 83.500 0.0000
25.000 83.750 0.0000
25.000 84.000 0.0000
25.000 84.250 0.0000
25.000 84.500 0.0000
25.000 84.750 0.0000
25.000 85.000 0.0000
25.000 85.250 0.0000
25.000 85.500 0.0000
25.000 85.750 0.0000
25.000 86.000 0.6740
25.000 86.250 3.7727
25.000 86.500 0.0000
25.000 86.750 0.0000
25.000 87.000 0.0000
25.000 87.250 0.0000
25.000 87.500 0.0000
25.000 87.750 0.0000
25.000 88.000 1.9253
25.000 88.250 6.4654
25.000 88.500 10.6191
25.000 88.750 18.4200
25.000 89.000 23.8606
25.000 89.250 27.8211
25.000 89.500 27.6202
25.000 89.750 3.1970
25.000 90.000 1.5273
25.000 90.250 2.7655
25.000 90.500 1.6348
25.000 90.750 0.0562
25.000 91.000 0.0000
25.000 91.250 0.0000
25.000 91.500 0.0000
25.000 91.750 0.5729
25.000 92.000 0.0000
25.000 92.250 0.0000
25.000 92.500 0.0000
25.000 92.750 0.0000
25.000 93.000 0.0000
25.250 81.000 0.0000
25.250 81.250 0.0000
25.250 81.500 0.0000
25.250 81.750 0.0000
25.250 82.000 0.0000
25.250 82.250 0.0000
25.250 82.500 0.0000
25.250 82.750 0.0000
25.250 83.000 0.0000
25.250 83.250 0.0000
25.250 83.500 0.0000
25.250 83.750 0.0000
25.250 84.000 0.0000
25.250 84.250 0.0000
25.250 84.500 0.0000
25.250 84.750 0.0000
25.250 85.000 0.0000
25.250 85.250 0.0000
25.250 85.500 0.0000
25.250 85.750 0.0000
25.250 86.000 2.4606
25.250 86.250 5.6687
25.250 86.500 5.2623
25.250 86.750 0.6120
25.250 87.000 0.4580
25.250 87.250 0.7092
25.250 87.500 0.0000
25.250 87.750 0.0000
25.250 88.000 0.0000
25.250 88.250 0.0000
25.250 88.500 3.9000
25.250 88.750 23.7553
25.250 89.000 18.7957
25.250 89.250 8.4121
25.250 89.500 3.7616
25.250 89.750 2.9547
25.250 90.000 1.4157
25.250 90.250 0.6391
25.250 90.500 0.2186
25.250 90.750 0.5867
25.250 91.000 2.6515
25.250 91.250 2.7777
25.250 91.500 1.7957
25.250 91.750 0.3232
25.250 92.000 0.0000
25.250 92.250 0.0000
25.250 92.500 0.0000
25.250 92.750 0.0000
25.250 93.000 0.0000
25.500 81.000 0.0000
25.500 81.250 0.0000
25.500 81.500 0.0000
25.500 81.750 0.0000
25.500 82.000 0.0000
25.500 82.250 0.0000
25.500 82.500 0.0000
25.500 82.750 0.0000
25.500 83.000 0.0000
25.500 83.250 0.0000
25.500 83.500 0.0000
25.500 83.750 0.0000
25.500 84.000 0.0000
25.500 84.250 0.0000
25.500 84.500 0.0000
25.500 84.750 0.0000
25.500 85.000 0.0000
25.500 85.250 0.0000
25.500 85.500 0.0000
25.500 85.750 0.0000
25.500 86.000 3.2335
25.500 86.250 2.7110
25.500 86.500 2.5816
25.500 86.750 0.9000
25.500 87.000 0.0000
25.500 87.250 1.2000
25.500 87.500 0.4504
25.500 87.750 0.0000
25.500 88.000 0.0000
25.500 88.250 0.0000
25.500 88.500 0.0000
25.500 88.750 0.0000
25.500 89.000 0.5967
25.500 89.250 0.0000
25.500 89.500 0.0000
25.500 89.750 0.0000
25.500 90.000 1.1404
25.500 90.250 2.8041
25.500 90.500 1.2789
25.500 90.750 0.7873
25.500 91.000 0.6205
25.500 91.250 0.0000
25.500 91.500 0.0000
25.500 91.750 0.0000
25.500 92.000 0.0000
25.500 92.250 0.0000
25.500 92.500 0.0000
25.500 92.750 0.0000
25.500 93.000 0.9389
25.750 81.000 0.0000
25.750 81.250 0.0000
25.750 81.500 0.0000
25.750 81.750 0.0000
25.750 82.000 0.0000
25.750 82.250 0.0000
25.750 82.500 0.0000
25.750 82.750 0.0000
25.750 83.000 0.0000
25.750 83.250 0.0000
25.750 83.500 0.0000
25.750 83.750 0.0000
25.750 84.000 0.0000
25.750 84.250 0.0000
25.750 84.500 0.0000
25.750 84.750 0.0000
25.750 85.000 0.0000
25.750 85.250 0.0000
25.750 85.500 0.0000
25.750 85.750 0.0000
25.750 86.000 0.6251
25.750 86.250 0.7778
25.750 86.500 0.0000
25.750 86.750 0.0000
25.750 87.000 0.0000
25.750 87.250 0.0000
25.750 87.500 0.0000
25.750 87.750 0.0000
25.750 88.000 0.0000
25.750 88.250 0.0000
25.750 88.500 0.0000
25.750 88.750 0.0000
25.750 89.000 0.0000
25.750 89.250 0.0000
25.750 89.500 0.0000
25.750 89.750 7.8566
25.750 90.000 7.2363
25.750 90.250 19.1341
25.750 90.500 3.0047
25.750 90.750 0.7375
25.750 91.000 0.5885
25.750 91.250 0.0000
25.750 91.500 0.0000
25.750 91.750 0.0000
25.750 92.000 0.0000
25.750 92.250 1.8258
25.750 92.500 3.1952
25.750 92.750 1.3504
25.750 93.000 0.0000
26.000 81.000 0.0000
26.000 81.250 0.0000
26.000 81.500 0.0000
26.000 81.750 0.0000
26.000 82.000 0.0000
26.000 82.250 0.0000
26.000 82.500 0.0000
26.000 82.750 0.0000
26.000 83.000 0.0000
26.000 83.250 0.0000
26.000 83.500 0.0000
26.000 83.750 0.0000
26.000 84.000 0.0000
26.000 84.250 0.0000
26.000 84.500 0.0000
26.000 84.750 0.0000
26.000 85.000 0.0000
26.000 85.250 0.0000
26.000 85.500 0.0000
26.000 85.750 0.0000
26.000 86.000 0.0000
26.000 86.250 0.0000
26.000 86.500 0.0000
26.000 86.750 0.0000
26.000 87.000 0.0000
26.000 87.250 0.0000
26.000 87.500 0.0000
26.000 87.750 0.0000
26.000 88.000 0.0000
26.000 88.250 0.0000
26.000 88.500 0.0000
26.000 88.750 0.0000
26.000 89.000 0.0000
26.000 89.250 0.0000
26.000 89.500 0.0000
26.000 89.750 1.5544
26.000 90.000 3.2695
26.000 90.250 2.6554
26.000 90.500 0.3053
26.000 90.750 0.0000
26.000 91.000 0.0000
26.000 91.250 0.0000
26.000 91.500 0.0000
26.000 91.750 0.0000
26.000 92.000 0.0000
26.000 92.250 5.0790
26.000 92.500 8.1001
26.000 92.750 0.0000
26.000 93.000 0.0000
26.250 81.000 0.0000
26.250 81.250 0.0000
26.250 81.500 0.0000
26.250 81.750 0.0000
26.250 82.000 0.0000
26.250 82.250 0.0000
26.250 82.500 0.0000
26.250 82.750 0.0000
26.250 83.000 0.0000
26.250 83.250 0.0000
26.250 83.500 0.0000
26.250 83.750 0.0000
26.250 84.000 0.0000
26.250 84.250 0.0000
26.250 84.500 0.0000
26.250 84.750 0.0000
26.250 85.000 0.0000
26.250 85.250 0.0000
26.250 85.500 0.0000
26.250 85.750 0.0000
26.250 86.000 0.0000
26.250 86.250 0.0000
26.250 86.500 0.0000
26.250 86.750 0.0000
26.250 87.000 0.0000
26.250 87.250 0.0000
26.250 87.500 0.0000
26.250 87.750 0.0000
26.250 88.000 0.0000
26.250 88.250 0.0000
26.250 88.500 0.0000
26.250 88.750 0.0000
26.250 89.000 0.0000
26.250 89.250 0.0000
26.250 89.500 1.2599
26.250 89.750 2.2836
26.250 90.000 25.4422
26.250 90.250 6.9100
26.250 90.500 1.7416
26.250 90.750 0.4511
26.250 91.000 1.1717
26.250 91.250 2.0048
26.250 91.500 0.0000
26.250 91.750 0.9111
26.250 92.000 0.6298
26.250 92.250 1.8069
26.250 92.500 1.5123
26.250 92.750 0.0994
26.250 93.000 0.0000
26.500 81.000 0.0000
26.500 81.250 0.0000
26.500 81.500 0.0000
26.500 81.750 0.0000
26.500 82.000 0.0000
26.500 82.250 0.0000
26.500 82.500 0.0000
26.500 82.750 0.0000
26.500 83.000 0.0000
26.500 83.250 0.0000
26.500 83.500 0.0000
26.500 83.750 0.0000
26.500 84.000 0.0000
26.500 84.250 0.0000
26.500 84.500 0.0000
26.500 84.750 0.0000
26.500 85.000 0.0000
26.500 85.250 0.0000
26.500 85.500 0.0000
26.500 85.750 0.0000
26.500 86.000 0.0000
26.500 86.250 0.0000
26.500 86.500 0.0000
26.500 86.750 0.0000
26.500 87.000 0.0000
26.500 87.250 0.0000
26.500 87.500 0.0000
26.500 87.750 0.0000
26.500 88.000 0.0000
26.500 88.250 0.0000
26.500 88.500 0.0000
26.500 88.750 0.0000
26.500 89.000 0.0000
26.500 89.250 0.0000
26.500 89.500 0.6693
26.500 89.750 1.6922
26.500 90.000 2.5320
26.500 90.250 13.9450
26.500 90.500 11.3046
26.500 90.750 0.4768
26.500 91.000 0.8766
26.500 91.250 5.5705
26.500 91.500 1.8254
26.500 91.750 1.2097
26.500 92.000 0.3286
26.500 92.250 4.5842
26.500 92.500 2.4961
26.500 92.750 0.1137
26.500 93.000 0.0000
26.750 81.000 0.0000
26.750 81.250 0.0000
26.750 81.500 0.0000
26.750 81.750 0.0000
26.750 82.000 0.0000
26.750 82.250 0.0000
26.750 82.500 0.0000
26.750 82.750 0.0000
26.750 83.000 0.0000
26.750 83.250 0.0000
26.750 83.500 0.0000
26.750 83.750 0.0000
26.750 84.000 0.0000
26.750 84.250 0.0000
26.750 84.500 0.0000
26.750 84.750 0.0000
26.750 85.000 0.0000
26.750 85.250 0.0000
26.750 85.500 0.0000
26.750 85.750 0.0000
26.750 86.000 0.0000
26.750 86.250 0.0000
26.750 86.500 0.0000
26.750 86.750 0.0000
26.750 87.000 0.0000
26.750 87.250 0.0000
26.750 87.500 0.0000
26.750 87.750 0.0000
26.750 88.000 0.0000
26.750 88.250 0.0000
26.750 88.500 0.0000
26.750 88.750 0.0000
26.750 89.000 0.0000
26.750 89.250 0.0000
26.750 89.500 0.0000
26.750 89.750 0.1314
26.750 90.000 0.2427
26.750 90.250 2.4675
26.750 90.500 3.4485
26.750 90.750 0.3351
26.750 91.000 0.0000
26.750 91.250 0.0000
26.750 91.500 0.0000
26.750 91.750 0.0000
26.750 92.000 0.0000
26.750 92.250 0.0000
26.750 92.500 0.0000
26.750 92.750 0.0000
26.750 93.000 0.0000
27.000 81.000 0.0000
27.000 81.250 0.0000
27.000 81.500 0.0000
27.000 81.750 0.0000
27.000 82.000 0.0000
27.000 82.250 0.0000
27.000 82.500 0.0000
27.000 82.750 0.0000
27.000 83.000 0.0000
27.000 83.250 0.0000
27.000 83.500 0.0000
27.000 83.750 0.0000
27.000 84.000 0.0000
27.000 84.250 0.0000
27.000 84.500 0.0000
27.000 84.750 0.0000
27.000 85.000 0.0000
27.000 85.250 0.0000
27.000 85.500 0.0000
27.000 85.750 0.0000
27.000 86.000 0.0000
27.000 86.250 0.0000
27.000 86.500 0.0000
27.000 86.750 0.0000
27.000 87.000 0.3511
27.000 87.250 0.4371
27.000 87.500 0.0000
27.000 87.750 0.1989
27.000 88.000 0.3567
27.000 88.250 0.3921
27.000 88.500 0.3459
27.000 88.750 1.1956
27.000 89.000 1.4188
27.000 89.250 2.0937
27.000 89.500 1.0448
27.000 89.750 0.4623
27.000 90.000 0.0000
27.000 90.250 0.0000
27.000 90.500 0.0000
27.000 90.750 0.0000
27.000 91.000 0.0000
27.000 91.250 0.0000
27.000 91.500 0.0000
27.000 91.750 0.0000
27.000 92.000 0.0000
27.000 92.250 0.7821
27.000 92.500 0.0000
27.000 92.750 0.0000
27.000 93.000 0.0000
27.250 81.000 0.0000
27.250 81.250 0.0000
27.250 81.500 0.0000
27.250 81.750 0.0000
27.250 82.000 0.0000
27.250 82.250 0.0000
27.250 82.500 0.0000
27.250 82.750 0.0000
27.250 83.000 0.0000
27.250 83.250 0.0000
27.250 83.500 0.0000
27.250 83.750 0.0000
27.250 84.000 0.0000
27.250 84.250 0.0000
27.250 84.500 0.0000
27.250 84.750 0.0000
27.250 85.000 0.0000
27.250 85.250 0.0000
27.250 85.500 0.0000
27.250 85.750 0.0000
27.250 86.000 0.0000
27.250 86.250 0.0000
27.250 86.500 0.0000
27.250 86.750 0.0000
27.250 87.000 1.4797
27.250 87.250 2.2486
27.250 87.500 3.6513
27.250 87.750 1.9698
27.250 88.000 1.0157
27.250 88.250 0.0000
27.250 88.500 0.5744
27.250 88.750 1.9637
27.250 89.000 1.4979
27.250 89.250 1.5428
27.250 89.500 1.1052
27.250 89.750 0.1601
27.250 90.000 0.0000
27.250 90.250 0.2092
27.250 90.500 0.0000
27.250 90.750 0.0000
27.250 91.000 0.0000
27.250 91.250 0.0000
27.250 91.500 0.0000
27.250 91.750 0.0000
27.250 92.000 0.0000
27.250 92.250 0.0000
27.250 92.500 0.0000
27.250 92.750 0.0000
27.250 93.000 0.0000
27.500 81.000 0.0000
27.500 81.250 0.0000
27.500 81.500 0.0000
27.500 81.750 0.0000
27.500 82.000 0.0000
27.500 82.250 0.0000
27.500 82.500 0.0000
27.500 82.750 0.0000
27.500 83.000 0.0000
27.500 83.250 0.0000
27.500 83.500 0.0000
27.500 83.750 0.0000
27.500 84.000 0.0000
27.500 84.250 0.0000
27.500 84.500 0.0000
27.500 84.750 0.0000
27.500 85.000 0.0000
27.500 85.250 0.0000
27.500 85.500 0.0000
27.500 85.750 0.0000
27.500 86.000 0.0000
27.500 86.250 0.0000
27.500 86.500 0.0000
27.500 86.750 0.3601
27.500 87.000 0.6042
27.500 87.250 0.0000
27.500 87.500 0.0000
27.500 87.750 0.0000
27.500 88.000 0.0000
27.500 88.250 0.0000
27.500 88.500 0.0000
27.500 88.750 0.0000
27.500 89.000 0.1947
27.500 89.250 0.2008
27.500 89.500 0.0000
27.500 89.750 0.0000
27.500 90.000 0.0000
27.500 90.250 0.0000
27.500 90.500 0.0000
27.500 90.750 0.0000
27.500 91.000 0.0000
27.500 91.250 0.0000
27.500 91.500 0.0000
27.500 91.750 0.0000
27.500 92.000 0.0000
27.500 92.250 0.0000
27.500 92.500 0.0000
27.500 92.750 0.0000
27.500 93.000 0.0000
27.750 81.000 0.0000
27.750 81.250 0.0000
27.750 81.500 0.0000
27.750 81.750 0.0000
27.750 82.000 0.0000
27.750 82.250 0.0000
27.750 82.500 0.0000
27.750 82.750 0.0000
27.750 83.000 0.0000
27.750 83.250 0.0000
27.750 83.500 0.0000
27.750 83.750 0.0000
27.750 84.000 0.0000
27.750 84.250 0.0000
27.750 84.500 0.0000
27.750 84.750 0.0000
27.750 85.000 0.0000
27.750 85.250 0.0000
27.750 85.500 0.0000
27.750 85.750 0.0000
27.750 86.000 0.0000
27.750 86.250 0.0000
27.750 86.500 0.0000
27.750 86.750 0.0000
27.750 87.000 0.0000
27.750 87.250 0.0000
27.750 87.500 0.0000
27.750 87.750 0.0000
27.750 88.000 0.0000
27.750 88.250 0.0000
27.750 88.500 0.0000
27.750 88.750 0.5924
27.750 89.000 0.0000
27.750 89.250 0.0000
27.750 89.500 0.0000
27.750 89.750 0.0000
27.750 90.000 0.0000
27.750 90.250 0.0000
27.750 90.500 0.0000
27.750 90.750 0.0000
27.750 91.000 0.0000
27.750 91.250 0.0000
27.750 91.500 0.0000
27.750 91.750 0.0000
27.750 92.000 0.0000
27.750 92.250 0.0000
27.750 92.500 0.0000
27.750 92.750 0.0000
27.750 93.000 3.5693
28.000 81.000 0.0000
28.000 81.250 0.0000
28.000 81.500 0.0000
28.000 81.750 0.0000
28.000 82.000 0.0000
28.000 82.250 0.0000
28.000 82.500 0.0000
28.000 82.750 0.1259
28.000 83.000 0.3984
28.000 83.250 0.0000
28.000 83.500 0.0000
28.000 83.750 0.0000
28.000 84.000 0.0000
28.000 84.250 0.0000
28.000 84.500 0.0000
28.000 84.750 0.0000
28.000 85.000 0.0000
28.000 85.250 0.2065
28.000 85.500 0.0715
28.000 85.750 0.0000
28.000 86.000 0.0000
28.000 86.250 0.0000
28.000 86.500 0.0000
28.000 86.750 0.0000
28.000 87.000 0.0000
28.000 87.250 0.0000
28.000 87.500 0.0000
28.000 87.750 0.0000
28.000 88.000 0.0000
28.000 88.250 0.0000
28.000 88.500 2.5652
28.000 88.750 0.0000
28.000 89.000 0.0000
28.000 89.250 0.0000
28.000 89.500 0.0000
28.000 89.750 0.0000
28.000 90.000 0.0000
28.000 90.250 0.0000
28.000 90.500 0.0000
28.000 90.750 0.4598
28.000 91.000 0.2349
28.000 91.250 0.9012
28.000 91.500 0.0000
28.000 91.750 0.0000
28.000 92.000 0.0000
28.000 92.250 0.0000
28.000 92.500 0.0000
28.000 92.750 0.7194
28.000 93.000 4.6987
------------------------------------------------
Subject: Error in converting TRMM 3-hr rainfall into MET netcdf format
From: John Halley Gotway
Time: Wed Oct 05 15:36:43 2011
Krishna,
OK, I figured it out. The problem was that for some reason, my
version of the R ncdf package was out of sync with what's currently
available. I've modified the trmm2nc.R script and it should now
work as expected. It required a very minor change.
Please see attached, and let me know if you continue to experience
problems.
Thanks,
John
On 10/05/2011 12:13 PM, krishna kishore via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>
> Thank you,
> I am attaching the same file with this mail.
> Please explain what is the error...
> Krishna.
>
> On Wed, Oct 5, 2011 at 11:40 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>> Krishna,
>>
>> It would probably be easiest for you just to send me the file
that's giving
>> the error. Can you please send me this file:
>> TRMM_rain_20070420-12_15UTC.txt
>>
>> I'll run it here and try to diagnose what's going on. You can
either send
>> it via email or post it to our anonymous ftp site, as described
here:
>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>> If you do post it to our ftp site, please write me back to tell me
when
>> it's there.
>>
>> Hopefully we can get this working for you.
>>
>> Thanks,
>> John
>>
>>
>>
>> On 10/05/2011 09:48 AM, krishna kishore via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>>>
>>> Hi John,
>>> Thank you very much..
>>> Now I tried to convert the TRMM data....Still giving error..
>>>
>>
**************************************************************************
>>> /home/osuri/Desktop/TRMM2NC_for_MET>Rscript trmm2nc.R
>>> TRMM_rain_20070420-12_15UTC.txt TRMM_rain_20070420-12_15UTC.nc
>>> Reading: TRMM_rain_20070420-12_15UTC.txt
>>> Error in dim.def.ncdf("lat", "degrees_north", lat, create_dimvar =
TRUE,
>> :
>>> unused argument(s) (longname = "latitude")
>>> Execution halted
>>>
>>
**************************************************************************
>>> Can you check what is the problem?
>>> Thanks,
>>> Krishna.
>>>
>>>
>>> On Wed, Oct 5, 2011 at 8:50 PM, John Halley Gotway via RT <
>> met_help at ucar.edu
>>>> wrote:
>>>
>>>> Krishna,
>>>>
>>>> Ah, I see the confusion. Those commands I sent should be run
from the
>>>> command line - NOT FROM WITHIN R. Please try running the
commands
>> again,
>>>> but do NOT call R to start up an R session first.
>>>>
>>>> Hopefully, that'll do it.
>>>>
>>>> Thanks,
>>>> John
>>>>
>>>> On 10/04/2011 04:23 PM, krishna kishore via RT wrote:
>>>>>
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>>>>>
>>>>> Hi John,
>>>>> I Installed the ncdf library in R.
>>>>> For your more information....I Installed R software in my
ubuntu linux
>>>>> machine (through sudo apt-get install r-base) ....and installed
the
>> ncdf
>>>>> library also.
>>>>> Then I tried with the commands
>>>>> R for R prompt and then try to run the trmm2nc.R script
>>>>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
>>>>> TRMM_rain_20070420-9_12UTC.nc
>>>>> Error: unexpected symbol in "Rscript trmm2nc.R"
>>>>>
>>>>>
>>>>>> R --version
>>>>> Error: object 'R' not found
>>>>>
>>>>>> which Rscript
>>>>> Error: unexpected symbol in "which Rscript"
>>>>>
>>>>>> library(ncdf)
>>>>> nothing is displaying....means installed the library correctly..
>>>>>
>>>>> But why the above errors are giving when i am executin Rscript
and R
>>>>> --version commands.
>>>>>
>>>>> Thats means, R software is not installed in my account properly?
>>>>> What I should do to pass through this problem..
>>>>>
>>>>> Thanks,
>>>>> Krishna.
>>>>>
>>>>>
>>>>> On Mon, Sep 26, 2011 at 11:50 PM, RAL HelpDesk {for John Halley
Gotway}
>> <
>>>>> met_help at ucar.edu> wrote:
>>>>>
>>>>>> Krishna,
>>>>>>
>>>>>> I used the trmm2nc.R script to convert the sample data file you
sent,
>>>> and
>>>>>> everything worked fine.
>>>>>>
>>>>>> So there is likely a problem in the environment somewhere. I
have
>> never
>>>>>> seen that sort of error message from R before. Let's try
looking how
>> R
>>>> is
>>>>>> installed on your system:
>>>>>>
>>>>>> Check what version of R you're running:
>>>>>> R --version
>>>>>> Check that you have the Rscript command in your path:
>>>>>> which Rscript
>>>>>> Try running the Rscript command:
>>>>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
>>>>>> TRMM_rain_20070420-9_12UTC.nc
>>>>>>
>>>>>> Please note that this script expects that the R "ncdf" package
is
>>>> installed
>>>>>> on your system. To check if it's there, just try executing the
>>>> following
>>>>>> commands:
>>>>>> R
>>>>>> library(ncdf)
>>>>>>
>>>>>> If it loads without error, everything is fine. If you get an
error,
>>>> you'll
>>>>>> need to install that "ncdf" package. To do so, execute the
following
>>>>>> commands:
>>>>>> R
>>>>>> install.packages("ncdf",
>>>>>> configure.args="--with-netcdf_include=/usr/local/netcdf/include
>>>>>> -with-netcdf_lib=/usr/local/netcdf/lib")
>>>>>>
>>>>>> Please let me know what you learn from these steps.
>>>>>>
>>>>>> Thanks,
>>>>>> John Halley Gotway
>>>>>> met_help at ucar.edu
>>>>>>
>>>>>> On 09/23/2011 08:01 PM, RAL HelpDesk {for krishna kishore}
wrote:
>>>>>>>
>>>>>>> Fri Sep 23 20:01:29 2011: Request 49784 was acted upon.
>>>>>>> Transaction: Ticket created by osurikishore at gmail.com
>>>>>>> Queue: met_help
>>>>>>> Subject: Error in converting TRMM 3-hr rainfall into MET
netcdf
>>>>>> format
>>>>>>> Owner: Nobody
>>>>>>> Requestors: osurikishore at gmail.com
>>>>>>> Status: new
>>>>>>> Ticket <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784>
>>>>>>>
>>>>>>>
>>>>>>> Dear Methelp,
>>>>>>> I tried to convert the TRMM 3-hrly accumulated data (from
>>>>>>>
http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml)
>> using
>>>>>>> trmm2nc.R file.
>>>>>>> I copied the ASCII rainfall data and attached for your
reference.
>>>>>>> When i am running, it is giving error that
>>>>>>> *Error: unexpected symbol in "Rscript trmm2nc.R"*
>>>>>>>
>>>>>>> Can you help me why it is like that.
>>>>>>> Is this with NETCDF library?
>>>>>>> For your information, I have netcdf library and already export
that
>>>> path
>>>>>> in
>>>>>>> my .bashrc file.
>>>>>>> Will the Rscript automatically takes or should I mentioned
somewhere
>>>> when
>>>>>>> using this Rscript.
>>>>>>>
>>>>>>> Can you suggest what went wrong?
>>>>>>> Thanks in advance,
>>>>>>> Krishna.
>>>>>>>
>>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>
>
------------------------------------------------
Subject: Error in converting TRMM 3-hr rainfall into MET netcdf format
From: John Halley Gotway
Time: Wed Oct 05 15:36:43 2011
########################################################################
##
## Name: trmm2nc.R
##
## Description:
## This script should be called via an Rscript command to
reformat
## ASCII TRMM precipitation pcp into a NetCDF format that MET
## can read. The input pcp for this script can be attained via:
## 24-hour accumulations:
##
http://disc2.nascom.nasa.gov/Giovanni/tovas/realtime.3B42RT_daily.shtml
## 03-hour accumulations:
##
http://disc2.nascom.nasa.gov/Giovanni/tovas/realtime.3B42RT.shtml
##
## Usage:
## Rscript trmm2nc.R
## trmm_file
## nc_file
## [-save]
##
## Arguments:
## "trmm_file" is an ASCII file containing TRMM accumulated
precipitation pcp.
## "nc_file" is the out NetCDF file to be written.
## "-save" to call save.image().
##
## Details:
##
## Examples:
## Rscript trmm2nc.R \
## trmm_sample.txt \
## trmm_sample.nc
##
## Author:
## John Halley Gotway (johnhg at ucar.edu), NCAR-RAL-DTC
## 09/02/2011
##
########################################################################
########################################################################
#
# Required libraries.
#
########################################################################
library(ncdf)
########################################################################
#
# Constants and command line arguments
#
########################################################################
hdr_len = 5 # Each pcp file begins with 5 header lines.
hdr_file = "trmm.hdr" # Temporary header file.
missing = -9999 # Missing pcp value to be used in MET.
save = FALSE # If set to TRUE, call save.image()
########################################################################
#
# Handle the arguments.
#
########################################################################
# Retreive the arguments
args = commandArgs(TRUE)
# Check the number of arguments
if(length(args) < 2) {
cat("Usage: trmm2nc.R\n")
cat(" trmm_file\n")
cat(" nc_file\n")
cat(" [-save]\n")
cat(" where \"trmm_file\" is an ASCII file containing TRMM
accumulated precipitation pcp.\n")
cat(" \"nc_file\" is the out NetCDF file to be
written.\n")
cat(" \"-save\" to call save.image().\n\n")
quit()
}
# Store the input file names
trmm_file = args[1]
nc_file = args[2]
# Parse optional arguments
for(i in 1:length(args)) {
if(args[i] == "-save") {
save = TRUE
}
}
########################################################################
#
# Process the header pcp.
#
########################################################################
# Read the header pcp
sys_cmd = paste("head -", hdr_len, " ", trmm_file, " > ", hdr_file,
sep="")
system(sys_cmd)
hdr = read.table(hdr_file, fill=TRUE)
sys_cmd = paste("rm -f", hdr_file)
system(sys_cmd)
# Determine the accumlation interval
if(as.character(hdr[1,3]) == "TRMM") {
acc_hr = 24
} else if(as.character(hdr[1,3]) == "3-hourly") {
acc_hr = 3
} else {
cat("\n\nERROR: Can\'t figure out the accumulation interval!\n\n")
exit(1)
}
acc_sec = acc_hr * 3600
acc_str = sprintf("%.2i", acc_hr)
# Parse the valid time and compute the initialization time.
if(acc_hr == 24) {
valid = as.POSIXct(strptime(unlist(strsplit(as.character(hdr[3,4]),
'-'))[2],
format="%d%b%Y)", tz="GMT"))
} else {
valid = as.POSIXct(strptime(unlist(strsplit(as.character(hdr[3,4]),
'-'))[2],
format="%HZ%d%b%Y)", tz="GMT"))
}
init = valid - acc_sec
########################################################################
#
# Process the accumulated precipitation pcp.
#
########################################################################
# Read the TRMM file
cat(paste("Reading:", trmm_file, "\n"))
trmm = read.table(trmm_file, header=TRUE, skip=hdr_len)
# Determine lat/lon dimensions
lat = unique(trmm$Latitude)
lon = unique(trmm$Longitude)
dlat = lat[2]-lat[1]
dlon = lon[2]-lon[1]
nlat = length(lat)
nlon = length(lon)
# Slice the 1D array into a 2D array
pcp = array(trmm$AccRain, dim=c(nlon, nlat), dimnames=c("lon", "lat"))
pcp[pcp < missing] = missing
########################################################################
#
# Create the NetCDF output file.
#
########################################################################
# Define dimensions
dim_lat = dim.def.ncdf("lat", "degrees_north", lat,
create_dimvar=TRUE)
dim_lon = dim.def.ncdf("lon", "degrees_east", lon,
create_dimvar=TRUE)
# Define variables
var_pcp = var.def.ncdf(paste("APCP_", acc_str, sep=''), "kg/m^2",
list(dim_lon, dim_lat), missing,
longname="Total precipitation", prec="single")
# Define file
nc = create.ncdf(nc_file, var_pcp)
# Add variable attributes
att.put.ncdf(nc, var_pcp, "name", "APCP")
att.put.ncdf(nc, var_pcp, "level", paste("A", acc_hr, sep=''))
att.put.ncdf(nc, var_pcp, "grib_code", 61)
att.put.ncdf(nc, var_pcp, "_FillValue", missing, prec="single")
att.put.ncdf(nc, var_pcp, "init_time", format(init, "%Y%m%d_%H%M%S",
tz="GMT"), prec="text")
att.put.ncdf(nc, var_pcp, "init_time_ut", as.numeric(init),
prec="double")
att.put.ncdf(nc, var_pcp, "valid_time", format(valid, "%Y%m%d_%H%M%S",
tz="GMT"), prec="text")
att.put.ncdf(nc, var_pcp, "valid_time_ut", as.numeric(valid),
prec="double")
att.put.ncdf(nc, var_pcp, "accum_time", paste(acc_str, "0000",
sep=''))
att.put.ncdf(nc, var_pcp, "accum_time_sec", acc_sec)
# Add global attributes
cur_time = Sys.time()
att.put.ncdf(nc, 0, "FileOrigins", paste("File", nc_file, "generated",
format(Sys.time(), "%Y%m%d_%H%M%S"),
"on host", Sys.info()[4], "by
the Rscript trmm2nc.R"))
att.put.ncdf(nc, 0, "MET_version", "V3.0.1")
att.put.ncdf(nc, 0, "Projection", "LatLon", prec="text")
att.put.ncdf(nc, 0, "lat_ll", paste(min(lat), "degrees_north"),
prec="text")
att.put.ncdf(nc, 0, "lon_ll", paste(min(lon), "degrees_east"),
prec="text")
att.put.ncdf(nc, 0, "delta_lat", paste(dlat, "degrees"), prec="text")
att.put.ncdf(nc, 0, "delta_lon", paste(dlon, "degrees"), prec="text")
att.put.ncdf(nc, 0, "Nlat", paste(nlat, "grid_points"), prec="text")
att.put.ncdf(nc, 0, "Nlon", paste(nlon, "grid_points"), prec="text")
# Add pcp to the file
put.var.ncdf(nc, var_pcp, pcp)
# Close the file
close.ncdf(nc)
cat(paste("Writing:", nc_file, "\n"))
########################################################################
#
# Finished.
#
########################################################################
# Optionally, save all of the pcp to an .RData file
if(save == TRUE) save.image()
------------------------------------------------
Subject: Error in converting TRMM 3-hr rainfall into MET netcdf format
From: krishna kishore
Time: Wed Oct 12 14:58:08 2011
Hi John,
Thanks, It is working now.
I created 3-hr model accumulated rainfall.
When I am running the grid_stat, it is giving '
************************************************************
*** Running Grid-Stat using NETCDF input for both forecast and
observation
***
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744071610141774
Forecast File: WRFPRS_d01_24_27.nc
Observation File: TRMM_rain_20070430-12_15UTC.nc
Configuration File: config/GridStatConfig_APCP_3
*NetCDF: Attribute not found*
*****************************************************************
I checked both the NETCDF files. TRMM file (.25 degree resolution) is
in
lat-lon projection and model file (2 km resolution) is in Mercator.
Because of different resolutions, the number of points are different.
The dimensions of each file is attached with this mail.
Can you help me to bring them onto same grid.
Thank you very much,
Krishna.
On Thu, Oct 6, 2011 at 3:06 AM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Krishna,
>
> OK, I figured it out. The problem was that for some reason, my
version of
> the R ncdf package was out of sync with what's currently available.
I've
> modified the trmm2nc.R script and it should now
> work as expected. It required a very minor change.
>
> Please see attached, and let me know if you continue to experience
> problems.
>
> Thanks,
> John
>
>
> On 10/05/2011 12:13 PM, krishna kishore via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
> >
> > Thank you,
> > I am attaching the same file with this mail.
> > Please explain what is the error...
> > Krishna.
> >
> > On Wed, Oct 5, 2011 at 11:40 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> >> Krishna,
> >>
> >> It would probably be easiest for you just to send me the file
that's
> giving
> >> the error. Can you please send me this file:
> >> TRMM_rain_20070420-12_15UTC.txt
> >>
> >> I'll run it here and try to diagnose what's going on. You can
either
> send
> >> it via email or post it to our anonymous ftp site, as described
here:
> >> http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >> If you do post it to our ftp site, please write me back to tell
me when
> >> it's there.
> >>
> >> Hopefully we can get this working for you.
> >>
> >> Thanks,
> >> John
> >>
> >>
> >>
> >> On 10/05/2011 09:48 AM, krishna kishore via RT wrote:
> >>>
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
> >>>
> >>> Hi John,
> >>> Thank you very much..
> >>> Now I tried to convert the TRMM data....Still giving error..
> >>>
> >>
>
**************************************************************************
> >>> /home/osuri/Desktop/TRMM2NC_for_MET>Rscript trmm2nc.R
> >>> TRMM_rain_20070420-12_15UTC.txt TRMM_rain_20070420-12_15UTC.nc
> >>> Reading: TRMM_rain_20070420-12_15UTC.txt
> >>> Error in dim.def.ncdf("lat", "degrees_north", lat, create_dimvar
=
> TRUE,
> >> :
> >>> unused argument(s) (longname = "latitude")
> >>> Execution halted
> >>>
> >>
>
**************************************************************************
> >>> Can you check what is the problem?
> >>> Thanks,
> >>> Krishna.
> >>>
> >>>
> >>> On Wed, Oct 5, 2011 at 8:50 PM, John Halley Gotway via RT <
> >> met_help at ucar.edu
> >>>> wrote:
> >>>
> >>>> Krishna,
> >>>>
> >>>> Ah, I see the confusion. Those commands I sent should be run
from the
> >>>> command line - NOT FROM WITHIN R. Please try running the
commands
> >> again,
> >>>> but do NOT call R to start up an R session first.
> >>>>
> >>>> Hopefully, that'll do it.
> >>>>
> >>>> Thanks,
> >>>> John
> >>>>
> >>>> On 10/04/2011 04:23 PM, krishna kishore via RT wrote:
> >>>>>
> >>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784
>
> >>>>>
> >>>>> Hi John,
> >>>>> I Installed the ncdf library in R.
> >>>>> For your more information....I Installed R software in my
ubuntu
> linux
> >>>>> machine (through sudo apt-get install r-base) ....and
installed the
> >> ncdf
> >>>>> library also.
> >>>>> Then I tried with the commands
> >>>>> R for R prompt and then try to run the trmm2nc.R script
> >>>>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
> >>>>> TRMM_rain_20070420-9_12UTC.nc
> >>>>> Error: unexpected symbol in "Rscript trmm2nc.R"
> >>>>>
> >>>>>
> >>>>>> R --version
> >>>>> Error: object 'R' not found
> >>>>>
> >>>>>> which Rscript
> >>>>> Error: unexpected symbol in "which Rscript"
> >>>>>
> >>>>>> library(ncdf)
> >>>>> nothing is displaying....means installed the library
correctly..
> >>>>>
> >>>>> But why the above errors are giving when i am executin Rscript
and R
> >>>>> --version commands.
> >>>>>
> >>>>> Thats means, R software is not installed in my account
properly?
> >>>>> What I should do to pass through this problem..
> >>>>>
> >>>>> Thanks,
> >>>>> Krishna.
> >>>>>
> >>>>>
> >>>>> On Mon, Sep 26, 2011 at 11:50 PM, RAL HelpDesk {for John
Halley
> Gotway}
> >> <
> >>>>> met_help at ucar.edu> wrote:
> >>>>>
> >>>>>> Krishna,
> >>>>>>
> >>>>>> I used the trmm2nc.R script to convert the sample data file
you
> sent,
> >>>> and
> >>>>>> everything worked fine.
> >>>>>>
> >>>>>> So there is likely a problem in the environment somewhere. I
have
> >> never
> >>>>>> seen that sort of error message from R before. Let's try
looking
> how
> >> R
> >>>> is
> >>>>>> installed on your system:
> >>>>>>
> >>>>>> Check what version of R you're running:
> >>>>>> R --version
> >>>>>> Check that you have the Rscript command in your path:
> >>>>>> which Rscript
> >>>>>> Try running the Rscript command:
> >>>>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
> >>>>>> TRMM_rain_20070420-9_12UTC.nc
> >>>>>>
> >>>>>> Please note that this script expects that the R "ncdf"
package is
> >>>> installed
> >>>>>> on your system. To check if it's there, just try executing
the
> >>>> following
> >>>>>> commands:
> >>>>>> R
> >>>>>> library(ncdf)
> >>>>>>
> >>>>>> If it loads without error, everything is fine. If you get an
error,
> >>>> you'll
> >>>>>> need to install that "ncdf" package. To do so, execute the
> following
> >>>>>> commands:
> >>>>>> R
> >>>>>> install.packages("ncdf",
> >>>>>> configure.args="--with-
netcdf_include=/usr/local/netcdf/include
> >>>>>> -with-netcdf_lib=/usr/local/netcdf/lib")
> >>>>>>
> >>>>>> Please let me know what you learn from these steps.
> >>>>>>
> >>>>>> Thanks,
> >>>>>> John Halley Gotway
> >>>>>> met_help at ucar.edu
> >>>>>>
> >>>>>> On 09/23/2011 08:01 PM, RAL HelpDesk {for krishna kishore}
wrote:
> >>>>>>>
> >>>>>>> Fri Sep 23 20:01:29 2011: Request 49784 was acted upon.
> >>>>>>> Transaction: Ticket created by osurikishore at gmail.com
> >>>>>>> Queue: met_help
> >>>>>>> Subject: Error in converting TRMM 3-hr rainfall into
MET
> netcdf
> >>>>>> format
> >>>>>>> Owner: Nobody
> >>>>>>> Requestors: osurikishore at gmail.com
> >>>>>>> Status: new
> >>>>>>> Ticket <URL:
> >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784>
> >>>>>>>
> >>>>>>>
> >>>>>>> Dear Methelp,
> >>>>>>> I tried to convert the TRMM 3-hrly accumulated data (from
> >>>>>>>
http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml)
> >> using
> >>>>>>> trmm2nc.R file.
> >>>>>>> I copied the ASCII rainfall data and attached for your
reference.
> >>>>>>> When i am running, it is giving error that
> >>>>>>> *Error: unexpected symbol in "Rscript trmm2nc.R"*
> >>>>>>>
> >>>>>>> Can you help me why it is like that.
> >>>>>>> Is this with NETCDF library?
> >>>>>>> For your information, I have netcdf library and already
export that
> >>>> path
> >>>>>> in
> >>>>>>> my .bashrc file.
> >>>>>>> Will the Rscript automatically takes or should I mentioned
> somewhere
> >>>> when
> >>>>>>> using this Rscript.
> >>>>>>>
> >>>>>>> Can you suggest what went wrong?
> >>>>>>> Thanks in advance,
> >>>>>>> Krishna.
> >>>>>>>
> >>>>>>>
> >>>>>>
> >>>>>>
> >>>>>
> >>>>>
> >>>>
> >>>>
> >>>
> >>>
> >>
> >>
> >
> >
>
>
>
########################################################################
> ##
> ## Name: trmm2nc.R
> ##
> ## Description:
> ## This script should be called via an Rscript command to
reformat
> ## ASCII TRMM precipitation pcp into a NetCDF format that MET
> ## can read. The input pcp for this script can be attained
via:
> ## 24-hour accumulations:
> ##
>
http://disc2.nascom.nasa.gov/Giovanni/tovas/realtime.3B42RT_daily.shtml
> ## 03-hour accumulations:
> ##
> http://disc2.nascom.nasa.gov/Giovanni/tovas/realtime.3B42RT.shtml
> ##
> ## Usage:
> ## Rscript trmm2nc.R
> ## trmm_file
> ## nc_file
> ## [-save]
> ##
> ## Arguments:
> ## "trmm_file" is an ASCII file containing TRMM accumulated
> precipitation pcp.
> ## "nc_file" is the out NetCDF file to be written.
> ## "-save" to call save.image().
> ##
> ## Details:
> ##
> ## Examples:
> ## Rscript trmm2nc.R \
> ## trmm_sample.txt \
> ## trmm_sample.nc
> ##
> ## Author:
> ## John Halley Gotway (johnhg at ucar.edu), NCAR-RAL-DTC
> ## 09/02/2011
> ##
>
########################################################################
>
>
########################################################################
> #
> # Required libraries.
> #
>
########################################################################
>
> library(ncdf)
>
>
########################################################################
> #
> # Constants and command line arguments
> #
>
########################################################################
>
> hdr_len = 5 # Each pcp file begins with 5 header lines.
> hdr_file = "trmm.hdr" # Temporary header file.
> missing = -9999 # Missing pcp value to be used in MET.
> save = FALSE # If set to TRUE, call save.image()
>
>
########################################################################
> #
> # Handle the arguments.
> #
>
########################################################################
>
> # Retreive the arguments
> args = commandArgs(TRUE)
>
> # Check the number of arguments
> if(length(args) < 2) {
> cat("Usage: trmm2nc.R\n")
> cat(" trmm_file\n")
> cat(" nc_file\n")
> cat(" [-save]\n")
> cat(" where \"trmm_file\" is an ASCII file containing TRMM
> accumulated precipitation pcp.\n")
> cat(" \"nc_file\" is the out NetCDF file to be
> written.\n")
> cat(" \"-save\" to call save.image().\n\n")
> quit()
> }
>
> # Store the input file names
> trmm_file = args[1]
> nc_file = args[2]
>
> # Parse optional arguments
> for(i in 1:length(args)) {
> if(args[i] == "-save") {
> save = TRUE
> }
> }
>
>
########################################################################
> #
> # Process the header pcp.
> #
>
########################################################################
>
> # Read the header pcp
> sys_cmd = paste("head -", hdr_len, " ", trmm_file, " > ", hdr_file,
sep="")
> system(sys_cmd)
> hdr = read.table(hdr_file, fill=TRUE)
> sys_cmd = paste("rm -f", hdr_file)
> system(sys_cmd)
>
> # Determine the accumlation interval
> if(as.character(hdr[1,3]) == "TRMM") {
> acc_hr = 24
> } else if(as.character(hdr[1,3]) == "3-hourly") {
> acc_hr = 3
> } else {
> cat("\n\nERROR: Can\'t figure out the accumulation interval!\n\n")
> exit(1)
> }
> acc_sec = acc_hr * 3600
> acc_str = sprintf("%.2i", acc_hr)
>
> # Parse the valid time and compute the initialization time.
> if(acc_hr == 24) {
> valid = as.POSIXct(strptime(unlist(strsplit(as.character(hdr[3,4]),
> '-'))[2],
> format="%d%b%Y)", tz="GMT"))
> } else {
> valid = as.POSIXct(strptime(unlist(strsplit(as.character(hdr[3,4]),
> '-'))[2],
> format="%HZ%d%b%Y)", tz="GMT"))
> }
> init = valid - acc_sec
>
>
########################################################################
> #
> # Process the accumulated precipitation pcp.
> #
>
########################################################################
>
> # Read the TRMM file
> cat(paste("Reading:", trmm_file, "\n"))
> trmm = read.table(trmm_file, header=TRUE, skip=hdr_len)
>
> # Determine lat/lon dimensions
> lat = unique(trmm$Latitude)
> lon = unique(trmm$Longitude)
> dlat = lat[2]-lat[1]
> dlon = lon[2]-lon[1]
> nlat = length(lat)
> nlon = length(lon)
>
> # Slice the 1D array into a 2D array
> pcp = array(trmm$AccRain, dim=c(nlon, nlat), dimnames=c("lon",
"lat"))
> pcp[pcp < missing] = missing
>
>
########################################################################
> #
> # Create the NetCDF output file.
> #
>
########################################################################
>
> # Define dimensions
> dim_lat = dim.def.ncdf("lat", "degrees_north", lat,
> create_dimvar=TRUE)
> dim_lon = dim.def.ncdf("lon", "degrees_east", lon,
> create_dimvar=TRUE)
>
> # Define variables
> var_pcp = var.def.ncdf(paste("APCP_", acc_str, sep=''), "kg/m^2",
> list(dim_lon, dim_lat), missing,
> longname="Total precipitation", prec="single")
>
> # Define file
> nc = create.ncdf(nc_file, var_pcp)
>
> # Add variable attributes
> att.put.ncdf(nc, var_pcp, "name", "APCP")
> att.put.ncdf(nc, var_pcp, "level", paste("A", acc_hr, sep=''))
> att.put.ncdf(nc, var_pcp, "grib_code", 61)
> att.put.ncdf(nc, var_pcp, "_FillValue", missing, prec="single")
> att.put.ncdf(nc, var_pcp, "init_time", format(init, "%Y%m%d_%H%M%S",
> tz="GMT"), prec="text")
> att.put.ncdf(nc, var_pcp, "init_time_ut", as.numeric(init),
prec="double")
> att.put.ncdf(nc, var_pcp, "valid_time", format(valid,
"%Y%m%d_%H%M%S",
> tz="GMT"), prec="text")
> att.put.ncdf(nc, var_pcp, "valid_time_ut", as.numeric(valid),
> prec="double")
> att.put.ncdf(nc, var_pcp, "accum_time", paste(acc_str, "0000",
sep=''))
> att.put.ncdf(nc, var_pcp, "accum_time_sec", acc_sec)
>
> # Add global attributes
> cur_time = Sys.time()
> att.put.ncdf(nc, 0, "FileOrigins", paste("File", nc_file,
"generated",
> format(Sys.time(), "%Y%m%d_%H%M%S"),
> "on host", Sys.info()[4],
"by the
> Rscript trmm2nc.R"))
> att.put.ncdf(nc, 0, "MET_version", "V3.0.1")
> att.put.ncdf(nc, 0, "Projection", "LatLon", prec="text")
> att.put.ncdf(nc, 0, "lat_ll", paste(min(lat), "degrees_north"),
> prec="text")
> att.put.ncdf(nc, 0, "lon_ll", paste(min(lon), "degrees_east"),
prec="text")
> att.put.ncdf(nc, 0, "delta_lat", paste(dlat, "degrees"),
prec="text")
> att.put.ncdf(nc, 0, "delta_lon", paste(dlon, "degrees"),
prec="text")
> att.put.ncdf(nc, 0, "Nlat", paste(nlat, "grid_points"), prec="text")
> att.put.ncdf(nc, 0, "Nlon", paste(nlon, "grid_points"), prec="text")
>
> # Add pcp to the file
> put.var.ncdf(nc, var_pcp, pcp)
>
> # Close the file
> close.ncdf(nc)
>
> cat(paste("Writing:", nc_file, "\n"))
>
>
########################################################################
> #
> # Finished.
> #
>
########################################################################
>
> # Optionally, save all of the pcp to an .RData file
> if(save == TRUE) save.image()
>
>
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
------------------------------------------------
Subject: Error in converting TRMM 3-hr rainfall into MET netcdf format
From: krishna kishore
Time: Wed Oct 12 14:58:08 2011
/home/ba01/u135/kosuri/met/met3/scripts>ncdump -h WRFPRS_IC-
2007042912_FCST24_27.nc
netcdf WRFPRS_IC-2007042912_FCST24_27 {
dimensions:
lat = 599 ;
lon = 699 ;
variables:
float lat(lat, lon) ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
lat:standard_name = "latitude" ;
float lon(lat, lon) ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
lon:standard_name = "longitude" ;
float APCP_03(lat, lon) ;
APCP_03:name = "APCP" ;
APCP_03:long_name = "Total precipitation" ;
APCP_03:level = "A3" ;
APCP_03:units = "kg/m^2" ;
APCP_03:grib_code = 61 ;
APCP_03:_FillValue = -9999.f ;
APCP_03:init_time = "20070429_120000" ;
APCP_03:init_time_ut = 1177848000 ;
APCP_03:valid_time = "20070430_150000" ;
APCP_03:valid_time_ut = 1177945200 ;
APCP_03:accum_time = "030000" ;
APCP_03:accum_time_sec = 10800 ;
// global attributes:
:FileOrigins = "File WRFPRS_d01_24_27.nc generated 20111005_222108
UTC on host coates-fe03.rcac.purdue.edu by the MET pcp_combine tool" ;
:MET_version = "V3.0" ;
:MET_tool = "pcp_combine" ;
:RunCommand = "Subtraction:
/scratch/lustreA/k/kosuri/DOMAINS/CZIL/2007Apr29/cntl_def/postprd/WRFPRS_d01.027
with accumulation of 270000 minus
/scratch/lustreA/k/kosuri/DOMAINS/CZIL/2007Apr29/cntl_def/postprd/WRFPRS_d01.024
with accumulation of 240000." ;
:Projection = "Mercator" ;
:lat_ll = "17.963000 degrees_north" ;
:lon_ll = "80.722000 degrees_east" ;
:lat_ur = "27.856000 degrees_north" ;
:lon_ur = "93.278000 degrees_east" ;
:nx = "699" ;
:ny = "599" ;
}
------------------------------------------------
Subject: Error in converting TRMM 3-hr rainfall into MET netcdf format
From: krishna kishore
Time: Wed Oct 12 14:58:08 2011
/home/ba01/u135/kosuri/met/met3/scripts>ncdump -h TRMM_rain_20070430-
12_15UTC.nc |more
netcdf TRMM_rain_20070430-12_15UTC {
dimensions:
lon = 49 ;
lat = 41 ;
variables:
double lon(lon) ;
lon:units = "degrees_east" ;
double lat(lat) ;
lat:units = "degrees_north" ;
float APCP_03(lat, lon) ;
APCP_03:units = "kg/m^2" ;
APCP_03:missing_value = -9999.f ;
APCP_03:long_name = "Total precipitation" ;
APCP_03:name = "APCP" ;
APCP_03:level = "A3" ;
APCP_03:grib_code = 61.f ;
APCP_03:_FillValue = -9999.f ;
APCP_03:init_time = "20070430_120000" ;
APCP_03:init_time_ut = 1177934400. ;
APCP_03:valid_time = "20070430_150000" ;
APCP_03:valid_time_ut = 1177945200. ;
APCP_03:accum_time = "030000" ;
APCP_03:accum_time_sec = 10800.f ;
// global attributes:
:FileOrigins = "File TRMM_rain_20070430-12_15UTC.nc generated
20111005_183646 on host coates-fe02.rcac.purdue.edu by the Rs
cript trmm2nc.R" ;
:MET_version = "V3.0.1" ;
:Projection = "LatLon" ;
:lat_ll = "18 degrees_north" ;
:lon_ll = "81 degrees_east" ;
:delta_lat = "0.25 degrees" ;
:delta_lon = "0.25 degrees" ;
:Nlat = "41 grid_points" ;
:Nlon = "49 grid_points" ;
}
------------------------------------------------
Subject: Re: Fwd: [rt.rap.ucar.edu #49784] Error in converting TRMM 3-hr rainfall into MET netcdf format
From: John Halley Gotway
Time: Wed Oct 12 15:30:54 2011
Krishna,
My suggestions are listed below:
(1) In your Grid-Stat config file, you need to define the masking
region over which verification scores should be computed. Please try
using the following:
mask_grid[] = [ "FULL" ];
mask_poly[] = [];
Specifying "FULL" for the mask_grid tells Grid-Stat to summarize
statistics over the entire domain.
You had "" in the mask_poly setting, but that will cause an error.
Just leave the array empty if you don't have any polyline regions to
apply.
(2) I see why you're getting the following error:
*NetCDF: Attribute not found*
This is due to some logic in parsing the grid specification. Here's
what it comes down to - the value of the "MET_version" attribute
listed in the NetCDF header info must match the version of the
code you're running. In your case, it appears that you're running
METv3.0 but the NetCDF global attribute is this:
:MET_version = "V3.0.1" ;
That silly little difference is causing the problem. The quick fix is
to modify the TRMM2NC.R script and replace "V3.0.1" with "V3.0".
Sorry for this problem - longer term, we'll try to figure out a
better solution.
(3) Once you get past the problem listed above, you'll get this error:
ERROR: process_scores() -> The forecast and observation grids do
not match for field 1.
As you mention, you need to get the model data onto the same grid as
the observations. Assuming that you have GRIB output from WPP or UPP,
you can regrid it using the copygb tool. Listed here are
examples of how you can set up copygb to regrid your data:
http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/index.php
In your case, you want to regrid to a LatLon grid, which is described
here:
http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/run2.php
Here is the grid definition you're trying to get to:
:Projection = "LatLon" ;
:lat_ll = "18 degrees_north" ;
:lon_ll = "81 degrees_east" ;
:delta_lat = "0.25 degrees" ;
:delta_lon = "0.25 degrees" ;
:Nlat = "41 grid_points" ;
:Nlon = "49 grid_points" ;
Following the website instructions, that gets translated into a grid
definition for copygb as:
"255 0 49 41 18000 81000 128 28250 93250 250 250 64"
So for each GRIB output file you have, run it through copygb like
this:
copygb -xg"255 0 49 41 18000 81000 128 28250 93250 250 250 64"
WRFPRS_in.grib WRFPRS_regrid.grib
Once it's all on the same grid as the observations, you can run it
through pcp_combine - or whatever you need to do to get the correct
accumulation interval.
Hope that helps.
John Halley Gotway
On 10/12/2011 02:58 PM, krishna kishore via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>
> Hi John,
> Thanks, It is working now.
> I created 3-hr model accumulated rainfall.
> When I am running the grid_stat, it is giving '
> ************************************************************
> *** Running Grid-Stat using NETCDF input for both forecast and
observation
> ***
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744071610141774
> Forecast File: WRFPRS_d01_24_27.nc
> Observation File: TRMM_rain_20070430-12_15UTC.nc
> Configuration File: config/GridStatConfig_APCP_3
> *NetCDF: Attribute not found*
> *****************************************************************
>
> I checked both the NETCDF files. TRMM file (.25 degree resolution)
is in
> lat-lon projection and model file (2 km resolution) is in Mercator.
> Because of different resolutions, the number of points are
different.
>
> The dimensions of each file is attached with this mail.
>
> Can you help me to bring them onto same grid.
>
> Thank you very much,
> Krishna.
>
>
>
>
>
>
> On Thu, Oct 6, 2011 at 3:06 AM, John Halley Gotway via RT
<met_help at ucar.edu
>> wrote:
>
>> Krishna,
>>
>> OK, I figured it out. The problem was that for some reason, my
version of
>> the R ncdf package was out of sync with what's currently available.
I've
>> modified the trmm2nc.R script and it should now
>> work as expected. It required a very minor change.
>>
>> Please see attached, and let me know if you continue to experience
>> problems.
>>
>> Thanks,
>> John
>>
>>
>> On 10/05/2011 12:13 PM, krishna kishore via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>>>
>>> Thank you,
>>> I am attaching the same file with this mail.
>>> Please explain what is the error...
>>> Krishna.
>>>
>>> On Wed, Oct 5, 2011 at 11:40 PM, John Halley Gotway via RT <
>>> met_help at ucar.edu> wrote:
>>>
>>>> Krishna,
>>>>
>>>> It would probably be easiest for you just to send me the file
that's
>> giving
>>>> the error. Can you please send me this file:
>>>> TRMM_rain_20070420-12_15UTC.txt
>>>>
>>>> I'll run it here and try to diagnose what's going on. You can
either
>> send
>>>> it via email or post it to our anonymous ftp site, as described
here:
>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>>> If you do post it to our ftp site, please write me back to tell
me when
>>>> it's there.
>>>>
>>>> Hopefully we can get this working for you.
>>>>
>>>> Thanks,
>>>> John
>>>>
>>>>
>>>>
>>>> On 10/05/2011 09:48 AM, krishna kishore via RT wrote:
>>>>>
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784 >
>>>>>
>>>>> Hi John,
>>>>> Thank you very much..
>>>>> Now I tried to convert the TRMM data....Still giving error..
>>>>>
>>>>
>>
**************************************************************************
>>>>> /home/osuri/Desktop/TRMM2NC_for_MET>Rscript trmm2nc.R
>>>>> TRMM_rain_20070420-12_15UTC.txt TRMM_rain_20070420-12_15UTC.nc
>>>>> Reading: TRMM_rain_20070420-12_15UTC.txt
>>>>> Error in dim.def.ncdf("lat", "degrees_north", lat, create_dimvar
=
>> TRUE,
>>>> :
>>>>> unused argument(s) (longname = "latitude")
>>>>> Execution halted
>>>>>
>>>>
>>
**************************************************************************
>>>>> Can you check what is the problem?
>>>>> Thanks,
>>>>> Krishna.
>>>>>
>>>>>
>>>>> On Wed, Oct 5, 2011 at 8:50 PM, John Halley Gotway via RT <
>>>> met_help at ucar.edu
>>>>>> wrote:
>>>>>
>>>>>> Krishna,
>>>>>>
>>>>>> Ah, I see the confusion. Those commands I sent should be run
from the
>>>>>> command line - NOT FROM WITHIN R. Please try running the
commands
>>>> again,
>>>>>> but do NOT call R to start up an R session first.
>>>>>>
>>>>>> Hopefully, that'll do it.
>>>>>>
>>>>>> Thanks,
>>>>>> John
>>>>>>
>>>>>> On 10/04/2011 04:23 PM, krishna kishore via RT wrote:
>>>>>>>
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784
>
>>>>>>>
>>>>>>> Hi John,
>>>>>>> I Installed the ncdf library in R.
>>>>>>> For your more information....I Installed R software in my
ubuntu
>> linux
>>>>>>> machine (through sudo apt-get install r-base) ....and
installed the
>>>> ncdf
>>>>>>> library also.
>>>>>>> Then I tried with the commands
>>>>>>> R for R prompt and then try to run the trmm2nc.R script
>>>>>>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
>>>>>>> TRMM_rain_20070420-9_12UTC.nc
>>>>>>> Error: unexpected symbol in "Rscript trmm2nc.R"
>>>>>>>
>>>>>>>
>>>>>>>> R --version
>>>>>>> Error: object 'R' not found
>>>>>>>
>>>>>>>> which Rscript
>>>>>>> Error: unexpected symbol in "which Rscript"
>>>>>>>
>>>>>>>> library(ncdf)
>>>>>>> nothing is displaying....means installed the library
correctly..
>>>>>>>
>>>>>>> But why the above errors are giving when i am executin Rscript
and R
>>>>>>> --version commands.
>>>>>>>
>>>>>>> Thats means, R software is not installed in my account
properly?
>>>>>>> What I should do to pass through this problem..
>>>>>>>
>>>>>>> Thanks,
>>>>>>> Krishna.
>>>>>>>
>>>>>>>
>>>>>>> On Mon, Sep 26, 2011 at 11:50 PM, RAL HelpDesk {for John
Halley
>> Gotway}
>>>> <
>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>
>>>>>>>> Krishna,
>>>>>>>>
>>>>>>>> I used the trmm2nc.R script to convert the sample data file
you
>> sent,
>>>>>> and
>>>>>>>> everything worked fine.
>>>>>>>>
>>>>>>>> So there is likely a problem in the environment somewhere. I
have
>>>> never
>>>>>>>> seen that sort of error message from R before. Let's try
looking
>> how
>>>> R
>>>>>> is
>>>>>>>> installed on your system:
>>>>>>>>
>>>>>>>> Check what version of R you're running:
>>>>>>>> R --version
>>>>>>>> Check that you have the Rscript command in your path:
>>>>>>>> which Rscript
>>>>>>>> Try running the Rscript command:
>>>>>>>> Rscript trmm2nc.R TRMM_rain_20070420-9_12UTC.txt
>>>>>>>> TRMM_rain_20070420-9_12UTC.nc
>>>>>>>>
>>>>>>>> Please note that this script expects that the R "ncdf"
package is
>>>>>> installed
>>>>>>>> on your system. To check if it's there, just try executing
the
>>>>>> following
>>>>>>>> commands:
>>>>>>>> R
>>>>>>>> library(ncdf)
>>>>>>>>
>>>>>>>> If it loads without error, everything is fine. If you get an
error,
>>>>>> you'll
>>>>>>>> need to install that "ncdf" package. To do so, execute the
>> following
>>>>>>>> commands:
>>>>>>>> R
>>>>>>>> install.packages("ncdf",
>>>>>>>> configure.args="--with-
netcdf_include=/usr/local/netcdf/include
>>>>>>>> -with-netcdf_lib=/usr/local/netcdf/lib")
>>>>>>>>
>>>>>>>> Please let me know what you learn from these steps.
>>>>>>>>
>>>>>>>> Thanks,
>>>>>>>> John Halley Gotway
>>>>>>>> met_help at ucar.edu
>>>>>>>>
>>>>>>>> On 09/23/2011 08:01 PM, RAL HelpDesk {for krishna kishore}
wrote:
>>>>>>>>>
>>>>>>>>> Fri Sep 23 20:01:29 2011: Request 49784 was acted upon.
>>>>>>>>> Transaction: Ticket created by osurikishore at gmail.com
>>>>>>>>> Queue: met_help
>>>>>>>>> Subject: Error in converting TRMM 3-hr rainfall into
MET
>> netcdf
>>>>>>>> format
>>>>>>>>> Owner: Nobody
>>>>>>>>> Requestors: osurikishore at gmail.com
>>>>>>>>> Status: new
>>>>>>>>> Ticket <URL:
>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49784>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Dear Methelp,
>>>>>>>>> I tried to convert the TRMM 3-hrly accumulated data (from
>>>>>>>>>
http://disc2.nascom.nasa.gov/Giovanni/tovas/TRMM_V6.3B42.shtml)
>>>> using
>>>>>>>>> trmm2nc.R file.
>>>>>>>>> I copied the ASCII rainfall data and attached for your
reference.
>>>>>>>>> When i am running, it is giving error that
>>>>>>>>> *Error: unexpected symbol in "Rscript trmm2nc.R"*
>>>>>>>>>
>>>>>>>>> Can you help me why it is like that.
>>>>>>>>> Is this with NETCDF library?
>>>>>>>>> For your information, I have netcdf library and already
export that
>>>>>> path
>>>>>>>> in
>>>>>>>>> my .bashrc file.
>>>>>>>>> Will the Rscript automatically takes or should I mentioned
>> somewhere
>>>>>> when
>>>>>>>>> using this Rscript.
>>>>>>>>>
>>>>>>>>> Can you suggest what went wrong?
>>>>>>>>> Thanks in advance,
>>>>>>>>> Krishna.
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>>
########################################################################
>> ##
>> ## Name: trmm2nc.R
>> ##
>> ## Description:
>> ## This script should be called via an Rscript command to
reformat
>> ## ASCII TRMM precipitation pcp into a NetCDF format that MET
>> ## can read. The input pcp for this script can be attained
via:
>> ## 24-hour accumulations:
>> ##
>>
http://disc2.nascom.nasa.gov/Giovanni/tovas/realtime.3B42RT_daily.shtml
>> ## 03-hour accumulations:
>> ##
>> http://disc2.nascom.nasa.gov/Giovanni/tovas/realtime.3B42RT.shtml
>> ##
>> ## Usage:
>> ## Rscript trmm2nc.R
>> ## trmm_file
>> ## nc_file
>> ## [-save]
>> ##
>> ## Arguments:
>> ## "trmm_file" is an ASCII file containing TRMM accumulated
>> precipitation pcp.
>> ## "nc_file" is the out NetCDF file to be written.
>> ## "-save" to call save.image().
>> ##
>> ## Details:
>> ##
>> ## Examples:
>> ## Rscript trmm2nc.R \
>> ## trmm_sample.txt \
>> ## trmm_sample.nc
>> ##
>> ## Author:
>> ## John Halley Gotway (johnhg at ucar.edu), NCAR-RAL-DTC
>> ## 09/02/2011
>> ##
>>
########################################################################
>>
>>
########################################################################
>> #
>> # Required libraries.
>> #
>>
########################################################################
>>
>> library(ncdf)
>>
>>
########################################################################
>> #
>> # Constants and command line arguments
>> #
>>
########################################################################
>>
>> hdr_len = 5 # Each pcp file begins with 5 header lines.
>> hdr_file = "trmm.hdr" # Temporary header file.
>> missing = -9999 # Missing pcp value to be used in MET.
>> save = FALSE # If set to TRUE, call save.image()
>>
>>
########################################################################
>> #
>> # Handle the arguments.
>> #
>>
########################################################################
>>
>> # Retreive the arguments
>> args = commandArgs(TRUE)
>>
>> # Check the number of arguments
>> if(length(args) < 2) {
>> cat("Usage: trmm2nc.R\n")
>> cat(" trmm_file\n")
>> cat(" nc_file\n")
>> cat(" [-save]\n")
>> cat(" where \"trmm_file\" is an ASCII file containing TRMM
>> accumulated precipitation pcp.\n")
>> cat(" \"nc_file\" is the out NetCDF file to be
>> written.\n")
>> cat(" \"-save\" to call save.image().\n\n")
>> quit()
>> }
>>
>> # Store the input file names
>> trmm_file = args[1]
>> nc_file = args[2]
>>
>> # Parse optional arguments
>> for(i in 1:length(args)) {
>> if(args[i] == "-save") {
>> save = TRUE
>> }
>> }
>>
>>
########################################################################
>> #
>> # Process the header pcp.
>> #
>>
########################################################################
>>
>> # Read the header pcp
>> sys_cmd = paste("head -", hdr_len, " ", trmm_file, " > ", hdr_file,
sep="")
>> system(sys_cmd)
>> hdr = read.table(hdr_file, fill=TRUE)
>> sys_cmd = paste("rm -f", hdr_file)
>> system(sys_cmd)
>>
>> # Determine the accumlation interval
>> if(as.character(hdr[1,3]) == "TRMM") {
>> acc_hr = 24
>> } else if(as.character(hdr[1,3]) == "3-hourly") {
>> acc_hr = 3
>> } else {
>> cat("\n\nERROR: Can\'t figure out the accumulation interval!\n\n")
>> exit(1)
>> }
>> acc_sec = acc_hr * 3600
>> acc_str = sprintf("%.2i", acc_hr)
>>
>> # Parse the valid time and compute the initialization time.
>> if(acc_hr == 24) {
>> valid =
as.POSIXct(strptime(unlist(strsplit(as.character(hdr[3,4]),
>> '-'))[2],
>> format="%d%b%Y)", tz="GMT"))
>> } else {
>> valid =
as.POSIXct(strptime(unlist(strsplit(as.character(hdr[3,4]),
>> '-'))[2],
>> format="%HZ%d%b%Y)", tz="GMT"))
>> }
>> init = valid - acc_sec
>>
>>
########################################################################
>> #
>> # Process the accumulated precipitation pcp.
>> #
>>
########################################################################
>>
>> # Read the TRMM file
>> cat(paste("Reading:", trmm_file, "\n"))
>> trmm = read.table(trmm_file, header=TRUE, skip=hdr_len)
>>
>> # Determine lat/lon dimensions
>> lat = unique(trmm$Latitude)
>> lon = unique(trmm$Longitude)
>> dlat = lat[2]-lat[1]
>> dlon = lon[2]-lon[1]
>> nlat = length(lat)
>> nlon = length(lon)
>>
>> # Slice the 1D array into a 2D array
>> pcp = array(trmm$AccRain, dim=c(nlon, nlat), dimnames=c("lon",
"lat"))
>> pcp[pcp < missing] = missing
>>
>>
########################################################################
>> #
>> # Create the NetCDF output file.
>> #
>>
########################################################################
>>
>> # Define dimensions
>> dim_lat = dim.def.ncdf("lat", "degrees_north", lat,
>> create_dimvar=TRUE)
>> dim_lon = dim.def.ncdf("lon", "degrees_east", lon,
>> create_dimvar=TRUE)
>>
>> # Define variables
>> var_pcp = var.def.ncdf(paste("APCP_", acc_str, sep=''), "kg/m^2",
>> list(dim_lon, dim_lat), missing,
>> longname="Total precipitation",
prec="single")
>>
>> # Define file
>> nc = create.ncdf(nc_file, var_pcp)
>>
>> # Add variable attributes
>> att.put.ncdf(nc, var_pcp, "name", "APCP")
>> att.put.ncdf(nc, var_pcp, "level", paste("A", acc_hr, sep=''))
>> att.put.ncdf(nc, var_pcp, "grib_code", 61)
>> att.put.ncdf(nc, var_pcp, "_FillValue", missing, prec="single")
>> att.put.ncdf(nc, var_pcp, "init_time", format(init,
"%Y%m%d_%H%M%S",
>> tz="GMT"), prec="text")
>> att.put.ncdf(nc, var_pcp, "init_time_ut", as.numeric(init),
prec="double")
>> att.put.ncdf(nc, var_pcp, "valid_time", format(valid,
"%Y%m%d_%H%M%S",
>> tz="GMT"), prec="text")
>> att.put.ncdf(nc, var_pcp, "valid_time_ut", as.numeric(valid),
>> prec="double")
>> att.put.ncdf(nc, var_pcp, "accum_time", paste(acc_str, "0000",
sep=''))
>> att.put.ncdf(nc, var_pcp, "accum_time_sec", acc_sec)
>>
>> # Add global attributes
>> cur_time = Sys.time()
>> att.put.ncdf(nc, 0, "FileOrigins", paste("File", nc_file,
"generated",
>> format(Sys.time(), "%Y%m%d_%H%M%S"),
>> "on host", Sys.info()[4],
"by the
>> Rscript trmm2nc.R"))
>> att.put.ncdf(nc, 0, "MET_version", "V3.0.1")
>> att.put.ncdf(nc, 0, "Projection", "LatLon", prec="text")
>> att.put.ncdf(nc, 0, "lat_ll", paste(min(lat), "degrees_north"),
>> prec="text")
>> att.put.ncdf(nc, 0, "lon_ll", paste(min(lon), "degrees_east"),
prec="text")
>> att.put.ncdf(nc, 0, "delta_lat", paste(dlat, "degrees"),
prec="text")
>> att.put.ncdf(nc, 0, "delta_lon", paste(dlon, "degrees"),
prec="text")
>> att.put.ncdf(nc, 0, "Nlat", paste(nlat, "grid_points"),
prec="text")
>> att.put.ncdf(nc, 0, "Nlon", paste(nlon, "grid_points"),
prec="text")
>>
>> # Add pcp to the file
>> put.var.ncdf(nc, var_pcp, pcp)
>>
>> # Close the file
>> close.ncdf(nc)
>>
>> cat(paste("Writing:", nc_file, "\n"))
>>
>>
########################################################################
>> #
>> # Finished.
>> #
>>
########################################################################
>>
>> # Optionally, save all of the pcp to an .RData file
>> if(save == TRUE) save.image()
>>
>>
>
>
------------------------------------------------
More information about the Met_help
mailing list