气候诊断及绘图实用技术讲座.ppt

第二步 计算相关 利用Fortran 第三步 CTL文件的编写 dset d:/2010/sample2/hgt09rain.grd undef 9999.99 title xgxs xdef 144 linear 0 2.5 ydef 73 linear -90 2.5 zdef 1 linear 1 1 tdef 1 linear jan2009 1mo vars 1 corr 0 99 xgxs endvars 第四步 gs文件的编写 reinit open d:/2010/sample2/corr.ctl‘ set mproj nps set mpdset mres set lon -90 270 set lat 0 90’ d corr’ enable print d:/2010/sample2/hgt09rain.gmf print disable print ; 增加显著性检验及图的修饰 set black -0.2732 0.2732 set gxout shaded set clevs -0.3541 -0.2732 0.2732 0.3541 set ccols 7 11 0 11 7 d corr (1)增加显著性检验区 (2) 图的修饰 两种图的叠加显示，线条的粗细、 线型、阴影区的色彩选择等。 set cthick 10 ‘ set gxout contour set cint 0.1 set cmin 0 d corr set cmax -0.00001 set cstyle 2 d corr 例3------合成场及差值场 第一步 典型年的选取 第二步 gs文件的编写 reinit sdfopen d:/2010/month/hgt.mon.mean.nc set lon 40 180 set lat 0 70‘ set lev 500‘ 计算某时间段平均(如1971-2000年) set t 282 define h5007106=hgt set t 294 define h5007206=hgt set t 306’ …….. define h50006=(h5007106+….+h5000006)/30 典型年合成 （多雨年) set t 6 define h5004806=hgt‘ define h500q=(h5004806+h5007906+ h5008806+h5009906)/4 （少雨年) set t 642 define h5000106=hgt‘ define h500r=(h5008706+h5009706+ h5009906+h5000106)/4 方差计算 （多雨年) define s5004806=(h5004806-h500q)*(h5004806-h500q) define s5007906=(h5007906-h500q)*(h5007906-h500q) define s5008806=(h5008806-h500q)*(h5008806-h500q) define s5009906=(h5009906-h500q)*(h5009906-h500q)‘ define s500q=(s5004806+s5007906+s5008806+s5009906)/4 （少雨年) define s5008706=(h5008706-h500r)*(h5008706-h500r) define s5009706=(h5009706-h500r)*(h5009706-h500r) define s5009906=(h5009906-h500r)*(h5009906-h500r) define s5000106=(h5000106-h500r)*(h5000106-h500r)‘ define s500r=(s5008706+s5009706+s5009906+s5000106)/4 T检验 魏凤英《现代气候统计诊断 与预测技术》构造两个总体的均值有无显著差异的统计量. define t1=h500q-h500r define t2=sqrt(((4-1)*s500q+(4-1)*s500r)/(8-2)) define t3=sqrt(1/4+1/4) define t=t1/(t2*t3) 修改部分语句，绘制不同的图 （1）差值场及t检验； （2）距平场; （3）平均场; 例4------剖面图 reinit ‘sdfopen j:/2010/day/vwnd.1999.nc