气象程序设计及绘图实习报告.pptxVIP

实习要求： 根据所提供的实习说明和资料，按要求完成实习内容； 将完成的实习内容按步骤写入实习报告中，包括资料说明、Fortran 程序、CTL 文 件、gs 文件、绘制的图形（要求在图形上方用 draw title 的方式添加标题，注明姓 名拼音及学号） 根据所绘图形及查找相关文献，说明本次实习所求物理量场的意义。 实习内容： Nino34 海温指数与海平面气压相关图 Nino34 海温指数与中国 160 站降水相关图 绘制 1 月份 Nino34 海温指数与 1 月份海平面气压相关图 资料说明 Nino34 指数：Nino34.txt 是 1951 年 1 月至 2013 年 12 月 Nino3.4 区（热带太平洋： 西经 170 度~西经 120 度、北纬 5 度~南纬 5 度）区域平均海温指数资料，共有 63 行数据，每行数据第一个数字为年份，后面 12 个数字为该年 1-12 月的海温指数； 全球海平面气压月平均数据：slp.jan.grd 为 1951 至 2013 年的 1 月全球海平面气压 场二进制数据，年数为 63 年，格点数 144*73，水平分辨率 2.5?*2.5?； 步骤： 1 用 Fortran 编写 corr.grid.f90 文件，计算 1951-2013 年 1 月 Nino34 海温指数与海平面气压 相关系数，计算结果保存于corr.slp.grd 文件中； corr.grid.f90 文件：;relation=sfenzi/sqrt(sfenmu1*sfenmu2) end function relation 保存于 corr.slp.grd 文件中。 书写 corr.slp.grd 的数据描述文件 1.ctl： dset e:\gradsimo\test1\corr.slp.grd title slp undef -9.99e33 xdef 144 linear -180 2.5 ydef 73 linear -90 2.5 zdef 1 levels 1000 tdef 63 linear jun1951 1yr vars 1 slp 0 0 sea level pressure endvars 编写 1.gs 文件： reinit open e:\gradsimo\test1\1.ctl set lev 1000 set lon -180 180 set lat -90 90 set t 1 set gxout shaded set grads off set black -0.2 0.2 d slp set gxout contour d slp draw title Lu Yu 20131301070 cbar 1.0 0 printim e:\gradsimo\test1\1.png white ; 绘制 1951-2013 年 1 月 Nino34 海温指数与海平面气压相关系数图：;绘制 1 月份 Nino34 海温指数与 7 月份中国 160 站降水相关图： 资料说明 Nino34 指数：Nino34.txt 是 1951 年 1 月至 2013 年 12 月 Nino3.4 区（热带太平洋： 西经 170 度~西经 120 度、北纬 5 度~南纬 5 度）区域平均海温指数资料，共有 63 行数据，每行数据第一个数字为年份，后面 12 个数字为该年 1-12 月的海温指数； 160 站降水资料：r1607.txt 是中国 160 站 1951~2013 年共 63 年 7 月降水资料，数 据是 1951 年 7 月各站降水量按照 id.txt 中各站顺序写入，该年 160 站都写入后再写 下 1952 年 7 月各站降水量，...................，直至 2013 年 7 月各站降水量按照 id.txt 中各站顺序写入； 160 站站点资料：id.txt 是 160 站站号和纬度、经度资料，共有 160 行。每行第一个 数据为 8 位字符型的站号，后面两个数字分别为该站点纬度和经度； 步骤： 用 Fortran 编写 corr.station.f90 文件，计算 1951-2013 年 1 月 Nino34 海温指数与 7 月 160 站 降水相关系数，计算结果保存于 corr.7.txt 和corr.7.grd 文件中；;relation=sfenzi/sqrt(sfenmu1*sfenmu2) end function relation;准备 GrADS 画站点图时插值时所需格点数据 Grid.grd 和及其描述文件 Grid161.ctl（已提 供）； 编写 2.gs 文