Changes between Version 4 and Version 5 of ProcessingOutput
- Timestamp:
- 09/10/13 10:40:36 (11 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
ProcessingOutput
v4 v5 1 1 = Processing DYNAMICO output = 2 2 3 DYNAMICO outputs data on the native icasohedral grid in NetCDF-CF format, one file per field. Most fields are associated to hexagons, but some are associated to triangles.4 Native files can be visualized with the VTKMapper tool developed by Patrick Brockmann. Recently ParaView has become able to visualize NetCDF-CF files too. For quantitative post-processing it is probably best to interpolate to Cartesian curvilinear coordinates such as lat-lon.3 DYNAMICO outputs data on the native icasohedral grid in NetCDF-CF format, one file per field. Each file contains the data and metadata describing the grid. Most fields are associated to hexagons, but some are associated to triangles. 4 Native files can be visualized with the VTKMapper tool developed by Patrick Brockmann. Recently !ParaView has become able to visualize NetCDF-CF files too. For quantitative post-processing it is probably best to interpolate to Cartesian curvilinear coordinates such as lat-lon. 5 5 6 6 == Visualization with VTKMapper == … … 24 24 }}} 25 25 26 At this point merged.nc contains the surface pressure data (ps.nc) and the metadata expected by NCO. Additional fieldsprovided they are associated hexagons. This includes potential temperature (theta.nc), surface geopotential (phis.nc), geopotential (phi.nc).26 At this point merged.nc contains the surface pressure data (ps.nc) and the metadata expected by CDO. Additional fields can be merged into the same file using 'ncks -A' provided they are associated hexagons. This includes potential temperature (theta.nc), surface geopotential (phis.nc), geopotential (phi.nc). 27 27 28 Then nco can be invoked to interpolate all fields to a lat-lon grid : 28 CDO needs a text file describing the target grid : 29 29 30 30 {{{ 31 cdo remapdis,r180x90 merged.nc latlon.nc 31 cat > mygrid << EOF 32 gridtype = lonlat 33 xsize = 360 34 ysize = 180 35 xfirst = −179.5 36 xinc = 1 37 yfirst = -89.5 38 yinc = 1 39 EOF 32 40 }}} 33 41 34 In this example a regular mesh with 2-degree resolution is used. latlon.nc can be visualized using ncview or read/processed by standard tools. 42 In the above example the text file mygrid describing regular mesh with 1-degree resolution is created. 43 Finally cdo can be invoked to interpolate all fields to a lat-lon grid : 44 45 {{{ 46 cdo remapdis,mygrid merged.nc latlon.nc 47 }}} 48 49 latlon.nc can now be visualized using ncview or read/processed by standard tools. 35 50 36 51 This procedure is an interim solution which should be replaced by built-in output to lat-lon grids in the near future.