GrADS绘图实习

阅读: 评论:0

实习 数据文件的转换及数据描述文件的建立
一、实习资料
现在“data”文件夹下有ASCII码数据资料文件:月平均风场 (u200.dat、u850.dat、v200.dat、v850.dat)mh370写给2014的一封信、月平均高度场 (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和d文件)的Fortran程序,给出相应的数据描述文件(uv.ctl和l文件),并利用GrADS基本操作命令显示该地区850hPa和200hpa uv风场,最后保存图形。
三、实习目的:
通过编写该程序,进一步熟悉数据GrADS的数据格式,熟练编写将ASCII码数据转换为二进制格式的Fortran程序以及建立相应的CTL文件。
四、实习步骤:
4.1 利用“写字板”程序打开*.dat数据文件,熟悉该数据资料;
4.2 按要求编写Fortran程序,将所给的ASCII码数据资料文件转换成二进制无格式直接存取文件,结果保存为d;
4.3 通过“写字板”或“记事本”程序编写相应的数据描述文件,保存为l;
丙烯酸乳液4.4 画出2002年1月850hpa风场图,与图1比较验证数据正确性;
4.5 利用GrADS基本操作命令(open,display,print等)显示2002年7月850hpa和200hpa uv风场图形;
4.5 完成实习报告
1)说明所用资料
  资料名称:月平均风场篮子鱼 (u200.dat雷电颂朗诵u850.datv200.datv850.dat)
月平均高度场 (hgt500.dat)
范围:60~150E0~40N.
层次:风场,850 hPa200hPa;高度场,500hPa
时段:2002.1~2005.1248个月。
变量名:u  v  hgt
分辨率:2.5*2.5,格点数:37*17
资料读取文件:data.f90
2)给出所编写的Fortran程序: 
parameter(nx=37,ny=17,nz=1,nt=48)
新河东狮吼下载
dimension u1(nx,ny,nz,nt),v1(nx,ny,nz,nt),u2(nx,ny,nz,nt),v2(nx,ny,nz,nt),hgt(nx,ny,nz,nt)
open(1,file='u200.dat')
open(2,file='v200.dat')
open(3,file='u850.dat')
open(4,file='v850.dat')
open(11,file='hgt500.dat')
open(12,file='uv.grd',form='binary')
open(13,file='d',form='binary')
do 100 it=1,nt许仕龙
  do 50 iz=1,nz
      read(1,*) ((u1(i,j,iz,it),i=1,nx),j=1,ny)
      read(3,*) ((u2(i,j,iz,it),i=1,nx),j=1,ny) 
  50 continue
  do 51 iz=1,nz
      read(2,*) ((v1(i,j,iz,it),i=1,nx),j=1,ny)
      read(4,*) ((v2(i,j,iz,it),i=1,nx),j=1,ny)
  51 continue
  do 52 iz=1,nz
      read(11,*) ((hgt(i,j,iz,it),i=1,nx),j=1,ny)
  52 continue
100  continue
do 200 it=1,nt
  do 30 iz=1,nz
      write(12) ((u2(i,j,iz,it),i=1,nx),j=1,ny) 
      write(12) ((u1(i,j,iz,it),i=1,nx),j=1,ny)     
  30 Continue

本文发布于:2023-07-06 12:57:46,感谢您对本站的认可!

本文链接:https://patent.en369.cn/xueshu/178461.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:文件   数据   程序   资料   实习   相应   数据资料
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2022 Comsenz Inc.Powered by © 369专利查询检索平台 豫ICP备2021025688号-20 网站地图