[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[netCDF #PSZ-497222]: Combination of nested NetCDF files
- Subject: [netCDF #PSZ-497222]: Combination of nested NetCDF files
- Date: Fri, 03 Sep 2010 13:45:33 -0600
Hi Brian,
> This is Brian Etherton, out in North Carolina. Hope all is well with
> you guys out in Colorado, and all in the netcdfgroup.
>
> I have a quandary, and the information presented to me on the Unidata
> web site is so great as to be overwhelming. Hence. My note.
>
> I have a triple nested WRF run, a 27km outer nest, a 9km middle nest,
> and a 3km inner nest. The 9km and 3km nests move every time step,
> following a hurricane (in this case, Earl).
>
> What I want to do is to combine the three NetCDF files valid at a
> moment in time into one file. In this new file, I want to use the 3km
> resolution data where I have it, then 9km where I have it, and the 27km
> where I have it.
>
> In essence, I want to remove the 27km data and replace it with the 9km
> data where I can, and remove the 9km data and replace it with the 3km
> data when I can. I want to maintain the native resolution of my data
> in all instances. I do not want to degrade everything to the 27km grid.
>
> My question: what tool would allow me to do this? I have looked at a few,
> but none has this exact feature.
>
> Thoughts?
I think you want the capabilities that will be provided by "libcf",
specifically
the implementation of the GRIDSPEC library, which is currently under development
at Unidata and several other institutions. The "libcf" library provides C and
Fortran (and eventually other language) programming interfaces for support of
the
Climate and Forecast (CF) metadata conventions, which include specifications of
coordinate systems and complex mosaic and other advanced grids that software
can
handle.
Here's the abstract of a GRIDSPEC talk at the upcoming AGU annual meeting in
December, by 9 authors oincluding Unidata's Ed Hartnett:
Increasingly earth system models perform computations on grids that
are not describable as simple, rectangular arrays (e.g. lon by lat),
instead requiring a mosaic of interacting, logically rectangular
tiles. Such grids are developed for a variety of reasons that
include removal of coordinate singularities that may degrade
numerical reliability in a region of interest (e.g. the north pole
in an ocean model) and increasing the uniformity of numerical
precision over the globe. Coupled earth system models, typically
characterized by independent coordinate reference systems for
modeling atmosphere, ocean, ice, and terrestrial processes, are
themselves examples of such mosaic grids. GRIDSPEC is a proposed set
of conventions to the Climate and Forecast library (LibCF)
describing data on mosaic grids developed by V. Balaji et
al. (Geophysical Fluid Dynamics Laboratory). It supports
unstructured assemblies of structured grids, including the
cubed-sphere and tripolar meshes. Here we review a GRIDSPEC NetCDF
format based on host, contact, grid, and data files. We will show
how mosaic grids can be created from the ground up using a C API and
the Python Climate Data Anaysis Tools (CDAT) for visualization. As
an application we use GRIDSPEC to regrid cubed-sphere data onto a
longitude-latitude grid.
GRIDSPEC includes nested grids, and an early implementation of GRIDSPEC
is currently available in the netCDF software, if configured with
--enable-libcf, but I don't know of any widely used analysis and
visualization applications have been adapted to it yet. Because it's
still under development, and because GRIDSPEC is only a proposed
extension to the CF standard that hasn't been approved yet, I expect
developers will want to wait until it's nailed down before they adapt
their applications to deal with grids specified by GRIDSPEC.
I'm Cc:ing Ed, in case he has anything to add, but he's out of email contact
until next week ...
--Russ
Russ Rew UCAR Unidata Program
address@hidden http://www.unidata.ucar.edu
Ticket Details
===================
Ticket ID: PSZ-497222
Department: Support netCDF
Priority: Normal
Status: Closed