chem/KPP/kpp/kpp-2.1/drv/main.f

   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