[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[IDV #ZMS-449881]: FW: construct a 3D grid from multifold 2D grids at discrete pressure-levels

Hi Tyn.

After consulting with Yuan, we are wondering if you may be going down the wrong 
path. Have you looked at NcML to do these sorts of aggregations?

http://www.unidata.ucar.edu/software/netcdf/ncml/v2.2/Aggregation.html

Alternatively, you may wish to ask at address@hidden.

-Julien

> Dear support,
> 
> Below my support request.
> 
> Cheers, Tyn
> 
> On 11/29/10 12:16 AM, "Valentijn Venus" <address@hidden> wrote:
> 
> >Hi,
> >
> >In Unidata's GridUtil package there is a method to slice a 3D grid to
> >obtain a 2D fieldimpl.
> >
> >Now, I basically want the reverse: construct a 3D grid from many 2D grids
> >with similar coordinate systems, notably 20 discrete (MODIS) vertical
> >levels (nlevs=20) corresponding to the following pressure-levels (in mb):
> > 5,  10,  20,  30,  50,  70,  100 , 150 ,  200 , 250 ,  300 ,   400  ,
> >500 , 620 ,  700 ,   780 ,  850 ,  920 ,  950, 1000
> >
> >For each level there is the retrieved temperature field and a geopotential
> >height field.
> >
> >So the idea is to take each temperature field and resample the
> >geopotential height field at each x,y,z point to end up with a field
> >temperature at pressure point (assuming a standard atmosphere).
> >This can be done, but the solution is not straight forward. As I have
> >limited understanding of the VisAD data model, could you help me in the
> >right direction?
> >
> >
> >Attached an IDV-bundle containing the input 2D grid, and below this
> >message a _very_ ugly attempt to get this going.
> >
> >Cheers, Tyn
> >
> >
> >def makeMODISSwathTo3DSet(a,b):
> >  from visad import FunctionType, FieldImpl, Gridded3DDoubleSet, Set, Unit
> >  from visad import RealType
> >  from visad import RealTupleType
> >  from visad import Gridded3DSet
> >  from visad import CartesianProductCoordinateSystem
> >
> >  from ucar.unidata.data.grid import GridUtil
> >  import ucar.visad.quantities.AirPressure as ap
> >
> >  timedom = a.getDomainSet()
> >  numTimes = timedom.getLength()
> >
> >  dom1 = GridUtil.getSpatialDomain(a)
> >  xyCs = dom1.getCoordinateSystem()
> >
> >  samples = dom1.getSamples(0)
> >  zCs = ap.getStandardAtmosphereCS()
> >  #print zCs
> >  print zCs.getLengths()
> >  newCS = CartesianProductCoordinateSystem(xyCs, zCs)
> >
> >  sdomain = Gridded3DSet(RealTupleType.SpatialEarth3DTuple, samples,
> >                        dom1.getX().getLength(), dom1.getY().getLength(),
> >                        None,
> >                        None,
> >                        dom1.getSetErrors(), 0, 0)
> >
> >  image = GridUtil.setSpatialDomain(a,sdomain)
> >  newTime = Gridded3DDoubleSet(RealType.Time, timedom.getSamples(0),
> >                               numTimes,
> >                               timedom.getCoordinateSystem(),
> >                               timedom.getSetUnits(),
> >                               timedom.getSetErrors(), 0)
> >  newType = FunctionType(RealTupleType.Time1DTuple, image[0].getType())
> >  timeImage = FieldImpl(newType, newTime)
> >  for i in range(numTimes):
> >    timeImage.setSample(i, image[i],0)
> >  return timeImage
> >
> >
> 
> 
> Faculty of Geo-Information Science and Earth Observation (ITC)
> University of Twente
> Chamber of Commerce: 501305360000
> 
> E-mail disclaimer
> The information in this e-mail, including any attachments, is intended for 
> the addressee only. If you are not the intended recipient, you are hereby 
> notified that any disclosure, copying, distribution or action in relation to 
> the content of this information is strictly prohibited. If you have received 
> this e-mail by mistake, please delete the message and any attachment and 
> inform the sender by return e-mail. ITC accepts no liability for any error or 
> omission in the message content or for damage of any kind that may arise as a 
> result of e-mail transmission.
> 
> 


Ticket Details
===================
Ticket ID: ZMS-449881
Department: Support IDV
Priority: Normal
Status: Closed