留言板

shangyuaermod 2019-5-27 15:20: 前辈，你好！我是浙大能院的小徐同学，最近在用空气质量模型AERMOD模拟大气中污染物扩散，现在会跑出一个基本的结果，但关于AERMOD还有很多疑问，感觉对这个模型整体掌握的还很粗糙。由于身边的师兄师姐同学都不是做这个的，缺乏同行的交流，所以请问前辈研究过这个模型吗？或者前辈认识什么人或浙大其他院的同学做这个的吗？感谢

方自之 2016-11-8 09:09: 兄台你好，家园里搜有关浙大考博的帖子，搜索到你，请问兄台是否在浙大，想简单咨询点情况。先简单自我介绍：13年南信大大气科学本科，16年中国科学院大学大气硕士，研究方向偏陆面过程和WRF。静候回复

dcx 2014-9-18 14:54: 你好，你曾经在气象家园上下载的wrf结果画台风路径：先用exwsd.gs提取WRF输出结果的所需变量，再用readminslp.f程序计算出台风中心经纬度和最大地面风速和最低海平面气压，再用6htrack.gs根据经纬度画图。可以给我发一下吗，谢谢了。邮箱295079276@qq.com

freekiller 2014-6-13 09:41: 你在前面定义了p ，即p = short2flt(f4 ->msl(:,{20.25},{110.25}))
然后又使用 p = conform(u(time|:,lev|:),pressure_levels,1)
所以两者不相等了，可以把 p = conform(u(time|:,lev|:),pressure_levels,1) 改为其他
比如 p1 = conform(u(time|:,lev|:),pressure_levels,1)

freekiller 2014-6-13 09:41: 你在前面定义了p ，即p = short2flt(f4 ->msl(:,{20.25},{110.25}))
然后又使用 p = conform(u(time|:,lev|:),pressure_levels,1)
所以两者不相等了，可以把 p = conform(u(time|:,lev|:),pressure_levels,1) 改为其他
比如 p1 = conform(u(time|:,lev|:),pressure_levels,1)

freekiller 2014-6-13 09:40: 你在前面定义了p ，即p = short2flt(f4 ->msl(:,{20.25},{110.25}))
然后又使用 p = conform(u(time|:,lev|:),pressure_levels,1)
所以两者不相等了，可以把 p = conform(u(time|:,lev|:),pressure_levels,1) 改为其他
比如 p1 = conform(u(time|:,lev|:),pressure_levels,1)

seafairy 2014-6-12 17:04: 你好，百度了一下ncl的一个相关问题，看到提到你的，想跟你请教一个问题。conform的应用，感觉没错，但是还是提示了“fatal:Number of dimensions on right hand side do not match number of dimension in left hand side”

begin
  diri = "/mnt/hgfs/E/WUTIP/dat/"
  fili1 = "ERA.temp.0920-1005.nc"
  fili2 = "ERA.uv.0920-1005.nc"
  fili3 = "ERA.rh.0920-1005.nc"
  fili4 = "ERA.surf.0920-1005.nc"

  diro = "/mnt/hgfs/E/WUTIP/dat/pic/"
  filo = "tse_lat-time"

  f1  = addfile(diri+fili1,"r")
  f2  = addfile(diri+fili2,"r")
  f3  = addfile(diri+fili3,"r")
  f4  = addfile(diri+fili4,"r")
  tk = short2flt(f1 ->t(:,{1000:200},{20.25},{110.25}))
  u = short2flt(f2 ->u(:,{1000:200},{20.25},{110.25}))
  v = short2flt(f2 ->v(:,{1000:200},{20.25},{110.25}))
  rh = short2flt(f3 ->r(:,{1000:200},{20.25},{110.25}))
  p = short2flt(f4 ->msl(:,{20.25},{110.25}))
  u!0 ="time"
  u!1 ="lev"
  pressure_levels = (/1000.,975.,950.,925.,900.,875.,850.,825.,800.,775.,750.,700.,650.,600.,550.,500.,\
            450.,400.,350.,300.,250.,225.,200./)
  nlevels       = dimsizes(pressure_levels)    ; number of pressure levels
;计算假相当位温
p = conform(u(time|:,lev|:),pressure_levels,1)

能帮我看一下是什么问题吗？谢谢！

freekiller的个人空间 http://bbs.06climate.com/?243 [收藏] [复制] [分享] [RSS]