WRF后处理得到某个高度的值

roc

WRF模式计算完成后 得到wrfout_d01_2012-11-29_00:00:00

内容如下：
netcdf wrfout_d01 {
        dimensions:
            Time = UNLIMITED ; // (25 currently)
                             DateStrLen = 19 ;
                             west_east = 9 ;
                             south_north = 9 ;
                             bottom_top = 33 ;
                             bottom_top_stag = 34 ;
                             soil_layers_stag = 4 ;
                             west_east_stag = 10 ;
                             force_layers = 8 ;
                             south_north_stag = 10 ;
variables:
                             float U(Time, bottom_top, south_north, west_east_stag) ;
                                                  U:FieldType = 104 ;
                                                  U:MemoryOrder = "XYZ" ;
                                                  U:description = "x-wind component" ;
                                                  U:units = "m s-1" ;
                                                  U:stagger = "X" ;
                                                  U:coordinates = "XLONG_U XLAT_U" ;
                ...

data:
                     U =
            1.229847, 1.296391, 1.247952, 1.141632, 1.09145, 1.057645, 1.027596,
            1.002366, 0.9899619, 0.9751331,
……











这样的文件后 我想得到某个高度的风速度 温度等等

不知道怎么样通过bottom_top得到我想要的高度值 比如 100m高度时某个时间的U风速值

bottom_top与高度的对应关系是什么
求大神们解惑

mofangbao

看你用什么工具啊，你连用什么工具都没说

roc

后处理是用NCL
想把某个高度的值 最后输出到文本文件中

roc

mofangbao 发表于 2012-12-12 18:40

                               
登录/注册后可看大图

看你用什么工具啊，你连用什么工具都没说

NCL的部分代码如下：
……

FILE_NC = addfile("./wrfout_temp.nc","r")

……
    U  = wrf_user_getvar(FILE_NC,"U",-1)     
   
    loc = wrf_user_ij_to_ll(FILE_NC,tmp,tmp,res)

    u = new((/nt,nz,ny,nx/),float)
        do k = 0,nt-1
    do j = ystart,ystart+ny-1
    do i = xstart,xstart+nx-1
    u(k,0:0,j-ystart,i-xstart)   = (/u10(k,j,i)/)
        u(k,3:3,j-ystart,i-xstart)   = (/U(k,5:5,j,i)/)
        end do
    end do

    end do  

    h = (/10,70/)
  npt   = -1
  do t=0,nt-1         
    do k=0,nz-1
      do j=0,ny-1         
        do i=0,nx-1

       npt  = npt + 1   
       data(npt) = sprintf("%10.2f", loc(0,i))  
       data(npt) = data(npt) + sprintf("%10.2f", loc(1,j))  
       data(npt) = data(npt) + sprinti("%5i", height(k))  
       data(npt) = data(npt) + sprinti("%10i", date)  
       data(npt) = data(npt) + sprinti("%5.4i", time(t))  
       data(npt) = data(npt) + sprintf("%10.2f ",u(t,k,j,i))
         end do
       end do
     end do
  end do

  asciiwrite (fName , data)

……

嗳づ左岸の兔

请问楼主问题解决了吗