Re: drawing vertical cross section charts from GRIB files

The data are alreay on pressure levels.
Why are you using vinth2p??
===============
http://www.ncl.ucar.edu/Applications/height_lat.shtml

EG Example 6
;***********************
; h_lat_6.ncl
;***********************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;***********************
begin
   LON = 100

   f = addfile ("foo.grb","r")
   x = f->X(:,:{LON}) ; read pressure and latitudes at LON
   printVarSummary(x)

;***********************
; Create Plot
;***********************
   wks = gsn_open_wks ("ps", "h_lat" ) ; open workstation
   gsn_define_colormap(wks,"amwg") ; choose colormap

   res = True ; plot mods desired

  :res@cnLevelSelectionMode = "ManualLevels" ; manually select levels
  ;res@cnLevelSpacingF = 5.0 ; contour spacing
  ;res@cnMinLevelValF = -50. ; min level
  ;res@cnMaxLevelValF = 50. ; max level
   res@cnLineLabelsOn = True ; turn on line labels
   res@cnFillOn = True ; turn on color fill
   res@gsnSpreadColors = True ; use full range of map
   res@lbLabelAutoStride = True ; optimal labels

   plot = gsn_csm_pres_hgt(wks, x, res ) ; plaace holder

end

On 3/30/11 10:05 PM, Leung Samuel wrote:
> Hi NCL Team,
>
> I would like to study the feasibility of plotting an image in pressure
> vs latitude coordinates (attachment is the example). I need to extract
> data from grib files, but I guess I have to "convert" the elements
> inside in advance before I can use "gsn_csm_pres_hgt" for the plot.
>
> The element extracted from the grib files is a 3d array
> Dimensions and sizes: [lv_ISBL2 | 16] x [g0_lat_0 | 348] x [g0_lon_1 |
> 568]
> Coordinates:
> lv_ISBL2: [100..1000]
> g0_lat_0: [42.1..7.4]
> g0_lon_1: [89.9..146.6]
>
> I tried to put the array to the function "vinth2p" but it looked like it
> is not a correct approach.
> /fatal:vinth2p: Surface pressure must have same number of [case],
> [time], lat and lon elements as first input argument; number of
> dimensions do not match./
>
> Any suggestion? Let me know if you want to take a look at the grib file
> I got.
>
> Thanks for your help!
>
> Samuel Leung
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Mar 30 22:49:31 2011