Commands to run this model:
1. Unzip the file: tar -xzvf iNRACM.zip
2. Go to the folder "simulate": cd simulate
3. Edit the file "date.input" for simulation of photolysis rates (see Note 1 for details)
4. Compile photo module: gfortran photoRACM.f
5. Simulate photolysis rates: ./a.out (then enter the input file "date.input")
6. Generate photolysis rates for 15NOy: python arrfactor.py (see Note 2 for details)
7. Edit the file "plume_input" for simulation of concentrations (see Note 3 for details)
8. Compile sbox module: gfortran 15Nsbox.f
9. Simulate concentrations: ./a.out (then enter the input file "plume_input")


Note 1: Variables in "date.input"
Cloud base (1-50): 20   #Do NOT change typically, the number is the layer order#
Cloud top (1-51): 30    #Do NOT change typically#
Date (yymmdd): 
Longitude (Real):      #Negative number for Western hemisphere
Latitude (Real):        #Negative number for Southern hemisphere
Starting time (GMT): 3.0     #GMT-Greenwich Mean Time
Delta time: 0.5             #Step length#
Stopping time (GMT): 21    

The starting time and stopping time are to constrain the start and ending time of daylight (but the real daytime length could still be a function of the latitude and date), which are set to be 3 and 21 for the Greenwich site, respectively. Therefore, when using the real longitude, MAKE SURE the starting and stopping time are converted from the local time to the GMT; while when using the longitude of zero, the starting and stopping time can be the default values of 3 and 21, respectively. To make things easy, for each run, you can just put the longitude as zero and stick to using the default values of 3 and 21 for the starting time and stopping time, respectively.   

Note 2: arrfactor.py
Isotope factor (showing the fractionation factor, Z can be considered as 15N)
ZO2: 1.0042  #Fractionation factor α for the NO2 photolysis#

Note 3: Variables in "plume_input"
START-TIME:   #The start time of the simulation#   
END-TIME:    
PRINT-INTERVALS:  #The same as the DETA TIME in Step 1

Initial concentrations ppb - the initial concentrations in ppb at the START-TIME
Inorganic Compounds-Z can be considered as 15N, 15NO concentration that can be derived from NO concentration by adopting the 15N/14N ratio in air or by measurements.
Organic Compounds-the concentrations of individual species can be derived from the total VOCs concentration by assuming the relative abundances or by measurements 

Emissions ppm min-emission data in the unit of ppm/min

Note 4: How to use results to determine the nitrate Δ17O

In the RESULT (PlumYYMMDD) spreadsheet, 
(1)	QNOX-the amount of NO2 produced from O3+NO reaction
(2)	RNOX-the amount of NO2 produced from NO+organic proxy radicals reactions
(3)	ANOX-the amount of nitrate produced from HO+NO2 reaction
(4)	BNOY-the amount of nitrate produced from NO3+organic carbons reactions
(5)	CNOY-the amount of nitrate produced from the N2O5 reaction: 

Finally, calculate the ratios QNOX/(QNOX+RNOX), ANOX/(ANOX+BNOY+CNOY), BNOY/(ANOX+BNOY+CNOY) and CNOY/(ANOX+BNOY+CNOY) to determine the nitrate Δ17O