Re: WRF data fro streamlines

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Tue, 23 Oct 2007 08:56:24 -0600

I f you look at the Example 7 carefully you will not that
on the right hand sude the variable ar capitalized [U,V]
while in the left they are small [u,v].

Given that you are already using small u,v ... I would merely
change the left hand side and plot the U,V

            U = 0.5*(u(level,:,0:mlonU-2)+u(level,:,1:mlonU-1))
           V = 0.5*(v(level,0:nlatV-2,:)+v( level,1:nlatV-1,:))

           copy_VarAtts (u,U)
           copy_varAtts (v,V)

           lat2d = a->XLAT(0,:,:)
           lon2d = a->XLONG(0,:,)

           U_at_lat2d = lat2d
           U_at_lon2d = lon2d
           V_at_lat2d = lat2d
           V_at_lon2d = lon2d

Erik Noble wrote:
> Hi. I rewrote the portion of the my code, following the the method
> outlined in example 7 of the WRF lambert conformal page:
>
> http://www.ncl.ucar.edu/Applications/wrflc.shtml as follows:
>
> do it = 324,324 ;;;;;Time step you want!!!!
> ;Variables
> p = wrf_user_getvar(a, "pressure",it) ; pressure is our vertical coordinate
> u = wrf_user_getvar(a,"ua",it) ; u averaged to mass points
> v = wrf_user_getvar(a,"va",it) ; v averaged to mass points
> z = wrf_user_getvar(a, "z",it) ; grid point height
> dimU = dimsizes(u) ; demo "dimsizes"
> mlonU = dimU(2) ; number of "west_east_stag" longitudes
> dimV = dimsizes(v)
> nlatV = dimV(2)
> printVarSummary(u)
> printVarSummary(v)
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> do level = 0,nlevels-1 ; LOOP OVER LEVELS
>
> pressure = pressure_levels(level)
>
> u = 0.5*(u(level,:,0:mlonU-2)+u(level,:,1:mlonU-1))
> v = 0.5*(v(level,0:nlatV-2,:)+v( level,1:nlatV-1,:))
>
> And this is the error that I am getting;
>
> noble:/Volumes/Data_and_Models/ncl_scripts enoble$ ncl wrf_Streamlines.ncl
> Copyright (C) 1995-2007 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 4.3.1
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> (0) Working on time: 2006-09-10_12:00:00
>
>
> Variable: u
> Type: float
> Total Size: 1636308 bytes
> 409077 values
> Number of Dimensions: 3
> Dimensions and sizes: [27] x [109] x [139]
> Coordinates:
> Number Of Attributes: 2
> description : u Velocity
> units : m/s
>
>
> Variable: v
> Type: float
> Total Size: 1636308 bytes
> 409077 values
> Number of Dimensions: 3
> Dimensions and sizes: [27] x [109] x [139]
> Coordinates:
> Number Of Attributes: 2
> description : v Velocity
> units : m/s
> fatal:Number of dimensions on right hand side do not match number of dimension
> n left hand side
> fatal:Execute: Error occurred at or near line 70 in file wrf_Streamlines.ncl
>
> noble:/Volumes/Data_and_Models/ncl_scripts enoble$
>
> How do I fix this?
> -Erik
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Oct 23 2007 - 08:56:24 MDT

This archive was generated by hypermail 2.2.0 : Wed Oct 24 2007 - 09:14:00 MDT