WRF、MEGAN以及CMAQ答疑帖

qq531884069

特纳波 发表于 2023-3-13 16:38
大佬你好，我安装cmaq的时候遇到了一些问题，在编译CCTM的时候有这样报错/home/changfy/WRF-CMAQ/CMAQ-5.2/ ...

检查你netcdf的编译   去按照github上CMAQ官方的安装方式逐步来就绝对不会有错

xiuyuanyang

请问，MEGAN3.2，运行物中转换>./run.mgn2mech.v3.csh，报错如下，该怎么解决？PS:预处理程序和排放因子计算以及物种转换之前步骤都跑了一遍，运行正常，生成对应文件
--------------------------错误---------------------------
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source            
mgn2mech           00000000004F96CA  Unknown               Unknown  Unknown
libpthread-2.17.s  00002B16EBD34630  Unknown               Unknown  Unknown
mgn2mech           000000000040858C  MAIN__                    685  mgn2mech.F
mgn2mech           0000000000404EA2  Unknown               Unknown  Unknown
libc-2.17.so       00002B16EC167555  __libc_start_main     Unknown  Unknown
mgn2mech           0000000000404DA9  Unknown               Unknown  Unknown

qq531884069

xiuyuanyang 发表于 2023-3-16 21:07
请问，MEGAN3.2，运行物中转换>./run.mgn2mech.v3.csh，报错如下，该怎么解决？PS:预处理程序和排放因子计 ...

仔细去看你的csh文件  你肯定是文件路径写错或者要求的文件写成了其他文件导致程序直接crash

xiuyuanyang

qq531884069 发表于 2023-3-17 01:08
仔细去看你的csh文件  你肯定是文件路径写错或者要求的文件写成了其他文件导致程序直接crash

嗯，找到原因了，是源代码mgn2mech.F第685行outer数组越界了，源代码里是按照开发人员测试的数组写的，改成数组大小就行了，感觉本身这句代码没啥用，注释掉也可以
-----------------------------------源代码----------------------------------
WRITE(MESG,*) 'conv_fac,outer',conv_fac(1),outer(1,101,101)
-----------------------------------改为--------------------------------------
WRITE(MESG,*) 'conv_fac,outer',conv_fac(1),
     &                   outer(1,size(outer, dim=2),size(outer, dim=3))

特纳波

qq531884069 发表于 2023-3-16 06:29
检查你netcdf的编译   去按照github上CMAQ官方的安装方式逐步来就绝对不会有错

重装了netcdf就好了 谢谢

Xin-Jiayi

qq531884069 发表于 2023-3-12 10:02
自己写个循环就行  你应该已经解决了吧   forum上我看Christin和我都回复你了

已解决，感谢大佬

xiuyuanyang

大佬，开启ISAM模块编译CCTM，报错如下，怎么解决呀
---------------------------------------错误---------------------------------------
mpifort -c -ffixed-form -ffixed-line-length-132 -funroll-loops -finit-character=32 -O3 -funroll-loops -finit-character=32 -Wtabs -Wsurprising -ftree-vectorize -ftree-loop-if-convert -finline-limit=512 -I /home/atmos/MODEL/CMAQv54/lib/x86_64/gcc/ioapi/lib -I /home/atmos/MODEL/CMAQv54/lib/x86_64/gcc/ioapi/include_files -I /home/atmos/MODEL/CMAQv54/lib/x86_64/gcc/mpi -I. -Disam -Dparallel -Dsens -Dm3dry_opt -DSUBST_BARRIER=SE_BARRIER -DSUBST_GLOBAL_MAX=SE_GLOBAL_MAX -DSUBST_GLOBAL_MIN=SE_GLOBAL_MIN -DSUBST_GLOBAL_MIN_DATA=SE_GLOBAL_MIN_DATA -DSUBST_GLOBAL_TO_LOCAL_COORD=SE_GLOBAL_TO_LOCAL_COORD -DSUBST_GLOBAL_SUM=SE_GLOBAL_SUM -DSUBST_GLOBAL_LOGICAL=SE_GLOBAL_LOGICAL -DSUBST_GLOBAL_GATHER=SE_GLOBAL_GATHER -DSUBST_GLOBAL_BCAST=SE_GLOBAL_BCAST -DSUBST_LOOP_INDEX=SE_LOOP_INDEX -DSUBST_SUBGRID_INDEX=SE_SUBGRID_INDEX -DSUBST_HI_LO_BND_PE=SE_HI_LO_BND_PE -DSUBST_SUM_CHK=SE_SUM_CHK -DSUBST_SE_INIT=SE_INIT -DSUBST_INIT_ARRAY=SE_INIT_ARRAY -DSUBST_COMM=SE_COMM -DSUBST_MY_REGION=SE_MY_REGION -DSUBST_SLICE=SE_SLICE -DSUBST_GATHER=SE_GATHER -DSUBST_DATA_COPY=SE_DATA_COPY -DSUBST_IN_SYN=SE_IN_SYN -DSUBST_PE_COMM=\"./PE_COMM.EXT\" -DSUBST_CONST=\"./CONST.EXT\" -DSUBST_FILES_ID=\"./FILES_CTM.EXT\" -DSUBST_EMISPRM=\"./EMISPRM.EXT\" -DSUBST_MPI=\"/home/atmos/MODEL/CMAQv54/lib/x86_64/gcc/mpi/include/mpif.h\" wr_sa.F
vppm.F:341.11:

      DO NP = 1, NPMAX                                                         
           1
错误： (1)处的符号‘np’没有隐式类型
vppm.F:341.22:

      DO NP = 1, NPMAX                                                         
                      1
错误： (1)处的符号‘npmax’没有隐式类型
vppm.F:169.46:

      REAL,    INTENT( INOUT ) :: SEN( :,:,: )                                 
                                              1
错误： (1)处的符号不是一个 DUMMY 变量
make: *** [vppm.o] 错误 1
make: *** 正在等待未完成的任务....
**ERROR** while running make command

if ( ! 0 ) then
mv Makefile Makefile.gcc
if ( -e Makefile.gcc && -e Makefile ) rm Makefile
ln -s Makefile.gcc Makefile
endif
if ( 0 != 0 ) then
if ( -e /home/atmos/MODEL/CMAQv54/CCTM/scripts/BLD_CCTM_v54_ISAM_gcc/CCTM_v54_ISAM.cfg ) then
mv CCTM_v54_ISAM.cfg.bld /home/atmos/MODEL/CMAQv54/CCTM/scripts/BLD_CCTM_v54_ISAM_gcc/CCTM_v54_ISAM.cfg
if ( 0 ) then
exit

特纳波

大佬你好，我在运行mcip的时候遇到一下问题，**********************************************************************
*** SUBROUTINE: RDWRFEM
***   ERROR RETRIEVING VARIABLE FROM WRF FILE
***   VARIABLE = LANDUSEF
***   RCODE = NetCDF: Start+count exceeds dimension bound
**********************************************************************
能解答一下吗

qq531884069

特纳波 发表于 2023-4-3 09:56
大佬你好，我在运行mcip的时候遇到一下问题，********************************************************** ...

用ncdump -h检查你mcip读取的文件是否网格全部一致

Landsay

楼主你好，我想请问一下，我在两台超算上安装了WRF-Chem，用MEGAN制作排放源，得到的wrfbiochemi数据是一样的，用的namelist.input也是一样的，同一区域同一时期模拟结果中的臭氧浓度也是基本一样的， 但是为什么一个里MEBIO_ISOP变量的值是正常逐小时变化的，另一个超算上wrfout里的MEBIO_ISOP变量值都是0呢？麻烦楼主了，谢谢