source: palm/trunk/UTIL/chemistry/gasphase_preproc/kpp/drv/main.f @ 3359

Last change on this file since 3359 was 2696, checked in by kanani, 7 years ago

Merge of branch palm4u into trunk

File size: 1.1 KB
Line 
1      PROGRAM driver
2
3      INCLUDE 'KPP_ROOT_Parameters.h'
4      INCLUDE 'KPP_ROOT_Global.h'
5
6      INTEGER i
7 
8C ---- TIME VARIABLES ------------------     
9
10      RTOLS = 1e-3
11      do i=1,NVAR
12        RTOL(i) = RTOLS
13        ATOL(i) = 1E-18
14      end do
15     
16      CALL Initialize()
17     
18      open(10, file='Extrapd.m')
19      write(10,*) 'ed=['
20       
21      TSTART = 3600*12
22      TEND = TSTART + 3600*24*5
23      STEPMIN = 0.01
24      STEPMAX = 900
25      DT = 3600.
26      TEMP = 236.21
27     
28C -- BELOW THIS LIMIT USE TAYLOR INSTEAD OF EXP ---
29C ********** TIME LOOP *************************
30      TIME = TSTART
31      do while (TIME .le. TEND)
32
33        write(6,991) (C(MONITOR(i))/CFACTOR, i=1,NMONITOR)
34
35        write(10,992) (TIME-TSTART)/3600.D0, 
36     *                (C(LOOKAT(i))/CFACTOR, i=1,NLOOKAT)
37
38        call Update_SUN() 
39        call Update_RCONST()
40
41        call INTEGRATE( TIME, TIME+DT )
42
43      end do
44C *********** END TIME LOOP ********
45
46      write(10,*) '];'
47      close(10)
48      STOP
49
50991   FORMAT('Monitor:',10(1X,E12.6))
51992   FORMAT(F6.1,100(1X,D24.16))
52      END
53
Note: See TracBrowser for help on using the repository browser.