Last change
on this file since 3698 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 | |
---|
8 | C ---- 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 | |
---|
28 | C -- BELOW THIS LIMIT USE TAYLOR INSTEAD OF EXP --- |
---|
29 | C ********** 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 |
---|
44 | C *********** END TIME LOOP ******** |
---|
45 | |
---|
46 | write(10,*) '];' |
---|
47 | close(10) |
---|
48 | STOP |
---|
49 | |
---|
50 | 991 FORMAT('Monitor:',10(1X,E12.6)) |
---|
51 | 992 FORMAT(F6.1,100(1X,D24.16)) |
---|
52 | END |
---|
53 | |
---|
Note: See
TracBrowser
for help on using the repository browser.