Re: uv2vr_cfd Help!.

patton wrote:
> Hello list,
>
> I'm using :
>
> plot = gsn_csm_contour_map_overlay(wks,vrt,Hgeo,res,resh)
>
> to ploting a overlay between "Relative Vorticity and Geopotential Height"
>
> I have:
> .
> .
> .
> tf = addfile("CANUHELPME?.grib2","r")
> U = tf->U_GRD_GDS0_ISBL_10({650},:,:) V =
> tf->V_GRD_GDS0_ISBL_10({650},:,:)
> Hgeo = tf->HGT_GDS0_ISBL_10({650},:,:)
>
> wks = gsn_open_wks("x11","HELPMEXD") ;
> gsn_define_colormap(wks,"BlWhRe")
> scale = 1.e05 vrt =
> U vrt = uv2vrG_Wrap(U,V)
> .
> .
> .
> and is perfect, but is not what i want.
>
> uv2vrG_Wrap(U,V) Computes the vorticity via spherical harmonics and i
> need
> Compute it using centered finite differences.
>
> For that i founded :
>
> uv2vr_cfd()
> cool!! =D:::....
>
> but i tried writting:
> .
> .
> .
> scale = 1.e05 vrt =
> U vrt = uv2vr_cfd(U,V,:,:,3)
> .
> .
> .
> and simply i can't plot it.
>
> I have been looking for examples and more documentation about that
> function
> but... i'm lost =D.
>
====================================

I speculate that the coordinate variables have not been associated with
the vr variable as coordinate variables. Try
      vr = *uv2vr_cfd* (u,v,lat,lon, 3)
      printVarSummary(vr)
      vr_at_long_name = "vorticity"
      vr_at_units = "1/s"
      copy_VarCoords(u,vr(
      printVarSummary(vr)
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jan 17 2008 - 14:19:04 MST