ncl 5> begin
ncl 6> ;-- define file name
ncl 7>   fili  = "NPR.MIRS.V7.SND.AAMH.NP.D18289.S0854.E1047.B4992829.NS.nc"
ncl 8> 
ncl 9> ;-- open file and read variable 
ncl 10>   f     =   addfile(fili, "r")  
ncl 11>   var   =   f->PVapor(0,:,:)    ;-- first time step, latitude=40N, longitude=0-60E.
ncl 12>   lev   =  f->Plevel
ncl 13>   lat2d = f->Latitude
ncl 14>  
ncl 15>   minlat = -87 ; or some val
ncl 16>   maxlat = 76 ; or some val
ncl 17>   minlev = 0
ncl 18>   maxlev = 100 ;  
ncl 19>   ii=ind(lat2d.ge.minlat.and.lat2d.le.maxlat.and.lev.ge.minlev.and.lev.le.maxlev)  
ncl 20> 
ncl 21> ;-- define workstation
ncl 22>   wks = gsn_open_wks("png","plot_slices")
ncl 23> 
ncl 24> ;-- set resources
ncl 25>   res                          =  True
ncl 26>   res@tiMainString             = "DKRZ NCL Example: Slice plot at 40N"   ;-- title string
ncl 27> 
ncl 28>   cmap = read_colormap_file("BlAqGrYeOrReVi200")      ;-- read the color map
ncl 29> 
ncl 30>   res@cnFillOn                 =  True       ;-- turn on color fill
ncl 31>   res@cnFillPalette            = cmap(::-1,:)
ncl 32>   res@cnLineLabelsOn           =  False      ;-- turns off contour line labels
ncl 33>   res@cnInfoLabelOn            =  False      ;-- turns off contour info label
ncl 34>   res@lbOrientation            = "vertical"  ;-- vertical label bar
ncl 35>   res@tiYAxisString            =  var@long_name+"  [hPa]"
ncl 36>   res@gsnAddCyclic             = False
ncl 37>   
ncl 38>   res@sfXArray                 =  lat2d(ii)     ;-- uses lon_t as plot x-axis
ncl 39>   res@sfYArray                 =  lev(ii)  ;-- uses lev_t in hPa as plot y-axis
ncl 40> 
ncl 41>   res@trYReverse               =  True       ;-- reverses y-axis
ncl 42> 
ncl 43> ;-- generate the plot
ncl 44>   plot = gsn_csm_contour(wks,var(ii),res)
ncl 45> 
ncl 46> end
fatal:And: Number of dimensions do not match, can't continue
fatal:["Execute.c":8638]:Execute: Error occurred at or near line 19

ncl 47> stop record
ncl 48>