文件后缀为.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