科学网 › 标签 › mload

标签: mload

相关帖子	版块	作者	回复/查看	最后发表

没有相关内容

相关日志

关于mload处理海洋非潮汐ECCO数据的一个问题: 热度 2 zhoufcumt 2012-3-27 18:53; 如果mload的版本为： Current version: v1.07 12/30/2009 The name of driving file is missing. runstring: mload drvfil 处理海洋非潮汐ECCO数据将出现问题， * ALGO_GPS lon= 281.928600 colat= 44.044200 2001 7 2 12 2452093.00********** 2873.3294********** -0.2133488E+00 2001 7 3 12 2452094.00********** 2873.3939********** -0.2133490E+00 2001 7 4 12 2452095.00********** 2873.4921********** -0.2133492E+00 2001 7 5 12 2452096.00********** 2873.3397********** -0.2133489E+00 2001 7 6 12 2452097.00********** 2873.2047********** -0.2133490E+00 2001 7 7 12 2452098.00********** 2873.0894********** -0.2133491E+00 2001 7 8 12 2452099.00********** 2872.8303********** -0.2133491E+00 2001 7 9 12 2452100.00********** 2872.7016********** -0.2133494E+00 2001 7 10 12 2452101.00********** 2872.9059********** -0.2133495E+00 2001 7 11 12 2452102.00********** 2872.9851********** -0.2133492E+00 2001 7 12 12 2452103.00********** 2873.1177********** -0.2133490E+00 2001 7 13 12 2452104.00********** 2873.1014********** -0.2133495E+00 2001 7 14 12 2452105.00********** 2873.0174********** -0.2133499E+00 2001 7 15 12 2452106.00********** 2872.8788********** -0.2133499E+00 2001 7 16 12 2452107.00********** 2872.6647********** -0.2133496E+00 2001 7 17 12 2452108.00********** 2872.5502********** -0.2133497E+00 2001 7 18 12 2452109.00********** 2872.5314********** -0.2133499E+00 2001 7 19 12 2452110.00********** 2872.4692********** -0.2133497E+00 2001 7 20 12 2452111.00********** 2872.6958********** -0.2133496E+00 2001 7 21 12 2452112.00********** 2872.7277********** -0.2133495E+00 2001 7 22 12 2452113.00********** 2872.8434********** -0.2133494E+00 2001 7 23 12 2452114.00********** 2872.9935********** -0.2133496E+00 2001 7 24 12 2452115.00********** 2872.8649********** -0.2133496E+00 2001 7 25 12 2452116.00********** 2872.8305********** -0.2133495E+00 2001 7 26 12 2452117.00********** 2872.8374********** -0.2133494E+00 2001 7 27 12 2452118.00********** 2872.8707********** -0.2133493E+00 2001 7 28 12 2452119.00********** 2873.0374********** -0.2133493E+00 2001 7 29 12 2452120.00********** 2873.2057********** -0.2133492E+00 2001 7 30 12 2452121.00********** 2873.0368********** -0.2133490E+00 2001 7 31 12 2452122.00********** 2873.0663********** -0.2133490E+00 但是上面的结果出现很多*号，明显不对，起初以为是数据的问题，后来在中国地震局田老师的帮助下，发现是mload软件版本的问题，当用 Current version: v1.05 06/18/2008 The name of driving file is missing. runstring: mload drvfil 这个版本的mload时，程序运行正常。 * ALGO_GPS lon= 281.928600 colat= 44.044200 2001 7 2 12 2452093.00 0.3185 0.1545 0.0621 0.1293137E-06 2001 7 3 12 2452094.00 0.3421 0.1377 -0.0176 0.1407710E-06 2001 7 4 12 2452095.00 0.3727 0.1342 0.0234 0.1591735E-06 2001 7 5 12 2452096.00 0.1318 0.0683 0.0734 0.1121511E-06 2001 7 6 12 2452097.00 0.1061 0.0920 0.0549 0.1474834E-06 2001 7 7 12 2452098.00 0.2334 0.1451 -0.0052 0.1753573E-06 2001 7 8 12 2452099.00 0.0614 0.0821 -0.0448 0.9426171E-07 2001 7 9 12 2452100.00 -0.1500 -0.0143 -0.0509 0.1597562E-07 2001 7 10 12 2452101.00 0.0114 -0.0047 0.0585 0.3456892E-07 2001 7 11 12 2452102.00 0.1186 0.0514 0.0243 0.9984900E-07 2001 7 12 12 2452103.00 0.1906 0.0742 -0.0379 0.1153199E-06 2001 7 13 12 2452104.00 0.1172 0.0663 -0.0625 0.1252564E-06 2001 7 14 12 2452105.00 0.1867 0.0919 -0.1185 0.1520511E-06 2001 7 15 12 2452106.00 0.2222 0.1137 -0.0265 0.1414194E-06 2001 7 16 12 2452107.00 0.1469 0.0652 -0.0104 0.6827845E-07 2001 7 17 12 2452108.00 0.2163 0.0494 -0.0507 0.4460967E-07 2001 7 18 12 2452109.00 0.2948 0.0558 0.0297 0.6096472E-07 2001 7 19 12 2452110.00 0.1338 0.0161 0.0841 0.8744006E-10 2001 7 20 12 2452111.00 0.1505 0.0341 0.0157 0.3220546E-07 2001 7 21 12 2452112.00 0.0634 0.0334 0.0310 0.2757080E-07 2001 7 22 12 2452113.00 0.0513 0.0406 0.0705 0.3372522E-07 2001 7 23 12 2452114.00 0.1440 0.0789 0.1073 0.7496755E-07 2001 7 24 12 2452115.00 0.1028 0.0413 0.0959 0.3255861E-07 2001 7 25 12 2452116.00 0.0732 0.0269 0.1139 0.3924788E-07 2001 7 26 12 2452117.00 0.1020 0.0518 0.0803 0.3385971E-07 2001 7 27 12 2452118.00 0.1059 0.0853 0.0102 0.6758766E-07 2001 7 28 12 2452119.00 0.1344 0.1131 -0.0449 0.1140044E-06 2001 7 29 12 2452120.00 0.0816 0.0853 0.0393 0.1060210E-06 2001 7 30 12 2452121.00 -0.1202 0.0226 0.1370 0.1010329E-07 2001 7 31 12 2452122.00 -0.0720 0.0483 0.1141 0.3242029E-07 数据格式应为： Year month day hour JD vertical north east gravity，单位为毫米和gal。非常感谢田老师的帮助，并附上 v1.05 06/18/2008版本的mload。 mload_linux; 个人分类: qoca|4346 次阅读|3 个评论

读取GLDAS数据出错问题的解决: 热度 4 zhoufcumt 2011-10-18 10:44; 用wgrib读取GLDAS的grb格式的数据时，会出现如下错误：Undefined parameter table (center 221-221 table 1), using NCEP-opn。解决方案如下： wgrib的安装从ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib/下载最新的安装包，目前最新的是 wgrib.tar.v1.8.1.2a ，解压后编辑makefile， CC=gcc make即可，将生成的可执行程序wgrib拷贝到/usr/local/bin下即可。终端输入wgrib显示如下： Portable Grib decoder for NCEP/NCAR Reanalysis etc. it slices, dices v1.8.1.0b (5-10) Wesley Ebisuzaki DWD-tables 2,201-205 (11-28-2005) Helmut P. Frank spectral: Luis Kornblueh (MPI) usage: wgrib Inventory/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 digits -min print minutes -ncep_ens ensemble info encoded in ncep format Decoding GRIB selection -d decode record number -p 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) -dwdgrib output dwd headers, grib (do not append) -H output will include PDS and GDS (-bin/-ieee only) -append append to output file -o output file name, 'dump' is default Misc -cmc use NCEP tables for CMC (dangerous) 即证明正确安装。用一例子验证下，wgrib GLDAS_NOAH10_M.A201001.001.grb 显示如下： 1:0:d=10010100:NSWRS:kpds5=111:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 2:35306:d=10010100:NLWRS:kpds5=112:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 3:68714:d=10010100:LHTFL:kpds5=121:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 4:104020:d=10010100:SHTFL:kpds5=122:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 Undefined parameter table (center 221-221 table 1), using NCEP-opn 5:139326:d=10010100:GFLUX:kpds5=155:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 6:168938:d=10010100:LFTX:kpds5=131:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 7:202346:d=10010100:4LFTX:kpds5=132:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 8:239550:d=10010100:EVP:kpds5=57:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 9:271060:d=10010100:SSRUN:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 10:304468:d=10010100:BGRUN:kpds5=234:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 11:339774:d=10010100:SNOM:kpds5=99:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 12:371284:d=10010100:BVF2:kpds5=138:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 13:404692:d=10010100:WEASD:kpds5=65:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 14:457078:d=10010100:TSOIL:kpds5=85:kpds6=112:kpds7=4:TR=0:P1=0:P2=0:TimeU=1:0-4 cm down:anl:NAve=12 15:490486:d=10010100:TSOIL:kpds5=85:kpds6=112:kpds7=3:TR=0:P1=0:P2=0:TimeU=1:0-3 cm down:anl:NAve=12 16:523894:d=10010100:TSOIL:kpds5=85:kpds6=112:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:0-2 cm down:anl:NAve=12 17:557302:d=10010100:TSOIL:kpds5=85:kpds6=112:kpds7=1:TR=0:P1=0:P2=0:TimeU=1:0-1 cm down:anl:NAve=12 18:588812:d=10010100:SOILM:kpds5=86:kpds6=112:kpds7=4:TR=0:P1=0:P2=0:TimeU=1:0-4 cm down:anl:NAve=12 19:626016:d=10010100:SOILM:kpds5=86:kpds6=112:kpds7=3:TR=0:P1=0:P2=0:TimeU=1:0-3 cm down:anl:NAve=12 20:663220:d=10010100:SOILM:kpds5=86:kpds6=112:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:0-2 cm down:anl:NAve=12 21:700424:d=10010100:SOILM:kpds5=86:kpds6=112:kpds7=1:TR=0:P1=0:P2=0:TimeU=1:0-1 cm down:anl:NAve=12 22:737628:d=10010100:TCDC:kpds5=71:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 23:769138:d=10010100:WIND:kpds5=32:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 24:794956:d=10010100:TMP:kpds5=11:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 25:826466:d=10010100:SPFH:kpds5=51:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 26:842794:d=10010100:PRES:kpds5=1:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 27:879998:d=10010100:DSWRF:kpds5=204:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 28:917202:d=10010100:DLWRF:kpds5=205:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 出现了这样的错误Undefined parameter table (center 221-221 table 1), using NCEP-opn 无论加 -necp_opn 还是 -necp_rean参数都出现上述错误。解决方案如下（以我的电脑设置为主，过程一致，你可以参照你的电脑设置，只要你掌握基本的linux知识）：http://disc.sci.gsfc.nasa.gov/hydrology/documentation/hydro_doc.shtml#NOAH下载GRIBTAB for Noah output 表文件 gribtab_GLDAS_NOAH.txt重命名为 gribtab_GLDAS_NOAH 放入root下的data目录，在.bashrc下设置 # GRIBTAB for Noah output export GRIBTAB=/root/data/gribtab_GLDAS_NOAH 这样，重新运行下 wgrib GLDAS_NOAH10_M.A201001.001.grb 显示： 1:0:d=10010100:SWnet:kpds5=111:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 2:35306:d=10010100:LWnet:kpds5=112:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 3:68714:d=10010100:Qle:kpds5=121:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 4:104020:d=10010100:Qh:kpds5=122:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 5:139326:d=10010100:Qg:kpds5=155:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 6:168938:d=10010100:Snowf:kpds5=131:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 7:202346:d=10010100:Rainf:kpds5=132:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 8:239550:d=10010100:Evap:kpds5=57:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 9:271060:d=10010100:Qs:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 10:304468:d=10010100:Qsb:kpds5=234:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 11:339774:d=10010100:Qsm:kpds5=99:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 12:371284:d=10010100:AvgSurfT:kpds5=138:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 13:404692:d=10010100:SWE:kpds5=65:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 14:457078:d=10010100:SoilTemp4:kpds5=85:kpds6=112:kpds7=4:TR=0:P1=0:P2=0:TimeU=1:0-4 cm down:anl:NAve=12 15:490486:d=10010100:SoilTemp4:kpds5=85:kpds6=112:kpds7=3:TR=0:P1=0:P2=0:TimeU=1:0-3 cm down:anl:NAve=12 16:523894:d=10010100:SoilTemp4:kpds5=85:kpds6=112:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:0-2 cm down:anl:NAve=12 17:557302:d=10010100:SoilTemp4:kpds5=85:kpds6=112:kpds7=1:TR=0:P1=0:P2=0:TimeU=1:0-1 cm down:anl:NAve=12 18:588812:d=10010100:SoilMoist4:kpds5=86:kpds6=112:kpds7=4:TR=0:P1=0:P2=0:TimeU=1:0-4 cm down:anl:NAve=12 19:626016:d=10010100:SoilMoist4:kpds5=86:kpds6=112:kpds7=3:TR=0:P1=0:P2=0:TimeU=1:0-3 cm down:anl:NAve=12 20:663220:d=10010100:SoilMoist4:kpds5=86:kpds6=112:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:0-2 cm down:anl:NAve=12 21:700424:d=10010100:SoilMoist4:kpds5=86:kpds6=112:kpds7=1:TR=0:P1=0:P2=0:TimeU=1:0-1 cm down:anl:NAve=12 22:737628:d=10010100:Canopint:kpds5=71:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 23:769138:d=10010100:Wind:kpds5=32:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 24:794956:d=10010100:Tair:kpds5=11:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 25:826466:d=10010100:Qair:kpds5=51:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 26:842794:d=10010100:PSurf:kpds5=1:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12 27:879998:d=10010100:SWdown:kpds5=204:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 28:917202:d=10010100:LWdown:kpds5=205:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12 错误解决，且NOAH提供的28个变量正确，正确解决这些问题，利用qoca下的模块mload将可以利用GLDAS下的数据做各项改正。wgrib的命令如何使用，可参照相关说明文档，在此不做说明！相关的网址链接： ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib/ http://www.cpc.ncep.noaa.gov/products/wesley/wgrib.html http://agdisc.gsfc.nasa.gov/dods/GLDAS_NOAH10_M http://disc.sci.gsfc.nasa.gov/hydrology/documentation/hydro_doc.shtml#NOAH http://www.nco.ncep.noaa.gov/pmb/docs/on388/; 个人分类: 空间大地测量相关软件|8383 次阅读|7 个评论

更多...

帐号		自动登录	找回密码
密码			注册

关闭 安全验证

标签: mload

相关帖子

相关日志

关闭安全验证