[ncl-talk] Fw: Fw:plot error

Elham Ghasemifar elham.ghasemifar at modares.ac.ir
Wed Sep 12 07:51:02 MDT 2018


dear colleges

any one can help me to correct the error,

thanks

________________________________
From: Elham Ghasemifar
Sent: Wednesday, September 12, 2018 1:44 PM
To: ncl-talk at ucar.edu
Subject: Fw:plot error




________________________________
From: Elham Ghasemifar
Sent: Wednesday, September 12, 2018 1:43 PM
To: ncl-talk at ucar.edu
Subject:


Dear colleges

i ran the script for cloudsat data , but in the final section i got the error, i put a script here,

Script started on mié 12 sep 2018 10:39:20 CEST
[?1034hsh-4.1$  module load ncl/6.4.0
sh-4.1$ ncl
 Copyright (C) 1995-2017 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 6.4.0
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.
ncl 0> load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
ncl 1> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
ncl 2> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
ncl 3> begTime = get_cpu_time()
  ncl 4> latS     = 35
 ncl 5> latN     = 45
ncl 6>   file_name = "2011040011644_25451_CS_2B-GEOPROF_GRANULE_P1_R05_E03_F00.hdf"
ncl 7>   eos_file=addfile(file_name+ ".he2","r")
ncl 8> print(eos_file)

Variable: eos_file
Type: file
filename:    2011040011644_25451_CS_2B-GEOPROF_GRANULE_P1_R05_E03_F00.hdf
path:    2011040011644_25451_CS_2B-GEOPROF_GRANULE_P1_R05_E03_F00.hdf
   file global attributes:
      sem_MDSignal_missop_2B_GEOPROF : ==
      sem_MDSignal_missing_2B_GEOPROF : 9999
      sem_MDSignal_valid_range_2B_GEOPROF : ( -35, -20 )
      sem_MDSignal_units_2B_GEOPROF : dBZe
      sem_MDSignal_long_name_2B_GEOPROF : minimum detectable signal
      sem_MDSignal_offset_2B_GEOPROF :  0
      sem_MDSignal_factor_2B_GEOPROF :  1
      Clutter_reduction_flag_long_name_2B_GEOPROF : Clutter_reduction_flag
      Clutter_reduction_flag_offset_2B_GEOPROF :  0
      Clutter_reduction_flag_factor_2B_GEOPROF :  1
      Navigation_land_sea_flag_valid_range_2B_GEOPROF : ( 1, 4 )
      Navigation_land_sea_flag_units_2B_GEOPROF : --
      Navigation_land_sea_flag_long_name_2B_GEOPROF : Land Sea Flag
      Navigation_land_sea_flag_offset_2B_GEOPROF :  0
      Navigation_land_sea_flag_factor_2B_GEOPROF :  1
      sem_NoiseGate_missop_2B_GEOPROF : ==
      sem_NoiseGate_missing_2B_GEOPROF : 0
--More--       sem_NoiseGate_long_name_2B_GEOPROF : Noise Floor and Noise Variance estima
--More-- te  based on window of a fixed size centered at this range bin.
--More--       sem_NoiseGate_offset_2B_GEOPROF :  0
--More--       sem_NoiseGate_factor_2B_GEOPROF :  1
--More--       sem_NoiseFloorVar_missop_2B_GEOPROF : ==
--More--       sem_NoiseFloorVar_missing_2B_GEOPROF :  0
--More--       sem_NoiseFloorVar_long_name_2B_GEOPROF : Noise Variance in measured return
--More--  (i.e. Echo) power
--More--       sem_NoiseFloorVar_offset_2B_GEOPROF :  0
--More--       sem_NoiseFloorVar_factor_2B_GEOPROF :  1
--More--       sem_NoiseFloor_missop_2B_GEOPROF : ==
--More--       sem_NoiseFloor_missing_2B_GEOPROF :  0
--More--       sem_NoiseFloor_long_name_2B_GEOPROF : Noise Floor calculated in Cloud_Mask
--More-- ing Routine
--More--       sem_NoiseFloor_offset_2B_GEOPROF :  0
--More--       sem_NoiseFloor_factor_2B_GEOPROF :  1
--More--       CPR_Echo_Top_missop_2B_GEOPROF : ==
--More--       CPR_Echo_Top_missing_2B_GEOPROF : -9
--More--       CPR_Echo_Top_valid_range_2B_GEOPROF : ( 0, 5 )
--More--       CPR_Echo_Top_long_name_2B_GEOPROF : CPR Echo Top Characterizations
--More--       CPR_Echo_Top_offset_2B_GEOPROF :  0
--More--       CPR_Echo_Top_factor_2B_GEOPROF :  1
--More--       MODIS_scene_var_missop_2B_GEOPROF : ==
--More--       MODIS_scene_var_missing_2B_GEOPROF : -9
--More--       MODIS_scene_var_valid_range_2B_GEOPROF : ( 0, 5 )
--More--       MODIS_scene_var_long_name_2B_GEOPROF : MODIS scene variability
--More--       MODIS_scene_var_offset_2B_GEOPROF :  0
--More--       MODIS_scene_var_factor_2B_GEOPROF :  1
--More--       MODIS_scene_char_missop_2B_GEOPROF : ==
--More--       MODIS_scene_char_missing_2B_GEOPROF : -9
--More--       MODIS_scene_char_valid_range_2B_GEOPROF : ( 0, 9 )
--More--       MODIS_scene_char_long_name_2B_GEOPROF : MODIS scene characterizations
--More--       MODIS_scene_char_offset_2B_GEOPROF :  0
--More--       MODIS_scene_char_factor_2B_GEOPROF :  1
--More--       MODIS_Cloud_Fraction_missop_2B_GEOPROF : ==
--More--       MODIS_Cloud_Fraction_missing_2B_GEOPROF : -99
--More--       MODIS_Cloud_Fraction_valid_range_2B_GEOPROF : ( 0, 100 )
--More--       MODIS_Cloud_Fraction_long_name_2B_GEOPROF : MODIS 250m Cloud Fraction
--More--       MODIS_Cloud_Fraction_offset_2B_GEOPROF :  0
--More--       MODIS_Cloud_Fraction_factor_2B_GEOPROF :  1
--More--       MODIS_cloud_flag_missop_2B_GEOPROF : ==
--More--       MODIS_cloud_flag_missing_2B_GEOPROF : 99
--More--       MODIS_cloud_flag_valid_range_2B_GEOPROF : ( 0, 3 )
--More--       MODIS_cloud_flag_units_2B_GEOPROF : None
--More--       MODIS_cloud_flag_long_name_2B_GEOPROF : MOD35_bit_2and3_cloud_flag
--More--       MODIS_cloud_flag_offset_2B_GEOPROF :  0
--More--       MODIS_cloud_flag_factor_2B_GEOPROF :  1
--More--       Sigma_Zero_missop_2B_GEOPROF : ==
--More--       Sigma_Zero_missing_2B_GEOPROF : -9999
--More--       Sigma_Zero_valid_range_2B_GEOPROF : ( -1000, 4000 )
--More--       Sigma_Zero_units_2B_GEOPROF : dB*100
--More--       Sigma_Zero_long_name_2B_GEOPROF : Sigma-Zero
--More--       Sigma_Zero_offset_2B_GEOPROF :  0
--More--       Sigma_Zero_factor_2B_GEOPROF :  1
--More--       Radar_Reflectivity_missop_2B_GEOPROF : ==
--More--       Radar_Reflectivity_missing_2B_GEOPROF : -8888
--More--       Radar_Reflectivity_valid_range_2B_GEOPROF : ( -4000, 5000 )
--More--       Radar_Reflectivity_units_2B_GEOPROF : dBZe
--More--       Radar_Reflectivity_long_name_2B_GEOPROF : Radar Reflectivity Factor
--More--       Radar_Reflectivity_offset_2B_GEOPROF :  0
--More--       Radar_Reflectivity_factor_2B_GEOPROF : 100
--More--       Gaseous_Attenuation_missop_2B_GEOPROF : ==
--More--       Gaseous_Attenuation_missing_2B_GEOPROF : -9999
--More--       Gaseous_Attenuation_valid_range_2B_GEOPROF : ( 0, 1000 )
--More--       Gaseous_Attenuation_units_2B_GEOPROF : dBZe
--More--       Gaseous_Attenuation_long_name_2B_GEOPROF : Gaseous_Attenuation
--More--       Gaseous_Attenuation_offset_2B_GEOPROF :  0
--More--       Gaseous_Attenuation_factor_2B_GEOPROF : 100
--More--       CPR_Cloud_mask_missop_2B_GEOPROF : ==
--More--       CPR_Cloud_mask_missing_2B_GEOPROF : -9
--More--       CPR_Cloud_mask_valid_range_2B_GEOPROF : ( 0, 40 )
--More--       CPR_Cloud_mask_long_name_2B_GEOPROF : CPR Cloud Mask
--More--       CPR_Cloud_mask_offset_2B_GEOPROF :  0
--More--       CPR_Cloud_mask_factor_2B_GEOPROF :  1
--More--       SurfaceHeightBin_fraction_missing_2B_GEOPROF :  0
--More--       SurfaceHeightBin_fraction_units_2B_GEOPROF : none
--More--       SurfaceHeightBin_fraction_long_name_2B_GEOPROF : SurfaceHeightBin_fraction
--More--       SurfaceHeightBin_fraction_offset_2B_GEOPROF :  0
--More--       SurfaceHeightBin_fraction_factor_2B_GEOPROF :  1
--More--       SurfaceHeightBin_missop_2B_GEOPROF : ==
--More--       SurfaceHeightBin_missing_2B_GEOPROF : -1
--More--       SurfaceHeightBin_valid_range_2B_GEOPROF : ( 1, 125 )
--More--       SurfaceHeightBin_long_name_2B_GEOPROF : Location of Surface Bin as determi
--More-- ned by 1B CPR algorithm.  The value here is shifted (as Height).
--More--       SurfaceHeightBin_offset_2B_GEOPROF :  0
--More--       SurfaceHeightBin_factor_2B_GEOPROF :  1
--More--       RayStatus_validity_valid_range_2B_GEOPROF : ( 0, -1 )
--More--       RayStatus_validity_units_2B_GEOPROF : --
--More--       RayStatus_validity_long_name_2B_GEOPROF : Ray Status Validity
--More--       RayStatus_validity_offset_2B_GEOPROF :  0
--More--       RayStatus_validity_factor_2B_GEOPROF :  1
--More--       Data_targetID_valid_range_2B_GEOPROF : ( 0, -1 )
--More--       Data_targetID_units_2B_GEOPROF : --
--More--       Data_targetID_long_name_2B_GEOPROF : CPR bus orientation (target ID)
--More--       Data_targetID_offset_2B_GEOPROF :  0
--More--       Data_targetID_factor_2B_GEOPROF :  1
--More--       Data_status_valid_range_2B_GEOPROF : ( 0, -1 )
--More--       Data_status_units_2B_GEOPROF : --
--More--       Data_status_long_name_2B_GEOPROF : Data status flags
--More--       Data_status_offset_2B_GEOPROF :  0
--More--       Data_status_factor_2B_GEOPROF :  1
--More--       Data_quality_valid_range_2B_GEOPROF : ( 0, -1 )
--More--       Data_quality_units_2B_GEOPROF : --
--More--       Data_quality_long_name_2B_GEOPROF : Data Quality
--More--       Data_quality_offset_2B_GEOPROF :  0
--More--       Data_quality_factor_2B_GEOPROF :  1
--More--       Roll_offset_valid_range_2B_GEOPROF : ( -90, 90 )
--More--       Roll_offset_units_2B_GEOPROF : degrees
--More--       Roll_offset_long_name_2B_GEOPROF : Nominal satellite roll angle offset fro
--More-- m nadir
--More--       Roll_offset_offset_2B_GEOPROF :  0
--More--       Roll_offset_factor_2B_GEOPROF :  1
--More--       Pitch_offset_valid_range_2B_GEOPROF : ( -90, 90 )
--More--       Pitch_offset_units_2B_GEOPROF : degrees
--More--       Pitch_offset_long_name_2B_GEOPROF : Nominal satellite pitch angle offset f
--More-- rom nadir
--More--       Pitch_offset_offset_2B_GEOPROF :  0
--More--       Pitch_offset_factor_2B_GEOPROF :  1
--More--       Vertical_binsize_missop_2B_GEOPROF : ==
--More--       Vertical_binsize_missing_2B_GEOPROF : -9999
--More--       Vertical_binsize_units_2B_GEOPROF : m
--More--       Vertical_binsize_offset_2B_GEOPROF :  0
--More--       Vertical_binsize_factor_2B_GEOPROF :  1
--More--       DEM_elevation_missop_2B_GEOPROF : ==
--More--       DEM_elevation_missing_2B_GEOPROF : 9999
--More--       DEM_elevation_valid_range_2B_GEOPROF : ( -9999, 8850 )
--More--       DEM_elevation_units_2B_GEOPROF : meters
--More--       DEM_elevation_long_name_2B_GEOPROF : Digital Elevation Map
--More--       DEM_elevation_offset_2B_GEOPROF :  0
--More--       DEM_elevation_factor_2B_GEOPROF :  1
--More--       Range_to_intercept_valid_range_2B_GEOPROF : ( 600, 800 )
--More--       Range_to_intercept_units_2B_GEOPROF : km
--More--       Range_to_intercept_long_name_2B_GEOPROF : Range from the spacecraft to the
--More--  CPR boresight intercept with the geoid
--More--       Range_to_intercept_offset_2B_GEOPROF :  0
--More--       Range_to_intercept_factor_2B_GEOPROF :  1
--More--       Height_missop_2B_GEOPROF : ==
--More--       Height_missing_2B_GEOPROF : -9999
--More--       Height_valid_range_2B_GEOPROF : ( -5000, 30000 )
--More--       Height_units_2B_GEOPROF : m
--More--       Height_long_name_2B_GEOPROF : Height of range bin in Reflectivity/Cloud Ma
--More-- sk above reference surface (~ mean sea level).
--More--       Height_offset_2B_GEOPROF :  0
--More--       Height_factor_2B_GEOPROF :  1
--More--       Longitude_valid_range_2B_GEOPROF : ( -180, 180 )
--More--       Longitude_units_2B_GEOPROF : degrees
--More--       Longitude_long_name_2B_GEOPROF : Spacecraft Longitude
--More--       Longitude_offset_2B_GEOPROF :  0
--More--       Longitude_factor_2B_GEOPROF :  1
--More--       Latitude_valid_range_2B_GEOPROF : ( -90, 90 )
--More--       Latitude_units_2B_GEOPROF : degrees
--More--       Latitude_long_name_2B_GEOPROF : Spacecraft Latitude
--More--       Latitude_offset_2B_GEOPROF :  0
--More--       Latitude_factor_2B_GEOPROF :  1
--More--       TAI_start_valid_range_2B_GEOPROF : (    0, 600000000 )
--More--       TAI_start_units_2B_GEOPROF : seconds
--More--       TAI_start_long_name_2B_GEOPROF : TAI time for the first profile.
--More--       TAI_start_offset_2B_GEOPROF :  0
--More--       TAI_start_factor_2B_GEOPROF :  1
--More--       UTC_start_valid_range_2B_GEOPROF : (  0, 86400 )
--More--       UTC_start_units_2B_GEOPROF : seconds
--More--       UTC_start_long_name_2B_GEOPROF : UTC seconds since 00:00 Z of the first pr
--More-- ofile
--More--       UTC_start_offset_2B_GEOPROF :  0
--More--       UTC_start_factor_2B_GEOPROF :  1
--More--       Profile_time_valid_range_2B_GEOPROF : (  0, 6000 )
--More--       Profile_time_units_2B_GEOPROF : seconds
--More--       Profile_time_long_name_2B_GEOPROF : Seconds since the start of the granule
--More-- .
--More--       Profile_time_offset_2B_GEOPROF :  0
--More--       Profile_time_factor_2B_GEOPROF :  1
--More--       end_time_2B_GEOPROF : 20110209025537
--More--       start_time_2B_GEOPROF : 20110209011644
--More--       product_version_2B_GEOPROF : P1_R0
--More--       algorithm_version_2B_GEOPROF : P1_R05.
--More--       algorithm_name_2B_GEOPROF : 2B-GEOPROF
--More--       ID_CREATED_2B_GEOPROF : Wed Jan 24 21:49:27 2018
--More--       ID_MACHINE_2B_GEOPROF : CSATLAS1
--More--       ID_CENTER_URL_2B_GEOPROF : http://cloudsat.cira.colostate.edu
--More--       ID_CENTER_2B_GEOPROF : CloudSat Data Processing Center
--More--       ID_SITE_2B_GEOPROF : Cooperative Institute for Research in the Atmosphere
--More--       granule_number_2B_GEOPROF : 25451
--More--    dimensions:
--More--       nray_2B_GEOPROF = 37082
--More--       scalar_2B_GEOPROF = 1
--More--       nbin_2B_GEOPROF = 125
--More--    variables:
--More--       float sem_MDSignal_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    sem_MDSignal
--More--          _FillValue :    9999
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    minimum detectable signal
--More--          units :    dBZe
--More--          valid_range :    -35
--More--          missing :    9999
--More--          missop :    ==
--More--  
--More--       byte Clutter_reduction_flag_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    Clutter_reduction_flag
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Clutter_reduction_flag
--More--  
--More--       byte Navigation_land_sea_flag_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    Navigation_land_sea_flag
--More--          unsigned :    True
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Land Sea Flag
--More--          units :    --
--More--          valid_range :    1
--More--  
--More--       byte sem_NoiseGate_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    sem_NoiseGate
--More--          _FillValue :    0
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Noise Floor and Noise Variance estimate  based on window
--More--  of a fixed size centered at this range bin.
--More--          missing :    0
--More--          missop :    ==
--More--  
--More--       float sem_NoiseFloorVar_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    sem_NoiseFloorVar
--More--          _FillValue :     0
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Noise Variance in measured return (i.e. Echo) power
--More--          missing :     0
--More--          missop :    ==
--More--  
--More--       float sem_NoiseFloor_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    sem_NoiseFloor
--More--          _FillValue :     0
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Noise Floor calculated in Cloud_Masking Routine
--More--          missing :     0
--More--          missop :    ==
--More--          ar_factor :     1
--More--          ar_offset :     0
--More--          ar_long_name :    Noise Variance in measured return (i.e. Echo) power
--More--          ar_missing :     0
--More--          ar_missop :    ==
--More--  
--More--       byte CPR_Echo_Top_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    CPR_Echo_Top
--More--          _FillValue :    -9
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    CPR Echo Top Characterizations
--More--          valid_range :    0
--More--          missing :    -9
--More--          missop :    ==
--More--  
--More--       byte MODIS_scene_var_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    MODIS_scene_var
--More--          _FillValue :    -9
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    MODIS scene variability
--More--          valid_range :    0
--More--          missing :    -9
--More--          missop :    ==
--More--  
--More--       byte MODIS_scene_char_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    MODIS_scene_char
--More--          _FillValue :    -9
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    MODIS scene characterizations
--More--          valid_range :    0
--More--          missing :    -9
--More--          missop :    ==
--More--  
--More--       byte MODIS_Cloud_Fraction_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    MODIS_Cloud_Fraction
--More--          _FillValue :    -99
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    MODIS 250m Cloud Fraction
--More--          valid_range :    0
--More--          missing :    -99
--More--          missop :    ==
--More--  
--More--       byte MODIS_cloud_flag_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    MODIS_cloud_flag
--More--          _FillValue :    99
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    MOD35_bit_2and3_cloud_flag
--More--          units :    None
--More--          valid_range :    0
--More--          missing :    99
--More--          missop :    ==
--More--  
--More--       short Sigma_Zero_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    Sigma-Zero
--More--          _FillValue :    -9999
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Sigma-Zero
--More--          units :    dB*100
--More--          valid_range :    -1000
--More--          missing :    -9999
--More--          missop :    ==
--More--  
--More--       short Radar_Reflectivity_2B_GEOPROF ( nray_2B_GEOPROF, nbin_2B_GEOPROF )
--More--          hdfeos_name :    Radar_Reflectivity
--More--          _FillValue :    -8192
--More--          factor :    100
--More--          offset :     0
--More--          long_name :    Radar Reflectivity Factor
--More--          units :    dBZe
--More--          valid_range :    -4000
--More--          missing :    -8888
--More--          missop :    ==
--More--          _FillValue :    -8192
--More--  
--More--       short Gaseous_Attenuation_2B_GEOPROF ( nray_2B_GEOPROF, nbin_2B_GEOPROF )
--More--          hdfeos_name :    Gaseous_Attenuation
--More--          _FillValue :    15360
--More--          factor :    100
--More--          offset :     0
--More--          long_name :    Gaseous_Attenuation
--More--          units :    dBZe
--More--          valid_range :    0
--More--          missing :    -9999
--More--          missop :    ==
--More--          _FillValue :    15360
--More--  
--More--       byte CPR_Cloud_mask_2B_GEOPROF ( nray_2B_GEOPROF, nbin_2B_GEOPROF )
--More--          hdfeos_name :    CPR_Cloud_mask
--More--          _FillValue :    -9
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    CPR Cloud Mask
--More--          valid_range :    0
--More--          missing :    -9
--More--          missop :    ==
--More--          _FillValue :    -9
--More--  
--More--       float SurfaceHeightBin_fraction_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    SurfaceHeightBin_fraction
--More--          _FillValue :     0
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    SurfaceHeightBin_fraction
--More--          units :    none
--More--          missing :     0
--More--  
--More--       byte SurfaceHeightBin_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    SurfaceHeightBin
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Location of Surface Bin as determined by 1B CPR algorith
--More-- m.  The value here is shifted (as Height).
--More--          valid_range :    1
--More--          missing :    -1
--More--          missop :    ==
--More--          fraction_factor :     1
--More--          fraction_offset :     0
--More--          fraction_long_name :    SurfaceHeightBin_fraction
--More--          fraction_units :    none
--More--          fraction_missing :     0
--More--  
--More--       byte RayStatus_validity_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    RayStatus_validity
--More--          unsigned :    True
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Ray Status Validity
--More--          units :    --
--More--          valid_range :    0
--More--  
--More--       byte Data_targetID_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    Data_targetID
--More--          unsigned :    True
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    CPR bus orientation (target ID)
--More--          units :    --
--More--          valid_range :    0
--More--  
--More--       short Data_status_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    Data_status
--More--          unsigned :    True
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Data status flags
--More--          units :    --
--More--          valid_range :    0
--More--  
--More--       byte Data_quality_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    Data_quality
--More--          unsigned :    True
--More--          factor :     1
--More--          offset :     0
--More--          long_name :    Data Quality
--More--          units :    --
--More--          valid_range :    0
--More--  
--More--       float Roll_offset_2B_GEOPROF ( scalar_2B_GEOPROF )
--More--          hdfeos_name :    Roll_offset
--More--  
--More--       float Pitch_offset_2B_GEOPROF ( scalar_2B_GEOPROF )
--More--          hdfeos_name :    Pitch_offset
--More--  
--More--       float Vertical_binsize_2B_GEOPROF ( scalar_2B_GEOPROF )
--More--          hdfeos_name :    Vertical_binsize
--More--  
--More--       short DEM_elevation_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    DEM_elevation
--More--          _FillValue :    9999
--More--  
--More--       float Range_to_intercept_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    Range_to_intercept
--More--  
--More--       short Height_2B_GEOPROF ( nray_2B_GEOPROF, nbin_2B_GEOPROF )
--More--          hdfeos_name :    Height
--More--          _FillValue :    -9999
--More--  
--More--       float Longitude_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          long_name :    longitude
--More--          units :    degrees_east
--More--          hdfeos_name :    Longitude
--More--  
--More--       float Latitude_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          long_name :    latitude
--More--          units :    degrees_north
--More--          hdfeos_name :    Latitude
--More--  
--More--       double TAI_start_2B_GEOPROF ( scalar_2B_GEOPROF )
--More--          hdfeos_name :    TAI_start
--More--  
--More--       float UTC_start_2B_GEOPROF ( scalar_2B_GEOPROF )
--More--          hdfeos_name :    UTC_start
--More--  
--More--       float Profile_time_2B_GEOPROF ( nray_2B_GEOPROF )
--More--          hdfeos_name :    Profile_time
--More--  
ncl 9>
ncl 10>
ncl 11>
ncl 12>
ncl 13>
ncl 14>
ncl 15>
ncl 16>
ncl 17>
ncl 18>
ncl 19>
ncl 20>
ncl 21>
ncl 22>
ncl 23>   LAT   = eos_file->Latitude_2B_GEOPROF
ncl 24>   iTrop = ind(LAT.ge.latS .and. LAT.le.latN)
ncl 25>   printVarSummary(iTrop)

Variable: iTrop
Type: integer
Total Size: 8364 bytes
            2091 values
Number of Dimensions: 1
Dimensions and sizes:    [2091]
Coordinates:
ncl 26>   data_raw = eos_file->Radar_Reflectivity_2B_GEOPROF(iTrop,:) ; (nray_2B_GEOPROF, nbin_2B_GEOPROF)  ; (N, 125)
ncl 27> printVarSummary(data_raw)

Variable: data_raw
Type: short
Total Size: 522750 bytes
            261375 values
Number of Dimensions: 2
Dimensions and sizes:    [nray_2B_GEOPROF | 2091] x [nbin_2B_GEOPROF | 125]
Coordinates:
Number Of Attributes: 9
  hdfeos_name :    Radar_Reflectivity
  _FillValue :    -8192
  factor :    100
  offset :     0
  long_name :    Radar Reflectivity Factor
  units :    dBZe
  valid_range :    -4000
  missing :    -8888
  missop :    ==
ncl 28>   lat = eos_file->Latitude_2B_GEOPROF(iTrop)
ncl 29>   lon = eos_file->Longitude_2B_GEOPROF(iTrop)
ncl 30>   lev = eos_file->Height_2B_GEOPROF(iTrop,:)
ncl 31> lev at long_name  ="COARDS & CF conformant "
ncl 32> lev at units      = eos_file at Height_units_2B_GEOPROF
ncl 32>  ncl 32>  ncl 32>  lev at long_name  ="COARDS & CF conformant " ncl 32>                                             lev = eos_file->Height_2B_GEOPROF(iTrop,:)  ncl 32>                                               lev at long_name  ="COARDS & CF conformant " ncl 32>                                            time           = eos_file->Profile_time_2B_GEOPROF(iTrop)
ncl 33> time at long_name = eos_file at Profile_time_long_name_2B_GEOPROF
ncl 34> time at units     = eos_file at Profile_time_units_2B_GEOPROF
ncl 35> hdf4_file=addfile(file_name,"r")
ncl 36>   print(hdf4_file)

Variable: hdf4_file
Type: file
filename:    2011040011644_25451_CS_2B-GEOPROF_GRANULE_P1_R05_E03_F00
path:    2011040011644_25451_CS_2B-GEOPROF_GRANULE_P1_R05_E03_F00.hdf
   file global attributes:
      HDFEOSVersion : HDFEOS_V2.5
      StructMetadata_0 : GROUP=SwathStructure
    GROUP=SWATH_1
        SwathName="2B-GEOPROF"
        GROUP=Dimension
            OBJECT=Dimension_1
                DimensionName="nray"
                Size=37082
            END_OBJECT=Dimension_1
            OBJECT=Dimension_2
                DimensionName="scalar"
                Size=1
            END_OBJECT=Dimension_2
            OBJECT=Dimension_3
                DimensionName="nbin"
                Size=125
            END_OBJECT=Dimension_3
--More--         END_GROUP=Dimension
--More--         GROUP=DimensionMap
--More--         END_GROUP=DimensionMap
--More--         GROUP=IndexDimensionMap
--More--         END_GROUP=IndexDimensionMap
--More--         GROUP=GeoField
--More--             OBJECT=GeoField_1
--More--                 GeoFieldName="Profile_time"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("nray")
--More--             END_OBJECT=GeoField_1
--More--             OBJECT=GeoField_2
--More--                 GeoFieldName="UTC_start"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("scalar")
--More--             END_OBJECT=GeoField_2
--More--             OBJECT=GeoField_3
--More--                 GeoFieldName="TAI_start"
--More--                 DataType=DFNT_FLOAT64
--More--                 DimList=("scalar")
--More--             END_OBJECT=GeoField_3
--More--             OBJECT=GeoField_4
--More--                 GeoFieldName="Latitude"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("nray")
--More--             END_OBJECT=GeoField_4
--More--             OBJECT=GeoField_5
--More--                 GeoFieldName="Longitude"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("nray")
--More--             END_OBJECT=GeoField_5
--More--             OBJECT=GeoField_6
--More--                 GeoFieldName="Height"
--More--                 DataType=DFNT_INT16
--More--                 DimList=("nray","nbin")
--More--                 CompressionType=HDFE_COMP_DEFLATE
--More--                 DeflateLevel=4
--More--             END_OBJECT=GeoField_6
--More--             OBJECT=GeoField_7
--More--                 GeoFieldName="Range_to_intercept"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("nray")
--More--             END_OBJECT=GeoField_7
--More--             OBJECT=GeoField_8
--More--                 GeoFieldName="DEM_elevation"
--More--                 DataType=DFNT_INT16
--More--                 DimList=("nray")
--More--             END_OBJECT=GeoField_8
--More--             OBJECT=GeoField_9
--More--                 GeoFieldName="Vertical_binsize"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("scalar")
--More--             END_OBJECT=GeoField_9
--More--             OBJECT=GeoField_10
--More--                 GeoFieldName="Pitch_offset"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("scalar")
--More--             END_OBJECT=GeoField_10
--More--             OBJECT=GeoField_11
--More--                 GeoFieldName="Roll_offset"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("scalar")
--More--             END_OBJECT=GeoField_11
--More--         END_GROUP=GeoField
--More--         GROUP=DataField
--More--             OBJECT=DataField_1
--More--                 DataFieldName="Data_quality"
--More--                 DataType=DFNT_UINT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_1
--More--             OBJECT=DataField_2
--More--                 DataFieldName="Data_status"
--More--                 DataType=DFNT_UINT16
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_2
--More--             OBJECT=DataField_3
--More--                 DataFieldName="Data_targetID"
--More--                 DataType=DFNT_UINT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_3
--More--             OBJECT=DataField_4
--More--                 DataFieldName="RayStatus_validity"
--More--                 DataType=DFNT_UINT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_4
--More--             OBJECT=DataField_5
--More--                 DataFieldName="SurfaceHeightBin"
--More--                 DataType=DFNT_INT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_5
--More--             OBJECT=DataField_6
--More--                 DataFieldName="SurfaceHeightBin_fraction"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_6
--More--             OBJECT=DataField_7
--More--                 DataFieldName="CPR_Cloud_mask"
--More--                 DataType=DFNT_INT8
--More--                 DimList=("nray","nbin")
--More--                 CompressionType=HDFE_COMP_DEFLATE
--More--                 DeflateLevel=4
--More--             END_OBJECT=DataField_7
--More--             OBJECT=DataField_8
--More--                 DataFieldName="Gaseous_Attenuation"
--More--                 DataType=DFNT_INT16
--More--                 DimList=("nray","nbin")
--More--                 CompressionType=HDFE_COMP_DEFLATE
--More--                 DeflateLevel=4
--More--             END_OBJECT=DataField_8
--More--             OBJECT=DataField_9
--More--                 DataFieldName="Radar_Reflectivity"
--More--                 DataType=DFNT_INT16
--More--                 DimList=("nray","nbin")
--More--                 CompressionType=HDFE_COMP_DEFLATE
--More--                 DeflateLevel=4
--More--             END_OBJECT=DataField_9
--More--             OBJECT=DataField_10
--More--                 DataFieldName="Sigma-Zero"
--More--                 DataType=DFNT_INT16
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_10
--More--             OBJECT=DataField_11
--More--                 DataFieldName="MODIS_cloud_flag"
--More--                 DataType=DFNT_INT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_11
--More--             OBJECT=DataField_12
--More--                 DataFieldName="MODIS_Cloud_Fraction"
--More--                 DataType=DFNT_INT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_12
--More--             OBJECT=DataField_13
--More--                 DataFieldName="MODIS_scene_char"
--More--                 DataType=DFNT_INT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_13
--More--             OBJECT=DataField_14
--More--                 DataFieldName="MODIS_scene_var"
--More--                 DataType=DFNT_INT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_14
--More--             OBJECT=DataField_15
--More--                 DataFieldName="CPR_Echo_Top"
--More--                 DataType=DFNT_INT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_15
--More--             OBJECT=DataField_16
--More--                 DataFieldName="sem_NoiseFloor"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_16
--More--             OBJECT=DataField_17
--More--                 DataFieldName="sem_NoiseFloorVar"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_17
--More--             OBJECT=DataField_18
--More--                 DataFieldName="sem_NoiseGate"
--More--                 DataType=DFNT_INT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_18
--More--             OBJECT=DataField_19
--More--                 DataFieldName="Navigation_land_sea_flag"
--More--                 DataType=DFNT_UINT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_19
--More--             OBJECT=DataField_20
--More--                 DataFieldName="Clutter_reduction_flag"
--More--                 DataType=DFNT_INT8
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_20
--More--             OBJECT=DataField_21
--More--                 DataFieldName="sem_MDSignal"
--More--                 DataType=DFNT_FLOAT32
--More--                 DimList=("nray")
--More--             END_OBJECT=DataField_21
--More--         END_GROUP=DataField
--More--         GROUP=MergedFields
--More--         END_GROUP=MergedFields
--More--     END_GROUP=SWATH_1
--More-- END_GROUP=SwathStructure
--More-- GROUP=GridStructure
--More-- END_GROUP=GridStructure
--More-- GROUP=PointStructure
--More-- END_GROUP=PointStructure
--More-- END
--More--    dimensions:
--More--       nray_2B_GEOPROF = 37082
--More--       nbin_2B_GEOPROF = 125
--More--    variables:
--More--       short Height ( nray_2B_GEOPROF, nbin_2B_GEOPROF )
--More--          _FillValue :    -9999
--More--          hdf_name :    Height
--More--  
--More--       byte CPR_Cloud_mask ( nray_2B_GEOPROF, nbin_2B_GEOPROF )
--More--          _FillValue :    -9
--More--          hdf_name :    CPR_Cloud_mask
--More--  
--More--       short Gaseous_Attenuation ( nray_2B_GEOPROF, nbin_2B_GEOPROF )
--More--          _FillValue :    15360
--More--          hdf_name :    Gaseous_Attenuation
--More--  
--More--       short Radar_Reflectivity ( nray_2B_GEOPROF, nbin_2B_GEOPROF )
--More--          _FillValue :    -8192
--More--          hdf_name :    Radar_Reflectivity
--More--  
ncl 37>
ncl 38>
ncl 39>
ncl 40>
ncl 41>
ncl 42>
ncl 43>
ncl 44>
ncl 45>
ncl 46>
ncl 47>
ncl 48>
ncl 49>
ncl 50>
ncl 51>
ncl 52>  data_hdf4 = hdf4_file->Radar_Reflectivity(iTrop,:)
ncl 53>  print("data_hdf4: min="+min(data_hdf4)+"  max="+max(data_hdf4))
(0)    data_hdf4: min=-8888  max=16874
ncl 54> data_hdf4 at valid_range=(/5000,-4000/)
ncl 55> data_hdf4 at factor =100
ncl 56> data  = where(data_raw.ge.data_hdf4 at valid_range(0) .and. data_raw.le.data_hdf4 at valid_range(1) ,data_raw, data_hdf4 at _FillValue)
ncl 57> data at _FillValue = data_hdf4 at _FillValue
ncl 58> dataf   = tofloat(data)
ncl 59> dataf   = dataf / data_hdf4 at factor
ncl 60>   dataf!0 = data_raw!0
ncl 61>   dataf!1 = data_raw!1
ncl 62>  printVarSummary(dataf)

Variable: dataf
Type: float
Total Size: 1045500 bytes
            261375 values
Number of Dimensions: 2
Dimensions and sizes:    [nray_2B_GEOPROF | 2091] x [nbin_2B_GEOPROF | 125]
Coordinates:
Number Of Attributes: 1
  _FillValue :    -8192
ncl 63> ;;dataf&nbin_2B_GEOPROF = lev(0,:)
ncl 64> dataf&nray_2B_GEOPROF = time
ncl 65> data_plt = dataf(nbin_2B_GEOPROF|:,nray_2B_GEOPROF|:)
ncl 66> data_plt       = where(data_plt.lt.-20, data_plt at _FillValue, data_plt) t                                       .-20, data_plt at _FillValue, data_plt)  ncl 66> data_plt       = where(data_plt.le.-20, data_plt at _FillValue, data_plt)  ncl 66> data_plt       = where(data_plt.le.-20                                   , data_plt at _FillValue, data_plt)  ncl 66> data_plt       = where(data_plt.le.-2                                   , data_plt at _FillValue, data_plt)  ncl 66> data_plt       = where(data_plt.le.-                                   , data_plt at _FillValue, data_plt)  ncl 66> data_plt       = where(data_plt.le.3, data_plt at _FillValue, data_plt)  ncl 66> data_plt       = where(data_plt.le.35, data_plt at _FillValue, data_plt)  ncl 66> data_plt       = where(data_plt.le.35
ncl 67> if (latS.lt.0) then
ncl 68> latLabelS = abs(latS)+"S"
ncl 69> else
ncl 70> latLabelS = latS+"N"
ncl 71> end if
ncl 72> if (latN.lt.0) then
ncl 73> latLabelN = abs(latN)+"S"
ncl 74> else
ncl 75> latLabelN = latN+"N"
ncl 76> end if
ncl 77> latLabel = latLabelS+"_"+latLabelN
ncl 78>   pltTitle = str_get_cols(file_name,0,18)+"_"+latLabel
ncl 79> xwks = gsn_open_wks("pdf", file_name)
ncl 80>   xwks = gsn_open_wks("pdf", pltTitle)
ncl 81>   gsn_define_colormap(xwks,"BlAqGrYeOrReVi200")
ncl 82>  res = True
ncl 83>   res at cnFillOn = True
ncl 84>   res at gsnMaximize = True
ncl 85>   res at gsnPaperOrientation = "portrait"
ncl 86>   res at cnLinesOn = False
ncl 87>   res at cnLineLabelsOn = False
ncl 88>   res at gsnSpreadColors = True
ncl 89>   res at cnFillMode = "RasterFill"
ncl 90>   res at lbOrientation = "vertical"
ncl 91>   res at cnMissingValFillPattern = 0
ncl 92>   res at cnMissingValFillColor = 0
ncl 93>   res at lbLabelAutoStride = True
ncl 94> res at tiMainString = file_name
ncl 95>  res at tiMainString = str_get_cols(file_name,0,18)
ncl 96>   res at tiXAxisString = time at long_name+" ("+time at units+")"
ncl 97> res at tiYAxisString = lev at long_name+" ("+lev at units+")"
warning:Attempt to reference attribute (units) which is undefined
 ncl 98> lev at units=height       eight ncl 98> lev at units=Height ncl 98> lev at units=HHeight  ncl 98> lev at units=Heigh  ncl 98> lev at units=Heig  ncl 98> lev at units=Hei  ncl 98> lev at units=He  ncl 98> lev at units=H  ncl 98> lev at units=meter" "meter"  ncl 98> lev at units=" meter"  ncl 98> lev at units="          meter"  ncl 98> lev at units="
ncl 99> res at tiYAxisString = lev at long_name+" ("+lev at units+")"
ncl 100> res at gsnLeftString=data_hdf4 at long_name
warning:Attempt to reference attribute (long_name) which is undefined
 ncl 101> lev at long_name  ="COARDS & CF conformant "
ncl 102>  ncl 102>  lev at long_name  ="COARDS & CF conformant " ncl 102>                                           res at gsnLeftString=data_hdf4 at long_name
warning:Attempt to reference attribute (long_name) which is undefined
 ncl 103> data_hdf4 at long_name  ="COARDS & CF conformant "
ncl 104> res at gsnRightString=data_hdf4 at units
warning:Attempt to reference attribute (units) which is undefined
 ncl 105> data_hdf4 at units=dbz""=dbz" ncl 105> data_hdf4 at units"       =dbz" ncl 105> data_hdf4 at units="dbz" ncl 105> data_hdf4 at units="
ncl 106> res at sfXArray   = conform(data_plt, time, 1)
ncl 107> res at sfYArray   = lev(nbin_2B_GEOPROF|:,nray_2B_GEOPROF|:)
ncl 108>   res at trGridType = "Triangularmesh"
ncl 109> plot=gsn_csm_contour(xwks, data_plt, res)
warning:ContourPlotInitialize: no valid values in scalar field; ContourPlot not possible:[errno=1101]
 ncl 110>
[1]+  Stopped(SIGTSTP)        ncl
sh-4.1

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180912/61255930/attachment.html>


More information about the ncl-talk mailing list