雷达基数据读取程序

雷达基数据读取程序-Fortran version

敏视达公期以来对国家巨额投入的雷达系统采取了极为保守 的数据策略，雷达基数据格式只能通过非正常的渠道获得，

阻碍了雷达在气象业务和科研上的应用。敏视达的软件终端产品PUP采取了极为无聊的注册码策略以保障其所谓的版权，真不知这种操作设置极为不便、性能低劣 的软件有何保护的必要。

WSR98D雷达数据完全雷同与WSR88D雷达数据格式，相关的数据格式在国外众多模式和雷达处理软件中都有描述，只不过大多数是C版本，鉴于

Fortran在数值计算中的优势，在此本人提供Fortran版本的数据接口，适用于读取体扫描的基数据文件，并可使用grads直接插值绘图，为大家

提供科研上的方便。其格式与以前国家气象中心网站上所描述的数据格式完全一致。

type radar_record character*14 unused1 integer*2 Message_Type character*2 channel character*10 unused2

integer*4 radical_collect_time !! 径向资料采集的GMT时间(毫秒)

integer*2 radical_collect_date !! 儒略日（Julian）表示，自1970年1月1日开始

integer*2 unambiguousRange !! 不模糊距离,单位:0.1Km integer*2 AzimuthAngle !!方位角（[数值/8.]*[180./4096.]=度） integer*2 DataNumber !! 当前仰角内径向数据序号 integer*2 DataStatus !! 径向数据状态

integer*2 ElevationAngle !!仰角

integer*2 ElevationNumber !!体扫内的仰角编号 integer*2 FirstGateRangeOfRef !!第一个强度库的距离(米) integer*2 FirstGateRangeOfDoppler !!第一个速度/谱宽库的距离(米) integer*2 ReflectivityGateSize !! 强度库长(米) integer*2 DopplerGateSize !!速度/谱宽库数 integer*2 ReflectivityGates !!强度库数 integer*2 DopplerGates !!速度/谱宽库数 integer*2 radicalnumber integer*4 coefofsys

integer*2 RefPointer !!从雷达数据头到强度数据开始的字节数 integer*2 VelPointer !!从雷达数据头到速度数据开始的字节数 integer*2 SWPointer !!从雷达数据头到谱宽数据开始的字节数 integer*2 VelResolution !!速度分辨率:2=0.5m/s;4=1.0m/s integer*2 VCP !!体扫VCP模式 (11,21,31,32)

character*14 unused3

integer*2 NyquistVelocity !!Nyquist速率(0.01m/s)不模糊速度 character*38 unused4

character*1 dbz(460) !! 回波 character*1 vel(920) !! 速度 character*1 sw(920) !! 谱宽 character*4 unused5 endtype

注： 单个record可理解为一束雷达射线所得到的数据记录。 dbz\\\\vel\\\\sw可通过ichar函数返回数值。

感谢您的阅读，祝您生活愉快。