文件后缀为.asc
文件头部分格式：
nrows 1232
xllcorner 484314.627259768022
yllcorner 4327241.646642711945
cellsize 9.994747084015
NODATA_value -9999
如果QGIS或者GDAL命令中导出来的asc文件格式是：
nrows 1232
xllcorner 484314.627259768022
yllcorner 4327241.646642711945
dx 9.994747084015
dy 9.994747084015
NODATA_value -9999
则需要手动进行修改一下。
如果要自动化执行填充，注意文件格式的不同。或者修改anuga的源代码
anuga.asc2dem()方法，但修改的话有风险，因为不知道这个方法是否也用在别处。
2)很奇怪的是第二种格式的asc，也就是anuga使用的asc格式，QGIS也能打开。
所以最佳的修改方式是anuga源码中修改的都兼容。（尚未修改20200222）
相关asc的投影prj仅仅支持下面的格式（源代码中默认的）
Projection UTM
Zone 55
Datum WGS84
Zunits NO
Units m
Spheroid WGS84
Xshift 500000
Yshift 10000000
Parameters
但在asc2dem中没有看到投影转换的地方。数据是直接从asc读出来放到dem（一个netcdf文件，不是地理的dem文件）中的。
这个投影文件的使用地方在set_quantity('elevation')的时候。看看到底有什么用。（20200222）
没有仔细研究。投影保持源代码默认的 或者改为如下，
Projection NO
Zone 50
Datum WGS84
Zunits NO
Units m
Spheroid WGS84
Xshift 500000
Yshift 0
Parameters
运行结果：sww数据投影将和输入的原DEM数据坐标一样。所以怀疑投影文件没有发挥作用。
2020-02-27日：用来计算的prj文件保持和源anuga一样的投影，并不影响计算结果。但结果计算出来后，在QGIS显示是，依然使用源tif文件的投影。jumahe案例的投影文件为：
+proj=tmerc +lat_0=0 +lon_0=116 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs
