GrADS实习报告5

合集下载

实习报告

实习报告

2011年12月实习一:求500hPa高度场气候场、距平场和均方差场1.资料介绍:500hPa高度场,资料文件名h500.dat范围:60~150E,0~40N.时段:1982.1~1985.12共48个月分辨率:2.5*2.5,格点数:37*172.实习要求:编fortran程序,求500hPa高度场的(1)气候场;(2)距平场;(3)均方差场。

H500.For给出了如何用fortran读取资料h500.dat.并能用Grads做出图形。

3.实习原理:1.气候场:所谓的气候场即气象资料的平均场,资料中给出4年的500hpa 高度场的资料,需要求的气候场即四年资料的每个月一个平均的高度场,最终得出一共有12个平均场,这12个高度场即12个月的气候场。

求气候场的计算公式也就是求平均值的公式:。

2.距平场:距平值反映数据偏离平均值的状况,也是通常所说的异常,距平值是用来表征某个气象要素的波动幅度。

距平更明显地看出异常,消除了季节的变化。

距平场的计算公式为:,。

我们需要计算的距平场是原始数据与平均值之差,4年中一共有48个月,每月一个距平场,共由48个距平场。

3.均方差场:均方差用来描述样本中资料与平均值差异的平均状况,反映变量围绕平均值的平均变化程度(离散程度)。

均方差的计算公式:,。

4.实习步骤①利用Fortran程序进行数据处理和计算,Fortran编写的程序如下:program mainparameter(ii=37,jj=17,tt=48)real var(ii,jj,tt)integer i,j,texternal averageopen(5,file='d:\shixi\1\h500.dat')open(11,file='d:\shixi\1\ave.grd',form='binary')open(12,file=' d:\shixi\1\h5001.grd',form='binary')open(13,file=' d:\shixi\1\h5002.grd',form='binary')open(14,file=' d:\shixi\1\h5003.grd',form='binary')DO t=1,48read(5,1000)read(5,3000) ((var(i,j,t),i=1,ii),j=1,jj) !读取数据end docall average(var)1000 format(2i7)2000 FORMAT(37F6.2)3000 format(37f8.1)4000 format(37f7.2)close(5)close(11)close(12)close(13)end program!求气候场subroutine average(var)integer rreal var(37,17,48),av(37,17,12),jp(37,17,48),jfc(37,17,12)real bc(37,17,12)av(37,17,12)=0.0jfc(37,17,12)=0.0do r=1,12do t=r,48,12do i=1,37do j=1,17av(i,j,r)=av(i,j,r)+var(i,j,t)/4end doend doend dowrite(11) ((av(i,j,r),i=1,37),j=1,17)end do!求距平场r=0do t=1,48r=mod((r+1),12)if(r/=0) thendo i=1,37do j=1,17jp(i,j,t)=var(i,j,t)-av(i,j,r)end doend doelsedo i=1,37do j=1,17jp(i,j,t)=var(i,j,t)-av(i,j,12)end doend doend ifwrite(12) ((jp(i,j,t),i=1,37),j=1,17)end do!求均方差do r=1,12do t=r,48,12do i=1,37do j=1,17jfc(i,j,r)=jfc(i,j,r)+(var(i,j,t)-av(i,j,r))**2/4end doend doend dodo i=1,37do j=1,17jfc(i,j,r)=sqrt(jfc(i,j,r))end doend dowrite(13) ((jfc(i,j,r),i=1,37),j=1,17)end do!求变差系数do r=1,12do i=1,37do j=1,17bc(i,j,r)=jfc(i,j,r)/av(i,j,r)end doend dowrite(14) ((bc(i,j,r),i=1,37),j=1,17)end doend subroutine average②利用grads进行绘图,gs文件及绘出图像如下:'open d:\shixi\1\aver.ctl''set lat 0 40''set lon 60 150''set lev 500''enable print d:\shixi\1\aver.gmf'i=1while(i<=12)'set t 'i'''d aver''print''c'i=i+1endwhile'disable print';'open d:\shixi\1\h5001.ctl''set lat 0 40''set lon 60 150''set lev 500''enable print d:\shixi\1\h5001.gmf' i=1while(i<=12)'set t 'i'''d sx''print''c'i=i+1endwhile'disable print''open d:\shixi\1\h5002.ctl''set lat 0 40''set lon 60 150''set lev 500''enable print d:\shixi\1\vp.gmf' i=1while(i<=12)'set t 'i'''d vp''print''c'i=i+1endwhile'disable print''open d:\shixi\h5003.ctl''set lat 0 40''set lon 60 150''set lev 500''enable print d:\shixi\1\xd.gmf'i=1while(i<=48)'set t 'i'''d xd''print''c'i=i+1endwhile'disable print'实习二1.实习要求:根据表中给出的年平均气温和冬季平均气温的等级数据进行下列计算:1)计算两个气温之间的相关系数。

GrADS学习资料:实习五各类参数设置综合练习 实习六

GrADS学习资料:实习五各类参数设置综合练习 实习六

it=1
while(it<=5)
‘set lev 500 ’
‘set lon 120 270’ ‘set lat –10 10’
(设置维数环境)
‘set t ‘it’’
‘d u ’
it=it+1
endwhile
‘disable fwrite’
(关闭文件u.grd)
实习六 练习描述语言的使用
ctl写法见第二章,注意数据的维度设置要 按照数据实际大小来设置。
缺测值一定要与原数据一致。
fortran转换数据,与第一个实习刚好相反, 自己编程.
3.实习要求: 2)根据所得数据资料,利用For十进制数据文件sst.txt。 3)根据所得sst.grd数据文件,编写相应的数
据描述文件sst.ctl,并绘制图形。请将所绘 图形与使用原来sst.mnmean.nc资料在相同维 数环境下所画图形进行比较,看是否一致。
1. 实习目的: 掌握GrADS软件数据提取方法,学会资料处
理。
实习六 练习描述语言的使用
2.实习资料:
“data”文件夹下有全球海表温度资料 sst.mnmean.nc。
注意:海温时间从1854年1月开始。
nc文件打开:sdfopen 路径\ sst.mnmean.nc
nc文件的维数设置、格距、缺测值的查询: q ctlinfo
实习五 各类参数设置综合练习
3. 实习要求: 利用所提供的数据文件,编写.gs文件:
绘制出2002年1-12月120 E、 0-40N 200hPa 纬向风的纬度-时间剖面图。
(1)纬向风为西风时填色,东风绘制等值线, 给出色标,0值线加粗;
(2)X轴标注为“time”,Y轴标注为“lat”,标 题标注为“u 120E”。

实习报告

实习报告

中尺度气象模拟实习报告摘要:此次实习内容是通过数值模拟的方法对0812号台风“鹦鹉”发展及登陆过程进行模拟,利用grads软件画出各变量的二维图,最终结合所学知识(天气学原理、气象动力学等)分析模式对天气特征的模拟性能及影响这次天气过程前后的天气形势。

结论表明,(1)“鹦鹉”强度大,移速快,影响范围大(2)“鹦鹉”移动路径比较稳定,以西北偏西方向移动为主,但移动过程中仍出现明显北抬的现象(3)登陆时眼区大、风雨分布不对称,造成眼区内风雨小、外围风雨大,强降雨区域与台风中心区域分离引言:台风是以一种发生在热带海洋上的强烈风暴,在海洋面温度超过26℃以上的热带或副热带海洋上,由于近洋面气温高,大量空气膨胀上升,使近洋面气压降低,外围空气源源不断地补充流入上升去。

受地转偏向力的影响,流入的空气旋转起来。

而上升空气膨胀变冷,其中的水汽冷却凝结形成水滴时,要放出热量,又促使低层空气不断上升。

这样近洋面气压下降得更低,空气旋转得更加猛烈,最后形成了台风。

台风在全球三个海区影响最大,即西太平洋、西北太平洋和孟加拉湾。

西北太平洋是全球热带气旋发生频率最高的海域,而中国又是世界上登陆台风频率最高的国家,是西北太平洋和南海台风袭击次数最多的国家。

此次实习研究的0812号台风“鹦鹉”是属于西太平洋生成的台风,于8月22日16时登陆香港,之后穿过珠江口于22 时再次登陆中山,它是今年以来登陆我省最强的台风,登陆时风速达33 m/s,同时由于“鹦鹉”发展的不对称性形成偏心台风,使其登陆后出现直径达100 公里的风眼区,随着台风的西北移,横扫我省中北部大部分地区,形成台风中心经过的地方无风无雨,只有台风外围西南侧才出现大风大雨的现象。

总体来说台风“鹦鹉”的移动路径比较稳定,向西北偏西方向移动为主,但移动过程中仍出现明显北抬的现象,而且其登陆点及强度变化都与某些气象因子有很好的对应关系。

国内对台风“鹦鹉”的研究不是很多,其中易建坤指出影响台风移动路径变化的因素较多,环境场对台风的引导作用、台风自身结构与强度的变化及台风与周围大气环境的相互作用等都会影响台风的移动路径。

天气学诊断分析实习报告

天气学诊断分析实习报告

一、实习目得熟练使用天气学原理与中国天气中得知识,诊断分析常用得动力与热力参量(涡度、散度、相对湿度、水汽通量、水汽通量散度等)。

1、学会使用Fortran程序读入有格式得数据文件,输出二进制数据文件;2、学会使用Fortran编程计算涡度、散度、相对湿度、水汽通量等动力与热力参量。

3、学会使用Grads读入二进制数据文件,绘制等值线图、矢量图等,并存储。

4、掌握如何分析得到得各物理量得图形。

二、实习要求完成读取数据及利用给出得资料数据计算各物理量得fortran程序,将已有基本量数据及计算出得物理量数据写入out、grd文件里。

利用GRADS软件画出具有代表性得高度上得各物理量得图形,结合所画出得图形及韩国分析图与辅助图分析2009年4月18-20日南方气旋发展过程得各要素场配置情况,对这次过程进行天气诊断分析。

三、实习内容(含数据资料与方法说明)实习资料:计算范围:东经32~160,北纬12~80,格距 4*4。

所给要素为高度场,风场,温度以及温度露点差,共11层。

实习内容方法说明:1、利用所给资料(格点)与绘图软件绘制等高线、等压线、等温线,了解当时环流形势(主要就是500与850hPa)。

2、计算涡度(实测风涡度与地转风涡度)与散度(主要就是300,500,850hPa):计算18日08时 ,18日20时,19日08时,19日20点 ,20日08时500hPa实测风涡度平流,分析并比较实测风涡度与地转风涡度不同得原因,分析南方气旋发展生成得原因,以及散度高低层得配置。

绘图: 19日20点850hPa实测风涡度与地转风涡度以及她们得差, 19日08时,19日20点 ,20日08时500hPa实测风涡度平流。

3、计算1000-300hPa各层水汽通量并绘图,分析水汽通量输送较大得层次及其水汽通量散度,指出水汽输送通道大致范围,并说明水汽辐合区域得位置,什么位置有利于降水得发生,原因就是什么。

GrADS绘图实习

GrADS绘图实习

实习一数据文件的转换及数据描述文件的建立一、实习资料:现在“data”文件夹下有ASCII码数据资料文件:月平均风场 (u200.dat、u850.dat、v200.dat、v850.dat)、月平均高度场 (hgt500.dat),2002.1-2005.12共48个月,资料的水平网格范围是M×N个网格点(M=37,N=17),分辨率为2.5°×2.5°,自西向东经度60-150°E,由南至北纬度0-40°N。

二、实习要求:编写出将ASCII码数据资料文件:月平均风场和月平均高度场分别转换成二进制无格式直接存取(GrADS数据格式)文件(uv.grd和hgt.grd文件)的Fortran程序,给出相应的数据描述文件(uv.ctl和hgt.ctl文件),并利用GrADS基本操作命令显示该地区850hPa和200hpa uv风场,最后保存图形。

三、实习目的:通过编写该程序,进一步熟悉数据GrADS的数据格式,熟练编写将ASCII码数据转换为二进制格式的Fortran程序以及建立相应的CTL文件。

四、实习步骤:4.1 利用“写字板”程序打开*.dat数据文件,熟悉该数据资料;4.2 按要求编写Fortran程序,将所给的ASCII码数据资料文件转换成二进制无格式直接存取文件,结果保存为mhy.grd;4.3 通过“写字板”或“记事本”程序编写相应的数据描述文件,保存为mhy.ctl;4.4 画出2002年1月850hpa风场图,与图1比较验证数据正确性;4.5 利用GrADS基本操作命令(open,display,print等)显示2002年7月850hpa和200hpa uv 风场图形;4.5 完成实习报告1)说明所用资料资料名称:月平均风场(u200.dat、u850.dat、v200.dat、v850.dat);月平均高度场(hgt500.dat)范围:60~150E,0~40N.层次:风场,850 hPa、200hPa;高度场,500hPa;时段:2002.1~2005.12共48个月。

GrADS实习报告

GrADS实习报告

南京信息工程大学滨江学院GrADS实习报告院系:大气与遥感系专业:大气科学年级:14级班级:1学号:***********姓名:**实习一第一章气象程序设计与绘图软件安装与运行1.1实验实习目的1.1.1了解并掌握FORTRAN 90软件开发环境及基本操作1.1.2了解并掌握GrADS2.0的软件环境及基本操作方法(1)GrADS的安装(2)GrADS的启动与退出(3)GrADS文件处理(4)GrADS的使用流程1.2实验实习内容11.2.1问题的描述现有北京1951—1980年1月月平均气温文件beijingT.dat,编写程序计算北京30年的平均气温,并输出结果。

北京1951—1980年1月平均气温为1.0,—5.3,—2.0,—5.7,—2.1,—0.6,—1.7,—3.6,—3.0,0.1,—2.6,—1.4,—3.9,—4.7,—6.0,—1.7,—3.4,—3.1,—3.8,—2.0,—1.7,—3.6,—2.7,—2.4,—0.9,—2.7,—1.6,—3.91.2.2算法设计通过对问题的分析,设计求解算法,并绘制流图1.2.3程序编写program maininteger,parameter:: n=30realtave,sum,temp(n)open(1,file='c:users\bjxy\desktop\grads\chap1\data\beijingT.TXT',FORM='formatted' )do i=1,nread(1,*)temp(i)end doclose(1)open(2,file='c:users\bjxy\desktop\grads\chap1\data\tave.dat',FORM='formatted') sum=0.0do i=1,nsum=sum+temp(i)end dotave=sum/nwrite(2,*)taveclose(2)End1.2.4实验实习要求(1)创建新工作空间test1,工作空间文件夹创建在D盘上(2)创建新工作空间test1内创建新项目xml,项目文件夹创建在工作空间文件夹内(3)在项目xml内创建源程序文件test1.f90,文件创建在项目文件夹内,编辑输入给定的源程序文本,在前三行“”处输入班级、姓名、日期信息(以后要求相同)(4)在项目xml内创建辅助文档文件miaoshu1.txt,文件创建夹内,在文件中编辑输入问题描述文本(5)在项目xml内创建辅助文件suanfal.doc,文件创建在项目文件夹内,在文件中编辑绘制图1.1所示的流程图(6)编译源程序test1.f90构建可执行程序xml.exe,运行可执行程序xml.exe (7)强输入数据和输出结果以注释形式编辑输入到源程序末尾1.2.5 实验实习步骤(1)Compaq Visual FORTRAN 6.6系统安装后,在桌面创建一个Developer Studio 图标,在“开始”→“程序”菜单中创建Compaq Visual FORTRAN子菜单项,通过桌面图标或“开始”菜单可快速启动Microsoft Developer Studio。

grads实习报告

南京信息工程大学实验(实习)报告实验(实习)名称基于NCAR/NCEP再分析资料的降水和气温的气候特征分析日期2016.10 得分指导教师孙晓娟系大气科学专业大气科学班次长望实验班姓名夏新露一、实习目的(1)了解NCEP/NCAR再分析数据记得数据结构。

(2)通过NCEP/NCAR再分析数据集,掌握FORTRAN有格式输入输出语句(Read、Print、Write)的读写四维数据的方法。

(3)掌握FORTRAN多维数组的定义和使用方法。

(4)掌握隐式do循环的使用方法。

(5)巩固FORTRAN顺序结构概念和顺序结构程序设计方法。

(6)掌握GrADS对“*.nc”文件的读取方法。

(7)掌握GrADS中fwrite写出二进制数据的方法。

(8)掌握GrADS中While循环的使用方法。

(9)掌握GrADS中变量的定义和使用方法。

(10)掌握四维“*.grd”数据格式建立数据描述文件(*.ctl)的方法。

(11)巩固GrADS的基本操作命令的使用方法。

(12)掌握GrADS对contour和shaded图形类型的设置方法以及两种绘图类型的图形要素设置方法。

(13)掌握GrADS中基础绘图命令的设置方法。

(14) 掌握GrADS 中printim 生成文件(*.gmf 、*gif 、*.png )的方法。

二、 实习内容问题描述利用1948-2010年NCAR/NCEP 月平均气温和降水再分析资料(nc 格式文件),分析1月气温和降水气候特征。

要求利用Fortran 提取1948-2010年1月数据并进行计算,以“*.grd ”格式保存1948-2010年1月气温、降水气候场数据,用GrADS 生成1948-2010年1月气温、降水气候二维等值线和二维填色图两种格式叠加的图像,书写标题,并在高温、低温中心标记“H ”和“L ”。

保存并分析之。

问题分析已知:1948-2010年NCAR/NCEP 月降水和气温再分析资料。

气象统计方法实习报告

气象统计方法实习报告标准化管理处编码[BBX968T-XBB8968-NNJ668-MM9N]目录实习一求500hPa高度场气候场、距平场和均方差场--------------------- 31、资料介绍---------------------------------------------------- 32.要求-------------------------------------------------------- 33、实习结果---------------------------------------------------- 31)、FORTRAN源程序--------------------------------------------- 3(2)、grads文件--------------------------------------------- 5(3)、实习结果-------------------------------------------- 6实习二计算给定数据资料的简单相关系数和自相关系数----------------- 111、资料介绍--------------------------------------------------- 112、要求------------------------------------------------------- 113、实习结果--------------------------------------------------- 12(1)、Fortran源程序-------------------------------------- 12(2)、程序运行结果:------------------------------------- 14实习三分析中国夏季降水线性趋势的分布特征------------------------- 151.资料介绍及要求:-------------------------------------------- 152.实习结果--------------------------------------------------- 15(1).matlab程序------------------------------------------ 15(2).程序运行结果---------------------------------------- 16实习四求给定数据的一元线性回归方程------------------------------- 171、资料介绍及要求--------------------------------------------- 172、实习结果--------------------------------------------------- 18(1)、MATLAB程序----------------------------------------- 18(2)、程序运行结果--------------------------------------- 18(3)、结果分析------------------------------------------- 19实习五对给定的海温数据进行EOF分析------------------------------- 211、资料介绍--------------------------------------------------- 212、要求------------------------------------------------------- 213、实习结果:------------------------------------------------- 21(1)、FORTRAN源程序-------------------------------------- 21(2)空间场和时间序列的ctl文件--------------------------- 23(3)运行结果--------------------------------------------- 24(4)分析------------------------------------------------- 26实习三(附加)计算给定数据的11年滑动平均和累积距平-------------- 281、资料介绍--------------------------------------------------- 282、要求------------------------------------------------------- 283、实习结果--------------------------------------------------- 28实习四(附加)求给定数据的多元线性回归方程----------------------- 311、说明------------------------------------------------------- 312、要求------------------------------------------------------- 313、实习结果:------------------------------------------------- 31(1)Matlab源程序----------------------------------------- 31(2)运行结果--------------------------------------------- 35实习一求500hPa高度场气候场、距平场和均方差场1、资料介绍有一500hPa高度场资料,文件名,范围:60~150E,0~40N.时段:~共48个月。

转炉炼钢第五周实习报告

一、实习背景自实习开始以来,我已经在转炉炼钢车间度过了四周的时间。

这四周里,我深入了解了转炉炼钢的基本原理、操作流程以及安全生产规范。

随着实习的深入,我对炼钢工艺的认识逐渐加深,也对炼钢工人的辛勤付出有了更直观的感受。

二、实习内容本周,我的实习内容主要包括以下几个方面:1. 理论学习:继续学习转炉炼钢的基本理论知识,包括转炉炼钢的原理、设备结构、操作流程等。

通过查阅资料、请教师傅,我对转炉炼钢的整个过程有了更加全面的认识。

2. 现场操作:本周,我在师傅的指导下,参与了转炉炼钢的现场操作。

具体包括:- 出钢操作:学习如何控制转炉出钢的速度和温度,确保出钢质量。

- 吹氧操作:学习如何根据钢水成分调整吹氧强度,以控制钢水温度和成分。

- 炉前操作:学习如何观察炉内钢水的变化,判断钢水成分,为吹氧操作提供依据。

3. 设备维护:学习转炉炼钢设备的日常维护保养知识,包括炉体、吹氧系统、冷却系统等。

4. 安全生产:深入学习安全生产规范,提高安全意识,确保实习过程中的安全。

三、实习收获1. 理论知识:通过对转炉炼钢理论的学习,我对炼钢工艺有了更加全面的认识,为今后的工作打下了坚实的基础。

2. 操作技能:通过现场操作,我掌握了转炉炼钢的基本操作技能,为将来独立操作转炉打下了基础。

3. 安全意识:通过学习安全生产规范,我提高了安全意识,为实习过程中的安全提供了保障。

4. 团队协作:在实习过程中,我学会了与同事相互配合,共同完成工作任务。

四、实习体会1. 实践出真知:通过实习,我深刻体会到理论知识与实践操作相结合的重要性。

只有将所学知识运用到实际工作中,才能真正掌握一门技能。

2. 严谨认真:转炉炼钢工作需要严谨认真,一丝不苟。

在实习过程中,我学会了如何对待工作,对待每一个细节。

3. 安全第一:安全生产是转炉炼钢工作的重中之重。

在实习过程中,我时刻保持安全意识,确保实习过程中的安全。

4. 团队协作:转炉炼钢工作需要团队合作,只有团结协作,才能完成工作任务。

2024年毕业实习报告范文(七篇)

2024年毕业实习报告范文一、实习成果概述经过实习,各个作业流程的理论基础得到了实际应用与验证。

通过实际操作,不仅提高了知识运用能力,还深入了解了第三方物流服务的特点、主要设备和作业流程。

对此进行深入分析,对今后的学习具有积极的促进作用。

同时,结合所学理论,提出了针对性的改进意见。

二、实习内容本次实习的核心内容为深入了解第三方物流服务的特点、主要设备和作业流程,并在此基础上进行分析,结合理论知识提出改进意见。

三、实习工作流程初入顺丰物流公司,被分配至快运部。

首个星期,部分人员在客服部实习,负责接听客户电话,处理客户查询及投诉,并熟练掌握了制单与录单工作。

随后,两组人员岗位互换,进入库房实习,负责货物的标签粘贴、装卸及配送等工作。

库房工作中,特别注重货物摆放的合理性及装卸的规范性,确保货物安全无损。

四、个人工作心得与建议1. 货物码放应严格遵守规定,避免倒置或相邻货物相互影响,确保货物安全。

2. 卸货后应仔细核对清单,发现货物异常及时记录并报告理赔部。

3. 加强客服部整体水平的提升,领导需以身作则,定期培训员工,提升客服水平,并加大客服投入。

4. 加大网络宣传力度,提升公司在网络信息化时代的竞争力,加强员工网络知识培训。

5. 司机在装卸货物时应严格遵守规范,避免野蛮操作导致货损货差。

6. 公司应建立员工意见反馈机制,鼓励员工提出改进意见,促进公司和谐发展。

此次实习让我深刻认识到理论与实践相结合的重要性,为今后的工作和学习打下了坚实基础。

同时,也深感公司对于员工意见的重视对于公司发展的重要性。

2024年毕业实习报告范文(二)在会计的本职工作中,我始终保持敬业精神,同时利用空闲时间向出纳学习,以拓宽知识领域。

许多人对出纳的理解仅限于银行往来,但事实上,这仅是出纳职责的一部分,且并非易事。

出纳需要熟知各种银行业务,同时具备坚韧不拔的精神,无论严冬酷暑,都需频繁往返银行。

此外,出纳还需管理日常现金库,处理现金收支及开具支票和操作税控机。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
计算 1948-2010 年 1 月气温、降水气候值 FORTRAN 程序如下: program ex05 parameter(it=144,jt=73,lt=63)
问题分析 已知: 1948-2010 年 NCAR/NCEP 月降水和气温再分析资料。 绘制: 1948-2010 年 1 月气温、降水气候图。 通过分析,首先用 Grads 将从 NCAR/NCEP 月气温和降水再分析资料( nc
格式文件)中提取气温和降水 1 月逐年数据,并保存为为 grd 格式文件,再用
Fortran 程序利用 grd 格式资料按照公式 x
1
1 n xi 求得气温和降水的气候值,用 ni1
Ϊ¡°Êý¡±£©£®ÊýѧÆðÔ´ÓÚÈËÀàÔçÆÚµÄÉú²ú»î¶¯£¬¹Å°Í±ÈÂ×ÈË´ÓÔ¶¹Åʱ´ú¿ªÊ¼ÒѾ»ýÀÛÁËÒ»¶¨µÄÊýѧ֪ʶ£¬²¢ÄÜÓ¦ÓÃʵ¼ÊÎÊÌ⣮´ÓÊýѧ±¾Éí¿´£¬ËûÃǵÄÊýѧ֪ʶҲֻÊǹ۲ìºÍ¾ÑéËùµÃ£¬Ã»ÓÐ×ۺϽáÂÛºÍÖ¤Ã÷£¬µ«Ò²Òª³ä·Ö¿Ï¶¨ËûÃǶÔÊýѧËù×ö³öµÄ¹±Ï×£®»ù´¡ÊýѧµÄ֪ʶÓëÔËÓÃÊǸöÈËÓëÍÅÌåÉú»îÖв»¿É»òȱµÄÒ»²¿·Ö¡£Æä»ù±¾¸ÅÄîµÄ¾«Á¶ÔçÔڹŰ£¼°¡¢ÃÀË÷²»´ïÃ×ÑǼ°¹ÅÓ¡¶ÈÄڵĹŴúÊýѧÎı¾ÄÚ±ã¿É¹Û¼û£®´ÓÄÇʱ¿ªÊ¼£¬Æä·¢Õ¹±ã³ÖÐø²»¶ÏµØÓÐС·ù¶ÈµÄ½øÕ¹£®µ«µ±Ê±µÄ´úÊýѧºÍ¼¸ºÎѧ³¤¾ÃÒÔÀ´ÈÔ´¦ÓÚ¶ÀÁ¢µÄ״̬£®´úÊýѧ¿ÉÒÔ˵ÊÇ×îΪÈËÃǹ㷺½ÓÊܵġ°Êýѧ¡±£®¿ÉÒÔ˵ÿһ¸öÈË´ÓСʱºò¿ªÊ¼Ñ§ÊýÊýÆð£¬×îÏȽӴ¥µ½µÄÊýѧ¾ÍÊÇ´úÊýѧ£®¶øÊýѧ×÷Ϊһ¸öÑо¿¡°Êý¡±µÄѧ¿Æ£¬´úÊýѧҲÊÇÊýѧ×îÖØÒªµÄ×é³É²¿·ÖÖ®Ò»£®
出语句( Read、Print 、Write)的读写四维数据的方法。 ( 3) 掌握 FORTRAN 多维数组的定义和使用方法。 ( 4) 掌握隐式 do 循环的使用方法。 ( 5) 巩固 FORTRAN 顺序结构概念和顺序结构程序设计方法。 ( 6) 掌握 GrADS 对“ *.nc ”文件的读取方法。 ( 7) 掌握 GrADS 中 fwrite 写出二进制数据的方法。 ( 8) 掌握 GrADS 中 While 循环的使用方法。 ( 9) 掌握 GrADS 中变量的定义和使用方法。 ( 10)掌握四维“ *.grd”数据格式建立数据描述文件( *.ctl )的方法。 ( 11)巩固 GrADS 的基本操作命令的使用方法。 ( 12)掌握 GrADS 对 contour 和 shaded图形类型的设置方法以及两种绘
2
Ϊ¡°Êý¡±£©£®ÊýѧÆðÔ´ÓÚÈËÀàÔçÆÚµÄÉú²ú»î¶¯£¬¹Å°Í±ÈÂ×ÈË´ÓÔ¶¹Åʱ´ú¿ªÊ¼ÒѾ»ýÀÛÁËÒ»¶¨µÄÊýѧ֪ʶ£¬²¢ÄÜÓ¦ÓÃʵ¼ÊÎÊÌ⣮´ÓÊýѧ±¾Éí¿´£¬ËûÃǵÄÊýѧ֪ʶҲֻÊǹ۲ìºÍ¾ÑéËùµÃ£¬Ã»ÓÐ×ۺϽáÂÛºÍÖ¤Ã÷£¬µ«Ò²Òª³ä·Ö¿Ï¶¨ËûÃǶÔÊýѧËù×ö³öµÄ¹±Ï×£®»ù´¡ÊýѧµÄ֪ʶÓëÔËÓÃÊǸöÈËÓëÍÅÌåÉú»îÖв»¿É»òȱµÄÒ»²¿·Ö¡£Æä»ù±¾¸ÅÄîµÄ¾«Á¶ÔçÔڹŰ£¼°¡¢ÃÀË÷²»´ïÃ×ÑǼ°¹ÅÓ¡¶ÈÄڵĹŴúÊýѧÎı¾ÄÚ±ã¿É¹Û¼û£®´ÓÄÇʱ¿ªÊ¼£¬Æä·¢Õ¹±ã³ÖÐø²»¶ÏµØÓÐС·ù¶ÈµÄ½øÕ¹£®µ«µ±Ê±µÄ´úÊýѧºÍ¼¸ºÎѧ³¤¾ÃÒÔÀ´ÈÔ´¦ÓÚ¶ÀÁ¢µÄ״̬£®´úÊýѧ¿ÉÒÔ˵ÊÇ×îΪÈËÃǹ㷺½ÓÊܵġ°Êýѧ¡±£®¿ÉÒÔ˵ÿһ¸öÈË´ÓСʱºò¿ªÊ¼Ñ§ÊýÊýÆð£¬×îÏȽӴ¥µ½µÄÊýѧ¾ÍÊÇ´úÊýѧ£®¶øÊýѧ×÷Ϊһ¸öÑо¿¡°Êý¡±µÄѧ¿Æ£¬´úÊýѧҲÊÇÊýѧ×îÖØÒªµÄ×é³É²¿·ÖÖ®Ò»£®
图类型的图形要素设置方法。 ( 13)掌握 GrADS 中基础绘图命令的设置方法。 ( 14)掌握 GrADS 中 printim 生成文件( *.gmf 、*gif 、 *.png)的方法。
二、 实习内容
问题描述 利用 1948-2010 年 NCAR/NCEP 月平均气温和降水再分析资料( nc 格式文
程序设计: 提取 NCAR/NCEP 再分析资料中 1 月气温、降水二进制数据 ( 1) air1948-2010-1.gs 'reinit' 'sdfopen E:\daqi\GrADS\gradsSHIXI\shixi5\air.mon.mean.nc' 'setgxoutfwrite' 'setfwrite E:\daqi\GrADS\gradsSHIXI\shixi5\air1.grd' 'set x 1 144' 'set y 1 73' i=1 while(i<=756) 'set t 'i'' 'd air' i=i+12 endwhile 'disablefwrite' ;(2)pre1948-2010-1.gs 'reinit'
GrADS 绘制二者的气候图,分析 1948-2010 年 1 月蒙古高压气候特征。
三、 实习步骤
① 分析问题,理清算法和程序,设计和编写程序。 ② 书写“ air1948-2010-1.gs”和“ pre1948-2010-1.gs”从 NCAR/NCEP 再分析资
料( nc 格式),利用 fwrite 提取 1 月逐年数据,生成气温、降水 1 月的二进 制数据“air1.grd”、“pre1.grd”,编写 1 月气温、降水资料“air1.grd”、“pre11.grd” 的数据描述文件“ air1.ctl ”、“pre1.ctl”。 ③ 启动软件开发环境 Microsoft Developer Studio。 ④ 在 D 盘上创建新工作区 shixi05。 ⑤ 在工作区 shixi2 内创建新项目 shixi05。 ⑥ 在项目 shixi21 内创建源程序文件“ shixi05.f90 ”,编辑输入源程序文本,打 开 1948-2010 年 1 月气温、降水资料“ air1.grd”、“pre1.grd”,根据求均值公 式求得 1948-2010 年 1 月气温、降水气候场值。 ⑦ 编译、构建、运行、调试 Fortran 程序,生成 1948-2010 年 1 月气温、降水气 候场数据“ tc1.grd”、“ pc1.grd”。 ⑧ 打开记事本,编写“tc1.grd”和“ pc1.grd”的数据描述文件“tc1.ctl”和“pc1.ctl”。 ⑨ 利用 GrADS 的基本绘图命令和 contour、shaded要素设置方法以及基本绘图 命令,利用 printim 和 enable print 两种输出图像方法及输出图像格式设置, 编写“ tc1.gs”和“ pc1.gs”。 ⑩ 启动 GrADS,运行、调试“ tc1.gs”和“pc1.gs”,保存气温、降水气候图“tc1.gmf”、 “pc1. gmf”或者“ tc1.gif”、“ pc1. gif”。 ? 分析 1948-2010 年 1 月气温、降水的气候特征。
南京信息工程大学 实验(实习)报告
实验(实习)名称基于 NCAR/NCEP 再分析资料的降水和气温的气候特征分析 日期 2016.10 得分指导 实习目的
( 1) 了解 NCEP/NCAR 再分析数据记得数据结构。 ( 2) 通过 NCEP/NCAR 再分析数据集,掌握 FORTRAN 有格式输入输
'sdfopen E:\daqi\GrADS\gradsSHIXI\shixi5\pr_wtr.eatm.mon.mean.nc' 'setgxoutfwrite' 'setfwrite E:\daqi\GrADS\gradsSHIXI\shixi5\pre1.grd' 'set x 1 144' 'set y 1 73' i=1 while(i<=756) 'set t 'i'' 'd pr_wtr' i=i+12 endwhile 'disablefwrite' ; 编写“ air1.grd”、“ pre1.grd”数据描述文件“ air1.ctl”和“ pre1.ctl” ( 1) air1.ctl dset E:\daqi\GrADS\gradsSHIXI\shixi5\air1.grd title air temperature of NCEP Reanalysi s in Jan undef -9.96921e+36 xdef 144 linear 0 2.5 ydef 73 linear -90 2.5 zdef 1 linear 0 1 tdef 63 linear 00Z01JAN1948 1yr vars 1 air 0 t,y,x Winter Air Temperature endvars
Ϊ¡°Êý¡±£©£®ÊýѧÆðÔ´ÓÚÈËÀàÔçÆÚµÄÉú²ú»î¶¯£¬¹Å°Í±ÈÂ×ÈË´ÓÔ¶¹Åʱ´ú¿ªÊ¼ÒѾ»ýÀÛÁËÒ»¶¨µÄÊýѧ֪ʶ£¬²¢ÄÜÓ¦ÓÃʵ¼ÊÎÊÌ⣮´ÓÊýѧ±¾Éí¿´£¬ËûÃǵÄÊýѧ֪ʶҲֻÊǹ۲ìºÍ¾ÑéËùµÃ£¬Ã»ÓÐ×ۺϽáÂÛºÍÖ¤Ã÷£¬µ«Ò²Òª³ä·Ö¿Ï¶¨ËûÃǶÔÊýѧËù×ö³öµÄ¹±Ï×£®»ù´¡ÊýѧµÄ֪ʶÓëÔËÓÃÊǸöÈËÓëÍÅÌåÉú»îÖв»¿É»òȱµÄÒ»²¿·Ö¡£Æä»ù±¾¸ÅÄîµÄ¾«Á¶ÔçÔڹŰ£¼°¡¢ÃÀË÷²»´ïÃ×ÑǼ°¹ÅÓ¡¶ÈÄڵĹŴúÊýѧÎı¾ÄÚ±ã¿É¹Û¼û£®´ÓÄÇʱ¿ªÊ¼£¬Æä·¢Õ¹±ã³ÖÐø²»¶ÏµØÓÐС·ù¶ÈµÄ½øÕ¹£®µ«µ±Ê±µÄ´úÊýѧºÍ¼¸ºÎѧ³¤¾ÃÒÔÀ´ÈÔ´¦ÓÚ¶ÀÁ¢µÄ״̬£®´úÊýѧ¿ÉÒÔ˵ÊÇ×îΪÈËÃǹ㷺½ÓÊܵġ°Êýѧ¡±£®¿ÉÒÔ˵ÿһ¸öÈË´ÓСʱºò¿ªÊ¼Ñ§ÊýÊýÆð£¬×îÏȽӴ¥µ½µÄÊýѧ¾ÍÊÇ´úÊýѧ£®¶øÊýѧ×÷Ϊһ¸öÑо¿¡°Êý¡±µÄѧ¿Æ£¬´úÊýѧҲÊÇÊýѧ×îÖØÒªµÄ×é³É²¿·ÖÖ®Ò»£®
相关文档
最新文档