[ncl-talk] number of subscripts

Elham Ghasemifar elham.ghasemifar at modares.ac.ir
Sun Sep 16 04:42:01 MDT 2018


dear colleges

hello

i want to plot cloudsat 2b geoprof data using ncl example, unfortunately my lat and lon have different column and i cant to relate latitude and longitude to the data, in the example script , the lat and lon related to data,

any researcher can help me how i can overcome to the error, i appreciate any help from the ncl-talk,

please see the complete script:


ncl 0> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/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> fili = "2011040011644_25451_CS_2B-GEOPROF_GRANULE_P1_R05_E03_F00.hdf"

ncl 4> varname = "CPR_Cloud_mask_2B_GEOPROF"

ncl 5> f = addfile (fili+".he2","r")

dbze = f->Radar_Reflectivity_2B_GEOPROF

ncl 6> alt40 = f->Height_2B_GEOPROF

ncl 7> kalt = dimsizes(alt40)

ncl 8> time = f->Profile_time_2B_GEOPROF

ncl 9> ntime = dimsizes(time)

ncl 10> lat = f->Latitude_2B_GEOPROF

ncl 11> nlat = dimsizes(lat)

ncl 12> lon = f->Longitude_2B_GEOPROF

ncl 13> nlon = dimsizes(lon)

ncl 14> time at units = "hours since 2011-02-10 00:00:0.0"

ncl 15> utc_date = cd_calendar(time, 0)

ncl 16> yyyymm = cd_calendar(time,-1)

ncl 17> yyyyfrac= yyyymm_to_yyyyfrac(yyyymm,0)

ncl 18> months = utc_date(:,1)

ncl 19> data = f->$varname$

ncl 20> lat at units = "degrees_north"

ncl 21> lat at units = "degrees_north"

ncl 22> data1

ncl 23> data1= dim_sum_n_Wrap(data ,1)

ncl 24> data2 = dim_avg_n_Wrap(data1,0)

ncl 25> printVarSummary(data1)

Variable: data1

Type: byte

Total Size: 37082 bytes

37082 values

Number of Dimensions: 1

Dimensions and sizes: [nray_2B_GEOPROF | 37082]

Coordinates:

Number Of Attributes: 9

_FillValue : -9

hdfeos_name : CPR_Cloud_mask

factor : 1

offset : 0

long_name : CPR Cloud Mask

valid_range : 0

missing : -9

missop : ==

sum_op_ncl : dim_sum_n over dimension(s): nbin_2B_GEOPROF

ncl 26> printVarSummary(data2)

Variable: data2

Type: float

Total Size: 4 bytes

1 values

Number of Dimensions: 1

Dimensions and sizes: [1]

Coordinates:

Number Of Attributes: 11

hdfeos_name : CPR_Cloud_mask

factor : 1

offset : 0

long_name : CPR Cloud Mask

valid_range : 0

missing : -9

missop : ==

sum_op_ncl : dim_sum_n over dimension(s): nbin_2B_GEOPROF

_FillValue_original : -9

_FillValue : -9

average_op_ncl : dim_avg_n over dimension(s): nray_2B_GEOPROF

ncl 27> pltType = "png"

ncl 28> wks = gsn_open_wks(pltType,"cloudsat")

ncl 29> plot = new(4,graphic)

ncl 30> res = True

ncl 31> res at gsnDraw = False

ncl 32> res at gsnFrame = False

ncl 33> res at lbLabelBarOn = False

ncl 34> res at mpFillOn = False

ncl 36> res at cnFillOn = True

ncl 37> res at cnFillPalette = "WhiteBlueGreenYellowRed"

ncl 38> res at cnLinesOn = False

ncl 39> res at cnLineLabelsOn = False

ncl 40> resP = True

ncl 41> resP at gsnMaximize = True

ncl 42> resP at gsnPanelLabelBar = True

ncl 43> res at cnLevelSelectionMode = "ManualLevels"

ncl 44> res at cnMinLevelValF = 0.025

ncl 45> res at cnMaxLevelValF = 0.60

ncl 46> res at cnLevelSpacingF = 0.025

ncl 47> lvl = (/ 1200, 5520, 10230, 14160 /)

ncl 48> res at gsnLeftString = "total cloud fraction"

ncl 49> do i=0,3

ncl 50> res at gsnRightString

ncl 51> res at gsnRightString= alt40({lvl(i)})+" "+alt40 at units

ncl 52> plot(i) = gsn_csm_contour_map(wks,data2({lvl(i)},:,:),res)

ncl 53> end do

fatal:Number of subscripts do not match number of dimensions of variable,(1) Subscripts used, (2) Subscripts expected

fatal:["Execute.c":8638]:Execute: Error occurred at or near line 66

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


More information about the ncl-talk mailing list