短期气候预测实习报告一.doc

南京信息工程大学实验（实习）报告实验（实习）名称 短期气候预测实习 日期10.29得分 指导教师 系 专业 年级 班次 姓名 Trichtu 学号 实习目的：1掌握大气环流基本状况的表征平均、距平、纬偏图2大气环流基本状况分析3熟悉了解：资料格式FORTRAN语言4 GRADS绘图（数据描述文件及执行程序实习内容：(1)计算1948-2007年（60年）1月的平均高度场，绘制环流平均图；(2)计算2008年1月的高度距平，绘制高度距平场图（相对于1948-2007年共60年的平均）；(3)计算2008年1月的高度场纬偏值，绘制环流纬偏图。实习资料：NCEP/NCAR 19482008年112月的500hPa月平均高度场资料范围（90S-90N,0-360E)网格距2.52.5,纬向格点数144，经向格点数73GRD格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放。实习方法：计算（1948-2007）1月平均高度场，绘平均图；(i,j)为经、纬格点，m为样本长度，k为年份计算2008年1月的高度距平，绘距平图；(i,j)为经、纬格点，m为纬圈格点数，k为年份计算2008年1月的高度场纬偏值，绘纬偏图(i,j)为经、纬格点，m为样本长度，k为年份实习步骤编写计算1948-2007年（共60年）1月500hPa高度平均程序；编写计算2008年1月500hPa高度场距平程序（相对于1948-2007年共60年的平均）程序；编写2008年1月500hPa高度场纬向平均和偏差程序；根据结果数据编写GRADS数据描述文件及执行文件，并绘图。 program EX1real a(144,73,12,61),ave1(144,73),ave7(144,73),a2008(144,73)real aalt(73),p2008(144,73) integer i,j,imo,itopen(2,file=hgt500.grd,form=binary) !补充正确路径open(4,file=ave60.grd,form=binary) !补充数据输出路径 do it=1,61do imo=1,12do j=1,73do i=1,144read(2)a(i,j,imo,it)enddo;enddo;enddo;enddocccccccccccccc 求平均do i=1,144do j=1,73do it=1,60ave1(i,j)=ave1(i,j)+a(i,j,1,it)ave7(i,j)=ave7(i,j)+a(i,j,7,it)enddoave1(i,j)=ave1(i,j)/60.0ave7(i,j)=ave7(i,j)/60.0enddo;enddodo j=1,73do i=1,144 write(4)ave1(i,j)enddo;enddo do j=1,73do i=1,144 write(4)ave7(i,j)enddo;enddoCCCCCC2008 距平 do j=1,73do i=1,144 a2008(i,j)=a(i,j,1,61)-ave1(i,j)end doend do open(5,file=a2008.grd,form=binary) do j=1,73do i=1,144 write(5) a2008(i,j)enddoenddo CCCCCCC偏距平 do j=1,73do i=1,144 aalt(j)=aalt(j)+a2008(i,j)/144.0end do end do do j=1,73do i=1,144partial(i,j)=a2008(i,j)-aalt(j) end do end do open(6,file=p2008.grd,form=binary) do j=1,73do i=1,144 write(6) p2008(i,j)enddoenddoendCtl文件：（其他的类似，不做重复）dset e:copya2008.grdundef -9.99E+33title monthly mean hgt sta valuexdef 144 linear 0.000 2.500ydef 73 linear -90.000 2.500zdef 1 levels 500tdef 1 linear jan1948 1movars 1ab 1 -999 hgt abnormal 2008endvars;画图gs文件：（其他类似，不做重复）reinitopen e:copyave60.ctlopen e:copya2008.ctlopen e:copypartial.ctlenable print e:copyave60.gmfset grads offset xlopts 1 1 0.15set ylopts 1 1 0.13set gxout shadedd ave1.1set gxout contourd ave1.1draw title Jan average 500hPa heightprintcset grads offset xlopts 1 1 0.15set ylopts 1 1 0.13set gxout shadedd ave7.1set gxout contourd ave7.1draw title Jul average 500hPa heightprintcset grads offset xlopts 1 1 0.15set ylopts 1 1 0.13set gxout shadedd ab.2set gxout contourd ab.2draw title 2008 Jan 500hPa anormal heightprintcset grads offset xlopts 1 1 0.15set ylopts 1 1 0.13set gxout shadedd pa.3set gxout contourd pa.3draw title 2008 Jan 500hPa partial heightprintcdisable printreinit;实习结果分析(1) 计算1948-2007年（60年）1月的平均高度场，绘制环流平均图60年冬季全球500hPa环流平均场从全球500hPa的高度场中看出，位势梯度在低纬地区较小，在中高纬地区较大。由于北半球的海陆分布特点，北半球的位势呈现冬季三槽三脊的特点可以明显的表现出来。而南半球受陆地较小，海洋热惯性较大的特点，分布较为平直。60年夏季全球500hPa环流平均场从全球位势梯度分布来看，南大北少，低纬小，中高纬大。最高的位势位于加勒比海地区。由于太阳直射北半球，北半球的受热增多，气块膨胀，平均位势高度均大于南半球。槽脊分布来看，北半球由弱的四槽四脊。(2) 计算2008年1月的高度距平，绘制高度距平场图（相对于1948-2007年共60年的平均）2008年1月500hPa距平场500hPa，北大西洋和欧洲大部分地区由强的正距平中心，东亚由一个较小的正距平中心。北大西洋和北极圈地区有强的副局评中心，青藏高原有弱的距平中心。而南半球表现为低纬地区的四个正距平中心和高纬