Re: ncl vertical interpolation

This is being forwarded to wrfhelp@ucar.edu



On Mon, Jul 7, 2014 at 2:08 AM, sk <sdk0084@yahoo.com> wrote:

> Hi All,
> I am using wrf_user_intrp3d to create cross sections through two points as
> follow:
> ....
> plane = (/ xs,ys, xe,ye /)
> .....
> u_plane = wrf_user_intrp3d( ua,z,"v",plane,0.,True)
> v_plane = wrf_user_intrp3d( va,z,"v",plane,0.,True)
> t_plane = wrf_user_intrp3d( ta,z,"v",plane,0.,True)
> zz = wrf_user_intrp3d( z,z,"v",plane,0., True )
> b = ind(zz(0,:,0) .gt. zmax )
> zmax_pos = b(0) - 1
> if ( abs(zz(0,zmax_pos,0)-zmax) .lt. abs(zz(0,zmax_pos+1,0)-zmax) ) then
> zspan = b(0) - 1
> else
> zspan = b(0)
> end if
> delete([/zz,b/])
> ter_plane = wrf_user_intrp2d(ter,plane,0.,True)
> x_plane = wrf_user_intrp2d(lon,plane,0.,True)
> y_plane = wrf_user_intrp2d(lat,plane,0.,True)
> ....
> plot1 = gsn_csm_contour(wks,vcross(0,0:zmax_pos,:),res) ;;cross section
> at 06 UTC
>
> however, interpolation doesn't follow the terrain (attached file, the red
> line) and in general the interpolation start much higher in the right part
> of the cross section about 200 m even that the first model level is at
> about 12m.
> In otherwise the contours looks correct. Any idea how to improve the plots
> and make them look better?
>
> Regards,
> SK
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

Received on Mon Jul 07 02:02:36 2014