爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 13170|回复: 27

wrf云顶温度

[复制链接]
发表于 2015-4-22 21:51:26 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
wrf得到的云顶温度为何都小于0
QQ图片20150422215219.png
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-4-22 21:53:58 | 显示全部楼层
我的程序是这样的
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"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"

  begin

  a = addfile("/public3/home/interior/careeri/huser037/WRF/WRFV3S/run/wrfout_d01_2010-06-29_00:00:00","r")

  wks = gsn_open_wks("pdf","test")

; lake=wrf_user_getvar(a,"LAKEMASK",1)
  lat = a->XLAT(1,:,:)   ;
  lon = a->XLONG(1,:,:)
; lon = wrf_user_getvar(a,"lon",1)
;time=ispan(16,264,1)
;print(time)
  actt = wrf_user_getvar(a,"ctt",(/48,49/))   ; U component of the wind at 10m
  printVarSummary(actt)
  avg_ctt=dim_avg_n_Wrap(actt,0)
; Plotting options for Wind Vectors
      res = True                        
      res@cnFillOn = True  
      res@pmLabelBarOrthogonalPosF = -0.1
    ;  res@ContourParameters = (/ 0.1 /)
    ;  res@cnFillColors = (/"White","Chartreuse"/)
      contour = wrf_contour(a,wks,avg_ctt,res)
      delete(res)
  
     pltres = True                              ; Set plot options
     mpres = True                               ; Set map options
     mpres@mpDataBaseVersion="MediumRes"
     mpres@mpDataSetName="Earth..4"
     mpres@mpOutlineSpecifiers=(/"China:states"/)

    ; MAKE PLOTS                                       
      plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)

frame(wks)


  end
密码修改失败请联系微信:mofangbao
发表于 2015-4-24 10:14:15 | 显示全部楼层
zhangxiao 发表于 2015-4-22 21:53
我的程序是这样的
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
  load "$NCARG_ROOT/lib ...

你好!请问wrf输出量中哪个是云顶温度?谢谢
密码修改失败请联系微信:mofangbao
发表于 2015-4-24 10:51:31 | 显示全部楼层
请问楼主这个云顶温度是修改后处理用ncl得到的?
密码修改失败请联系微信:mofangbao
发表于 2015-4-24 19:16:21 | 显示全部楼层
同问?ncl里用wrf_user_getvar,ctt就是云顶亮温?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-4-24 21:54:22 | 显示全部楼层
也不是很清楚 https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml 提取的"ctt"难道不是云顶温度吗
密码修改失败请联系微信:mofangbao
发表于 2015-4-27 08:36:20 | 显示全部楼层
zhangxiao 发表于 2015-4-24 21:54
也不是很清楚 https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml 提取的"ctt"难 ...

云顶温度的缩写肯定是这样的,是的
密码修改失败请联系微信:mofangbao
发表于 2015-4-28 20:57:52 | 显示全部楼层
请问楼主,这个脚本使用必须要用NCL6.3.0的吗
密码修改失败请联系微信:mofangbao
发表于 2015-4-28 22:29:28 | 显示全部楼层
请问楼主知道为什么‘云顶温度’都是负值了么?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-5-3 12:24:16 | 显示全部楼层
不知道,有知道的请讲解,不胜感激
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表