source: palm/trunk/UTIL/chemistry/gasphase_preproc/kpp/util/util.f @ 3639

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

Merge of branch palm4u into trunk

File size: 2.8 KB
Line 
1C ****************************************************************
2C                           
3C InitSaveData - Opens the data file for writing
4C   Parameters :                                                 
5C
6C ****************************************************************
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
17C End of InitSaveData function
18C ****************************************************************
19
20C ****************************************************************
21C                           
22C SaveData - Write LOOKAT species in the data file 
23C   Parameters :                                                 
24C
25C ****************************************************************
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)
36999   FORMAT(E24.16,100(1X,E24.16))
37
38      RETURN
39      END
40
41C End of SaveData function
42C ****************************************************************
43
44C ****************************************************************
45C                           
46C CloseSaveData - Close the data file 
47C   Parameters :                                                 
48C
49C ****************************************************************
50
51      SUBROUTINE CloseSaveData ()
52
53      INCLUDE 'KPP_ROOT_Parameters.h'
54
55      CLOSE(10)
56
57      RETURN
58      END
59
60C End of CloseSaveData function
61C ****************************************************************
62
63C ****************************************************************
64C                           
65C GenerateMatlab - Generates MATLAB file to load the data file 
66C   Parameters : 
67C                It will have a character string to prefix each 
68C                species name with.                                                 
69C
70C ****************************************************************
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
83990   FORMAT(A1,'c = KPP_ROOT;')
84      write(20,*) 'clear KPP_ROOT;'
85      write(20,991) PREFIX, PREFIX
86991   FORMAT(A1,'t=',A1,'c(:,1);')
87      write(20,992) PREFIX
88992   FORMAT(A1,'c(:,1)=[];')
89
90      do i=1,NLOOKAT
91        write(20,993) PREFIX, SPC_NAMES(LOOKAT(i)), PREFIX, i
92993     FORMAT(A1,A6,' = ',A1,'c(:,',I2,');')
93      end do
94     
95      CLOSE(20)
96
97      RETURN
98      END
99
100C End of GenerateMatlab function
101C ****************************************************************
102
Note: See TracBrowser for help on using the repository browser.