******************************************************************************
CAABA
******************************************************************************

- split CHANGELOG into caaba and mecca
- test KHET in caaba boxmodel
- test montecarlo with DifferentMaxTime = false
- define aqueous-phase type "ocean" "aerosol" or "rain"
  (then use if_ocean instead of if_a03)
- problem: meccalit.bib doesn't have the refs for mim1
- combine init_j and init_j_index into one struct (also for init_spec?)
- rewrite a subroutine that only prints the rate coefficients (get
  rconst(:) out of kp4)
- jval and sappho output: loop over all species (skip PHOTON)
- caaba_io_*.inc should not have STOP statements
- convert $C2/nc_parser.bash to tcsh (?) and integrate into xcaaba (?)
- define error subroutine also for caaba
- in SUBROUTINE open_output_file in caaba_io_netcdf.inc, check that
  size of specid(10000) is at least as big as MAX_EQN (max_spc?)
- add scenario "MYTEST" that users can change
- use ncap2 instead of ncclamp in multirun.tcsh:
  #ncclamp $murun/$ncfile time 0 $basemurun gt
  ncap2 -O -s 'time > 0 ? "$basemurun" : 0000' $murun/$ncfile $murun/$ncfile
  see ~/messy/mecca/delange-telis/caaba_3.0_arno.zip
- is seed for mc okay? (PJ: use status of previous run instead)

******************************************************************************
MECCA
******************************************************************************

- rename testfac -> fast_eq, then mention in meccanism.tex and cite ref2848
- write script with more checks to test if *.eqn syntax is okay, based
  on check_rxn_numbers.py
- look at Annas bcary rxn mech, test big mcm output with ferret (>2000
  variables)
- in aqueous.eqn, the ion balance is disturbed because of rxn <A9605>
  even though the charge balance is okay in that reaction. Why?
- better error message if reaction number is used twice (currently
  only nc error message)
- put mecca info also into the *.nc files (netcdf attributes: "mecca_wanted,
  mecca_decomp, ...)
- good ideas for mecca labels: L1...L9 complexity level
- add MAM, MBM and MOM from Domenico (alpha- and beta- pinene mechanisms)
- use ref2559 (Pratte and Rossi, 2006) for alpha(HOCl)?
- iodine:
  - add diurnal cycle to CH3I emissions? (according to results from Doug
    Wallace at HITT in Heidelberg)
- KPP: add f90 comment with rxn number at end of each line in update_rconst
- add link test (similar to $E2/messy/util/mecca_link.sh) to xcaaba
- for boxmodel, user can either supply "runtime" or "model_end"
- add nml-switch l_check_range (as in MECCA1)
- LaTeX:
  - meccanism.tex: setting N=1 in k_3rd_iupac produces k_3rd (check this first!)
- messy_mecca_box.f90-*: search for "value from MOCCA" and insert
  correct reference
- /modeldata/.../mecca: change unit in t63 file from kg/kg to mol/mol
- pKa(HOCl) and pKa(HOBr) from 1802 and 1803 or from [1039]? (no HOBr T-dep)
- MECCA/KPP problem ERR=0 if initial values of species = 0 (ask Andrea)
- search for better KH(HOBr)
- pKa(CH3COOH) from ref1790 for EQ42 (ignore T-dep)
- add BrCl^- to eqn file according to rates from [1707]
- cleanup MECCA structure (after update of mecca_mbl_cpl_sulfvi_before):
  zmr, zmrbc, zmrac: warum 3 verschiedene?
- update_physc: define tcorr=1./temp-1./T0 and use it in mecca (henry_func)
- PJ: use configure for boxmodel
- mv CALL get_tracer(status, GPTRSTR, 'H2O',idx=zidt_H2O) etc to init_coupling
- replace REQ_AEROSOL by l_aero and delete l_aero from CTRL
- add DMS+O3 in aq phase? (suggestion from RvG in Nice talk)
- add MSA- + OH -> S(VI)? (suggestion from RvG in Nice talk) 
- is "C(:)=0" at correct place, or move to mr2conc?
- (vmr)y2conc, molmr2conc in messy_main_tools.f90 massmr2conc
- #G3101 (N2+O1D) *79/78 to account for Ar? ask J. Crowley or G. Moortgat
- A9300, A9303: problem with negative NO2 as product?
- what about O3 in #G9401: DMS+NO3 rxn? Is it consumed?
- add homogeneous reaction of N2O5?
  <G32nn>  N2O5 + H2O = 2 HNO3 : {%TrG} 2.5E-22+C(ind_H2O)*1.8E-39; {RvK40}
- MECCA_KHET:
  - in MECCA: N2O5 + aerosol produces DUMMY or 0.3 O3
  - gamma hno3 (for KHET) = ?
  - gamma n2o5 (for KHET) = 0.02 ref2238?
  - gamma n2o5 (for KHET) = <0.0016...0.017 ref1797?
- graphviz:
  - debug: what if "-" sign is before a product?
  - aq-phase: correct treatment of "##"
- try new rxns from ref2314? (however, see also email from 04 Jul 2008
  from Terry Dillon)
  INO3 thermal diss.:       k = 2.1E15 exp(-13670/T)
  I + INO3 -> I2 + NO3      k = 9.1E-11 exp(-146/T)
- aerosol "phase" (apn) or "mode" or "bin"?
- after renaming temp, press, cair to T, p, ca, work on KPP: do not
  produce temp variable
- what about O3 in #G9401: DMS+NO3 rxn?
- update Kw to misc/kw-ref2810.txt
- tdep for KH*KA for hno3 should be in Kh not in KA

******************************************************************************
