nc数据处理

chongzika

用write_matrix可以么？

longlivehj

http://www.ncl.ucar.edu/Applications/o-netcdf.shtml
里面有两种方法，都可以！

chongzika

longlivehj 发表于 2014-5-21 11:00
http://www.ncl.ucar.edu/Applications/o-netcdf.shtml
里面有两种方法，都可以！

float SWDNBC ( Time, south_north, west_east )
         FieldType :        104
         MemoryOrder :        XY
         description :        INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM
         units :        W m-2
         stagger :        
         coordinates :        XLONG XLAT
我的变量大都是这样的3维数据，好像直接用第一个例子不行吧，我这不是经纬度的坐标。需要重新定义么？
Variable: f
Type: file
filename:        wrfout_d02_2001-04-07_00:00:00
path:        wrfout_d02_2001-04-07_00:00:00
   file global attributes:
      TITLE :  OUTPUT FROM WRF V3.4 MODEL
      START_DATE : 2001-04-07_00:00:00
      SIMULATION_START_DATE : 2001-04-07_00:00:00
      WEST-EAST_GRID_DIMENSION : 211
      SOUTH-NORTH_GRID_DIMENSION : 211
      BOTTOM-TOP_GRID_DIMENSION : 28
      DX : 18000
      DY : 18000

longlivehj

chongzika 发表于 2014-5-21 21:17
float SWDNBC ( Time, south_north, west_east )
         FieldType :        104
         MemoryOrd ...

第一种方法可以的。在把变量写入nc的同时，把坐标变量也一并读出并写入nc，不需要重新定义的。

chongzika

longlivehj 发表于 2014-5-22 08:40
第一种方法可以的。在把变量写入nc的同时，把坐标变量也一并读出并写入nc，不需要重新定义的。

像GLW变量这样的维数怎么定义的呢？
time=f->time
south_north=f->south_north
west_east=f->west_east
我这样不能读呢？错误显示：

dimensions:

     Time = 73  // unlimited

     DateStrLen = 19

     west_east = 210

     south_north = 210

     bottom_top = 27

     bottom_top_stag = 28

     soil_layers_stag = 4

     west_east_stag = 211

     south_north_stag = 211

longlivehj

chongzika 发表于 2014-5-22 14:27
像GLW变量这样的维数怎么定义的呢？
time=f->time
south_north=f->south_north

是 XLONG和XLAT！