短期气候预测实习报告一

南京信息工程大学实验 实习 报告南京信息工程大学实验 实习 报告 实验 实习 名称 短期气候预测实习 日期 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 1948 2008 年 1 12 月的 500hPa 月平均高度场资料 范围 90 S 90 N 0 360 E 网格距 2 5 2 5 纬向格点数 144 经向格点数 73 GRD 格式 资料从南到北 自西向东排列 每月为一个记录 按年逐月排放 实习方法 实习方法 计算 1948 2007 1 月平均高度场 绘平均图 i j 为经 纬格点 m 为样本长度 k 为年份 计算 2008 年 1 月的高度距平 绘距平图 i j 为经 纬格点 m 为纬圈格点数 k 为年份 计算 2008 年 1 月的高度场纬偏值 绘纬偏图 1 1 jiH m jiH m k k 1 1 jiH m jiH m k k 1 1 m kkk k HA i jHi jH i j m i j 为经 纬格点 m 为样本长度 k 为年份 实习步骤实习步骤 编写计算 1948 2007 年 共 60 年 1 月 500hPa 高度平均程序 编写计算 2008 年 1 月 500hPa 高度场距平程序 相对于 1948 2007 年共 60 年的 平均 程序 编写 2008 年 1 月 500hPa 高度场纬向平均和偏差程序 根据 结果数据编写 GRADS 数据描述文件及执行文件 并绘图 program EX1 real a 144 73 12 61 ave1 144 73 ave7 144 73 a2008 144 73 real aalt 73 p2008 144 73 integer i j imo it open 2 file hgt500 grd form binary 补充正确路径 open 4 file ave60 grd form binary 补充数据输出路径 do it 1 61 do imo 1 12 do j 1 73 do i 1 144 read 2 a i j imo it enddo enddo enddo enddo cccccccccccccc 求平均 do i 1 144 do j 1 73 do it 1 60 ave1 i j ave1 i j a i j 1 it ave7 i j ave7 i j a i j 7 it enddo ave1 i j ave1 i j 60 0 ave7 i j ave7 i j 60 0 enddo enddo do j 1 73 do i 1 144 write 4 ave1 i j enddo enddo do j 1 73 do i 1 144 write 4 ave7 i j enddo enddo 1 1 jiH m jiHjiHLA m i kkk CCCCCC2008 距平 do j 1 73 do i 1 144 a2008 i j a i j 1 61 ave1 i j end do end do open 5 file a2008 grd form binary do j 1 73 do i 1 144 write 5 a2008 i j enddo enddo CCCCCCC 偏距平 do j 1 73 do i 1 144 aalt j aalt j a2008 i j 144 0 end do end do do j 1 73 do i 1 144 partial i j a2008 i j aalt j end do end do open 6 file p2008 grd form binary do j 1 73 do i 1 144 write 6 p2008 i j enddo enddo end Ctl 文件 其他的类似 不做重复 文件 其他的类似 不做重复 dset e copy a2008 grd undef 9 99E 33 title monthly mean hgt sta value xdef 144 linear 0 000 2 500 ydef 73 linear 90 000 2 500 zdef 1 levels 500 tdef 1 linear jan1948 1mo vars 1 ab 1 999 hgt abnormal 2008 endvars 画图 gs 文件 其他类似 不做重复 reinit open e copy ave60 ctl open e copy a2008 ctl open e copy partial ctl enable print e copy ave60 gmf set grads off set xlopts 1 1 0 15 set ylopts 1 1 0 13 set gxout shaded d ave1 1 set gxout contour d ave1 1 draw title Jan average 500hPa height print c set grads off set xlopts 1 1 0 15 set ylopts 1 1 0 13 set gxout shaded d ave7 1 set gxout contour d ave7 1 draw title Jul average 500hPa height print c set grads off set xlopts 1 1 0 15 set ylopts 1 1 0 13 set gxout shaded d ab 2 set gxout contour d ab 2 draw title 2008 Jan 500hPa anormal height print c set grads off set xlopts 1 1 0 15 set ylopts 1 1 0 13 set gxout shaded d pa 3 set gxout contour d pa 3 draw title 2008 Jan 500hPa partial height print c disable print reinit 实习结果分析 1 计算 1948 2007 年 60 年 1 月的平均高度场 绘制环流平均图 60 年冬季全球 500hPa 环流平均场 从全球 500hPa 的高度场中看出 位势梯度在低纬地区较小 在中高纬地区较大 由于北半 球的海陆分布特点 北半球的位势呈现冬季三槽三脊的特点可以明显的表现出来 而南半球 受陆地较小 海洋热惯性较大的特点 分布较为平直 60 年夏季全球 500hPa 环流平均场 从全球位势梯度分布来看 南大北少 低纬小 中高纬大 最高的位势位于加勒比海地区 由于太阳直射北半球 北半球的受热增多 气块膨胀 平均位势高度均大于南半球 槽脊分 布来看 北半球由弱的四槽四脊 2 计算 2008 年 1 月的高度距平 绘制高度距平场图 相对于 1948 2007 年共 60 年的平均 2008 年 1 月 500hPa 距平场 500hPa 北大西洋和欧洲大部分地区由强的正距平中心 东亚由一个较小的正距平中心 北 大西洋和北极圈地区有强的副局评中心 青藏高原有弱的距平中心 而南半球表现为低纬地 区的四个