fwrite提取文件--将两个文件写到一个文件中

芳芳芳 · 发表于 2017-4-7 11:14:04

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

数据为1小时的10m风场，需要将研究区域从2月1日01时（第一个文件，t=2）到3月1日00时（第二个文件，t=1）写入到一个dat文件中。
gs文件如下：

'reinit'
'open D:\e\201602.ctl'
'open D:\e\201603.ctl'
'set gxout fwrite'
'set fwrite d:\ee\201602.dat'
'set lon 117.614 122.523'
'set lat 37.104 40.375'
t=2
while(t<=696)
'set t 't''
'd UGRD10m'
'd VGRD10m'
t = t+1
endwhile

'set lon 117.614 122.523'
'set lat 37.104 40.375'
'set t 1'
'd UGRD10m.2'
'd VGRD10m.2'
'disable fwrite'

但运行后报错，显示语法错误，即gs红色下划线所在行。

请各位大神帮忙看看gs有什么错误，谢谢大家！！！

芳芳芳 · 发表于 2017-4-7 11:15:10

求帮助，谢谢大家~~~~~

四叶草 · 发表于 2017-4-7 11:54:29

是不是格点精度问题？建议楼主贴出报错的信息

芳芳芳 · 发表于 2017-4-7 12:23:51

四叶草发表于 2017-4-7 11:54
是不是格点精度问题？建议楼主贴出报错的信息

file:///C:/Documents%20and%20Settings/fang/桌面/图片1.png

芳芳芳 · 发表于 2017-4-7 12:24:57

四叶草发表于 2017-4-7 11:54
是不是格点精度问题？建议楼主贴出报错的信息

这个是201602.ctl

dset d:\e\wnd10m.gdas.201602.grib2
index d:\e\wnd10m.gdas.201602.grib2.idx
undef 9.999E+20
title d:\e\wnd10m.gdas.201602.grib2
*  produced by g2ctl v0.0.3d
* griddef=1:0:(1760 x 880):grid_template=40: Gaussian grid: (1760 x 880) units 1e-006 input WE:NS output WE:SN number of latitudes between pole-equator=440 #points=1548800 lat 89.844000 to -89.844000 lon 0.000000 to 359.795000
dtype grib2
xdef 1760 linear 0.000000 0.204545454545455
ydef 880 levels
-89.844 -89.641 -89.437 -89.233 -89.028 -88.824 -88.620 -88.415 -88.211 -88.007 -87.802 -87.598 -87.393 -87.189
-86.984 -86.780 -86.576 -86.371 -86.167 -85.962 -85.758 -85.554 -85.349 -85.145 -84.940 -84.736 -84.531 -84.327
-84.123 -83.918 -83.714 -83.509 -83.305 -83.100 -82.896 -82.692 -82.487 -82.283 -82.078 -81.874 -81.669 -81.465
-81.261 -81.056 -80.852 -80.647 -80.443 -80.238 -80.034 -79.830 -79.625 -79.421 -79.216 -79.012 -78.807 -78.603
-78.399 -78.194 -77.990 -77.785 -77.581 -77.376 -77.172 -76.968 -76.763 -76.559 -76.354 -76.150 -75.945 -75.741
-75.537 -75.332 -75.128 -74.923 -74.719 -74.514 -74.310 -74.106 -73.901 -73.697 -73.492 -73.288 -73.083 -72.879
-72.675 -72.470 -72.266 -72.061 -71.857 -71.652 -71.448 -71.244 -71.039 -70.835 -70.630 -70.426 -70.221 -70.017
-69.813 -69.608 -69.404 -69.199 -68.995 -68.790 -68.586 -68.382 -68.177 -67.973 -67.768 -67.564 -67.359 -67.155
-66.951 -66.746 -66.542 -66.337 -66.133 -65.928 -65.724 -65.520 -65.315 -65.111 -64.906 -64.702 -64.497 -64.293
-64.089 -63.884 -63.680 -63.475 -63.271 -63.066 -62.862 -62.658 -62.453 -62.249 -62.044 -61.840 -61.635 -61.431
-61.227 -61.022 -60.818 -60.613 -60.409 -60.204 -60.000 -59.796 -59.591 -59.387 -59.182 -58.978 -58.773 -58.569
-58.365 -58.160 -57.956 -57.751 -57.547 -57.342 -57.138 -56.934 -56.729 -56.525 -56.320 -56.116 -55.911 -55.707
-55.503 -55.298 -55.094 -54.889 -54.685 -54.480 -54.276 -54.072 -53.867 -53.663 -53.458 -53.254 -53.049 -52.845
-52.641 -52.436 -52.232 -52.027 -51.823 -51.618 -51.414 -51.210 -51.005 -50.801 -50.596 -50.392 -50.187 -49.983
-49.779 -49.574 -49.370 -49.165 -48.961 -48.756 -48.552 -48.348 -48.143 -47.939 -47.734 -47.530 -47.325 -47.121
-46.917 -46.712 -46.508 -46.303 -46.099 -45.894 -45.690 -45.486 -45.281 -45.077 -44.872 -44.668 -44.463 -44.259
-44.055 -43.850 -43.646 -43.441 -43.237 -43.032 -42.828 -42.624 -42.419 -42.215 -42.010 -41.806 -41.601 -41.397
-41.192 -40.988 -40.784 -40.579 -40.375 -40.170 -39.966 -39.761 -39.557 -39.353 -39.148 -38.944 -38.739 -38.535
-38.330 -38.126 -37.922 -37.717 -37.513 -37.308 -37.104 -36.899 -36.695 -36.491 -36.286 -36.082 -35.877 -35.673
-35.468 -35.264 -35.060 -34.855 -34.651 -34.446 -34.242 -34.037 -33.833 -33.629 -33.424 -33.220 -33.015 -32.811
-32.606 -32.402 -32.198 -31.993 -31.789 -31.584 -31.380 -31.175 -30.971 -30.767 -30.562 -30.358 -30.153 -29.949
-29.744 -29.540 -29.336 -29.131 -28.927 -28.722 -28.518 -28.313 -28.109 -27.905 -27.700 -27.496 -27.291 -27.087
-26.882 -26.678 -26.474 -26.269 -26.065 -25.860 -25.656 -25.451 -25.247 -25.043 -24.838 -24.634 -24.429 -24.225
-24.020 -23.816 -23.612 -23.407 -23.203 -22.998 -22.794 -22.589 -22.385 -22.181 -21.976 -21.772 -21.567 -21.363
-21.158 -20.954 -20.750 -20.545 -20.341 -20.136 -19.932 -19.727 -19.523 -19.319 -19.114 -18.910 -18.705 -18.501
-18.296 -18.092 -17.888 -17.683 -17.479 -17.274 -17.070 -16.865 -16.661 -16.457 -16.252 -16.048 -15.843 -15.639
-15.434 -15.230 -15.026 -14.821 -14.617 -14.412 -14.208 -14.003 -13.799 -13.595 -13.390 -13.186 -12.981 -12.777
-12.572 -12.368 -12.164 -11.959 -11.755 -11.550 -11.346 -11.141 -10.937 -10.733 -10.528 -10.324 -10.119  -9.915
  -9.710  -9.506  -9.302  -9.097  -8.893  -8.688  -8.484  -8.279  -8.075  -7.871  -7.666  -7.462  -7.257  -7.053
  -6.848  -6.644  -6.440  -6.235  -6.031  -5.826  -5.622  -5.417  -5.213  -5.009  -4.804  -4.600  -4.395  -4.191
  -3.986  -3.782  -3.578  -3.373  -3.169  -2.964  -2.760  -2.555  -2.351  -2.147  -1.942  -1.738  -1.533  -1.329
  -1.124  -0.920  -0.716  -0.511  -0.307  -0.102 0.102 0.307 0.511 0.716 0.920 1.124 1.329 1.533
1.738 1.942 2.147 2.351 2.555 2.760 2.964 3.169 3.373 3.578 3.782 3.986 4.191 4.395
4.600 4.804 5.009 5.213 5.417 5.622 5.826 6.031 6.235 6.440 6.644 6.848 7.053 7.257
7.462 7.666 7.871 8.075 8.279 8.484 8.688 8.893 9.097 9.302 9.506 9.710 9.915  10.119
  10.324  10.528  10.733  10.937  11.141  11.346  11.550  11.755  11.959  12.164  12.368  12.572  12.777  12.981
  13.186  13.390  13.595  13.799  14.003  14.208  14.412  14.617  14.821  15.026  15.230  15.434  15.639  15.843
  16.048  16.252  16.457  16.661  16.865  17.070  17.274  17.479  17.683  17.888  18.092  18.296  18.501  18.705
  18.910  19.114  19.319  19.523  19.727  19.932  20.136  20.341  20.545  20.750  20.954  21.158  21.363  21.567
  21.772  21.976  22.181  22.385  22.589  22.794  22.998  23.203  23.407  23.612  23.816  24.020  24.225  24.429
  24.634  24.838  25.043  25.247  25.451  25.656  25.860  26.065  26.269  26.474  26.678  26.882  27.087  27.291
  27.496  27.700  27.905  28.109  28.313  28.518  28.722  28.927  29.131  29.336  29.540  29.744  29.949  30.153
  30.358  30.562  30.767  30.971  31.175  31.380  31.584  31.789  31.993  32.198  32.402  32.606  32.811  33.015
  33.220  33.424  33.629  33.833  34.037  34.242  34.446  34.651  34.855  35.060  35.264  35.468  35.673  35.877
  36.082  36.286  36.491  36.695  36.899  37.104  37.308  37.513  37.717  37.922  38.126  38.330  38.535  38.739
  38.944  39.148  39.353  39.557  39.761  39.966  40.170  40.375  40.579  40.784  40.988  41.192  41.397  41.601
  41.806  42.010  42.215  42.419  42.624  42.828  43.032  43.237  43.441  43.646  43.850  44.055  44.259  44.463
  44.668  44.872  45.077  45.281  45.486  45.690  45.894  46.099  46.303  46.508  46.712  46.917  47.121  47.325
  47.530  47.734  47.939  48.143  48.348  48.552  48.756  48.961  49.165  49.370  49.574  49.779  49.983  50.187
  50.392  50.596  50.801  51.005  51.210  51.414  51.618  51.823  52.027  52.232  52.436  52.641  52.845  53.049
  53.254  53.458  53.663  53.867  54.072  54.276  54.480  54.685  54.889  55.094  55.298  55.503  55.707  55.911
  56.116  56.320  56.525  56.729  56.934  57.138  57.342  57.547  57.751  57.956  58.160  58.365  58.569  58.773
  58.978  59.182  59.387  59.591  59.796  60.000  60.204  60.409  60.613  60.818  61.022  61.227  61.431  61.635
  61.840  62.044  62.249  62.453  62.658  62.862  63.066  63.271  63.475  63.680  63.884  64.089  64.293  64.497
  64.702  64.906  65.111  65.315  65.520  65.724  65.928  66.133  66.337  66.542  66.746  66.951  67.155  67.359
  67.564  67.768  67.973  68.177  68.382  68.586  68.790  68.995  69.199  69.404  69.608  69.813  70.017  70.221
  70.426  70.630  70.835  71.039  71.244  71.448  71.652  71.857  72.061  72.266  72.470  72.675  72.879  73.083
  73.288  73.492  73.697  73.901  74.106  74.310  74.514  74.719  74.923  75.128  75.332  75.537  75.741  75.945
  76.150  76.354  76.559  76.763  76.968  77.172  77.376  77.581  77.785  77.990  78.194  78.399  78.603  78.807
  79.012  79.216  79.421  79.625  79.830  80.034  80.238  80.443  80.647  80.852  81.056  81.261  81.465  81.669
  81.874  82.078  82.283  82.487  82.692  82.896  83.100  83.305  83.509  83.714  83.918  84.123  84.327  84.531
  84.736  84.940  85.145  85.349  85.554  85.758  85.962  86.167  86.371  86.576  86.780  86.984  87.189  87.393
  87.598  87.802  88.007  88.211  88.415  88.620  88.824  89.028  89.233  89.437  89.641  89.844
tdef 696 linear 00Z01feb2016 1hr
zdef 1 linear 1 1
vars 2
UGRD10m  0,103,10 0,2,2 ** 10 m above ground none U-Component of Wind [m/s]
VGRD10m  0,103,10 0,2,3 ** 10 m above ground none V-Component of Wind [m/s]
ENDVARS

芳芳芳 · 发表于 2017-4-7 12:26:12

这个是201603.ctl
这两个文件的精度是一样的。

dset d:\e\wnd10m.gdas.201603.grib2
index d:\e\wnd10m.gdas.201603.grib2.idx
undef 9.999E+20
title d:\e\wnd10m.gdas.201603.grib2
*  produced by g2ctl v0.0.3d
* griddef=1:0:(1760 x 880):grid_template=40: Gaussian grid: (1760 x 880) units 1e-006 input WE:NS output WE:SN number of latitudes between pole-equator=440 #points=1548800 lat 89.844000 to -89.844000 lon 0.000000 to 359.795000
dtype grib2
xdef 1760 linear 0.000000 0.204545454545455
ydef 880 levels
-89.844 -89.641 -89.437 -89.233 -89.028 -88.824 -88.620 -88.415 -88.211 -88.007 -87.802 -87.598 -87.393 -87.189
-86.984 -86.780 -86.576 -86.371 -86.167 -85.962 -85.758 -85.554 -85.349 -85.145 -84.940 -84.736 -84.531 -84.327
-84.123 -83.918 -83.714 -83.509 -83.305 -83.100 -82.896 -82.692 -82.487 -82.283 -82.078 -81.874 -81.669 -81.465
-81.261 -81.056 -80.852 -80.647 -80.443 -80.238 -80.034 -79.830 -79.625 -79.421 -79.216 -79.012 -78.807 -78.603
-78.399 -78.194 -77.990 -77.785 -77.581 -77.376 -77.172 -76.968 -76.763 -76.559 -76.354 -76.150 -75.945 -75.741
-75.537 -75.332 -75.128 -74.923 -74.719 -74.514 -74.310 -74.106 -73.901 -73.697 -73.492 -73.288 -73.083 -72.879
-72.675 -72.470 -72.266 -72.061 -71.857 -71.652 -71.448 -71.244 -71.039 -70.835 -70.630 -70.426 -70.221 -70.017
-69.813 -69.608 -69.404 -69.199 -68.995 -68.790 -68.586 -68.382 -68.177 -67.973 -67.768 -67.564 -67.359 -67.155
-66.951 -66.746 -66.542 -66.337 -66.133 -65.928 -65.724 -65.520 -65.315 -65.111 -64.906 -64.702 -64.497 -64.293
-64.089 -63.884 -63.680 -63.475 -63.271 -63.066 -62.862 -62.658 -62.453 -62.249 -62.044 -61.840 -61.635 -61.431
-61.227 -61.022 -60.818 -60.613 -60.409 -60.204 -60.000 -59.796 -59.591 -59.387 -59.182 -58.978 -58.773 -58.569
-58.365 -58.160 -57.956 -57.751 -57.547 -57.342 -57.138 -56.934 -56.729 -56.525 -56.320 -56.116 -55.911 -55.707
-55.503 -55.298 -55.094 -54.889 -54.685 -54.480 -54.276 -54.072 -53.867 -53.663 -53.458 -53.254 -53.049 -52.845
-52.641 -52.436 -52.232 -52.027 -51.823 -51.618 -51.414 -51.210 -51.005 -50.801 -50.596 -50.392 -50.187 -49.983
-49.779 -49.574 -49.370 -49.165 -48.961 -48.756 -48.552 -48.348 -48.143 -47.939 -47.734 -47.530 -47.325 -47.121
-46.917 -46.712 -46.508 -46.303 -46.099 -45.894 -45.690 -45.486 -45.281 -45.077 -44.872 -44.668 -44.463 -44.259
-44.055 -43.850 -43.646 -43.441 -43.237 -43.032 -42.828 -42.624 -42.419 -42.215 -42.010 -41.806 -41.601 -41.397
-41.192 -40.988 -40.784 -40.579 -40.375 -40.170 -39.966 -39.761 -39.557 -39.353 -39.148 -38.944 -38.739 -38.535
-38.330 -38.126 -37.922 -37.717 -37.513 -37.308 -37.104 -36.899 -36.695 -36.491 -36.286 -36.082 -35.877 -35.673
-35.468 -35.264 -35.060 -34.855 -34.651 -34.446 -34.242 -34.037 -33.833 -33.629 -33.424 -33.220 -33.015 -32.811
-32.606 -32.402 -32.198 -31.993 -31.789 -31.584 -31.380 -31.175 -30.971 -30.767 -30.562 -30.358 -30.153 -29.949
-29.744 -29.540 -29.336 -29.131 -28.927 -28.722 -28.518 -28.313 -28.109 -27.905 -27.700 -27.496 -27.291 -27.087
-26.882 -26.678 -26.474 -26.269 -26.065 -25.860 -25.656 -25.451 -25.247 -25.043 -24.838 -24.634 -24.429 -24.225
-24.020 -23.816 -23.612 -23.407 -23.203 -22.998 -22.794 -22.589 -22.385 -22.181 -21.976 -21.772 -21.567 -21.363
-21.158 -20.954 -20.750 -20.545 -20.341 -20.136 -19.932 -19.727 -19.523 -19.319 -19.114 -18.910 -18.705 -18.501
-18.296 -18.092 -17.888 -17.683 -17.479 -17.274 -17.070 -16.865 -16.661 -16.457 -16.252 -16.048 -15.843 -15.639
-15.434 -15.230 -15.026 -14.821 -14.617 -14.412 -14.208 -14.003 -13.799 -13.595 -13.390 -13.186 -12.981 -12.777
-12.572 -12.368 -12.164 -11.959 -11.755 -11.550 -11.346 -11.141 -10.937 -10.733 -10.528 -10.324 -10.119  -9.915
  -9.710  -9.506  -9.302  -9.097  -8.893  -8.688  -8.484  -8.279  -8.075  -7.871  -7.666  -7.462  -7.257  -7.053
  -6.848  -6.644  -6.440  -6.235  -6.031  -5.826  -5.622  -5.417  -5.213  -5.009  -4.804  -4.600  -4.395  -4.191
  -3.986  -3.782  -3.578  -3.373  -3.169  -2.964  -2.760  -2.555  -2.351  -2.147  -1.942  -1.738  -1.533  -1.329
  -1.124  -0.920  -0.716  -0.511  -0.307  -0.102 0.102 0.307 0.511 0.716 0.920 1.124 1.329 1.533
1.738 1.942 2.147 2.351 2.555 2.760 2.964 3.169 3.373 3.578 3.782 3.986 4.191 4.395
4.600 4.804 5.009 5.213 5.417 5.622 5.826 6.031 6.235 6.440 6.644 6.848 7.053 7.257
7.462 7.666 7.871 8.075 8.279 8.484 8.688 8.893 9.097 9.302 9.506 9.710 9.915  10.119
  10.324  10.528  10.733  10.937  11.141  11.346  11.550  11.755  11.959  12.164  12.368  12.572  12.777  12.981
  13.186  13.390  13.595  13.799  14.003  14.208  14.412  14.617  14.821  15.026  15.230  15.434  15.639  15.843
  16.048  16.252  16.457  16.661  16.865  17.070  17.274  17.479  17.683  17.888  18.092  18.296  18.501  18.705
  18.910  19.114  19.319  19.523  19.727  19.932  20.136  20.341  20.545  20.750  20.954  21.158  21.363  21.567
  21.772  21.976  22.181  22.385  22.589  22.794  22.998  23.203  23.407  23.612  23.816  24.020  24.225  24.429
  24.634  24.838  25.043  25.247  25.451  25.656  25.860  26.065  26.269  26.474  26.678  26.882  27.087  27.291
  27.496  27.700  27.905  28.109  28.313  28.518  28.722  28.927  29.131  29.336  29.540  29.744  29.949  30.153
  30.358  30.562  30.767  30.971  31.175  31.380  31.584  31.789  31.993  32.198  32.402  32.606  32.811  33.015
  33.220  33.424  33.629  33.833  34.037  34.242  34.446  34.651  34.855  35.060  35.264  35.468  35.673  35.877
  36.082  36.286  36.491  36.695  36.899  37.104  37.308  37.513  37.717  37.922  38.126  38.330  38.535  38.739
  38.944  39.148  39.353  39.557  39.761  39.966  40.170  40.375  40.579  40.784  40.988  41.192  41.397  41.601
  41.806  42.010  42.215  42.419  42.624  42.828  43.032  43.237  43.441  43.646  43.850  44.055  44.259  44.463
  44.668  44.872  45.077  45.281  45.486  45.690  45.894  46.099  46.303  46.508  46.712  46.917  47.121  47.325
  47.530  47.734  47.939  48.143  48.348  48.552  48.756  48.961  49.165  49.370  49.574  49.779  49.983  50.187
  50.392  50.596  50.801  51.005  51.210  51.414  51.618  51.823  52.027  52.232  52.436  52.641  52.845  53.049
  53.254  53.458  53.663  53.867  54.072  54.276  54.480  54.685  54.889  55.094  55.298  55.503  55.707  55.911
  56.116  56.320  56.525  56.729  56.934  57.138  57.342  57.547  57.751  57.956  58.160  58.365  58.569  58.773
  58.978  59.182  59.387  59.591  59.796  60.000  60.204  60.409  60.613  60.818  61.022  61.227  61.431  61.635
  61.840  62.044  62.249  62.453  62.658  62.862  63.066  63.271  63.475  63.680  63.884  64.089  64.293  64.497
  64.702  64.906  65.111  65.315  65.520  65.724  65.928  66.133  66.337  66.542  66.746  66.951  67.155  67.359
  67.564  67.768  67.973  68.177  68.382  68.586  68.790  68.995  69.199  69.404  69.608  69.813  70.017  70.221
  70.426  70.630  70.835  71.039  71.244  71.448  71.652  71.857  72.061  72.266  72.470  72.675  72.879  73.083
  73.288  73.492  73.697  73.901  74.106  74.310  74.514  74.719  74.923  75.128  75.332  75.537  75.741  75.945
  76.150  76.354  76.559  76.763  76.968  77.172  77.376  77.581  77.785  77.990  78.194  78.399  78.603  78.807
  79.012  79.216  79.421  79.625  79.830  80.034  80.238  80.443  80.647  80.852  81.056  81.261  81.465  81.669
  81.874  82.078  82.283  82.487  82.692  82.896  83.100  83.305  83.509  83.714  83.918  84.123  84.327  84.531
  84.736  84.940  85.145  85.349  85.554  85.758  85.962  86.167  86.371  86.576  86.780  86.984  87.189  87.393
  87.598  87.802  88.007  88.211  88.415  88.620  88.824  89.028  89.233  89.437  89.641  89.844
tdef 745 linear 00Z01mar2016 1hr
zdef 1 linear 1 1
vars 2
UGRD10m  0,103,10 0,2,2 ** 10 m above ground none U-Component of Wind [m/s]
VGRD10m  0,103,10 0,2,3 ** 10 m above ground none V-Component of Wind [m/s]
ENDVARS

lgz · 发表于 2017-4-7 13:05:58

建议楼主贴出报错的信息

四叶草 · 发表于 2017-4-7 13:43:03

芳芳芳发表于 2017-4-7 12:23

遇到这种情况，我一般都是直接把报错的语句注释，之后再一一排除。

芳芳芳 · 发表于 2017-4-7 13:47:30

四叶草发表于 2017-4-7 13:43
遇到这种情况，我一般都是直接把报错的语句注释，之后再一一排除。

我把第二段的经纬度设置都去掉之后还是报错，同样是语法错误，错误出现在'set t 1'的上一行，而且这行是空行。。

芳芳芳 · 发表于 2017-4-7 13:47:55

lgz 发表于 2017-4-7 13:05
建议楼主贴出报错的信息

已经贴出来啦

		自动登录	找回密码
密码			立即注册

[脚本编辑] fwrite提取文件--将两个文件写到一个文件中

登录后查看更多精彩内容~

回帖奖励 +1 金钱