资源描述:
《短期气候预测实习报告一》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、南京信息工程大学实验(实习)报告实验(实习)名称短期气候预测实习日期10.29得分指导教师系~专业~年级~班次~姓名Trichtu学号~实习目的:1掌握大气环流基本状况的表征--平均、距平、纬偏图2大气环流基本状况分析3熟悉了解:资料格式FORTRAN语言4GRADS绘图(数据描述文件及执行程序实习内容:(1)计算1948-2007年(60年)1月的平均高度场,绘制环流平均图;(2)计算2008年1月的高度距平,绘制高度距平场图(相对于1948-2007年共60年的平均);(3)计算2008年1月的高度场纬偏值,绘制环流纬偏图。实习资料:NCEP/N
2、CAR1948-2008年1~12月的500hPa月平均高度场资料范围(90°S-90°N,0-360°E)网格距2.5°×2.5°,纬向格点数144,经向格点数73GRD格式,资料从南到北、自西向东排列,每月为一个记录,按年逐月排放。实习方法:①计算(1948-2007)1月平均高度场,绘平均图;(i,j)为经、纬格点,m为样本长度,k为年份②计算2008年1月的高度距平,绘距平图;(i,j)为经、纬格点,m为纬圈格点数,k为年份③计算2008年1月的高度场纬偏值,绘纬偏图(i,j)为经、纬格点,m为样本长度,k为年份实习步骤①编写计算1948-2
3、007年(共60年)1月500hPa高度平均程序;②编写计算2008年1月500hPa高度场距平程序(相对于1948-2007年共60年的平均)程序;③编写2008年1月500hPa高度场纬向平均和偏差程序;④根据①②③结果数据编写GRADS数据描述文件及执行文件,并绘图。programEX1reala(144,73,12,61),ave1(144,73),ave7(144,73),a2008(144,73)realaalt(73),p2008(144,73)integeri,j,imo,itopen(2,file='hgt500.grd',form
4、='binary')!补充正确路径open(4,file='ave60.grd',form='binary')!补充数据输出路径doit=1,61doimo=1,12doj=1,73doi=1,144read(2)a(i,j,imo,it)enddo;enddo;enddo;enddocccccccccccccc求平均doi=1,144doj=1,73doit=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.0a
5、ve7(i,j)=ave7(i,j)/60.0`enddo;enddodoj=1,73doi=1,144write(4)ave1(i,j)enddo;enddodoj=1,73doi=1,144write(4)ave7(i,j)enddo;enddoCCCCCC2008距平doj=1,73doi=1,144a2008(i,j)=a(i,j,1,61)-ave1(i,j)enddoenddoopen(5,file='a2008.grd',form='binary')doj=1,73doi=1,144write(5)a2008(i,j)enddoendd
6、oCCCCCCC偏距平doj=1,73doi=1,144aalt(j)=aalt(j)+a2008(i,j)/144.0enddoenddodoj=1,73doi=1,144partial(i,j)=a2008(i,j)-aalt(j)enddoenddoopen(6,file='p2008.grd',form='binary')doj=1,73doi=1,144write(6)p2008(i,j)enddoenddoendCtl文件:(其他的类似,不做重复)dsete:copya2008.grdundef-9.99E+33titlemonthl
7、ymeanhgtstavaluexdef144linear0.0002.500ydef73linear-90.0002.500zdef1levels500tdef1linearjan19481movars1ab1-999hgtabnormal2008endvars;画图gs文件:(其他类似,不做重复)'reinit''opene:copyave60.ctl''opene:copya2008.ctl''opene:copypartial.ctl''enableprinte:copyave60.gmf''setgradsoff''setx
8、lopts110.15''setylopts110.13''setgxoutshaded''dave1.1'