/------------------------------------------------------------------------------------
/ Run : simpleCO
/ Purpose: case with CO chemistry for benchmark testing box model
/------------------------------------------------------------------------------------
/--flags--
DEPO  0  / deposition 
NOEM  0  / NO emission from tropical rain forest (Yienger/Levy)
/
/------------------------------------------------------------------------------------
NBOX 1 / fix number of boxes (use '1')
SCAP 0 / thermodynamic icalc for inorg aerosol (0=  no calc USE THIS; 1=calcul thermo)
TSTR 0 / start time (s): 21600=6am; 32400=9am - at midnight
TSTP 86400 / 172800 /259200 / stop time (s): 72000=8pm 
TLEN 60 / requested timestep size in seconds 
SKIP 15 /15 / # of timesteps per binary output point
RTOL 1.0E-2 / relative tolerance (solver convergence parameter)
ATOL 1.0E+2 / absolute tolerance (solver convergence parameter)
PSSF 1 / do we print the final timestep concentrations?
/-------------------------------------------
/ PHOTOLYSIS: SZAF = fixed; PPHO = varying
/-------------------------------------------
/SZAF 45.
PPHO  -3.0  0  0.0  2014  3  13
/------------------------------------------------------------------------------------
/ METEO PARAMS
/------------------------------------------------------------------------------------
/ The variation de l'altitude de la boite 1
/ est donne par le mot cle HBOX. 
/ (search for 'height' in spforcage3 and boxmod)
HTOP  2.0E5 / BOX HEIGHT (cm). HTOP fixes altitude of top of box 2. 
HBOX 2
0     1.0E5 
86400 1.0E5 
/
/ MIXING RATE (advective and/or diffusive) (s-1)
/ Mixes air in Box 1 with background air in horizontal or vertical direction.
/MIX 2
/0 3.06534e-05
/86400 3.06534e-05
/
AERO 0.0E-5 /Not used here : aerosol surface area for heterogeneous reactions (cm2/cm3)
/Invariant number density of non-volatile, unreactive aerosol (molec/cm3)
/Alma used = 1ug/m3 at 200 g/mol : 1e-6 / 200 * 1e-6 * 6.023* 1e23 -> 3.1e9 molecules/cc
SEED 3.0E+10
TEMP / temperature (temp): les parametres sont a fixer pour chaque couche
 284.9 0.    00000.
/ "saison" : cette variable a pour unique but de changer la vitesse
/ des resistances pour le depot sec. Deux "saisons" sont tabulees :
/ SEAS =1 : automne avance, apres le gel mais pas de neige au
/           sol (correspond a la categorie 3 du papier de Wesely, 89)
/ SEAS =2 : milieu de l'ete avec une vegetation luxuriante (correspond
/           a la categorie 1 du papier de Wesely, 89)
SEAS  2
/
/ humidite relative (rh): les parametres sont a fixer pour chaque couche
/ (autant de ligne que de niveau). rh est parametre comme une fonction
/ sinusoidale du temps. Le premier parametre donne rh moyen (en %),
/ le second donne l'amplitude, le troisieme le temps pour lequel le
/ maximum de rh est atteint (en s)
/
RHUM
  75.0 00.0 00000.
WATR  6.40E+17 / water content (overrides RHUM!)
NDEN  2.57E+19 / number density of atmosphere
/
/ lecture des proportions des differents environnements types pour la
/ definition du scenario. 4 environements sont possibles dans la structure
/ actuelle du code : urbain (URB1), terre cultivee(CULT), foret de
/ feuillu (FLEA) et foret de conifere (FCON). Le nombre apres le mot
/ cle SURF donne le nombre de donnee tabule en temps. l'ordre des donnees
/ est ensuite : temps, URB1,CULT,FLEA,FCON. Aucune ligne de commentaire
/ n'est autorise entre les donnees. Le dernier temps doit etre plus long
/ que le temps d'arret de la simulation. La surface est calculee par
/ interpolation lineaire entre les points
SURF 3
0       0.00 0.00 1.00 0.00
21600   0.00 0.00 1.00 0.00
2592000 0.00 0.00 1.00 0.00
/
/wind speed : for deposition
WIND 4.0  0.0  0.0
/
REAC GO3  9.80E+11  0.00E+12
REAC GCO  2.45E+12  0.00E+12
REAC GCH4 4.50E+13  0.00E+13
REAC GNO  2.45E+09  0.00E+09
/
CONS BOT 4
NOx 2
0.     1.00E+11
86400. 1.00E+11
GO3 2
0.      9.80E+11
86400.  9.80E+11
GC02000 2
0.     1.00E+11
86400. 1.00E+11
GC50001 2
0.     1.00E+12
86400. 1.00E+12
PREV 0
END