爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 15842|回复: 6

请问bio_emiss 怎么制作LAIv等相关输入数据

[复制链接]

新浪微博达人勋

发表于 2021-6-12 20:53:10 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
我准备生成中国地区的LAIv的nc数据文件。
我已经试验过示例数据可以运行成功。我编辑bio_emiss.f90文件,使其读取新的laiv数据文件。
这个nc文件是matlab创建的,报错信息以及nc文件如下:
请问各位有什么指导建议吗?
报错.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2021-6-15 04:02:18 | 显示全部楼层
仔细对比你的数据和示例数据之间的差别吧
可以注意一下示例数据是GIS做的,ncdump即可看到相关信息。
另外,LAI产品不同的话时间也不一样,是要改源码的。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-6-15 11:14:28 | 显示全部楼层
qq531884069 发表于 2021-6-15 04:02
仔细对比你的数据和示例数据之间的差别吧
可以注意一下示例数据是GIS做的,ncdump即可看到相关信息。
另 ...

谢谢,感谢指导!
我的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
mydata.png
example.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-6-15 15:03:06 | 显示全部楼层
报错写的很清楚了吧
是area_interp这个函数读取计算错误
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-6-17 11:19:52 | 显示全部楼层
qq531884069 发表于 2021-6-15 15:03
报错写的很清楚了吧
是area_interp这个函数读取计算错误

感谢指导,我查看了area_interp的函数,已经解决了问题。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-6-14 14:56:50 | 显示全部楼层
请问楼主,area_interp函数具体怎么修改
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-9-16 15:49:24 | 显示全部楼层
muziliyang 发表于 2021-6-15 11:14
谢谢,感谢指导!
我的LAIv数据是mapgis做的,但是导不出nc格式文件(我尝试过,没有成功)。导出格式为 ...

楼主感谢你的问题,不知道问题 解决了没,我的数据也是插值没有成功,全是0您,想请教一下您当时是怎么解决的。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表