chem/KPP/kpp/kpp-2.1/util/util.f

   1 C ****************************************************************
   2 C
   3 C InitSaveData - Opens the data file for writing
   4 C   Parameters :
   5 C
   6 C ****************************************************************
   7
   8       SUBROUTINE InitSaveData ()
   9
  10       INCLUDE 'KPP_ROOT_Parameters.h'
  11
  12       open(10, file='KPP_ROOT.dat')
  13
  14       RETURN
  15       END
  16
  17 C End of InitSaveData function
  18 C ****************************************************************
  19
  20 C ****************************************************************
  21 C
  22 C SaveData - Write LOOKAT species in the data file
  23 C   Parameters :
  24 C
  25 C ****************************************************************
  26
  27       SUBROUTINE SaveData ()
  28
  29       INCLUDE 'KPP_ROOT_Parameters.h'
  30       INCLUDE 'KPP_ROOT_Global.h'
  31
  32       INTEGER i
  33
  34       WRITE(10,999) (TIME-TSTART)/3600.D0,
  35      *              (C(LOOKAT(i))/CFACTOR, i=1,NLOOKAT)
  36 999   FORMAT(E24.16,100(1X,E24.16))
  37
  38       RETURN
  39       END
  40
  41 C End of SaveData function
  42 C ****************************************************************
  43
  44 C ****************************************************************
  45 C
  46 C CloseSaveData - Close the data file
  47 C   Parameters :
  48 C
  49 C ****************************************************************
  50
  51       SUBROUTINE CloseSaveData ()
  52
  53       INCLUDE 'KPP_ROOT_Parameters.h'
  54
  55       CLOSE(10)
  56
  57       RETURN
  58       END
  59
  60 C End of CloseSaveData function
  61 C ****************************************************************
  62
  63 C ****************************************************************
  64 C
  65 C GenerateMatlab - Generates MATLAB file to load the data file
  66 C   Parameters :
  67 C                It will have a character string to prefix each
  68 C                species name with.
  69 C
  70 C ****************************************************************
  71
  72       SUBROUTINE GenerateMatlab ( PREFIX )
  73
  74       INCLUDE 'KPP_ROOT_Parameters.h'
  75       INCLUDE 'KPP_ROOT_Global.h'
  76
  77       CHARACTER*8 PREFIX
  78       INTEGER i
  79
  80       open(20, file='KPP_ROOT.m')
  81       write(20,*) 'load KPP_ROOT.dat;'
  82       write(20,990) PREFIX
  83 990   FORMAT(A1,'c = KPP_ROOT;')
  84       write(20,*) 'clear KPP_ROOT;'
  85       write(20,991) PREFIX, PREFIX
  86 991   FORMAT(A1,'t=',A1,'c(:,1);')
  87       write(20,992) PREFIX
  88 992   FORMAT(A1,'c(:,1)=[];')
  89
  90       do i=1,NLOOKAT
  91         write(20,993) PREFIX, SPC_NAMES(LOOKAT(i)), PREFIX, i
  92 993     FORMAT(A1,A6,' = ',A1,'c(:,',I2,');')
  93       end do
  94
  95       CLOSE(20)
  96
  97       RETURN
  98       END
  99
 100 C End of GenerateMatlab function
 101 C ****************************************************************
 102