求助关于cctm的收敛问题

HarryChang

经过好久的努力，终于cctm开始计算了，然而，在计算的过程中还是报错了，报错信息为：
ERROR: Stopping because of EBI convergence failures
欧拉反迭代收敛失败
不太明白是怎么回事儿，这种问题有人遇到过吗，需要怎样解决，谢谢。

HarryChang

计算部分的日志如下：
Calculating emissions point source layer fractions for 000300
     >>> Processing Saturday June 15, 2019
         at time 0:03:00
     Value for KZMIN:  N returning FALSE

     This run uses Kz0UT, *NOT* KZMIN in subroutine edyintb.



     Maximum eddy diffusivity of:   322.25     (m**2/sec)
     at col, row, layer:  43,  17,   4
     corresponding to a free tropospheric wind shear of:  2.98389E-03 (/sec),
                       a bulk Richardson Number of:  -59.018    ,
     and pot. temps. in layer and layer+1:   301.57       299.07


WARNING: output file already exists: FLOOR_000.zjk_real_20190615
FLOOR_000.zjk_real_20190615


     File "FLOOR_000.zjk_real_20190615" opened for output on unit:  95
     FLOOR_000.zjk_real_20190615

after     VDIFF G  1.0507157E+00 A  1.4657550E+10 N  4.1565890E-03
after COUPLE_WR G  9.4741289E+03 A  1.3028908E+14 N  3.7483597E+01

     No BC's in file BNDY_GASC_1 for the following adv species: Set to 1.00E-30
          ALDX
          MEOH
          CAT1
          CRON
          CRPX
          OPAN
          ISPD
          TERP
          ETOH
          CL2
          HOCL
          FMCL
          HCL
          SESQ

     No BC's in file BNDY_AERO_1 for the following adv species: Set to 1.00E-30
          ANH4J
          ANH4I
          ANO3J
          ANO3I
          AORGCJ
          APOCJ
          APOCI
          APNCOMJ
          APNCOMI
          AECJ
          AECI
          AOTHRJ
          AOTHRI
AFEJ
          AALJ
          ASIJ
          ATIJ
          ACAJ
          AMGJ
          AKJ
          AMNJ
          AH2OJ
          AH2OI
          ANAJ
          ANAI
          ACLJ
          ACLI
          ASEACAT
          ACLK
          ASO4K
          ANH4K
          ANO3K
          AH2OK

     Reading U-windfield from MET_DOT_3D for variable: UWINDC

     Reading V-windfield from MET_DOT_3D for variable: VWINDC

     layer    S (X3FACE_GD) Delta S
       1     0.0070000     0.0070000     1.0000000
       2     0.0170000     0.0100000     1.0000000
       3     0.0300000     0.0130000     1.0000000
       4     0.0460000     0.0160000     1.0000000
       5     0.0660000     0.0200000     1.0000000
       6     0.0910000     0.0250000     1.0000000
       7     0.1200000     0.0290000     1.0000000
       8     0.1624703     0.0424703     1.0000000
       9     0.2049406     0.0424702     1.0000000
      10     0.2474109     0.0424703     1.0000000
      11     0.2898812     0.0424703     1.0000000
      12     0.3647836     0.0749025     1.0000000
      13     0.4335882     0.0688046     1.0000000
      14     0.4967090     0.0631208     1.0000000
      15     0.5545365     0.0578275     1.0000000
      16     0.6074384     0.0529019     1.0000000
      17     0.6557611     0.0483227     1.0000000
      18     0.6998305     0.0440694     1.0000000
19     0.7399533     0.0401227     1.0000000
      20     0.7764176     0.0364643     1.0000000
      21     0.8094943     0.0330767     1.0000000
      22     0.8394376     0.0299433     1.0000000
      23     0.8664861     0.0270485     1.0000000
      24     0.8908622     0.0243761     1.0000000
      25     0.9120463     0.0211840     1.0000000
      26     0.9304562     0.0184099     1.0000000
      27     0.9464553     0.0159991     1.0000000
      28     0.9603593     0.0139040     1.0000000
      29     0.9724425     0.0120832     1.0000000
      30     0.9829434     0.0105009     1.0000000
      31     0.9920692     0.0091258     1.0000000
      32     1.0000000     0.0079308     1.0000000

after       ADV G  9.4656973E+03 A  1.3018348E+14 N  3.7453659E+01

     H-eddy DT & integration steps:   3.6000000E+02       1
after     HDIFF G  9.4656973E+03 A  1.3018348E+14 N  3.7453659E+01
after DECOUPLE_ G  1.0502857E+00 A  1.4640354E+10 N  4.1556996E-03
     Value for CLD_DIAG:  N returning FALSE
     CTM_WET_DEP_1   :/usr/local/CMAQ-5.0.2/data/wrf36/cctm/CCTM_D502a_Linux3_x86_64gcc.WETDEP1.zjk_real_20190615

     >>--->> WARNING in subroutine OPEN3
     File not available.

     Could not open CTM_WET_DEP_1    file for update - try to open new

     "CTM_WET_DEP_1" opened as NEW(READ-WRITE )
     File name "/usr/local/CMAQ-5.0.2/data/wrf36/cctm/CCTM_D502a_Linux3_x86_64gcc.WETDEP1.zjk_real_20190615"
     File type GRDDED3
     Execution ID "CCTM_D502a_Linux3_x86_64gcc"
     Grid name "zjk_grid_real"
     Dimensions: 61 rows, 46 cols, 1 lays, 80 vbles
     NetCDF ID:   1245184  opened as VOLATILE READWRITE
     Starting date and time  2019166:010000 (1:00:00   June 15, 2019)
     Timestep                          010000 (1:00:00 hh:mm:ss)
     Maximum current record number         0

     >>--->> WARNING in subroutine RESCLD on PE 000
     Parameter QI (cloud ice) was not found on file MET_CRO_3D
     M3WARN:  DTBUF 0:00:00   June 15, 2019 (2019166:000000)

          YOU SHOULD VERIFY that the cloud microphysics scheme used
in the Meteorological Model did not include ice/snow.  If
          it did, then you need to reprocess the meteorological data
          through MCIP and pass QI to file  MET_CRO_3D  to avoid
          errors in the wet deposition.

     Processing will continue with QI set to ZERO.  <<---<<


     >>--->> WARNING in subroutine RESCLD on PE 000
     Parameter QS (snow) was not found on file MET_CRO_3D
     M3WARN:  DTBUF 0:00:00   June 15, 2019 (2019166:000000)

          YOU SHOULD VERIFY that the cloud microphysics scheme used
          in the Meteorological Model did not include ice/snow.  If
          it did, then you need to reprocess the meteorological data
          through MCIP and pass QS to file  MET_CRO_3D  to avoid
          errors in the wet deposition.

     Processing will continue with QS set to ZERO.  <<--<<

after   CLDPROC G  1.0502857E+00 A  1.4640354E+10 N  4.1556996E-03

          Euler Backward Iterative Parameters -
          Chemistry Integration Time Interval (min):      6.0000
          EBI maximum time step (min):                    2.5000


          Species convergence tolerances:
          NO2                   1.00E-03
          NO                    1.00E-03
          O                     1.00E+00
          O3                    1.00E-03
          NO3                   1.00E-03
          O1D                   1.00E+00
          OH                    1.00E-03
          HO2                   1.00E-03
          N2O5                  1.00E-03
          HNO3                  1.00E-03
          HONO                  1.00E-03
          PNA                   1.00E-03
          H2O2                  1.00E-03
          XO2                   1.00E-03
          XO2N                  1.00E-03
          NTR                   1.00E+00
ROOH                  1.00E-03
          FORM                  1.00E-03
          ALD2                  1.00E-03
          ALDX                  1.00E-03
          PAR                   1.00E-03
          CO                    1.00E-03
          MEO2                  1.00E-03
          MEPX                  1.00E-03
          MEOH                  1.00E-03
          HCO3                  1.00E+00
          FACD                  1.00E-03
          C2O3                  1.00E-03
          PAN                   1.00E-03
          PACD                  1.00E-03
          AACD                  1.00E-03
          CXO3                  1.00E-03
          PANX                  1.00E-03
          ROR                   1.00E-03
          OLE                   1.00E-03
          ETH                   1.00E-03
          IOLE                  1.00E-03
          TOL                   1.00E-03
          CRES                  1.00E-03
          TO2                   1.00E-03
          TOLRO2                1.00E-03
          OPEN                  1.00E-03
          MGLY                  1.00E-03
          CRO                   1.00E-03
          CAT1                  1.00E-03
          CRON                  1.00E-03
          CRNO                  1.00E-03
          CRN2                  1.00E-03
          CRPX                  1.00E-03
          OPO3                  1.00E-03
          CAO2                  1.00E-03
          OPAN                  1.00E-03
          XYL                   1.00E-03
          XYLRO2                1.00E-03
          ISOP                  1.00E-03
          ISPD                  1.00E-03
          ISOPRXN               1.00E+00
          TERP                  1.00E-03
          TRPRXN                1.00E+00
          SO2                   1.00E-03
SULF                  1.00E+00
          SULRXN                1.00E+00
          ETOH                  1.00E-03
          ETHA                  1.00E-03
          CL2                   1.00E-03
          CL                    1.00E-03
          HOCL                  1.00E-03
          CLO                   1.00E-03
          FMCL                  1.00E-03
          HCL                   1.00E-03
          TOLNRXN               1.00E+00
          TOLHRXN               1.00E+00
          XYLNRXN               1.00E+00
          XYLHRXN               1.00E+00
          BENZENE               1.00E-03
          BENZRO2               1.00E-03
          BNZNRXN               1.00E+00
          BNZHRXN               1.00E+00
          SESQ                  1.00E-03
          SESQRXN               1.00E+00
     Value for CTM_PHOTDIAG:  N returning FALSE

     File "CSQY_DATA" opened for input on unit:  94
     /usr/local/CMAQ-5.0.2/scripts/cctm/BLD_D502a/CSQY_DATA_cb05tucl_ae6_aq

Sucessfully Loaded JTABLE
     PHOT: Identified USGS24 land use scheme for surface albedo used by inline photolysis calculation.

     File "OMI" opened for input on unit:  93
     /usr/local/CMAQ-5.0.2/data/wrf36/raw/phot/OMI.dat


     >>--->> WARNING in subroutine O3TOTCOL on PE 000
     Requested date is beyond available data on OMI file:  <0:00:00   July 28, 2013
     M3WARN:  DTBUF 0:00:00   June 15, 2019 (2019166:000000)
     Total column ozone will be estimated from the corresponding Julian Day
     of the last available year on the OMI input file:0:00:00   June 15, 2013 <<---<<

application called MPI_Abort(MPI_COMM_WORLD, 0) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=0
:
system msg for write_line failure : Bad file descriptor
WARNING: EBI Euler convergence failure
         Reducing EBI time step because of MAXPRED convergence failure for
         Cell ( 14,   9,   1) and species CO  Back-up number 1
WARNING: EBI Euler convergence failure
         Reducing EBI time step because of MAXPRED convergence failure for
         Cell ( 14,   9,   1) and species CO  Back-up number 2
WARNING: EBI Euler convergence failure
         Reducing EBI time step because of MAXPRED convergence failure for
         Cell ( 14,   9,   1) and species CO  Back-up number 3
WARNING: EBI Euler convergence failure
         Reducing EBI time step because of MAXPRED convergence failure for
         Cell ( 14,   9,   1) and species CO  Back-up number 4
WARNING: EBI Euler convergence failure
         Reducing EBI time step because of MAXPRED convergence failure for
         Cell ( 14,   9,   1) and species CO  Back-up number 5
ERROR: Max number of EBI time step reductions exceeded
      Convergence failure for cell ( 14,   9,   1)
      Convergence failure for the following species:
          NO2
          NO
          O3
          NO3
          OH
          HO2
          N2O5
          PNA
          XO2
          XO2N


     *** ERROR ABORT in subroutine HRSOLVER on PE 000
     ERROR: Stopping because of EBI convergence failures
     Date and time 0:00:00   June 15, 2019  (2019166:000000)
date
Wed Aug 14 17:08:04 CST 2019
exit

HarryChang

欧拉反向迭代的迭代时间限制和判断收敛的条件在哪里修改呢？请问有哪位知道吗？

HarryChang

HarryChang 发表于 2019-8-14 17:33
欧拉反向迭代的迭代时间限制和判断收敛的条件在哪里修改呢？请问有哪位知道吗？

找到判断收敛的初始条件位置（5.0.2）为CMAQ/models/CCTM/gas/ebi_cb05tucl/hrinit.F

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c  Set species tolerances
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      RTOL( NO2     ) = 1.0E-03
      RTOL( NO      ) = 1.0E-03
      RTOL( O       ) = 1.0E+00
      RTOL( O3      ) = 1.0E-03
      RTOL( NO3     ) = 1.0E-03
      RTOL( O1D     ) = 1.0E+00
      RTOL( OH      ) = 1.0E-03
      RTOL( HO2     ) = 1.0E-03
      RTOL( N2O5    ) = 1.0E-03
      RTOL( HNO3    ) = 1.0E-03
      RTOL( HONO    ) = 1.0E-03
................

找到日志中报错的species，将收敛条件修改到容易满足，再重新编译CCTM，或许可以解决这一问题。
但猜测本质上还是由于我的排放清单存在问题导致的报错（图省事只做了线分配，没有做面源）。