Lats4d with RAMS output

[Arlindo da Silva]

On Sat, Nov 15, 2008 at 4:17 PM, Advait Godbole <agodbole at purdue.edu> wrote:

> Hi Arlindo,
>
> I am working on visualizing the CO2 concentrations over the US, using the
> Visualization Toolkit. I have a gridded binary output file which came from a
> RAMS model run, in the form of a " .gra " file. I can read it in through
> GrADS and also obtain plots, but to use VTK,


Have you considered using VTK through Python? In this case, use PyGrADS.


> I need to convert it into the netcdf format, parse it to a more manageable
> size ( the  current file is 4 GB) and then convert it to a format that VTK
> can read in.
>
> I tried to use lats4d but it gives some errors:
>
> **** Error when opening the .ctl file ****
>
> ga-> lats4d -i j:\VulcanWork\mjja-co2-proj-1.ctl


Try using lats4d from the command line.  Are you doing this on Windows?  You
may want to start one of the included shells (sh or tcsh) and issue
lats4d.sh from the command line  (using the shell script). Before actually
writing any files, try

C:\> sh lats4d.sh -i yourfile.ctl -format stats

to see whether the whole file can be read fine. I am attaching the latest
lats4d.sh and lats4d.gs, lats4d.sh just in case.

   Arlindo




No hardcopy metafile open
All files closed; all defined objects released
All GrADS attributes have been reinitialized
LATS warning (lats_parmtab): Closing the previously opened Parameter Table:
.grads.lats.table
LATS error (lats_grid): LATS (latsint.c) --> Duplicate grid name: grid-1
lats4d: error creating coards file grads.lats.nc


*** metadate of the file that I want to read in and convert (the binary file
is a ".gra" file) ********

ga-> open j:\VulcanWork\mjja-co2-proj-1.
>
> ctl
>
> Scanning description file:  j:\VulcanWork\mjja-co2-proj-1.ctl
>
> Data file j:\VulcanWork\mjja-co2-proj-1.gra is open as file 1
>
> LON set to -142.874 -52.826
> LAT set to 23.494 60.214
> LEV set to 34.7 34.7
> Time values set: 2004:5:1:0 2004:5:1:0
>
> ga-> query file
> File 1 : RAMS 5.04 Output
>   Descriptor: j:\VulcanWork\mjja-co2-proj-1.ctl
>   Binary: j:\VulcanWork\mjja-co2-proj-1.gra
>   Type = Gridded
>   Xsize = 193  Ysize = 109  Zsize = 45  Tsize = 990
>   Number of Variables = 1
>     co2 45 99 - RAMS : CO2 Concentration                       [ppm     ]
>
>
> Note: Even though gradsnc gave me errors, it did output a grads.lats.ncfile, but it contains stuff like " missing coordinates" and is only 988KB,
> which I am not sure should be the case, considering the extent of the data.
>
> I am attaching the file hele if thats going to help.
>
> I searched for forums but could not find any where I could sign up - and
> that is why the mail. I hope that is not a bother! I would appreciate any
> help that you can provide.
>
> Thanks,
>
> Advait Godbole, BSME
> Doctoral Student, Earth and Atmospheric Sciences
> Global Carbon Cycle Group
> Purdue University
> West Lafayette, IN
> voice: +1 317 730 5235
>
>

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20081117/d276d26c/attachment.html