Last change
on this file since 2938 was
2696,
checked in by kanani, 7 years ago
|
Merge of branch palm4u into trunk
|
File size:
1.1 KB
|
Rev | Line | |
---|
[2696] | 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.