气象程序设计及绘图实习报告
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
*12 c(63)
*7 a(63)
(5'D:\\\''')
(6'D:\\\34')
1,63
(6,'(A12)') c(i)
1,63
(c(i),'(A75.2)') a(i)(i)
1,63
(5) ((p()1,144)1,73)
*
1,73
1,144
1,63
z(k)()
r()(,63)
*
(7'D:\\\')
' -90 90'
' t 1'
' '
' '
' -0.2 0.2'
'd '
' '
'd '
' 20131301070'
' 1.0 0'
' e:\\20150706\1\1 '
;
绘制1951-2013年1月34海温指数与海平面气压相关系数图：
1951-2013年1月34海温指数与海平面气压相关系数图
2绘制1月份34海温指数与7月份中国160站降水相关图：
*12 c(63)
*7 a(63)
*8 (160)
(6'E:\\20150706\2\34')
(5'E:\\20150706\2\r1607')
(10:\\20150706\1\')
1,160
(10,*) (i)(i)(i)
*
1,63
(5,*) (p()1,160)
1,63
(6,'(A12)') c(i)
2.1资料说明
34指数：34是1951年1月至2013年12月3.4区（热带太平洋：西经170度~西经120度、北纬5度~南纬5度）区域平均海温指数资料，共有63行数据，每行数据第一个数字为年份，后面12个数字为该年1-12月的海温指数；
160站降水资料：r1607是中国160站1951~2013年共63年7月降水资料，数据是1951年7月各站降水量按照中各站顺序写入，该年160站都写入后再写下1952年7月各站降水量，...................，直至2013年7月各站降水量按照中各站顺序写入；
全球海平面气压月平均数据：为1951至2013年的1月全球海平面气压场二进制数据，年数为63年，格点数144*73，水平分辨率2.5︒*2.5︒；
1.2步骤：
1用编写90文件，计算1951-2013年1月34海温指数与海平面气压相关系数，计算结果保存于文件中；
90文件：
(63)(73,144,63)(73,144)(63)
''
' '
1951-2013年1月34海温指数与7月我国160站降水相关系数图
分析：
图1是1951-2013年1月34海温指数与海平面气压相关系数图，从图中可以看出，赤道中东太平洋地区海温指数和海平面气压成明显正相关，最大相关系数已超过0.8；而赤道锡太平洋地区海温指数和海平面气压则成负相关，最小相关系数低于-0.6。由此可以解释当发生厄尔尼诺事件时，由于赤道中东太平洋和秘鲁沿岸等地区海温上升，使沃克环流减弱，对中东太平洋地区的下沉气流产生负的影响，使中东太平洋地区海平面气压降低；与之相反，西太平洋地区上升气流减弱，海平面气压升高。
''
‘ E:\\20150ቤተ መጻሕፍቲ ባይዱ06\2\161'
' E:\\20150706\2\2'
' 3 54'
' 70 135'
' (.2)'
' a1(0.5)'
' 9(a1)'
' '
' '
' 1 1 1'
' E:\\20150706\1\7'
' '
' 0'
'd '
' '
' 2'
'd '
' 20131301070'
(1*2)
书写.7的数据描述文件2；
E:\\20150706\2\.7
E:\\20150706\2\7
2.5E33
1 1951 1
1
r 0 99
准备.7的站点映射文件7；
准备画站点图时插值时所需格点数据和及其描述文件161（已提供）；
编写2文件，绘制1951-2013年1月34海温指数与7月中国160站降水相关系数图，保存于.7文件中。
1,63
(c(i),'(A75.2)') a(i)(i)
1,160
1,63
s(j)()
r(i)(,63)
*
0
1
1
(7'E:\\20150706\2\.7')
(8'E:\\20150706\1\.7''')
1,160
(7,*) (i)(i)(i)(i)
1,160
(8)(i)(i)(i)(i)
0
(8)(160-1)(160-1)(160-1)(160-1)
•34海温指数与中国160站降水相关图
1.绘制1月份34海温指数与1月份海平面气压相关图
1.1资料说明
34指数：34是1951年1月至2013年12月3.4区（热带太平洋：西经170度~西经120度、北纬5度~南纬5度）区域平均海温指数资料，共有63行数据，每行数据第一个数字为年份，后面12个数字为该年1-12月的海温指数；
160站站点资料：是160站站号和纬度、经度资料，共有160行。每行第一个数据为8位字符型的站号，后面两个数字分别为该站点纬度和经度；
2.2步骤：
用编写90文件，计算1951-2013年1月34海温指数与7月160站降水相关系数，计算结果保存于.7和.7文件中；
p(160,63)
s(63)(160)(63)(160)(160)
(8'D:\\\''')
(7,'(144F5.2)') ((r()1,144)1,73)
(8) ((r()1,144)1,73)
()
()
()(n)(n)
!返回的相关系数
!循环控制变量
12 !加法器
向量的平均值
!计算平均值
0.
1
(i)
0.
1
(i)
!计算相关系数
0.
1=0.
2=0.
1
(a(i))*(b(i))
11+(a(i))**2
22+(b(i))**2
(1*2)
保存于文件中。
书写的数据描述文件1：
e:\\20150706\1\
-9.99e33
144 -180 2.5
73 -90 2.5
1 1000
63 1951 1
1
0 0
编写1文件：
''
' e:\\20150706\1\1'
' 1000'
' -180 180'
()
!本程序计算两列向量的相关系数
分别是待计算的向量
是向量的长度，要求两列向量等长
()
()(n)(n)
!返回的相关系数
!循环控制变量
12 !加法器
向量的平均值
!计算平均值
0.
1
(i)
0.
1
(i)
!计算相关系数
0.
1=0.
2=0.
1
(a(i))*(b(i))
11+(a(i))**2
22+(b(i))**2
实习要求：
•根据所提供的实习说明和资料，按要求完成实习内容；
•将完成的实习内容按步骤写入实习报告中，包括资料说明、程序、文件、文件、绘制的图形（要求在图形上方用的方式添加标题，注明姓名拼音及学号）
•根据所绘图形及查找相关文献，说明本次实习所求物理量场的意义。
实习内容：
•34海温指数与海平面气压相关图
*7 a(63)
(5'D:\\\''')
(6'D:\\\34')
1,63
(6,'(A12)') c(i)
1,63
(c(i),'(A75.2)') a(i)(i)
1,63
(5) ((p()1,144)1,73)
*
1,73
1,144
1,63
z(k)()
r()(,63)
*
(7'D:\\\')
' -90 90'
' t 1'
' '
' '
' -0.2 0.2'
'd '
' '
'd '
' 20131301070'
' 1.0 0'
' e:\\20150706\1\1 '
;
绘制1951-2013年1月34海温指数与海平面气压相关系数图：
1951-2013年1月34海温指数与海平面气压相关系数图
2绘制1月份34海温指数与7月份中国160站降水相关图：
*12 c(63)
*7 a(63)
*8 (160)
(6'E:\\20150706\2\34')
(5'E:\\20150706\2\r1607')
(10:\\20150706\1\')
1,160
(10,*) (i)(i)(i)
*
1,63
(5,*) (p()1,160)
1,63
(6,'(A12)') c(i)
2.1资料说明
34指数：34是1951年1月至2013年12月3.4区（热带太平洋：西经170度~西经120度、北纬5度~南纬5度）区域平均海温指数资料，共有63行数据，每行数据第一个数字为年份，后面12个数字为该年1-12月的海温指数；
160站降水资料：r1607是中国160站1951~2013年共63年7月降水资料，数据是1951年7月各站降水量按照中各站顺序写入，该年160站都写入后再写下1952年7月各站降水量，...................，直至2013年7月各站降水量按照中各站顺序写入；
全球海平面气压月平均数据：为1951至2013年的1月全球海平面气压场二进制数据，年数为63年，格点数144*73，水平分辨率2.5︒*2.5︒；
1.2步骤：
1用编写90文件，计算1951-2013年1月34海温指数与海平面气压相关系数，计算结果保存于文件中；
90文件：
(63)(73,144,63)(73,144)(63)
''
' '
1951-2013年1月34海温指数与7月我国160站降水相关系数图
分析：
图1是1951-2013年1月34海温指数与海平面气压相关系数图，从图中可以看出，赤道中东太平洋地区海温指数和海平面气压成明显正相关，最大相关系数已超过0.8；而赤道锡太平洋地区海温指数和海平面气压则成负相关，最小相关系数低于-0.6。由此可以解释当发生厄尔尼诺事件时，由于赤道中东太平洋和秘鲁沿岸等地区海温上升，使沃克环流减弱，对中东太平洋地区的下沉气流产生负的影响，使中东太平洋地区海平面气压降低；与之相反，西太平洋地区上升气流减弱，海平面气压升高。
''
‘ E:\\20150ቤተ መጻሕፍቲ ባይዱ06\2\161'
' E:\\20150706\2\2'
' 3 54'
' 70 135'
' (.2)'
' a1(0.5)'
' 9(a1)'
' '
' '
' 1 1 1'
' E:\\20150706\1\7'
' '
' 0'
'd '
' '
' 2'
'd '
' 20131301070'
(1*2)
书写.7的数据描述文件2；
E:\\20150706\2\.7
E:\\20150706\2\7
2.5E33
1 1951 1
1
r 0 99
准备.7的站点映射文件7；
准备画站点图时插值时所需格点数据和及其描述文件161（已提供）；
编写2文件，绘制1951-2013年1月34海温指数与7月中国160站降水相关系数图，保存于.7文件中。
1,63
(c(i),'(A75.2)') a(i)(i)
1,160
1,63
s(j)()
r(i)(,63)
*
0
1
1
(7'E:\\20150706\2\.7')
(8'E:\\20150706\1\.7''')
1,160
(7,*) (i)(i)(i)(i)
1,160
(8)(i)(i)(i)(i)
0
(8)(160-1)(160-1)(160-1)(160-1)
•34海温指数与中国160站降水相关图
1.绘制1月份34海温指数与1月份海平面气压相关图
1.1资料说明
34指数：34是1951年1月至2013年12月3.4区（热带太平洋：西经170度~西经120度、北纬5度~南纬5度）区域平均海温指数资料，共有63行数据，每行数据第一个数字为年份，后面12个数字为该年1-12月的海温指数；
160站站点资料：是160站站号和纬度、经度资料，共有160行。每行第一个数据为8位字符型的站号，后面两个数字分别为该站点纬度和经度；
2.2步骤：
用编写90文件，计算1951-2013年1月34海温指数与7月160站降水相关系数，计算结果保存于.7和.7文件中；
p(160,63)
s(63)(160)(63)(160)(160)
(8'D:\\\''')
(7,'(144F5.2)') ((r()1,144)1,73)
(8) ((r()1,144)1,73)
()
()
()(n)(n)
!返回的相关系数
!循环控制变量
12 !加法器
向量的平均值
!计算平均值
0.
1
(i)
0.
1
(i)
!计算相关系数
0.
1=0.
2=0.
1
(a(i))*(b(i))
11+(a(i))**2
22+(b(i))**2
(1*2)
保存于文件中。
书写的数据描述文件1：
e:\\20150706\1\
-9.99e33
144 -180 2.5
73 -90 2.5
1 1000
63 1951 1
1
0 0
编写1文件：
''
' e:\\20150706\1\1'
' 1000'
' -180 180'
()
!本程序计算两列向量的相关系数
分别是待计算的向量
是向量的长度，要求两列向量等长
()
()(n)(n)
!返回的相关系数
!循环控制变量
12 !加法器
向量的平均值
!计算平均值
0.
1
(i)
0.
1
(i)
!计算相关系数
0.
1=0.
2=0.
1
(a(i))*(b(i))
11+(a(i))**2
22+(b(i))**2
实习要求：
•根据所提供的实习说明和资料，按要求完成实习内容；
•将完成的实习内容按步骤写入实习报告中，包括资料说明、程序、文件、文件、绘制的图形（要求在图形上方用的方式添加标题，注明姓名拼音及学号）
•根据所绘图形及查找相关文献，说明本次实习所求物理量场的意义。
实习内容：
•34海温指数与海平面气压相关图