[Paraview] Trying to get sensible axes

Moreland, Kenneth kmorel at sandia.gov
Fri Aug 24 16:18:47 EDT 2012


Scott,

It's been a while since I've really looked at the CF convention, so apologies if I say something that is not quite right.

It looks like you have not specified your coordinate variables correctly.  You have variables lat, lon, x_disp, y_disp, and z_disp that look like they are supposed to be coordinate variables, but they are not according the CF convention.  A coordinate variable must have the exact same name as the dimension it describes.  Assuming you like the dimension names nx, ny, and nz, you will need to change your coordinate variables to be named nx, ny, and nz also.  (Personally, I would name the dimensions lat, lon, and height for clarity.)

It also looks like you are trying to simultaneously define Cartesian and spherical coordinates.  As far as I know, you cannot do that in the CF convention.  If it were me, I would define them as spherical coordinates (with lat/lon) and use the "Spherical Coordinates" checkbox of the ParaView reader to toggle between the two.

Hope that helps.

-Ken


Kenneth Moreland
Sandia National Laboratories
email: kmorel at sandia.gov
phone: (505) 844-8919
web:   http://www.sandia.gov/~kmorel
________________________________
From: paraview-bounces at paraview.org [paraview-bounces at paraview.org] on behalf of Scott Collis [scollis.acrf at gmail.com]
Sent: Friday, August 24, 2012 10:08 AM
To: paraview at paraview.org
Subject: [EXTERNAL] [Paraview] Trying to get sensible axes

G'day Paraview list!
I am using paraview 3.14 to read in some gridded netCDF files I have created that -semi- conforms to the CF conventions (there are still radar variables that need to be added)

I can read in and do all kinds of nice things.... but now I am getting picky :)

The reader seems to ignore my X,Y and height axes (or lat/lon)

Is my metadata wrong? (see attached header) what variable (or standard name or metadata attribute) does paraview look for... Main thing I want is nice looking cube axes..
Any advice gratefully recieved..


netcdf sgpcsapr_innermmcg_super_inneri7.c0.20110524.010010 {
dimensions:
        time = UNLIMITED ; // (1 currently)
        nz = 37 ;
        ny = 101 ;
        nx = 101 ;
variables:
        double time(time) ;
                time:units = "seconds since 1970-01-01 00:00:00.0" ;
                time:calendar = "gregorian" ;
                time:standard_name = "time" ;
        float z_disp(nz) ;
                z_disp:long_name = "height" ;
                z_disp:standard_name = "height" ;
                z_disp:units = "meter" ;
                z_disp:positive = "up" ;
        float y_disp(ny) ;
                y_disp:axis = "Y" ;
                y_disp:comment = "Y Displacement from the central facility" ;
                y_disp:long_name = "y-coordinate in Cartesian system" ;
                y_disp:standard_name = "y-coordinate in Cartesian system" ;
                y_disp:units = "m" ;
        float x_disp(nx) ;
                x_disp:axis = "X" ;
                x_disp:long_name = "x-coordinate in Cartesian system" ;
                x_disp:standard_name = "x-coordinate in Cartesian system" ;
                x_disp:units = "m" ;
                x_disp:comment = "X Displacement from the central facility" ;
        float lat(nx) ;
                lat:long_name = "latitude" ;
                lat:standard_name = "latitude" ;
                lat:units = "degrees_north" ;
        float lon(ny) ;
                lon:long_name = "longitude" ;
                lon:standard_name = "longitude" ;
                lon:units = "degrees_east" ;
        double attenuation_corrected_reflectivity_horizontal(time, nz, ny, nx) ;
                attenuation_corrected_reflectivity_horizontal:_FillValue = -9999. ;
                attenuation_corrected_reflectivity_horizontal:units = "dBZ" ;
                attenuation_corrected_reflectivity_horizontal:long_name = "equivalent_reflectivity_factor" ;
                attenuation_corrected_reflectivity_horizontal:standard_name = "equivalent_reflectivity_factor" ;
                attenuation_corrected_reflectivity_horizontal:valid_max = 80. ;
                attenuation_corrected_reflectivity_horizontal:valid_min = -45. ;
        double corrected_mean_doppler_velocity(time, nz, ny, nx) ;
                corrected_mean_doppler_velocity:_FillValue = -9999. ;
                corrected_mean_doppler_velocity:units = "m/s" ;
                corrected_mean_doppler_velocity:long_name = "radial_velocity_of_scatterers_away_from_instrument" ;
                corrected_mean_doppler_velocity:standard_name = "radial_velocity_of_scatterers_away_from_instrument" ;
                corrected_mean_doppler_velocity:valid_max = 45. ;
                corrected_mean_doppler_velocity:valid_min = -45. ;
        double diff_phase(time, nz, ny, nx) ;
                diff_phase:_FillValue = -9999. ;
                diff_phase:units = "degrees/km" ;
                diff_phase:long_name = "specific_differential_phase_hv" ;
                diff_phase:standard_name = "specific_differential_phase_hv" ;
                diff_phase:valid_max = 20. ;
                diff_phase:valid_min = -10. ;
        double doppler_spectral_width(time, nz, ny, nx) ;
                doppler_spectral_width:_FillValue = -9999. ;
                doppler_spectral_width:units = "m/s" ;
                doppler_spectral_width:long_name = "spectrum_width" ;
                doppler_spectral_width:standard_name = "spectrum_width" ;
                doppler_spectral_width:valid_max = 45. ;
                doppler_spectral_width:valid_min = 0. ;
        double unfolded_dp_phase_shift(time, nz, ny, nx) ;
                unfolded_dp_phase_shift:_FillValue = -9999. ;
                unfolded_dp_phase_shift:units = "degrees" ;
                unfolded_dp_phase_shift:long_name = "differential_phase_hv" ;
                unfolded_dp_phase_shift:standard_name = "differential_phase_hv" ;
                unfolded_dp_phase_shift:valid_max = 480. ;
                unfolded_dp_phase_shift:valid_min = 0. ;
        double reflectivity_horizontal(time, nz, ny, nx) ;
                reflectivity_horizontal:_FillValue = -9999. ;
                reflectivity_horizontal:units = "dBZ" ;
                reflectivity_horizontal:long_name = "equivalent_reflectivity_factor" ;
                reflectivity_horizontal:standard_name = "equivalent_reflectivity_factor" ;
                reflectivity_horizontal:valid_max = 80. ;
                reflectivity_horizontal:valid_min = -45. ;
        double diff_reflectivity(time, nz, ny, nx) ;
                diff_reflectivity:_FillValue = -9999. ;
                diff_reflectivity:units = "dB" ;
                diff_reflectivity:long_name = "log_differential_reflectivity_hv" ;
                diff_reflectivity:standard_name = "log_differential_reflectivity_hv" ;
                diff_reflectivity:valid_max = 8. ;
                diff_reflectivity:valid_min = -6. ;
        double copol_coeff(time, nz, ny, nx) ;
                copol_coeff:_FillValue = -9999. ;
                copol_coeff:units = "ratio" ;
                copol_coeff:long_name = "cross_correlation_ratio_hv" ;
                copol_coeff:standard_name = "cross_correlation_ratio_hv" ;
                copol_coeff:valid_max = 1. ;
                copol_coeff:valid_min = 0. ;
        double recalculated_diff_phase(time, nz, ny, nx) ;
                recalculated_diff_phase:_FillValue = -9999. ;
                recalculated_diff_phase:units = "degrees/km" ;
                recalculated_diff_phase:long_name = "specific_differential_phase_hv" ;
                recalculated_diff_phase:standard_name = "specific_differential_phase_hv" ;
                recalculated_diff_phase:valid_max = 20. ;
                recalculated_diff_phase:valid_min = -1. ;
        double norm_coherent_power(time, nz, ny, nx) ;
                norm_coherent_power:_FillValue = -9999. ;
                norm_coherent_power:comment = "Also know as Normalized Coherent Power" ;
                norm_coherent_power:valid_min = 0. ;
                norm_coherent_power:long_name = "signal_quality" ;
                norm_coherent_power:standard_name = "signal_quality" ;
                norm_coherent_power:units = "ratio" ;
                norm_coherent_power:valid_max = 1. ;
        double mean_doppler_velocity(time, nz, ny, nx) ;
                mean_doppler_velocity:_FillValue = -9999. ;
                mean_doppler_velocity:units = "m/s" ;
                mean_doppler_velocity:long_name = "radial_velocity_of_scatterers_away_from_instrument" ;
                mean_doppler_velocity:standard_name = "radial_velocity_of_scatterers_away_from_instrument" ;
                mean_doppler_velocity:valid_max = 95. ;
                mean_doppler_velocity:valid_min = -95. ;

// global attributes:
                :process_version = "Evaluation 0.6" ;
                :conventions = "CF 1.5" ;
                :command_line = "./map_sgp_csapr_super_inner.py /data-in/radar/sgp/sgpcsaprsur/sur/20110524/010544.mdv /data-in/radar/sgp/sgpsondewnpnC1/ " ;
                :dod_version = "0.5" ;
                :site_id = "I7: Nardin, Oklahoma" ;
                :facility_id = "sgp" ;
                :pyart_procs = " Dealised with the University of Washington 4DD code. Phase unfolded and attenation corrected for using the Z-PHI like Py-ART C-Band module. Mapped to a cartesian grid using the Ballsy mapper" ;
                :input_datastreams_num = "2" ;
                :input_datastreams = "SGPC-SAPR : 1 : 20110524.010010-20110524.010544;\n",
                        "sgpsondewnpn : ingest-sonde-10.6-0.sol5_10 : 20110523.232800-20110524.005842;\n",
                        "" ;
                :history = "created by user sc8 on bds2-vm1 at 24-Aug-2012,8:36:09 using ./map_sgp_csapr_super_inner.py" ;
}

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.paraview.org/pipermail/paraview/attachments/20120824/cfc366d7/attachment-0001.htm>


More information about the ParaView mailing list