source: palm/trunk/UTIL/chemistry/gasphase_preproc/kpp/drv/general.f90 @ 4725

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

Merge of branch palm4u into trunk

File size: 1.3 KB
Line 
1PROGRAM KPP_ROOT_Driver
2
3  USE KPP_ROOT_Model
4  USE KPP_ROOT_Initialize, ONLY: Initialize
5
6      KPP_REAL :: T, DVAL(NSPEC)
7      KPP_REAL :: RSTATE(20)
8      INTEGER :: i
9 
10!~~~> Initialization
11
12      STEPMIN = 0.0d0
13      STEPMAX = 0.0d0
14
15      DO i=1,NVAR
16        RTOL(i) = 1.0d-4
17        ATOL(i) = 1.0d-3
18      END DO
19     
20      CALL Initialize()
21      CALL InitSaveData()
22
23!~~~> Time loop
24      T = TSTART
25kron: DO WHILE (T < TEND)
26
27        TIME = T
28        CALL GetMass( C, DVAL )
29        WRITE(6,991) (T-TSTART)/(TEND-TSTART)*100, T,       &
30                   ( TRIM(SPC_NAMES(MONITOR(i))),           &
31                     C(MONITOR(i))/CFACTOR, i=1,NMONITOR )
32        CALL SaveData()
33        CALL Update_SUN() 
34        CALL Update_RCONST()
35
36        CALL INTEGRATE( TIN = T, TOUT = T+DT, RSTATUS_U = RSTATE, &
37        ICNTRL_U = (/ 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 /) )
38        T = RSTATE(1)
39
40      END DO kron
41!~~~> End Time loop
42
43      CALL GetMass( C, DVAL )
44      WRITE(6,991) (T-TSTART)/(TEND-TSTART)*100, T,     &
45               ( TRIM(SPC_NAMES(MONITOR(i))),           &
46                 C(MONITOR(i))/CFACTOR, i=1,NMONITOR ) 
47      TIME = T
48      CALL SaveData()
49      CALL CloseSaveData()
50
51991   FORMAT(F6.1,'%. T=',E9.3,2X,200(A,'=',E11.4,'; '))
52
53END PROGRAM KPP_ROOT_Driver
54
Note: See TracBrowser for help on using the repository browser.