第四章蒙古高压特征分析

第四章:蒙古高压特征分析实验实习目的掌握气象要素气候及异常特征的基本研究方法。进一步巩固FORTRAN顺序结构和循环结构的程序方法和变量、函数、数组的使用方法，重点掌握外部文件对数据进行输入输出操作的方法。⑶了解并掌握GrADS的使用流程。学习并掌握数据描述文件(*ctl文件)的写方法。学习书写简单的"*.gs”文件。掌握GrADS中open、reinit>d、c、quit、q、run等基本命令的使用方法。掌握GrADS维数环境设置方法。⑻掌握GrADS绘制"line”、“bar”图形类型的方法，绘制1951-2010年蒙古高压强度异常的时间序列图。(9)掌握保存GrADS绘制图形的一种方法。实验实习内容问题描述已知1951-2010年一月蒙古高压强度、而积、位置指数序列，计算蒙古高压各指数的气候值、变率和距平，绘制蒙古高压的强度、而积、位置指数序列的时间序列图，分析冬季蒙古高压的异常规律。问题分析已知：1951-2010年一月蒙古高压强度、而积、经度、纬度指数序列资料、、、。计算：蒙古高压各指数的气候、变率和距平值。绘制：1951-2010年一月蒙古高压强度、而积、位置指数序列图。通过分析，根据公式()、公式()，求得蒙古高压一月环流指数气候及异常值。根据GrADS中Iine和bar两种图形绘制方法，绘制蒙古高压环流指数距平的时间序列图。实验实习步骤蒙古高压环流指数的气候和异常值计算分析问题，理清算法和程序，设计程序流程图并编写程序。启动软件开发环境MicrosoftDeveloperStudioo在I盘上创建新工作区shixi04o在工作区shixi04内创建新项目shixi04。在项目shixiO4内创建源程序文件“”，编辑输入源程序文本。在源程序文本中打开数据文件“”、"”、“”、“”，并将其值读入到相应的数组中。编写计算均值、变率和距平的子程序。调用子程序分别计算强度、面积、经度、纬度环流指数的均值、变率和距平。将蒙古高压环流指数的均值和变率写入到“”和"”两个文件中。将蒙古高压环流指数的距平值写入到"”和"”两个文件中。(10)编译、构建、运行、调试FORTRAN程序。蒙古高压环流指数距平时间序列图绘制GrADS是对数据进行分析、处理和显示的软件，所以在使用该软件进行绘图时,必须有相应的数据文件，并且这些数据必须满＜GrADS所要求的数扌居格式，否则就要进行数据处理；另外，GrADS并不是直接使用这些数据文件进行操作，而是通过一个对应的“数据描述文件”间接使用数据文件；对数据的处理和显示是通i±GrADS系统提供的命令完成的；操作完成可以在图形输出窗口观察所绘图形，并将图形保存和输出。具体使用流程如所述，此处不再赘述。按所述的使用流程，完成蒙古高压环流指数距平时间序列图绘制，需要按以下步骤进行：为蒙古高压坏流指数距平数据文件“”书写数据描述文件“”，在此文件中定艾四个变量pa,sa,Iona,lata.编写"”可执行文件，利用GrADS基本操作命令和line、bar两种绘图类型的绘图要素设置，以不同颜色和线性显示蒙古高压强度和面积距平时间序列图(曲线)，分别以不同颜色显示蒙古高压经度和纬度距平时间序列图(柱状)。将蒙古高压强度和面积距平时间序列图(曲线)保存到“”，将蒙古高压经度和纬度距平时间序列图(柱状)分别保存到"”和"”中。启动GrADS,调试、执行“”。分析蒙古高压气候及其异常特征。实验实习关键技术及方法•v——I某数据资料时间序列的距平X为数据资料无与其平均值=之差■X=召一1=1,2,…,“()某数据资料的变率b为其均方差，反应变量围绕平均值的平均变化程度，其计算公式为实验实习程序编写FORTRAN程序编写以下程序用于蒙古高压环流指数气候及异常值计算。PROGRAMmhIMPLICITnoneinteger,parameter::ny=60!p(ny)、pa(ny)、pav和pd分别为强度指数原序列、距平序列、均值和差,其他参数量类似定艾reaIp(ny),s(ny),Ion(ny),Iat(ny),pa(ny),sa(ny),Iona(ny),lata(ny),pav,sav,Ionav,latav,pd,sd,Iond,latddo\7\J/\7\71234✓(X/(X/(X/(Xeessoointegeri,j,kopen(1,fiIe='i:\shixi\')open(2,fiIe='i:\shixi\')open(3,fiIe='i:\shixi\')open(4,fiIe='i:\shixi\')doi=1,nyread(1,*)p(i)read(2,*)s(i)read(3,*)Ion(i)read(4,*)lat(i)caIIcha(ny,p,pa,pav,pd)caIIcha(ny,s,sa,sav,sd)caIIcha(ny,Ion,Iona,Ionav,Iond)caIIcha(ny,Iat.lata,Iatav,latd)open(5,file=write(5,'*)write(5,*')write(5,*')write(5,'')close(5)i:\shixi\‘)pav,pdsav,sdIonav,IondIatav,latdopen(6,fiIe二’i:\shixi\‘，forirF*binary1)write(6)pav,pdwrite(6)sav,sdwrite(6)Ionav,Iondwrite(6)Iatav,latdcIose(6)!下面两种书写数据的数据结构不同，注意区分。open(7,fiIe=1i:\shixi\')write(7,'')write(7,'')write(7,1*)write(7,'')close(7)(pa(i)ti=1,ny)(sa(i),i=1,ny)(Iona(i),i=1,ny)(lata(i)ti=1,ny)open(8,fiIe二'i:\shixi\‘，forirF*binary1)doi=1,nywrite(8)pa(i)write(8)sa(i)write(8)Iona(i)write(8)lata(i)enddoclose(8)end!!!!!!!!!!!!!!!!!!!!!!!求特征值subroutinecha(ny,x,xa,xav,xd)integer::nyintegerireaI::x(ny),xa(ny),xav,xd,sumsum=0doi=1,nysum=sum+x(i)enddoxav=sum/nyxd=0doi=1,nyxa(i)=x(i)-xavxd=xa(i)*xa(i)+xdenddoxd二sqrt(xd/ny)returnendGrADS程序编写蒙古高压环流指数距平数据文件"”书写数据描述文件"”.dseti:\shixi\undef+33titleMongoIiahighpressurecirculationindexdeparturexdef1Iinear11ydef1Iinear11zdef1Iinear11tdef60IinearJan19511yrvars4pa10sa10Iona10lata10endvars绘制蒙古高压异常时间序列图，书写可执行文件“”。I•一•丄・remitopeni:\shixi\*enableprinti:\shixi\'setlatr*set1onr1set1ev1sett160**setgXOLitline'*setCCOIor2**setcstyler'setcthick4'*setcmark2'*dpa**setcco1or4'1setcstyle3'*setcthick6'*setcmarkT*d100*sa‘'print''disabIeprintenableprinti:\shixi\setgxoutbar1setbarbaseO'setbargap0*setccoIor2*dIona*print*disabIeprinc'enableprinti:\shixi\setgxoutbar1setbarbase0*setbargap0*setccoIor4*dlata'print*disableprin实验实习结果232015IOSo-s-1O-IS-20-25图11951-2010年1月蒙古高压强度与而积距平吋间序列图图21951-2010年1月蒙古高压经度距平时间序列图图31951-2010年1月蒙古高压纬度距平时间序列图实验实习结果分析从图1可以看出，1月蒙古高压强度与而积呈现显著的正相关关系。存在年际及年代际的周期震荡，震荡周期大约为3~4年以及20年，在1956年、1977年以及2008年出现较高值，20世纪60年代中期开始下降，大约1975年达到最低值，之后快速上升。从图2、3可以看出，1月蒙古高压的位置存在年代际的变化，周期大约为30~40年。在20世纪50年代初由南转北，由东转西，20世纪70年代末80年代初，蒙古高压位置由北转南。