用NCL已获取了10m高度处定点风速时程图。如何获取10米高度处风向时程图？？

YJane

用NCL已获取了10m高度处定点风速时程图。如何获取10米高度处风向时程图？？

opeveu

楼主这个问题解决了吗

YJane

opeveu 发表于 2019-12-2 15:35
楼主这个问题解决了吗

解决了，后来发现用matlab读取就好了，很方便

opeveu

YJane 发表于 2019-12-12 16:55
解决了，后来发现用matlab读取就好了，很方便

楼主提取wrfout里面的10m高度处的风速风向的脚本还有存留吗？望赐教

YJane

opeveu 发表于 2019-12-13 08:07
楼主提取wrfout里面的10m高度处的风速风向的脚本还有存留吗？望赐教

这个脚本可能需要你修改一下，你试试吧。真心建议见用matlab处理u10,v10；方便的很。


load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

a  = addfile("/home/yjane/B/0812/YW-YS-KF/wrfout_d01_2008-08-19_12:00:00.nc","r")


  type = "pdf"

  wks = gsn_open_wks(type,"A")


  FirstTime = True
  times  = wrf_user_list_times(a)  ; get times in the file
  ntimes = dimsizes(times)         ; number of times in the file


   it = 8
;  do it = 0,ntimes-1,2             ; TIME LOOP

    print("Working on time: " + times(it) )



    u  = wrf_user_getvar(a,"U",it)        ; u averaged to mass points
    v  = wrf_user_getvar(a,"V",it)        ; v averaged to mass points
    w  = wrf_user_getvar(a,"W",it)
    p  = wrf_user_getvar(a, "pressure",it) ; pressure is our vertical coordinate
    z  = wrf_user_getvar(a, "z",it)        ; grid point height




      ip_lats =   22.59

      ip_lons =  114.04





        loc = wrf_user_ll_to_ij(a, ip_lons, ip_lats, True)
        locX = loc(0)
        locY = loc(1)


     zz = z(0:35,locY,locX)
     uu = u(0:35,locY,locX)
     vv = v(0:35,locY,locX)
     ww = w(0:36,locY,locX)
     
     www = sqrt(uu*uu+vv*vv)



npts = 35
i = ispan(0,npts,1)
x = zz
y = www
write_table("114.04_22.59_VO.txt","w",[/i,x,y/],"%4.0i,%8.3f,%8.3f")
printVarSummary(w)
printVarSummary(v)
printVarSummary(u)
                                    
;     plot                       = gsn_csm_xy(wks,ww,zz,True)





end

opeveu

YJane 发表于 2019-12-22 21:17
这个脚本可能需要你修改一下，你试试吧。真心建议见用matlab处理u10,v10；方便的很。

万分感谢楼主的分享

shenyi

请问楼主，我用你这个脚本，只改了经纬度，为何显示这呢
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading z
fatal:["Execute.c":8637]:Execute: Error occurred at or near line 50 in file wrf_meteo_42.ncl

zly4814624

shenyi 发表于 2020-1-14 15:37
请问楼主，我用你这个脚本，只改了经纬度，为何显示这呢
fatal:Subscript out of range, error in subscri ...

错误提示很明显，维数超出了范围。

shenyi

zly4814624 发表于 2020-1-14 15:50
错误提示很明显，维数超出了范围。

我用两个不同范围大小的数据都提示这个，维数超过范围怎么办呢，就是重新选范围跑数据，怎么确定不超范围呢

zly4814624

shenyi 发表于 2020-1-14 16:38
我用两个不同范围大小的数据都提示这个，维数超过范围怎么办呢，就是重新选范围跑数据，怎么确定不超范围 ...

你可以把数组的各个维的维数读取一下，就知道每一维有多少了。