MeteoInfoLab脚本示例：wrfout转arl

李舒文

MeteoInfo 发表于 2017-5-17 20:23
参考这里：https://wenzhang.baidu.com/page/view?key=2caabb2fed30699d-1426755110

谢谢王老师，但是这个网页好像打不开~

MeteoInfo

李舒文 发表于 2017-5-17 20:37
谢谢王老师，但是这个网页好像打不开~

最近单位共享的数据里有T213模式的产品数据，是GRIB1格式。不过数据文件非常零碎，一个文件只包含了一个格点场（不知道为什么处理成这样）。好在GRIB数据是由一个个Message组成的，每个Message都是自描述的，没有统一的文件头什么的，因此数据文件的合并就非常简单了。在Windows系统中可以用DOS命令Copy来合并文件，例如：copy /b file1.grib+file2.grib newfile.grib，会将file1.grib和file2.grib合并为newfile.grib文件，/b参数表明数据是二进制的。文件较多的情况下可以用通配符，如 copy /b *.grib newfile.grib。Linux系统中用cat命令，如：cat file1.grib file2.grib > newfile.grib。



这种合并方式对于ARL格式数据也适用，因为每个网格数据也是自描述的。但对于NetCDF数据就不适用了。

李舒文

MeteoInfo 发表于 2017-5-17 20:42
最近单位共享的数据里有T213模式的产品数据，是GRIB1格式。不过数据文件非常零碎，一个文件只包含了一个 ...

好的，非常感谢王老师，我先去试试~

李舒文

李舒文 发表于 2017-5-17 20:46
好的，非常感谢王老师，我先去试试~

王老师您好，我发现用meteoinfo将wrfout转成arl格式后还是有点问题，但是用其他的程序转成arl是可以用的。错误提示为：*fatal error* main:all source points off all meteo grids! See MESSAGE file for more information.

MeteoInfo

李舒文 发表于 2017-5-17 21:27
王老师您好，我发现用meteoinfo将wrfout转成arl格式后还是有点问题，但是用其他的程序转成arl是可以用的 ...

你是怎么转的？

李舒文

MeteoInfo 发表于 2017-5-17 21:42
你是怎么转的？

http://bbs.06climate.com/forum.p ... p;extra=&page=1
我用您的这个脚本，就是改了一下路径，一个时次的wrfout转成一个arl文件，然后合并成一个arl文件。

李舒文

这个是我用来对比的方法，这个是可以运行通的，但是处理比较大的wrfout时会出错。http://bbs.06climate.com/forum.p ... p;extra=&page=3
我分别顺便贴出两个时次处理完后的arl文件。
第一个是用meteoinfo处理的，第二个是参考这个方法处理的。

*****************************************************************************
12062000000099INDX   0 0.0000000E+00 0.0000000E+00AWRF  0 090.00000.0000035.4999110.00040.50000.0000045.69001.000001.0000035.4999110.0000.00000150179 30 118281.000011SHGT 62 PRSS200 PBLH 90 USTR 90 DSWF 90 SHTF 90 LHTF 90 T02M207 U10M184 V10M 15 TPPA 90 0.9965 6PRES173 TEMP254 UWND169 VWND 45 WWND 64 SPHU250 0.9879 6PRES104 TEMP246 UWND182 VWND124 WWND 71 SPHU 41 0.9765 6PRES 17 TEMP157 UWND 37 VWND168 WWND 78 SPHU170 0.9620 6PRES155 TEMP132 UWND 66 VWND237 WWND 90 SPHU 69 0.9440 6PRES  3 TEMP127 UWND145 VWND119 WWND103 SPHU196 0.9214 6PRES 89 TEMP 46 UWND 74 VWND 16 WWND123 SPHU144 0.8945 6PRES119 TEMP130 UWND 82 VWND 85 WWND143 SPHU 46 0.8558 6PRES 66 TEMP 21 UWND  1 VWND 61 WWND167 SPHU241 0.8076 6PRES171 TEMP140 UWND223 VWND132 WWND 46 SPHU100 0.7593 6PRES 30 TEMP126 UWND196 VWND 47 WWND 94 SPHU251 0.7111 6PRES237 TEMP139 UWND 30 VWND153 WWND135 SPHU234 0.6452 6PRES197 TEMP147 UWND 22 VWND 88 WWND 42 SPHU 63 0.5657 6PRES214 TEMP166 UWND197 VWND170 WWND138 SPHU112 0.4936 6PRES 83 TEMP160 UWND113 VWND 99 WWND 61 SPHU 53 0.4285 6PRES 55 TEMP 90 UWND102 VWND128 WWND148 SPHU254 0.3698 6PRES178 TEMP105 UWND110 VWND134 WWND 83 SPHU190 0.3169 6PRES 29 TEMP250 UWND 53 VWND229 WWND177 SPHU103 0.2693 6PRES 35 TEMP194 UWND 59 VWND202 WWND124 SPHU 85 0.2267 6PRES209 TEMP103 UWND101 VWND250 WWND 75 SPHU137 0.1886 6PRES 15 TEMP210 UWND144 VWND179 WWND  8 SPHU238 0.1547 6PRES 53 TEMP173 UWND141 VWND144 WWND151 SPHU 73 0.1245 6PRES149 TEMP218 UWND211 VWND 93 WWND124 SPHU 99 0.0981 6PRES253 TEMP 99 UWND 92 VWND142 WWND102 SPHU177 0.0756 6PRES 75 TEMP 23 UWND131 VWND 40 WWND 99 SPHU192 0.0565 6PRES191 TEMP  8 UWND141 VWND 94 WWND117 SPHU241 0.0402 6PRES194 TEMP187 UWND 17 VWND180 WWND157 SPHU211 0.0264 6PRES231 TEMP105 UWND221 VWND117 WWND106 SPHU188 0.0146 6PRES188 TEMP144 UWND225 VWND107 WWND179 SPHU 90 0.0046 6PRES 94 TEMP 35 UWND 63 VWND 19 WWND 90 SPHU 90



****************************************************************************************
12 620 0 0 099INDX   0 0.0000000E+00 0.0000000E+00AWRF  0 045.6900110.00045.6900110.00040.5000.00000045.690075.500090.000035.5000110.000.000000150179 30 120681.000012SHGT 62 PRSS200 TPP1 90 DIFR 90 PBLH 90 USTR 90 DSWF 90 SHTF 90 LHTF 90 T02M207 U10M184 V10M 15 .99650 7PRES173 TEMP254 UWND169 VWND 45 WWND 67 DIFW156 SPHU250 .98800 7PRES104 TEMP246 UWND182 VWND124 WWND 72 DIFW 64 SPHU 41 .97650 7PRES 17 TEMP157 UWND 37 VWND168 WWND 77 DIFW248 SPHU170 .96200 7PRES155 TEMP132 UWND 66 VWND237 WWND 98 DIFW126 SPHU 69 .94400 7PRES  3 TEMP127 UWND145 VWND119 WWND111 DIFW148 SPHU196 .92150 7PRES 89 TEMP 46 UWND 74 VWND 16 WWND126 DIFW 60 SPHU144 .89450 7PRES119 TEMP130 UWND 82 VWND 85 WWND150 DIFW 56 SPHU 46 .85588 7PRES 66 TEMP 21 UWND  1 VWND 61 WWND178 DIFW 64 SPHU241 .80763 7PRES171 TEMP140 UWND223 VWND132 WWND 61 DIFW186 SPHU100 .75938 7PRES 30 TEMP126 UWND196 VWND 47 WWND106 DIFW 64 SPHU251 .71113 7PRES237 TEMP139 UWND 30 VWND153 WWND148 DIFW207 SPHU234 .64528 7PRES197 TEMP147 UWND 22 VWND 88 WWND 71 DIFW 30 SPHU 63 .56573 7PRES214 TEMP166 UWND197 VWND170 WWND157 DIFW139 SPHU112 .49370 7PRES 83 TEMP160 UWND113 VWND 99 WWND 84 DIFW119 SPHU 53 .42858 7PRES 55 TEMP 90 UWND102 VWND128 WWND171 DIFW236 SPHU254 .36982 7PRES178 TEMP105 UWND110 VWND134 WWND110 DIFW240 SPHU190 .31692 7PRES 29 TEMP250 UWND 53 VWND229 WWND 35 DIFW 35 SPHU103 .26939 7PRES 35 TEMP194 UWND 59 VWND202 WWND152 DIFW232 SPHU 85 .22679 7PRES209 TEMP103 UWND101 VWND250 WWND101 DIFW249 SPHU137 .18870 7PRES 15 TEMP210 UWND144 VWND179 WWND 46 DIFW152 SPHU238 .15473 7PRES 53 TEMP173 UWND141 VWND144 WWND168 DIFW 16 SPHU 73 .12452 7PRES149 TEMP218 UWND211 VWND 93 WWND149 DIFW158 SPHU 99 .09817 7PRES253 TEMP 99 UWND 92 VWND142 WWND135 DIFW 27 SPHU177 .07569 7PRES 75 TEMP 23 UWND131 VWND 40 WWND133 DIFW169 SPHU192 .05656 7PRES191 TEMP  7 UWND141 VWND 94 WWND147 DIFW134 SPHU241 .04028 7PRES194 TEMP187 UWND 17 VWND180 WWND 20 DIFW 10 SPHU211 .02643 7PRES231 TEMP105 UWND221 VWND117 WWND140 DIFW  4 SPHU188 .01464 7PRES188 TEMP144 UWND225 VWND107 WWND 50 DIFW 58 SPHU 90 .00461 7PRES 94 TEMP 36 UWND 63 VWND 19 WWND 90 DIFW 90 SPHU 90

xuzheyong

vwnd = inf['V'][t,lidx,:,:]   老师 请问这种是什么语法？如果是将list切片的话 应该是【：：】这样的 但是带，的就想不明白了

MeteoInfo

xuzheyong 发表于 2017-5-19 12:36
vwnd = inf['V'][t,lidx,:,:]   老师 请问这种是什么语法？如果是将list切片的话 应该是【：：】这样的 但 ...

参考这里：http://www.meteothinker.com/docs ... erguide/arrays.html

xuzheyong

嗯 谢谢老师