计算大气垂直运动引起的绝热变化的问题

LLL船猫 · 发表于 2024-6-6 21:00:33

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

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

x

本帖最后由 LLL船猫于 2024-6-6 21:04 编辑

公式如图，下面是我的脚本
; ;=============================计算Ver的数据=====================
f4 = addfile("h:/data1/LFCA/era5_t_w_800hPa_3level_1940-2022_q_1.nc", "r")
t=short2flt(f4->t(:,:,:,:))
w=short2flt(f4->w(:,:,:,:));pa/s
printVarSummary(w)
t_jja=new((/83,3,181,360/),"float")
w_jja=new((/83,3,181,360/),"float")
do j=0,82
      t_jja(j,:,:,:)=dim_avg_n(t(j*12+5:j*12+7,:,:,:), 0)；夏季平均
      w_jja(j,:,:,:)=dim_avg_n(w(j*12+5:j*12+7,:,:,:), 0)
end do

copy_VarMeta(t(0,:,:,:), t_jja(0,:,:,:))
copy_VarMeta(w(0,:,:,:), w_jja(0,:,:,:))

;===================干空气温度递减率========================
;===================RT/CpP==================================
dTdP = new((/83,3,181,360/),typeof(t_jja),t_jja@_FillValue)；3level  825hPa 850hPa 875hPa
Cp=1004.0;J/kg K
P=85000.0;Pa
R=287.0;J/K kg
dP=2500.0;Pa
dTdP(:,:,:,:) = center_finite_diff_n(t_jja(:,:,:,:),dP,False,0,1)
printVarSummary(dTdP)

Ver=new((/83,181,360/),"float")
do i=0,180
do j=0,360-1
Ver(:,i,j)=w_jja(:,1,i,j)*((R*t_jja(:,1,i,j)/(Cp*P))-dTdP(:,1,i,j))
end do
end do
copy_VarMeta(t(0,0,:,:), Ver(0,:,:))
printVarSummary(Ver)

  fnc       = "./era5_850hPa垂直运动的绝热加热_1940-2022_q_1.nc"
  system("/bin/rm -f "+fnc)             ; rm any pre-existing file
  fout       = addfile(fnc, "c") ; new netCDF file
  fout->Ver = Ver
  exit

然后我用这个数据画出的图比较怪，我看文献中有大面积的正异常和负异常，我这个就一小块一小块的，有没有大佬帮忙看看脚本哪里错了
我画的图如下