NCL输出模式1000hpa的温度数据，为什么陆地上都没有数据呢

shushamei

NCL输出模式1000hpa的温度数据，为什么陆地上都没有数据呢 ，850hpa没有数据的区域主要是青藏高原地区，用什么方法解决呢

  tc = wrf_user_getvar(a,"tc",it)        ; T in C
    u  = wrf_user_getvar(a,"ua",it)        ; u averaged to mass points
    v  = wrf_user_getvar(a,"va",it)        ; v averaged to mass points
    p  = wrf_user_getvar(a, "pressure",it) ; pressure is our vertical coordinate
    z  = wrf_user_getvar(a, "z",it)        ; grid point height
    rh = wrf_user_getvar(a,"rh",it)        ; relative humidity



tc_plane = wrf_user_intrp3d(tc,p,"h",pressure,0.,False)
      z_plane  = wrf_user_intrp3d( z,p,"h",pressure,0.,False)
      rh_plane = wrf_user_intrp3d(rh,p,"h",pressure,0.,False)
      u_plane  = wrf_user_intrp3d( u,p,"h",pressure,0.,False)
      v_plane  = wrf_user_intrp3d( v,p,"h",pressure,0.,False)
      spd     = (u_plane*u_plane + v_plane*v_plane)^(0.5) ; m/sec
      spd@description = "Wind Speed"
      spd@units = "m/s"
      u_plane = u_plane*1.94386     ; kts
      v_plane = v_plane*1.94386     ; kts
      u_plane@units = "kts"
      v_plane@units = "kts"

hzwjy

脚本贴出来看看。

jacci_zou

你确定是全球的数据？而不是海温的？咋觉着这图看着那么像海温图勒，，但是在俄罗斯中部和阿留申一带又出现非海温数据！！！

freekiller

这就对了，1000hPa 在陆地上很多都是在地下的，所以陆地区域大部分是缺测值，画图自然是空白

shushamei

freekiller 发表于 2012-10-6 12:29

                               
登录/注册后可看大图

这就对了，1000hPa 在陆地上很多都是在地下的，所以陆地区域大部分是缺测值，画图自然是空白

这样子啊 那地面的数据怎么获取呢

shushamei

jacci_zou 发表于 2012-10-5 06:38

                               
登录/注册后可看大图

你确定是全球的数据？而不是海温的？咋觉着这图看着那么像海温图勒，，但是在俄罗斯中部和阿留申一带又出现 ...

你看看我贴的脚本吧

freekiller

shushamei 发表于 2012-10-15 12:10

                               
登录/注册后可看大图

这样子啊 那地面的数据怎么获取呢

画sigma层               

shushamei

freekiller 发表于 2012-10-15 20:55

                               
登录/注册后可看大图

画sigma层

怎么画呢  能说的详细些吗

freekiller

wrf模式对应sigma值，画模式最下面那一层

shushamei

freekiller 发表于 2012-11-21 20:17

                               
登录/注册后可看大图

wrf模式对应sigma值，画模式最下面那一层

你的意思是指eta=0的那一层吗 就是bottom_top=0的时候那一层吗