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

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