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

1、前几次实习总结,1、同时打开几个ctl时，要求ctl中初始时间设置相同； 2、同时打开几个ctl文件时，不同文件的物理量使用时要在物理量之后加 .filenumber来区分，例如：u.2； 3、当写一个二进制数据的ctl时，要按照所写数据的维度来写，不能随便更改。 格点数：所写范围最大值与最小值之差除以格距再加1。,实习五各类参数设置综合练习,1实习目的： 掌握特殊图形处理方法 掌握阴影图形画图方法 学会处理特别等值线与特殊字符的方法,实习五 各类参数设置综合练习,2实习资料： “data”文件夹下二进制月平均风场文件：uv.grd 数据描述文件：uv.ctl 时间范围：2002.1-2005

2、.12 共48个月 水平范围：3717个网格点，分辨率为2.52.5 经纬度范围：60-150E， 0-40N 垂直层次：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. 实习步骤： 在写字板或记事本程序

3、中按要求编写*.gs文件，注意文件编写格式； 启动GrADS绘图软件，运行*.gs文件； 完成实习报告： 1）说明所用资料； 2）给出所编写的*.gs文件； 3）给出所绘图形。,实习五 各类参数设置综合练习,实习要求解读：,实习五 各类参数设置综合练习,(1)纬向风为西风时填色，东风绘制等值线，给出色标，0值线加粗； 提示： 西风：纬向风（u）大于0， 东风：纬向风（u）小于0。,色标： cbarn 1 0 4.5 0.35,实习五 各类参数设置综合练习,内部函数，画出色标,色标尺寸，可更改色标大小，例如可改为0.8,色标是横放还是竖放，0为横放，1为竖放,色标中心坐标（x，y）,色标 cba

4、rn 1 0 4.5 0.35,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,实习

5、六 练习描述语言的使用,3实习要求： 1）根据所提供的数据资料，利用fwrite命令将范围（纬度：-20-60N，经度：90E-200E，时间：1971年1月至2000年12月）海温场资料写入名为sst.grd的二进制数据文件中。 注意确定写出数据的格点数。,实习六 练习描述语言的使用,3实习要求： 2）根据所得数据资料，利用Fortran程序将 sst.grd文件转换为十进制数据文件sst.txt。 3）根据所得sst.grd数据文件，编写相应的数据描述文件sst.ctl，并绘制图形。请将所绘图形与使用原来sst.mnmean.nc资料在相同维数环境下所画图形进行比较，看是否一致。,实习六

6、练习描述语言的使用,4. 实习步骤： 4.1 在写字板或记事本程序中按要求1）编写*.gs文件，进行数据资料的提取； 4.2 根据实习要求2），利用Fortran程序进行二进制与十进制数据文件的转换，并查看sst.txt的数据信息； 4.3 在写字板或记事本程序中编写数据描述文件sst.ctl，并保存； 4.4 利用sst.grd和sst.ctl文件画图，完成实习内容。,实习六 练习描述语言的使用,格点数如何确定： x方向格点数目：-20-60N 60-（-20）=80 格距：2度*2度 格点数：80/2+1=41 同理可知：y方向格点数目 时间： 1971年1月至2000年12月： 30年，

7、一年12个月，则：30*12=360个月,实习六 练习描述语言的使用,fwrite: 将输出结果(数据)存入一个由“set fwrite 文件名”所指定的文件中 open c:pcgradsmodelmodel.le.ctl set gxout fwrite （设置输出类型） set fwrite d:resultu.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 （关闭文件u.grd）,实习六 练习描述语言的使用,ctl写