- 积分
- 274
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2020-9-22
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2021-6-15 11:14:28
|
显示全部楼层
谢谢,感谢指导!
我的LAIv数据是mapgis做的,但是导不出nc格式文件(我尝试过,没有成功)。导出格式为tif,所以我才用matlab提取数据去创建nc文件(ncl 应该更加合理)。改了matlab代码可以读取成功,但bio_megan没有插值成功,您可知这会是什么导致的呢?下面是数据的一个对比,以及运行时的提示信息。请你帮帮看一下。谢谢!
bio_emiss.f90的源代码我已经改过了。
Reading megan2 bio emiss file china_laiv201605.nc
megan2_bioemiss: nlon_megan, nlat_megan = 1234 945
megan2_bioemiss: data_dx,dx,has_area_map = -5559.21289 27000.0000 T
megan2_bioemiss: file china_laiv201605.nc has a new grid
megan2_bioemiss: xcen_megan(1,2) = 73.4499969482422 73.5000000000000
megan2_bioemiss: xedge_megan(1,2) = 73.4249954223633 73.4749984741211
megan2_bioemiss: dx = 19
megan2_bioemiss: nlon_megan,nlat_megan = 1234 945
megan2_bioemiss: ycen_megan = 6.40999984741211 6.36000013351440
megan2_bioemiss: yedge_megan = 6.38500022888184 6.33500027656555
area_interp: diagnostics
area_interp: min model lat ndx = ******
area_interp: max model lat ndx = ******
area_interp: mcells outside data domain = 17731
area_interp: wght size = 0.00000000000000
area_interp: total dcell cnt = 0
area_interp: active dcell cnt = 0
area_interp: interior dcell cnt = 0
area_interp: partial dcell cnt = 0
Finished megan2 bio emiss dataset china_laiv201605.nc |
-
-
|