WRF、MEGAN以及CMAQ答疑帖

loftydongshi

qq531884069 发表于 2021-3-26 10:47
这个比较看你Linux基础还有对空气质量模型的了解吧  有人回答问题的话还是很快的   最好的教程就是github ...

感谢感谢，linux还行，就是差人指导

qq531884069

loftydongshi 发表于 2021-3-26 13:48
感谢感谢，linux还行，就是差人指导

CMAQ的编译安装还有运行教程家园里数不胜数了吧哈哈    我就不班门弄斧了   有问题欢迎留言

累死的绿树

qq531884069 发表于 2021-3-26 10:41
这个问题其实很好解决   你下载一个CMAQv5.3.2   最新的里面更新到了2019年的OMI  替换即可

ok，我去下载看看，感谢

累死的绿树

qq531884069 发表于 2021-3-26 10:41
这个问题其实很好解决   你下载一个CMAQv5.3.2   最新的里面更新到了2019年的OMI  替换即可

CMAQ网站无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助），是不是对中国IP有限制？你那里有下载好的吗？

loftydongshi

qq531884069 发表于 2021-3-26 16:56
CMAQ的编译安装还有运行教程家园里数不胜数了吧哈哈    我就不班门弄斧了   有问题欢迎留言

楼主没别的，就是优秀！

qq531884069

累死的绿树 发表于 2021-3-26 17:51
CMAQ网站无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助），是不是对中国IP有限制？你 ...

CMAQ从5.3开始代码维护就在github上面了   到github下载即可
https://github.com/USEPA/CMAQ

xenos_an

qq531884069 发表于 2021-3-26 10:44
虽然我几个版本都用。。但是就国内而言我个人不建议你用MEGAN3.1，因为stress这块的数据都是当初Alex教授 ...

非常感谢，我先尝试解决。
另外不知道PMF有没有了解？

yycdyy

您好，我想请问下载编译bio_emiss的时候报错：
compilation aborted for bio_emiss.f90 (code 1)
make: *** [Makefile:20: bio_emiss.o] Error 1
Failed to build bio_emiss
这个应该怎么解决吗，看了家园的帖子，有说改Makefile的，改了之后也还是报错这个。求解答

yycdyy

接上，编译信息如下：

./make_util megan_bio_emiss
NETCDF_DIR = /mnt/hzc/model/megan/envForMegan/netcdf

ifort -g -c -I/mnt/hzc/model/megan/envForMegan/netcdf/include misc_definitions_modul
e.f90ifort -g -c -I/mnt/hzc/model/megan/envForMegan/netcdf/include constants_module.f90
ifort -g -c -I/mnt/hzc/model/megan/envForMegan/netcdf/include bio_types.f90
ifort -g -c -I/mnt/hzc/model/megan/envForMegan/netcdf/include area_mapper.f90
area_mapper.f90(185): remark #8291: Recommended relationship between field width 'W'
and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.   write(*,'('' sw corner @ ('',1p,g14.8,'','',g14.8,'')'')') lon(1,1),lat(1,1)
----------------------------------^
area_mapper.f90(185): remark #8291: Recommended relationship between field width 'W'
and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.   write(*,'('' sw corner @ ('',1p,g14.8,'','',g14.8,'')'')') lon(1,1),lat(1,1)
--------------------------------------------^
area_mapper.f90(186): remark #8291: Recommended relationship between field width 'W'
and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.   write(*,'('' se corner @ ('',1p,g14.8,'','',g14.8,'')'')') lon(ide,1),lat(ide,1)
----------------------------------^
area_mapper.f90(186): remark #8291: Recommended relationship between field width 'W'
and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.   write(*,'('' se corner @ ('',1p,g14.8,'','',g14.8,'')'')') lon(ide,1),lat(ide,1)
--------------------------------------------^
area_mapper.f90(187): remark #8291: Recommended relationship between field width 'W'
and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.   write(*,'('' ne corner @ ('',1p,g14.8,'','',g14.8,'')'')') lon(ide,jde),lat(ide,j
de)----------------------------------^
area_mapper.f90(187): remark #8291: Recommended relationship between field width 'W'
and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.   write(*,'('' ne corner @ ('',1p,g14.8,'','',g14.8,'')'')') lon(ide,jde),lat(ide,j
de)--------------------------------------------^
area_mapper.f90(188): remark #8291: Recommended relationship between field width 'W'
and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.   write(*,'('' nw corner @ ('',1p,g14.8,'','',g14.8,'')'')') lon(1,jde),lat(1,jde)
----------------------------------^
area_mapper.f90(188): remark #8291: Recommended relationship between field width 'W'
and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.   write(*,'('' nw corner @ ('',1p,g14.8,'','',g14.8,'')'')') lon(1,jde),lat(1,jde)
--------------------------------------------^
area_mapper.f90(1084): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') model_x(m),model_y(m)
-------------------------^
area_mapper.f90(1084): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') model_x(m),model_y(m)
-----------------------------------^
area_mapper.f90(1089): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') mvtx_lon(m),mvtx_lat(m)
-------------------------^
area_mapper.f90(1089): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') mvtx_lon(m),mvtx_lat(m)
-----------------------------------^
area_mapper.f90(1094): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') x(m),y(m)
-------------------------^
area_mapper.f90(1094): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') x(m),y(m)
-----------------------------------^
area_mapper.f90(1099): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') data_lon(m),data_lat(m)
-------------------------^
area_mapper.f90(1099): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') data_lon(m),data_lat(m)
-----------------------------------^
area_mapper.f90(1109): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') vtx_x(m),vtx_y(m)
-------------------------^
area_mapper.f90(1109): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.         write(*,'('' ('',g14.8,'','',g14.8,'')'')') vtx_x(m),vtx_y(m)
-----------------------------------^
area_mapper.f90(1142): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') model_x(m),model_y(m)
---------------------------^
area_mapper.f90(1142): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') model_x(m),model_y(m)
-------------------------------------^
area_mapper.f90(1147): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') mvtx_lon(m),mvtx_lat(m)
---------------------------^
area_mapper.f90(1147): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') mvtx_lon(m),mvtx_lat(m)
-------------------------------------^
area_mapper.f90(1152): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') x(m),y(m)
---------------------------^
area_mapper.f90(1152): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') x(m),y(m)
-------------------------------------^
area_mapper.f90(1157): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') data_lon(m),data_lat(m)
---------------------------^
area_mapper.f90(1157): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') data_lon(m),data_lat(m)
-------------------------------------^
area_mapper.f90(1167): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') vtx_x(m),vtx_y(m)
---------------------------^
area_mapper.f90(1167): remark #8291: Recommended relationship between field width 'W
' and the number of fractional digits 'D' in this edit descriptor is 'W>=D+7'.           write(*,'('' ('',g14.8,'','',g14.8,'')'')') vtx_x(m),vtx_y(m)
-------------------------------------^
ifort -g -c -I/mnt/hzc/model/megan/envForMegan/netcdf/include bio_emiss.f90
bio_emiss.f90(1885): error #5082: Syntax error, found ',' when expecting one of: ( :
% [ . = =>      ll_id,ilon-11,ilat-9,XLAT(ilon,ilat),XLONG(ilon,ilat), &
-----------^
bio_emiss.f90(1894): error #5082: Syntax error, found END-OF-STATEMENT when expectin
g one of: :: ) , :                                        MSEBIO_STRESS(ilon,ilat),MSEBIO_OTHER(ilon,i
lat)                            ------------------------------------------------------------------------------------
--------------------------------^bio_emiss.f90(1934): error #5082: Syntax error, found IDENTIFIER 'AT' when expecting
one of: :: ) ( , : * . % + - [ . = ** / // .LT. < .LE. <= .EQ. == ...      write (99, '(3(I,","), 46(F10.4,","))')cell_id,ilon-11,ilat-9,xlai(ilon,il   
at,:)------------------------------------------------------------------------------------
-^bio_emiss.f90(1884): error #6404: This name does not have a type, and must have an e
xplicit type.   [CE]      write (99, '(3(I5,","), 2(F10.4,","),2(I,","),19(F10.4,","),1(F10.4))')ce
-----------------------------------------------------------------------------^
bio_emiss.f90(1934): error #6404: This name does not have a type, and must have an e
xplicit type.   [IL]      write (99, '(3(I,","), 46(F10.4,","))')cell_id,ilon-11,ilat-9,xlai(ilon,il   
at,:)------------------------------------------------------------------------------^
bio_emiss.f90(1934): error #6404: This name does not have a type, and must have an e
xplicit type.   [AT]      write (99, '(3(I,","), 46(F10.4,","))')cell_id,ilon-11,ilat-9,xlai(ilon,il   
at,:)------------------------------------------------------------------------------------
-^compilation aborted for bio_emiss.f90 (code 1)
make: *** [Makefile:20: bio_emiss.o] Error 1
Failed to build bio_emiss

zhangyw

WRF模拟国内各个地区，风速总是偏大很多，不知道是不是下垫面的问题，如果是，应该怎么改呢？求指导！