[RegCNET] chem_icbc problem in CBMZ run
Athanasios Tsikerdekis
atsikerdekis at gmail.com
Wed Oct 29 12:08:19 CET 2014
Greetings RegCMers,
I would like to run an experiment with the gas chemistry option (CBMZ).
After fixing some minors problems I had I am left with a problem in the
chemical boundaries conditions. Everything else in the input files seems
to be created normally (terrain, sst, icbc, emcre_grid, inter_emissions
(using IIASA)) but in the last step where I run the script of chem_icbc
something is not created correctly because when I check the netcdf file
in the input folder (*CHBC*.nc) I can see that it has zero timesteps.
> ncdump -h input/CBMZrun_CHBC.2003090100.nc
> netcdf CBMZrun_CHBC.2003090100 {
> dimensions:
> jx = 136 ;
> iy = 124 ;
> kz = 18 ;
> * time = UNLIMITED ; // (0 currently)*
And when I try to run the model I get the following error as expected.
>
> HVREAD: TUVGRID
> READ CHEM EMISSION for 2003090002
> Opening ch. emission file input/CBMZrun_CHEMISS.nc
> CHE_EMISS: Reading record 165
> * Time var in ICBC has zero dim.*
> -------------- FATAL CALLED ---------------
> Fatal in file: mod_che_ncio.F90 at line: 745
> ICBC READ
> -------------------------------------------
> mod_che_ncio.F90 : 745: 1
> Abort called by computing node 0 at 2014-10-29
> 12:53:02.572 +0200
> application called MPI_Abort(comm=0x84000000, 1) - process 0[cli_0]:
> aborting job:
> application called MPI_Abort(comm=0x84000000, 1) - process 0
> -------------- FATAL CALLED ---------------
> Fatal in file: mod_mppparam.F90 at line: 677
> mpi_recv error.
> -------------------------------------------
>
> ===================================================================================
> = BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
> = EXIT CODE: 1
> = CLEANING UP REMAINING PROCESSES
> = YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
> ===================================================================================
I have downloaded again the average mean mz4 global files from the
OXIGLOB folder ( mz4_avg_1999-2009_**.nc) without any luck. Am I missing
something? Does anyone encounter a similar problem? I am attaching the
*.in file in case anyone wants to experiment on this. Thank you very
much in advance!
Best wishes,
Thanos
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20141029/2aa39754/attachment.html>
-------------- next part --------------
! QUADIEEMS domain
&dimparam
iy = 124,
jx = 136,
kz = 18,
nsg = 1,
/
&geoparam
iproj = 'LAMCON',
ds = 50.0,
ptop = 5.0,
clat = 47,
clon = 5,
plat = 45.39,
plon = 13.48,
truelatl = 30.0,
truelath = 60.0,
i_band = 0,
/
&terrainparam
domname = 'CBMZrun',
smthbdy = .false.,
lakedpth = .false.,
fudge_lnd = .false.,
fudge_lnd_s = .false.,
fudge_tex = .false.,
fudge_tex_s = .false.,
ltexture = .true.,
dirter = 'input',
inpter = 'DATASETS',
/
&ioparam
ibyte = 4,
/
&debugparam
debug_level = 1,
dbgfrq = 3,
/
&boundaryparam
nspgx = 12,
nspgd = 12,
/
&globdatparam
ibdyfrq = 6,
ssttyp = 'ERSST',
dattyp = 'EIN15',
gdate1 = 2003090100,
gdate2 = 2003090400,
dirglob = 'input',
inpglob = 'DATASETS',
/
&globwindow
lat0 = 0.0
lat1 = 0.0
lon0 = 0.0
lon1 = 0.0
/
&restartparam
ifrest = .false.,
mdate0 = 2003090200,
mdate1 = 2003090200,
mdate2 = 2003090300,
/
&timeparam
dtrad = 30.,
dtabem = 18.,
dtsrf = 600.,
dt = 120.,
/
&outparam
ifsave = .true. ,
savfrq = 7200.,
ifatm = .true. ,
atmfrq = 6.,
ifrad = .true. ,
radfrq = 6.,
ifsrf = .true. ,
ifsub = .true. ,
srffrq = 3.,
ifchem = .true.,
ifopt = .true.,
chemfrq = 6.,
dirout='output'
/
&physicsparam
iboudy = 5,
ibltyp = 1,
icup_lnd = 2,
icup_ocn = 4,
igcc = 1,
ipptls = 1,
iocnflx = 2,
ipgf = 0,
iemiss = 1,
lakemod = 0,
ichem = 1,
scenario = 'A1B',
idcsst = 0,
iseaice = 0,
idesseas = 1,
iconvlwp = 1,
/
&subexparam
qck1land = .250E-03,
qck1oce = .250E-03,
cevaplnd = .100E-02,
cevapoce = .100E-02,
caccrlnd = 3.000,
caccroce = 3.000,
cftotmax = 0.75,
/
&grellparam
/
&emanparam
/
&chemparam
chemsimtype = 'CBMZ', ! Which chemical tracers to be activated.
! One in :
! DUST : Activate 4 dust bins scheme
! SSLT : Activate 2 bins Sea salt scheme
! CARB : Activate 4 species black/anthropic
! carbon simulations
! SULF : Activate SO2 and SO4 tracers
! SUCA : Activate both SUKF and CARB
! AERO : Activate all DUST, SSLT, CARB and SULF
! CBMZ : Activate gas phase and sulfate
ichsolver = 1, ! Activate the chemical solver
idirect = 1, ! enable or not aerosol feedbacks on radiation and
! dynamics (aerosol direct and semi direct effcts):
! 0 = no coupling. Aerosol are only transported and
! don't interact with radiation scheme.
! 1 = no coupling to dynamic and thermodynamic. However
! the clear sky surface and top of atmosphere
! aerosol radiative forcings are diagnosed.
! 2 = allows aerosol feedbacks on radiative,
! thermodynamic and dynamic fields.
ichdrdepo = 1, ! 1 = enable tracer surface dry deposition. For dust,
! it is calculated by a size settling and dry
! deposition scheme. For other aerosol,a dry
! deposition velocity is simply prescribed further.
! Next release will include an improved aerosol dry
! deposition scheme for non dust aerosols.
ichcumtra = 1, ! 1 = enable tracer convective transport and mixing.
ichsursrc=1, ! Enable the emissions fluxes.
ichremlsc=1, ! 1 = enable tracer rainout
ichdustemd=1, ! Choice for parametrisation of dust emission size distribution
! 1 = use the standard scheme (Alfaro et al., Zakey et al.)
! 2 = use a revised soil granulometry ref +
! Kok et al emission size distribution : Menut et al.,2012; + Kok et al., 2011
rdstemfac = 1, ! Dust emission adjustment factor ( soil erodibility) : linearly reduce or increase the dust flux
/
&clmparam
dirclm = '/set/this/to/where/your/input/clm/data/are',
imask = 1,
clmfrq = 12.,
/
More information about the RegCNET
mailing list