GrADS绘图与编程实习五实习六.pptVIP

前几次实习总结 1、同时打开几个ctl时，要求ctl中初始时间设置相同； 2、同时打开几个ctl文件时，不同文件的物理量使用时要在物理量之后加 .filenumber来区分，例如：u.2； 3、当写一个二进制数据的ctl时，要按照所写数据的维度来写，不能随便更改。 格点数：所写范围最大值与最小值之差除以格距再加1。 实习五各类参数设置综合练习 1．实习目的： 掌握特殊图形处理方法 掌握阴影图形画图方法 学会处理特别等值线与特殊字符的方法 2．实习资料： “data”文件夹下二进制月平均风场文件：uv.grd 数据描述文件：uv.ctl 时间范围：2002.1-2005.12 共48个月 水平范围：37×17个网格点，分辨率为2.5°×2.5° 经纬度范围：60-150°E， 0-40°N 垂直层次：850hPa、200hPa 物理量：u，v 3. 实习要求： 利用所提供的数据文件，编写.gs文件： 绘制出2002年1-12月120 E、 0-40N 200hPa 纬向风的纬度-时间剖面图。 （1）纬向风为西风时填色，东风绘制等值线，给出色标，0值线加粗； （2）X轴标注为“time”,Y轴标注为“lat”,标题标注为“u 120E”。 （3）最终将图形保存为gmf格式。 4. 实习步骤： 在写字板或记事本程序中按要求编写*.gs文件，注意文件编写格式； 启动GrADS绘图软件，运行*.gs文件； 完成实习报告： 1）说明所用资料； 2）给出所编写的*.gs文件； 3）给出所绘图形。 实习要求解读： X轴标注为“time”,Y轴标注为“lat” draw xlab time draw ylab lat 要求x轴为时间，y轴为纬度，需要用到垂直坐标和水平坐标相互交换，在书中第四章4.5中找命令。 1. 实习目的： 掌握GrADS软件数据提取方法，学会资料处理。 2．实习资料： “data”文件夹下有全球海表温度资料sst.mnmean.nc。 注意：海温时间从1854年1月开始。 nc文件打开：sdfopen 路径\ sst.mnmean.nc nc文件的维数设置、格距、缺测值的查询： q ctlinfo 3．实习要求： 1）根据所提供的数据资料，利用fwrite命令将范围（纬度：-20-60oN，经度：90oE-200Eo，时间：1971年1月至2000年12月）海温场资料写入名为sst.grd的二进制数据文件中。 注意确定写出数据的格点数。 3．实习要求： 2）根据所得数据资料，利用Fortran程序将 sst.grd文件转换为十进制数据文件sst.txt。 3）根据所得sst.grd数据文件，编写相应的数据描述文件sst.ctl，并绘制图形。请将所绘图形与使用原来sst.mnmean.nc资料在相同维数环境下所画图形进行比较，看是否一致。 4. 实习步骤： 4.1 在写字板或记事本程序中按要求1）编写*.gs文件，进行数据资料的提取； 4.2 根据实习要求2），利用Fortran程序进行二进制与十进制数据文件的转换，并查看sst.txt的数据信息； 4.3 在写字板或记事本程序中编写数据描述文件sst.ctl，并保存； 4.4 利用sst.grd和sst.ctl文件画图，完成实习内容。 实习六 练习描述语言的使用 格点数如何确定： x方向格点数目：-20-60oN 60-（-20）=80 格距：2度*2度 格点数：80/2+1=41 同理可知：y方向格点数目 时间： 1971年1月至2000年12月： 30年，一年12个月，则：30*12=360个月 实习六 练习描述语言的使用 fwrite: 将输出结果(数据)存入一个由“set fwrite 文件名”所指定的文件中 ‘open c:\pcgrads\model\model.le.ctl’ ‘set gxout fwrite’ （设置输出类型） ‘set fwrite d:\result\u.grd’ 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’