ncl vertical interpolation

Hi All, I am using wrf_user_intrp3d to create cross sections through two =
points as follow: .... plane = (/ xs,ys, xe,ye /) .....  =

ane = wrf_user_intrp3d( va,z,"v",plane,0.,True)   t_plane = wrf_use=
r_intrp3d( ta,z,"v",plane,0.,True)   zz = wrf_user_intrp3d( z,z,"v",p=
lane,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)-zm=
ax) ) then             zspan = b(0) - 1   =
        else             zspan = b(=
0)           end if   delete([/zz,b/])   ter_pl=
ane = wrf_user_intrp2d(ter,plane,0.,True)   x_plane   = wrf_use=
r_intrp2d(lon,plane,0.,True)   y_plane   = wrf_user_intrp2d(lat,p=
lane,0.,True) ....  plot1 = gsn_csm_contour(wks,vcross(0,0:zmax_pos=
,:),res)   ;;cross section at 06 UTC   however, interpolation doe=
sn't follow the terrain (attached file, the red line) and in general the in=
terpolation 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 c=
ontours 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