load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" begin fixhgt = (/30.,60.,90./) nlvmax = 3 ;----------------------------------------------------------------------------------- ; Read WRF Files ;---------------------------------------------------------------------------------- f = addfile("/simulation_data/ktm1month/april/wrfout_d03_2012-04-01_00:00:00.nc","r") UV = wrf_user_getvar(f,"uvmet" ,-1) ;m/s Horizontal Z = wrf_user_getvar(f,"height", -1) ;m Model Height Zter = wrf_user_getvar(f,"ter" , -1) ;m Terrain Height wtim = wrf_user_getvar(f,"times" , -1) ;Time nvert = dimsizes(fixhgt) u=UV(0,:,:,:,:) v=UV(1,:,:,:,:) printVarSummary(u) printVarSummary(v) ;;;;;;;;;;;;;;;;;;;;;;;;defining new variables for interpolation or u and v;;;;;;;;;;;;;;;;;;;;;;;; ndim_u = dimsizes(u) u_int=new((/ndim_u(0),nvert,ndim_u(2),ndim_u(3)/),typeof(u)) ndim_v = dimsizes(v) v_int=new((/ndim_v(0),nvert,ndim_v(2),ndim_v(3)/),typeof(v)) ;u_int=u ;v_int=v ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;Change Reference Height Zcat = Z(:,0:nlvmax-1,:,:) do kk = 0, nlvmax-1 Zcat(:,kk,:,:) = Zcat(:,kk,:,:) - Zter end do ; print(Zcat(0,:,25,46)) ;--------------------interpolating u and v component--------------------------------------------------------------- printVarSummary(u_int) printVarSummary(v_int) do kk = 0, nlvmax-1 u_int(:,kk,:,:) = wrf_user_intrp3d( u(:,kk,:,:),Zcat,"h", fixhgt(kk),0.,False) v_int(:,kk,:,:) = wrf_user_intrp3d( v(:,kk,:,:),Zcat,"h", fixhgt(kk),0.,False) end do ; printVarSummary(u_int) ; printVarSummary(v_int) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;calculation of wind;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; wspeed=sqrt(u^2+v^2) wspeed_int=sqrt(u_int^2+v_int^2) ;------------------------------------------interpolation of wind speed----------------------------------------------- ; int_wspeed=wspeed ; do kk=0, nlvmax-1 ; int_wspeed(:,kk,:,:)=wrf_user_intrp3d(wspeed,Zcat,"h",fixhgt(kk),0.,False) ; end do ; printVarSummary(wspeed) ; printVarSummary(int_wspeed) ; printVarSummary(wspeed_int) ;print(" U-component interpolated-U ") ; print (sprintf("%15.4f",int_wspeed(0:5,1,23,24)) +" "+sprintf("%15.4f", wspeed_int(0:5,1,23,24))) end