求助！画整层水汽通量及散度出错

supermann

我想画1961-2016年的全国整层水汽通量和散度的分布图，以下是部分程序，比湿只有8层，所以人为添加了基几层然后设置为零，但是还是出现错误，积分函数那里出错：ier=-999，不知哪位大神能指点迷津，困扰很久了，求解答





itimes = 196101
itimee = 201612

f1 = addfile("/cygdrive/D/paperxu/data/pres.sfc.mon.mean.nc","r")
f2 = addfile("/cygdrive/D/paperxu/data/uwnd.mon.mean.nc","r")
f3 = addfile("/cygdrive/D/paperxu/data/vwnd.mon.mean.nc","r")
f4 = addfile("/cygdrive/D/paperxu/data/shum.mon.mean.nc","r")

lat = f2->lat
lon = f2->lon
lev = f2->level
nlat = dimsizes(lat)
nlon = dimsizes(lon)
nlev = dimsizes(lev)

time = f1->time
YYYYMM = cd_calendar(time,-1)
ist = ind(itimes.eq.YYYYMM)
ied = ind(itimee.eq.YYYYMM)

ps = f1->pres(ist:ied,:,:)
ps = ps * 0.01
u = f2->uwnd(ist:ied,:,:,:)
v = f3->vwnd(ist:ied,:,:,:)
q = f4->shum(ist:ied,:,:,:)
q = q * 1000.0
;printVarSummary(u)
;printVarSummary(v)

ps1 = month_to_season(ps,"SON")
u1 = month_to_season(u,"SON")
v1 = month_to_season(v,"SON")
q1 = month_to_season(q,"SON")
;printVarSummary(u1)
;printVarSummary(v1)


;*********************************************
    linlog = 1
    p  = (/ 1000.,925.,850.,700.,600.,500., 400.,300.,250.,200.,150.,100., 70.,50.,30.,20.,10./)

    pbot   = 1100.                                             
    ptop   = 300.

    qu=q1*u1(:,{1000:300},:,:)
    qv=q1*v1(:,{1000:300},:,:)
    copy_VarCoords(u1, qu)
    copy_VarCoords(v1, qv)
    nyear = 56
    qu_new = new((/nyear,nlev,nlat,nlon/),float)
    qv_new = new((/nyear,nlev,nlat,nlon/),float)

    qu_new(:,0:7,:,:) = qu
    qu_new(:,8:16,:,:) = 0.0

    qv_new(:,0:7,:,:) = qv
    qv_new(:,8:16,:,:) = 0.0
    copy_VarCoords(u1, qu_new)
    copy_VarCoords(u1, qv_new)
    uu = vibeta (p,qu_new(time|:,lat|:,lon|:,level|:),linlog,ps1,pbot,ptop)/9.8  ; retLIuns u(time,lat,lon)
    vv = vibeta (p,qv_new(time|:,lat|:,lon|:,level|:),linlog,ps1,pbot,ptop)/9.8  

    uuclimate = dim_avg_n_Wrap(uu, 0)
    vvclimate = dim_avg_n_Wrap(vv, 0)
    copy_VarCoords(u1(0,0,:,:), uuclimate)
    copy_VarCoords(u1(0,0,:,:), vvclimate)

    Div = uv2dv_cfd(uu, vv, lat, lon, 2)
    Div = dim_rmvmean_n_Wrap(Div, 0)
    Div = Div*10e5
    copy_VarCoords(u1(:,0,:,:), Div)

言墨

积分时缺测值太多了，同一个点，垂直方向上至少有三个数据不是缺测，你可以把缺测值设为0，反之积分时不影响

Airy佐

请问楼主画出了整层水汽通量吗？
想问一下为什么要q = q * 1000.0和ps = ps * 0.01，
还有算整层的水汽通量一定要算这个嘛？
ps1 = month_to_season(ps,"SON")
u1 = month_to_season(u,"SON")
v1 = month_to_season(v,"SON")
q1 = month_to_season(q,"SON"){:eb303:}