GrADS绘图软件实用手3.doc

第三章 GrADS数据格式每一组GrADS数据应至少包括两组数据文件，数据描述文件ASCII码和数据文件二进制，数据的真正存放地。数据文件中只是用户数据的有序排放，而关于数据种类、排放次序等是单独放在一个文件中的称数据描述文件。而象GRIB和NETCDF等通用数据格式，以上两者是存于同一个文件的或称为自定义/自解释格式数据。但考虑到GrADS传统，对这类自定义格式数据仍将生成相应的数据描述文件。GrADS并不直接使用“数据文件”，而是通过“描述文件”间接使用“数据文件”。上一章中我们已使用过了这样的一组数据。以此为例，介绍用户如何按GrADS的格式，将自己的数据生成相应的数据文件和数据描述文件。1. 格点数据描述文件model.le.ctl文件清单：以*开始的行为注解行。1.1 数据描述文件各项解释1. DSET 数据文件名定义与此数据描述文件相对应的数据文件名。若两者位于同一目录，前面的路经可以省略或以“”开始，代表两者位于同一目录。若不在同一目录下，应给出路经参数。如：c:/pcgrads/sample/model.le.dat注意路经的给法与DOS不同，而与UNIX环境一致，便于移植！或c:pcgradssamplemodel.le.dat,两种都行。2. TITLE 数据文件说明文字串。3. UNDEF vaule定义缺测值。一般给一很大的正/负值，表示，当取值超过这一正值/低于定义的负值，认为该值无效。（GrADS采用跳过或用周围有效点的值处理。）4. OPTIONS 这里定义了与二进制存储有关的选项，二进制存储的一大特点是可移植性差，因此通过keywords项来增加可移植性。若keywords省略，则OPTIONS也可省略。可取：sequential: 顺序无格式方式。yrev: Y维与YDEF定义相反方式存放。zrev: Z维与ZDEF定义相反方式存放。big_endian: 如数据是在sun,sgi,hp cray机器上生成的，而目前不在此类机器上使用。little_endian: 如数据是在iX86,dec机器上生成的，而目前不在此类机器上使用。byteswapped：反序位存放。5. XDEF number LINEAR X_Start incrementnumber（=1）给定X方向格点数，其后LINEAR参数指明X方向是等间隔分布格点，X_Start起点坐标，increment网格间距。XDEF LEVELS value_listLEVELS参数指明X方向是不等间隔分布格点，因此，其后要给出具体每个格点的坐标值（以空格分开）。6. YDEF number LINEAR Y_Start incrementYDEF LEVELS value_list7. ZDEF number LINEAR Z_Start incrementZDEF LEVELS value_list8. TDEF number LINEAR T_Start incrementTDEF LEVELS value_listT_Start和value_list给出时间格式如下：hh:mmZddmmmyyyyhh要以两位数代表小时（缺省00）；“：mm”以两位数代表分钟（缺省：00）；dd以两位或一位数代表日期（缺省1）；mmm为英文月分三个字符的缩写；yyyy以两位或四位数代表年。如：12Z1JAN90（省：mm） 14：45Z22JAN1987 JUN1960时间increment格式： vvkkvv以两位数字代表；kk取：mn-增量以分为单位；hr-以小时为单位；dy-以天为单位；mo以月为单位；yr以年为单位。9. VARS numberV_abrev levs units descriptionENDVARS以VARS指示开始定义变量，number代表变量总的个数；以ENDVARS 指示结束变量定义。从VARS以下，每行定义一个变量。V_abrev变量名称，GrADS将用到；levs，数字，代表变量层数，0表示只有一层；units单位，为GRIB预留，给99；description，对变量的文字描述。1.2生成model.le.dat和model.le.ctl文件的程序代码片段program mainReal ps(72,46,7,5),ts(72,46,7,5),z(72,46,7,5),u(72,46,7,5),v(72,46,7,5)Real t(72,46,7,5),q(72,46,7,5),p(72,46,7,5)C 第一部分， 读入model.le.dat数据。Open (7,file=model.le.dat, form=unformatted,access=direct,recl=72*46*4)C 在微机上，要乘4；用VisualFortran6.5 不要乘4；sgi 上不要乘4，具体情况可能要试。 IREC = 1 设记录累加器，初值为1。DO it=1,5 最外层时间循环15天。read(7,rec=IREC) (ps(j,i,it),j=1,72),i=1,46) IREC = IREC+1DO k=1,7 多层变量，每一个水平层作为一个记录（72X46个点）。 read(7,rec=IREC) (u(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDODO k=1,7 read(7,rec=IREC) (v(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDODO k=1,7 read(7,rec=IREC) (z(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDODO k=1,7 read(7,rec=IREC) (t(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDODO k=1,5 水气只有5层 read(,rec=IREC) (q(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDO read(,rec=IREC) (ts(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1read(7,rec=IREC) (p(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDOC 第二部分， 写出mmodel.le.dat数据。Open (8,file=mmodel.le.dat,form=unformatted,access=direct,recl=72*46*4)Open (9,file=mmodel.le.ctl)IREC = 1 设记录累加器，初值为1。DO it=1,5 最外层时间循环15天。Write(8,rec=IREC) (ps(j,I,it),j=1,72),i=1,46) 按在CTL文件中的先后次序输出。 IREC = IREC+1DO k=1,7 多层变量，每一个水平层作为一个记录（72X46个点）。 Write(8,rec=IREC) (u(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDODO k=1,7 Write(8,rec=IREC) (v(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDODO k=1,7 Write(8,rec=IREC) (z(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDODO k=1,7 Write(8,rec=IREC) (t(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDODO k=1,5 水气只有5层 Write(8,rec=IREC) (q(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDO Write(8,rec=IREC) (ts(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1Write(8,rec=IREC) (p(j,i,k,it),j=1,72),i=1,46) IREC = IREC+1ENDDOC write out the correspond the “ctl file”Write(9,(a) DSET model.le.datWrite(9, (a) OPTIONS little_endian cray_32bit_ieeeWrite(9, (a) TITLE 5 Days of Sample Model OutputWrite(9, (a) UNDEF , -2.56e33Write(9, (a) XDEF , 72, Linear , 0.0, 5.0Write(9, (a) YDEF , 46, Linear , -90.0, 4.0Write(9, (a) ZDEF , 7 , levels ,1000,850,700,500,300,200,100Write(9, (a) TDEF , 5 , linear 02JAN1987 1DYWrite(9, (a) VARS, 8Write(9, (a) ps, 0, 99, Surface PressureWrite(9, (a) u , 7, 99, U windsWrite(9, (a) v , 7, 99, V windsWrite(9, (a) z , 7, 99, Geopotential HeightsWrite(9, (a) t , 7, 99, TemperatureWrite(9, (a) q , 5, 99, Specific HumidityWrite(9, (a) ts, 0, 99, Surface TemperatureWrite(9, (a) p , 0, 99, PrecipitationWrite(9, (a) ENDVARS一定要在第一列。STOPEND2. 站点数据的格式例1 1980年1月和2月4个站的降水量rain.ch文件清单：1980 1 qqq 34.3 -85.5 123.31980 1 rrr 44.2 -84.5 99.11980 1 sss 22.4 -83.5 412.81980 1 ttt 30.4 -82.5 23.31980 2 qqq 34.3 -85.5 145.11980 2 rrr 44.2 -84.5 871.41980 2 sss 22.4 -83.5 223.11980 2 ttt 30.4 -82.5 45.5年 月 站名 纬度 经度 降水量rain.c把上述文件转为GrADS数据格式（rain.dat）。rain.ctlrain.map数据的生成：dosstnmap i rain.ctl 结果生成rain.map数据。例2 多个站点，多层数据的生成rsond.ch1991 7 5 0 时间53.75 91.40 台站的纬度、经度（台站号丢失） 1005.92 256.0 288.64 286.04 270.00 2.00地面的气压、海拔、 T、 Td、 风向、 风速。850 1433.00 290.84 278.84 230.00 7.00高空层（hPa）、位势、T、 Td、 风向、 风速。700 3064.00 279.84 268.84 335.00 9.00500 5730.00 262.14 243.14 315.00 4.00400 7400.00 249.34 238.34 345.00 3.00300 9440.00 232.34 224.34 330.00 22.00250 10660.00 223.54 216.54 340.00 2.00200 12090.00 216.34 210.34 40.00 2.00150 13920.00 219.94 99999.99 305.00 5.00100 16550.00 220.74 99999.99 310.00 4.0035.75 111.40 54511 台站的纬度、经度、台站号1005.92 256.0 288.64 286.04 270.00 2.00850 1433.00 290.84 278.84 230.00 7.00700 3064.00 279.84 268.84 335.00 9.00 rsond.c生成rsond.datrsond.ctldset rsond.datoptions big_endian dtype stationstnmap rsond.mapundef 99999.0title raidio sondtdef 1 linear jan1980 1movars 11ps 0 99 prouserh0 0 99 heightts 0 99 temperaturetds 0 99 dewpointsus 0 99 U windsvs 0 99 V windshg 1 99 height 高层资料层数为1。t 1 99 tempstd 1 99 dewpointsu 1 99 U windsv 1 99 V windsendvarsDosstnmap i rsond.ctl 结果生成rsond.map数据。3. 台站资料的显示rain.gsopen rain.ctlset digsize 0.2set lon -140 -40set lat 15 85set t 1*set stid ond p;rsond.gsopen rsond.ctlset dignum 1set digsize 0.19set stid onset lon 50 140set lat 0 80set lev 300*set gxout valueset gxout modeld u;v;t-273.16;td-273.16;0.0 ; 0.0 ; 21*d us;vs;ts-273.16;tds-273.16; ps; 0.0,21set gxout findstnd t ; 6.4 ; 3.2 ; say resultd td ;5.6 ;2.3 ; say resultd t ;3.9 ; 4.0 ; say result4. Cressman.gs客观分析方法rain.ch：是一组观测，用rain.c生成GrADS格式的数据：rain.dat; rain.ctl，最后用cressman.gs 生成分析结果。open rain.ctlopen cmodel.ctl 与model.le.ctl一样，只是对数据的网格作了重新定义。 set digsize 0.13set dignum 3set lon 80 106set lat 24 46set t 1*set gxout valueset mproj nps*d p*set gxout contourd oacres(ts.2 , p.1*100 ) 将把站点资料p分析到利用cmodel.ctl定义的网格上。fprint *以下将分析结果保存到文件中。set gxout fwriteset fwrite dummy.datd oacres(ts.2 , p.1*100 )set dfile 2 设置当前文件为cmodel.ctlq dim 显示当前文件维数设置，按此结果写dummy.ctl文件。disable fwrite;Rain.ctldset rain.datoptions big_endiandtype stationstnmap rain.mapundef -999.0title Rainfalltdef 1 linear jan1980 1movars 1p 0 99 rainfallendvarscmodel.ctldset model.le.dat options little_endian cray_32bit_ieeeUNDEF -2.56E33 TITLE 5 Days of Sample Model OutputXDEF 180 LINEAR 0.0 1. 用户定义网格。可调整网格的范围和分辨率。 YDEF 90 LINEAR 0.0 1. ZDEF 1 LEVELS 1000 TDEF 1 linear jan1980 1mo 此时间要与rain.ctl中的一致vars 1 ts 0 99 Surface TemperatureENDVARS q dim 命令显示的结果Default file is: 2 X is varying lon = 80 to 106 x = 41 to 54Y is varying lat = 24 to 46 y = 13 to 24Z is fixed dummy.ctldset dummy.dat options little_endian cray_32bit_ieeeUNDEF -2.56E33 TITLE Sample OutputXDEF 14 LINEAR 80.0 2. (14=54-41+1) 2 =(106-80)/(14-1)YDEF 12 LINEAR 24.0 2. (12=24-13+1) 2 =(46-24)/(12-1)ZDEF 1 LEVELS 1000 TDEF 1 LINEAR 02JAN1987 1DY 时间与垂直层要与rain.ctl中的一致vars 1 p 0 99 sampleENDVARS 5. GRIB格式数据处理GRIB格式数据是WMO规定的一种通用的气象数据格式，它不像二进制数据，是与机器无关的，因此可以在各类机器上交换而不受限制，并且具有高压缩比。GRADS中数据文件与之相应的数据格式说明文件是分开的，而GRIB数据文件中包含了数据格式说明，是一种自定义或自解释压缩数据文件。但在GRADS中使用GRIB数据时，关键还是要生成一个与数据文件分离的CTL文件。GRADS演示数据中就有一套GRIB数据。因此要用到以下一些外部功能：grib2ctl.pl grib数据文件名 控制文件名.ctl该命令并不在GrADS的支持系统里，但可以从网上下载，下载后放在/grads/bin目录下，该命令是用perl语言写的，可以在RedHat Linux下运行。控制文件生成后还要运行：Dosgribmap i 控制文件名.ctl生成index文件就可以开始画图了。另外GrADS还提供了两个外部命令gribscan（旧版本）和wgrib（新版，grib2ctl.pl也是依赖wgrib）可以扫描GRIB数据或提取数据。例1：NCEP grib数据：grib20030731181) grib2ctl.pl grib2003073118 grib.ctl2) gribmap i grib.ctlgrib.ctldset grib2003073118index grib2003073118.idxundef 9.999E+20title grib2003073118 *produced by grib2ctl vp33edtype grib 3options yrevydef 181 linear -90.000000 1xdef 360 linear 0.000000 1.000000tdef 1 linear 18Z31jul2003 1mozdef 26 levels1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10 vars 93LFTXsfc 0 131,1,0 * surface Surface lifted index KO3MRprs 6 154,100,0 * Ozone mixing ratio kg/kg 6层臭氧混合率的信息是错误的50：500hPa以上变量名有点复杂，用户可自行简化。位势2：12UTC例2: T106/T213 模式数据：CPGHYI50.08208：8日（2001年9月8日）I：预报长度48小时1. grib2ctl.pl CPGHYI50.082 cpp.ctl2. gribmap -0 E i cpp.ctlcpp.ctldset CPGHYI50.082index CPGHYI50.082.idxundef 9.999E+20纬度间隔为1度title CPGHYI50.082* produced by grib2ctl vp33edtype grib 255options yrev y方向是倒置的，即从北极向南极放置经度间隔应为1度xdef 180 linear 0.000000 2ydef 91 levels -88.494 -86.543 -84.581 -82.616 -80.651 -78.684 -76.718 -74.751 -72.784 -70.817 -68.850 -66.883 -64.916 -62.949 -60.982 -59.015 -57.048 -55.081 -53.114 -51.146 -49.179 -47.212 -45.245 -43.278 -41.311 -39.344 -37.376 -35.409 -33.442 -31.475 -29.508 -27.541 -25.573 -23.606 -21.639 -19.672 -17.705 -15.737 -13.770 -11.803 -9.836 -7.869 -5.902 -3.934 -1.967 0.000 1.967 3.934 5.902 7.869 9.836 11.803 13.770 15.737 17.705 19.672 21.639 23.606 25.573 27.541 29.508 31.475 33.442 35.409 37.376 39.344 41.311 43.278 45.245 47.212 49.179 51.146 53.114 55.081 57.048 59.015 60.982 62.949 64.916 66.883 68.850 70.817 72.784 74.751 76.718 78.684 80.651 82.616 84.581 86.543时间错误。预报从2001年9月8日12UTC开始，48小时预报，正确时间是2001年9月10日12UTC。但先不要改，在执行完gribmap后再改。 88.494tdef 1 linear 12Z10sep1901 1mozdef 1 linear 1 1vars 1z 0 7,100,500 * Geopotential height gpmENDVARS用wgrib和gribscan分析数据，建立ctl文件wgrib v -GDS10 CPGHYI50.082 显示以下结果：gribscan d gv gd i CPGHYI50.082显示以下结果：Grib数据定义48小时预报变量名和单位wgrib1:0:D=1901090812:HGT:500 mb:kpds=7,100,500:48hr fcst:Geopotential height gpm:GDS10= 0 0 32 0 0 4 0 180 0 91 1 95 144 0 0 0 128 0 0 0 2 187 56 3 232 3 232 0 0 0 0 048小时预报预报开始时刻记录序号gribscanPDS,1,255,7,100,500,1,244,1901,9,8,12,0,1,48,0,0,1,hgt,Geopotential height,gpm,GDS,4,180,91,0.000 ,179.000,1000.000,90.000,0.000 ,1000.000,0,0,0网格定义：180x91纬度从90到0，间隔1度经度从0到179，间隔1度根据以上数据建立cpn.ctldset CPGHYI50.082title t106 model 500hpa heightundef -9999.00先设为预报开始时刻，运行完gribmap后再改为：12Z10sep2001options yrev y方向是倒置的，即从北极向南极放置dtype gribindex CPGHYI50.082.gmpxdef 180 linear 0.0 1.0ydef 91 linear 0.0 1.0zdef 1 linear 500 0tdef 1 linear 12Z08sep1901 1movars 1 z 0 7,100,500 Geopotential heightendvarsgribmap -0 -e -i cpn.ctl注：在unix下gribmap不能处理dos格式的ctl文件，要用dos2unix处理一下。如果不用0选项，tdef要定义成“预报开始时刻预报长度”，如上例：12Z10sep1901，即使这样，时间也是错误的！并 grib2ctl.pl需要perl环境支持，本手册实例都是在linux环境完成的。6. T213 Grib 数据转换目前中央台T213模式结果是按grib格式向下发送的，以下一组工具可以将T213，按要素、按等压面存放的单个grib数据文件组合成GrADS要求的4维数据（空间3维加时间维）并生成相应的ctl文件。工具组是在unix或linux运行环境支持下的shell编程工具。具体流程如下：Autftp工具： 用户自定义的一组文件名，通过ftp工具下载你所需的T213预报结果。 grib解码工具： 执行由c语言编写的解码程序把T213grib码数据转换成ASCII码文件，此时要素还是按要素、按等压面存放在一个个单独的文件中。 GrADS数据生成工具： 执行由Fortran语言编写的程序，把上述ASCII码文件，转换成GrADS二进制数据，并生成相应的ctl文件。用wgrib解码grib数据在dos命令行窗口运行wgrib命令，结果显示wgrib的用法doswgribPortable Grib decoder for NCEP/NCAR Reanalysis etc. it slices, dices v (8-5-99) Wesley Ebisuzaki usage: wgrib grib file optionsInventory/diagnostic-output selections -s/-v short/verbose inventory -V diagnostic output (not inventory) (none) regular inventory Options -PDS/-PDS10 print PDS in hex/decimal -GDS/-GDS10 print GDS in hex/decimal -verf print forecast verification time -ncep_opn/-ncep_rean default T62 NCEP grib table -4yr print year using 4 digitsDecoding GRIB selection -d record number|all decode record number -p byte position decode record at byte position -i decode controlled by stdin (inventory list) (none) no decoding Options -text/-ieee/-grib/-bin convert to text/ieee/grib/bin (default) -nh/-h output will have no headers/headers (default) -H output will include PDS and GDS (-bin/-ieee only) -append append to output file -o file output file name, dump is default用wgrib扫描grib数据Doswgrib grib20030731181:0:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=10002:114114:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=9753:228228:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=9504:342342:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=9255:456456:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=9006:570570:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=8507:684684:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=8008:798798:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=7509:912912:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=70010:1027026:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=65011:1141140:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=60012:1255254:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=55013:1369368:d=03073118:HGT:kpds5=7:kpds6=100:kpds7=500. （grib.ctl中溴氧的描述是错误的）121:10411198:d=03073118:O3MR:kpds5=154:kpds6=100:kpds7=100 100hPa122:10509022:d=03073118:O3MR:kpds5=154:kpds6=100:kpds7=70 70hPa123:10606846:d=03073118:O3MR:kpds5=154:kpds6=100:kpds7=50 50hPa124:10712814:d=03073118:O3MR:kpds5=154:kpds6=100:kpds7=30 30hPa125:10818782:d=03073118:O3MR:kpds5=154:kpds6=100:kpds7=20 20hPa126:10932896:d=03073118:O3MR:kpds5=154:kpds6=100:kpds7=10 10hPa.280:23351278:d=03073118:UGRD:kpds5=33:kpds6=105:kpds7=10281:23424666:d=03073118:VGRD:kpds5=34:kpds6=105:kpds7=10282:23498054:d=03073118:GPA:kpds5=27:kpds6=100:kpds7=1000283:23612168:d=03073118:GPA:kpds5=27:kpds6=100:kpds7=500284:23726282:d=03073118:5WAVA:kpds5=230:kpds6=100:kpds7=500总共有284个记录，每个记录代表一个水平层。例如解码第13号记录，500hPa位势高度ACSII码Doswgrib grib2003073118 d 13 -text -o dump.dat解码后的数据按ASCII码存放在dump.dat文件中。dump.dat文件360 181 x方向360个格点；y方向181个格点5274.4 数据按先x从1到360；再y从1到181。 5274.4 问题是y方向第一个点是南极还是北极？5274.4 要回到此问题，需要查看由grib2ctl.pl 命令生成的*.ctl文件中，5274.4 如果是： options yrev y方向是倒置的，即从北极向南极放置5274.45274.45274.45274.4例如解码第13号记录，500hPa位势高度binary格式Doswgrib grib2003073118 d 13 -nh -o dump.datDump.dat相应的ctl文件格式参照grib.ctl构造。dset dump.datundef 9.999E+20title grib2003073118 options yrevydef 181 linear -90.000000 1xdef 360 linear 0.000000 1.000000tdef 1 linear 18Z31jul2003 1mozdef 26 levels1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10 vars 1HGTprs 0 7,100,500 * Geopotential height gpmENDVARS如果要画溴氧混合率需要修改grib.ctl文件dset grib2003073118index grib2003073118.idxundef 9.999E+20title grib2003073118 *produced by grib2ctl vp33edtype grib 3options yrevydef 181 linear -90.000000 1xdef 360 linear 0.000000 1.000000tdef 1 linear 18Z31jul2003 1mozdef 26 levels1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10 vars 98.O3MRprs1 0 154,100,100 * Ozone mixing ratio kg/kgO3MRprs2 0 154,100,70 * Ozone mixing ratio kg/kgO3MRprs3 0 154,100,50 * Ozone mixing ratio kg/kgO3MRprs4 0 154,100,30 * Ozone mixing ratio kg/kgO3MRprs5 0 154,100,20 * Ozone mixing ratio kg/kgO3MRprs6 0 154,100,10 * Ozone mixing ratio kg/kg.再运行gribmap -i grib.ctl 生成新的 *.idx文件即可。7. WRF模式数据处理WRF模式MM5模式都是目前从网上可以下载的气象软件，因此在国内经常可以见到。但这两种模式的数据特点数据的水平网格都不是标准的经纬度网格。需要在ctl文件中加入PDEF定义说明把这种非标准的数据经过GrADS内部的计算转换成标准的经纬度网格数据使用。wrfout_d01_000000.nc.bin_map.ctl文件dset wrfout_d01_000000.nc.bin.ieeeoptions sequential big_endian undef -999999 pdef 130 120 lcc 31.