绘制HCHO3个分布图

最新推荐文章于 2024-08-21 15:43:29 发布

阿水水水水水

最新推荐文章于 2024-08-21 15:43:29 发布

阅读量382

点赞数

分类专栏：脚本积累文章标签：脚本语言

本文链接：https://blog.youkuaiyun.com/qq_40256777/article/details/124504712

版权

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/csm/contributed.ncl"
f=addfile("/data1/public/loucx/hcho_resample_out.nc","r")
m=addfile("/data1/public/loucx/outtest.nc","r")
lat1=m->XLAT ;cmaq
lon1=m->XLONG
lat2=f->latitude(0,:,:) ;tropomi
lon2=f->longitude(0,:,:)
hcho1=m->HCHO
hcho2=f->HCHO_column(0,:,:)  ;(/174,168/)
var_0=hcho1
hcho=hcho2*6.02214*10000
var_1=new((/174,168/),float)
var_2=new((/174,168/),float)
do i=0,173
 do j=0,167
 temp1=hcho(i,j)
if (temp1 .gt.0) then
  var_1(i,j)=hcho(i,j)
  var_2(i,j)=var_0(i,j)-var_1(i,j)
end if
end do
end do
;--------plot
 txid_tr = new(3,graphic)
 amid_tr = new(3,graphic)
 txres                       = True
 txres@txPerimOn             = False ;True
 txres@txFontHeightF         = 0.028
 amres_tr                  = True
 amres_tr@amParallelPosF   =  0.48    ; This is the right edge of the plot.
 amres_tr@amOrthogonalP