GAS-GRAIN MODEL

A.
-	cp jcddata* jcddata

-	f90 -dp jeqwrt99.16.f

-	a.out<jcddata

-	look up DISK7 for errors

-	look up in DISK23 for how many FCHM routines are called
	#rats routines = # called
	rename DISK23 jscdfr99.16.f

-	remove all DISK*

-	edit Cd.com - jscdfr99.16.f

-	sh Cd.com - produces jscdfr99.16.o

-	if the Temperature is outside 10-15K open jchnsmc.00.m1.mod1.f and modify 
	it. Currently it is set to handle T only within the 10-15 range.

-	cp jchnsmc.00.m1.mod1.f jchnsmc.f

-	look up in jchnsmc.f for rats and fchm routines called

-	edit C.com - to compile jchnsmc.f

-	sh C.com

-	edit L.com - link jchnsmc.o jscdfr99.16.o sgear.o
	jchnsmc.e - main program

-	sh L.com

-	edit B.com - use jchnsmc.e as a main program
	change the directory if necessary
	
- 	make sure no B.log, DISK* or jcdoutput exist

- 	sh A.com
_________________________________________________________________________________
B. How to change the density

- 	open the jcddata file and change the number of the density. It scales 
  	nicely, so no need to do anything else

-	go to A.
_________________________________________________________________________________

C. How to change the temperature

- 	surf.00.m1.e (in the sur directory which is subdire of gas-grain)

	-enter input name: surf.00.m1.in
	-enter output name: surf.00.m1.11k
	-do you want jcddata output made -yes
	-enter the jcddata name: jcd.m1.11k (!be consistent with the names)
	-option to choose model: currently m1 - so "1"
		otherwise m2 - "2" and other "3"
	-do you want to change the T - yes
	-enter the T (currently 11)
	-do you want a jmod.m - yes
	-enter the name: jmod.m1.11k
	
-	edit jmod.m1.11k

	- C + NO ->2 - correct the rates (the first number only), the number 
	should be twice as little
	- CH + NH -> 3 - correct rates are three times as little
	- NH + NH -> 2 - correct rates are twice as little
	- O + HCO -> 2 - correct rates are twice as little
	
	-for the reactions 14 263-271 - all the rates are 0.00 (all the numbers)
	
-	put back jmod.m1.11k into gas-grain/input directory, where the res of the 
	jmod.m* files are!!!

-	pick a data file - jcddata.00.1

	-rename it jcddata.00.1.11k (be consistent)
	-remove all the reactions of type 14,15,16
	-copy the reactions 14 , 1-271 from jmod.m1.11k in the place of the 
	already deleted reactions in jcddata.00.1.11k
	-copy all the reaction 15,16 from jmod.m1.11k into jcddata.00.1.11k
	-in the end - in the comments - correct the name of the file
	-change the temperature : TEMPO = 1.1e1K
	
-	go to A.
_________________________________________________________________________