1、GrADS教学内容400-500元 13 26%调研要解决的问题：图1-1大学生月生活费分布我们从小学、中学到大学，学的知识总是限制在一定范围内，缺乏在商业统计、会计，理财税收等方面的知识；也无法把自己的创意准确而清晰地表达出来，缺少个性化的信息传递。对目标市场和竞争对手情况缺乏了解，分析时采用的数据经不起推敲，没有说服力等。这些都反映出我们大学生创业知识的缺乏；（2） 文化优势功能性手工艺品。不同的玉石具有不同的功效，比如石榴石可以促进血液循环，改善风湿和关节炎；白水晶则可以增强记忆力；茶晶能够帮助镇定情绪，缓解失眠、头昏等症状。顾客可以根据自己的需要和喜好自行搭配，每一件都独一无二、与众不
2、同。五、创业机会和对策分析十字绣 编制类 银饰制品类 串珠首饰类此次调查以女生为主，男生只占很少比例，调查发现58的学生月生活费基本在400元左右，其具体分布如（图1-1）现在是个飞速发展的时代，与时俱进的大学生当然也不会闲着，在装扮上也不俱一格，那么对作为必备道具的饰品多样性的要求也就可想而知了。实习一 数据文件的转换及数据描述文件的建立1实习目的：(1)通过编写程序，熟悉数据GrADS的数据格式；(2)熟练编写将ASCII码数据转换为二进制格式的Fortran程序；(3)熟练建立相应的CTL文件。2. 实习要求：(1)编写Fortran程序，将ASCII码数据资料文件转换成二进制文件：月平
3、均风场（uv.grd） 月平均高度场（hgt.grd）(2)写出二进制文件相应的数据描述文件： uv.ctl，hgt.ctl(3)利用GrADS基本操作命令显示该地区850hPa和200hPa uv风场，最后保存图nm形。3 实习步骤：(1) 利用“写字板”程序打开*.dat数据文件，熟悉该数据资料；(2 )按要求编写Fortran程序，将所给的ASCII码数据资料文件转换成二进制无格式直接存取文件，结果保存为*.grd；(3) 通过“写字板”或“记事本”程序编写相应的数据描述文件，保存为*.ctl；(4) 画出2002年1月850hPa风场图，与图1比较验证数据正确性；4.实习报告1）所用资
4、料“data”文件夹下有ASCII码数据资料文件：月平均风场：u200.dat、u850.dat、 v200.dat、v850.dat月平均高度场：hgt500.dat经纬度范围：60-150E， 0-40N时间范围：2002.1-2005.12 共48个月水平范围：3717个网格点，分辨率为2.52.52）所编写的Fortran程序 parameter(nx=37,ny=17,nt=48) dimension u200(nx,ny,nt),v200(nx,ny,nt),h500(nx,ny,nt), u850(nx,ny,nt),v850(nx,ny,nt) open(1,file=e:dw
5、ju200.dat) open(2,file=e:dwjv200.dat) open(3,file=e:dwju850.dat) open(4,file=e:dwjv850.dat) open(5,file=e:dwjhgt500.dat) do it=1,nt read(1,*)(u200(i,j,it),i=1,nx),j=1,ny) read(2,*)(v200(i,j,it),i=1,nx),j=1,ny) read(3,*)(u850(i,j,it),i=1,nx),j=1,ny) read(4,*)(v850(i,j,it),i=1,nx),j=1,ny) read(5,*)(h5
6、00(i,j,it),i=1,nx),j=1,ny) enddoclose(1);close(2);close(3);close(4);close(5)open(12,file=e:dwjuv0.grd,form=binary) do it=1,nt write(12) (u850(i,j,it),i=1,nx),j=1,ny) write(12) (u200(i,j,it),i=1,nx),j=1,ny) write(12) (v850(i,j,it),i=1,nx),j=1,ny) write(12) (v200(i,j,it),i=1,nx),j=1,ny) enddoclose(12)
7、 open(12,file=e:dwjhgt0.grd,form=binary) do it=1,nt write(12) (h500(i,j,it),i=1,nx),j=1,ny) enddoclose(12)end3）所编写的数据描述文件reinitopen F:1uv.ctlset lon 60 150set lat 0 40set t 1d u;vprintim F:1shuang.jpg;4）所绘制的图形实习二 基本命令的上机使用1. 实习目的：熟悉GrADS的工作环境及使用流程，掌握基本操作命令，练习编写和使用gs文件。2实习要求：利用所提供的数据文件，使用基本命令编辑gs文件绘制
8、图形一维单线图和两维平面图。3. 实习步骤：（1）熟悉GrADS绘图软件操作界面命令窗口和图形显示窗口。利用所提供的数据文件model.le.dat及相应的数据描述文件model.le.ctl，在命令窗口输入基本绘图指令，在图形显示窗口观察相应图形（2）练习编写和使用gs文件。利用文本编辑器（如记事本，或者写字板）将所需要的基本命令编辑成gs文件，然后使用run命令运行该文件；（3）根据所学过的图形要素设置命令“set”， 编辑gs文件绘制北半球200hPa纬向风一维单线图（经度固定于120E）和500hPa高度场平面图，并用“enable print”命令将所画图形分别存放。4实习报告：1）
9、说明所用资料“data”文件夹下有GrADS数据格式的二进制资料文件model.le.dat以及对应的数据描述文件model.le.ctl（这是一组模式输出的全球5天数值预报结果，包括了多个要素、多层，按经纬度网格存放的数据，具体见数据描述文件model.le.ctl）2）所编写的数据描述文件A北半球200hPa纬向风一维单线图（经度固定于120E）reinitopen F:1model.le.ctlset lon 120set lat 0 90set lev 200d uprintim F:1shuang2.1.jpg;B.500hPa高度场平面图reinitopen F:1model.le
10、.ctlset lon 0 360set lat 0 90set lev 500d zprintim F:1shuang2.2.jpg;3）所绘制的图形A北半球200hPa纬向风一维单线图（经度固定于120E）B.500hPa高度场平面图实习三 set命令的使用练习 1实习目的：(1)熟练使用set命令的各类参数设置(2)掌握.gs文件的编写格式和运行(3)学会保存图形文件2. 实习要求：(1)利用所提供的数据文件，绘制出第7时刻500hPa高度场图，要求设置等值线间隔单位为40位势米，颜色，标记方式，并且给出相应标题（请包含姓名拼音与学号），最终将图形保存。所有命令编写于.gs文件中。(2)
11、利用所提供的数据文件，绘制出第7时刻850hPa风场图，要求以三种方式(风矢量、流线、风向杆)显示风场，并将三幅图上下排列于同一图形窗口，并最终将图形保存。所有命令编写于.gs文件中。3. 实习步骤：在写字板或记事本程序中按要求编写*.gs文件，注意文件编写格式；启动GrADS绘图软件，运行*.gs文件。4.5 完成实习报告1）所用资料 “data”文件夹下二进制月平均风场文件：uv.grd数据描述文件：uv.ctl时间范围：2002.1-2005.12 共48个月水平范围：3717个网格点，分辨率为2.52.5经纬度范围：60-150E， 0-40N垂直层次：850hPa、200hPa物理量
12、：u，v“data”文件夹下二进制月平均高度场文件：hgt.grd数据描述文件：hgt.ctl时间范围：2002.1-2005.12 共48个月水平范围：3717个网格点，分辨率为2.52.5经纬度范围：60-150E， 0-40N垂直层次：500hPa物理量：hgt2）给出所编写的数据描述文件A第7时刻500hPa高度场图reinitopen F:1hgt.ctlset t 7set cint 40set ccolor 5set clab ond hgtdraw title fengshuang20112348046printim F:1shuang3.png whiteB第7时刻850hP
13、a风场图reinitnoopen F:1uv.ctlset vpage 0 8.5 7 10.5set t 7set gxout vectordisplay u;vset vpage 0 8.5 3.5 7set t 7set gxout vectordisplay u;vset vpage 0 8.5 0 3.5set t 7set gxout vectordisplay u;v3）所绘制的图形A第7时刻500hPa高度场图B第7时刻850hPa风场图实习四 各类参数设置综合练习 1. 实习目的：（1）掌握阴影图形与等值线图形的叠加画图方法（2）练习基础绘图指令draw命令的使用2实习要求：
14、（1）绘制出2003年7月60-150E、0-40N区域内850hPa流线图，且地形高度场超过1500米以上用黑色阴影显示，并且给出相应标题（请包含姓名拼音与学号），最终将图形保存。（2）所有命令编写于.gs文件中。3. 实习步骤：（1）在写字板或记事本程序中按要求编写*.gs文件，注意文件编写格式（2）启动GrADS绘图软件，运行*.gs文件4.实习报告1）所用资料地形高度资料：dxgd.grddxgd.ctl范围：2.5-360 E，-90-90N（全球）层次：1层，地面；时间：2002年1月变量名：h分辨率：2.5*2.5，格点数：144*73风场资料：uv.grduv.ctl范围：60
15、-150E，0-40N层次：风场，850 hPa、200hPa；时段：2002.12005.12共48个月变量名：u、v分辨率：2.5*2.5，格点数：37*172）所编写的*.gs文件reinitopen F:1uv.ctlopen F:1dxgd.ctlset grid offset t 1set gxout shadedset clevs 1500set ccols 0 1d h.2(z=1)set lev 850set t 19set gxout streamd u;vdraw title fengshuang20112348046printim F:1shuang4.png whit
16、e3）所绘图形.实习五 各类参数设置综合练习 1实习目的：（1）掌握特殊图形处理方法（2）掌握阴影图形画图方法（3）学会处理特别等值线与特殊字符的方法2. 实习要求：利用所提供的数据文件，编写.gs文件： 绘制出2002年1-12月120 E、 0-40N 200hPa纬向风的纬度-时间剖面图。（1）纬向风为西风时填色，东风绘制等值线，给出色标，0值线加粗；（2）X轴标注为“time”,Y轴标注为“lat”,标题标注为“u 120E”。（3）最终将图形保存为gmf格式。3. 实习步骤：（1）在写字板或记事本程序中按要求编写*.gs文件，注意文件编写格式；（2）启动GrADS绘图软件，运行*.g
17、s文件；完成实习报告：1） 所用资料；“data”文件夹下二进制月平均风场文件：uv.grd数据描述文件：uv.ctl时间范围：2002.1-2005.12 共48个月水平范围：3717个网格点，分辨率为2.52.5经纬度范围：60-150E， 0-40N垂直层次：850hPa、200hPa物理量：u，v2）所编写的*.gs文件；reinitopen F:1uv.ctlset lon 120set lev 200set t 1 12set xyrev onset gxout shadedset cmin 0d uset clab onset gxout contourset cmax 0d u
18、cbarn 1 0 4.5 0.35set clevs 0set cthick 11d udraw title u 120Edraw xlab timedraw ylab latprintim F:1shuang5.gmf;3）所绘图形。实习六 练习描述语言的使用 1. 实习目的：掌握GrADS软件数据提取方法，学会资料处理。2实习要求：(1)根据所提供的数据资料，利用fwrite命令将范围（纬度：-20-60N，经度：90E-200E，时间：1971年1月至2000年12月）海温场资料写入名为sst.grd的二进制数据文件中。注意确定写出数据的格点数。(2)根据所得数据资料，利用Fortra
19、n程序将 sst.grd文件转换为十进制数据文件sst.txt。(3)根据所得sst.grd数据文件，编写相应的数据描述文件sst.ctl，并绘制图形。请将所绘图形与使用原来sst.mnmean.nc资料在相同维数环境下所画图形进行比较，看是否一致。3. 实习步骤：(1)在写字板或记事本程序中按要求1）编写*.gs文件，进行数据资料的提取；(2)根据实习要求(2)，利用Fortran程序进行二进制与十进制数据文件的转换，并查看sst.txt的数据信息；在写字板或记事本程序中编写数据描述文件sst.ctl，并保存；(3)利用sst.grd和sst.ctl文件画图，完成实习内容。4. 实习报告：1
20、）所用资料；“data”文件夹下有全球海表温度资料sst.mnmean.nc。注意：海温时间从1854年1月开始。nc文件打开：sdfopen 路径 sst.mnmean.ncnc文件的维数设置、格距、缺测值的查询：q ctlinfo2）所编写的件；A） *.ctl文件dset F:1sst.grdtitle NOAA Extended Reconstructed SST V3undef 32767xdef 56 linear 90 2ydef 41 linear -20 2zdef 1 linear 1 1tdef 360 linear 00Z01JAN1971 1movars 1sst 0
21、 -999 Monthly Means of Sea Surface TemperatureendvarsB） *.gs文件reinitsdfopen F:1sst.mnmean.ncset gxout fwriteset fwrite F:1sst.grdit=1405while(it=1764)set lon 90 200set lat -20 60set t itd sstit=it+1endwhiledisable fwrite;实习七 变量与函数的设置 1实习目的：(1)掌握函数的使用(2)进一步练习阴影图形与等值线图形的叠加画图方法。2. 实习要求：(1)所有命令编写于.gs文件中
22、。(2)根据所提供的数据资料，计算1-12月多年平均（1971年-2000年）海温场，并绘出1月海温多年平均图。 (注意原始nc数据自1854年1月开始)(3进一步计算1971年至2000年1月海温距平场，绘制1998年1月和1999年1月海温距平场，海温正距平区填色，画色标，并分别写标题Jan1998和Jan1999。要求设置虚页，将两幅图并排放在一起。最终将图形保存。3. 实习步骤： (1字板或记事本程序中按要求编写*.gs文件，注意文件编写格式；(2GrADS绘图软件，运行*.gs文件；4. 实习报告：1）所用资料； “data”文件夹下有全球海表温度资料sst.mnmean.nc。nc
23、文件打开：sdfopen 路径 sst.mnmean.ncnc文件的维数设置、格距、缺测值的查询：q ctlinfo2）所编写的*.gs文件；reinitreinitsdfopen F:1sst.mnmean.ncset t 1 12define avesst=ave(sst,t+1404+0,t+1764,12)modify avesst seasonalset t 1405 1764define anom=sst-avesstset t 1729set vpage 0 5.5 2 8d anomset gxout shadedset cmin 0d anomset clab onset g
24、xout contourset cmax 0d anomcbarn 0.8 0 4 0.5draw title JAN 1998set t 1741set vpage 5.5 11 2 8d anomset gxout shadedset cmin 0d anomset clab onset gxout contourset cmax 0d anomcbarn 0.8 0 4 0.5draw title JAN 1999printim F:1shuang7.png white3）所绘图形。实习八 描述语言的使用 1. 实习目的：掌握循环的使用，进一步练习函数的使用。2. 实习要求：（1）计算夏
25、季海温距平场；（2）循环绘出1971年至2000年夏季海温距平图。要求海温正距平区填色，画色标，并将年份作为标题。最终将图形保存。（3）所有命令编写于.gs文件中。3. 实习步骤：(1)在写字板或记事本程序中按要求编写*.gs文件；(2)启动GrADS绘图软件，运行*.gs文件4. 实习报告：1）所用资料 “data”文件夹下有全球海表温度资料sst.mnmean.nc。nc文件打开：sdfopen 路径 sst.mnmean.nc nc文件的维数设置、格距、缺测值的查询：q ctlinfo2）所编写的*.gs文件reinitsdfopen F:1sst.mnmean.ncset gxout
26、fwriteset fwrite e:dwj8.6.grdit=1410while(it=1758)set lev 0set lon 90 200set lat -20 60set t itd sstit=it+12endwhile disable fwrite;*.ctl文件dset F:18.6.grdtitle NOAA Extended Reconstructsd SST V3undef 32767xdef 56 linear 90 2ydef 41 linear -20 2zdef 1 linear 1 1 tdef 30 linear 00Z01Jun1971 1yrvars 1s
27、st 0 -999 Monthly Means of Sea Surface Temperature ENDVARs7月：*.gs文件：reinitsdfopen F:1 sst.mnmean.ncset gxout fwriteset fwrite e:dwj 8.7.grdit=1411while(it=1759)set lev 0set lon 90 200set lat -20 60set t itd sstit=it+12endwhile disable fwrite;*.ctl文件dset F:1 8.7.grdtitle NOAA Extended Reconstructsd S
28、ST V3undef 32767xdef 56 linear 90 2ydef 41 linear -20 2zdef 1 linear 1 1 tdef 30 linear 00Z01Jun1971 1yrvars 1sst 0 -999 Monthly Means of Sea Surface Temperature ENDVARs8月：*.gs文件：reinitsdfopen F:1 sst.mnmean.ncset gxout fwriteset fwrite F:18.8.grdit=1412while(it=1760)set lev 0set lon 90 200set lat -20 60set t itd sstit=it+12endwhile disable fwrite;*.ctl文件dset F:18.8.grdtitle NOAA Extended Reconstructsd SST V3undef 32767xdef 56 linear 90 2ydef 41 linear -20 2zdef 1 linear 1 1 tdef 30 linear 00Z01Jun1971 1yrvar
