CESM模式资料后处理

麦田_smile

各位好，想请教一下cesm模式资料后处理问题，我的脚本如下，主要是采用了vinth2p函数

    fils     = systemfunc("ls " + ncdir + "ctrl.cam2.h0.*.nc")   
    f        = addfiles (fils, "r")   
    y       = f[:]->U
    printVarSummary(y)

    ff          = addfile(ncdir+"ctrl.cam2.h0.0001-01.nc","r")   
    pnew    = (/1000.,925.,850.,700.,600.,500.,400.,300.,250.,200.,150.,100.,70.,50.,30.,20.,10./)
    pnew@long_name = "pressure"
    pnew@units           = "hPa"
    hyam    = ff->hyam                     ; read to memory [optional]
    hybm    = ff->hybm
    P0mb    = ff->P0*0.01
    PS0       = ff->PS
    PS         = new((/240,96,144/),"float")
    do    i     = 0,239
      PS(i,:,:)= PS0(0,:,:)
    end do
    time       = ispan(0,239,1)
    time@units = "months since 1990-01-01"   
    PS!0    = "time"
    PS!1    = "lat"
    PS!2    = "lon"
    PS&time = time
    PS&lat  = y&lat
    PS&lon  = y&lon
    y&time  = time

    ynew      = vinth2p(y,hyam,hybm,pnew,PS,1,P0mb,1,True)
    ynew&time = time
    printVarSummary(ynew)
   
    fname = "uwnd.nc"
    system("/bin/rm -f "+ncdir+"output/"+fname)
    a = addfile(ncdir+"output/"+fname,"c")
    a->uwnd = ynew

但是得到的结果总是看着有点不太正常，尤其是青藏高原那里，(貌似高层就正常了)，难道是地形的影响？下图就是grads显示的850hPa纬向风场、位势高度场。请问这是什么原因？谢谢各位！

talkd

不仅是青藏高原，格陵兰岛也是，应该是地形因素，850hPa那里没有数据，模式输出的是虚假的值。

麦田_smile

talkd 发表于 2015-9-12 09:38
不仅是青藏高原，格陵兰岛也是，应该是地形因素，850hPa那里没有数据，模式输出的是虚假的值。

谢谢，您的意思就是本来就是如此？这样转化没有问题吧？

66kqkq

什么版本？我记得我用的模式结果会自动mask地形的

麦田_smile

66kqkq 发表于 2015-12-8 15:29
什么版本？我记得我用的模式结果会自动mask地形的

1.0的呢 版本比较老

rpf

我想请问楼主一个问题，这里面只需要将p0转换成hPa的单位，ps不需要转换吗？我的模式数据的ps单位也是Pa。

zhaofei

谢谢分享！！！

zhaofei

谢谢分享！！！