source: palm/trunk/UTIL/chemistry/gasphase_preproc/kpp/util/util.f90 @ 4116

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

Merge of branch palm4u into trunk

File size: 2.8 KB
Line 
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)
35999   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
82990   FORMAT(A1,'c = KPP_ROOT;')
83      write(20,*) 'clear KPP_ROOT;'
84      write(20,991) PREFIX, PREFIX
85991   FORMAT(A1,'t=',A1,'c(:,1);')
86      write(20,992) PREFIX
87992   FORMAT(A1,'c(:,1)=[];')
88
89      do i=1,NLOOKAT
90        write(20,993) PREFIX, SPC_NAMES(LOOKAT(i)), PREFIX, i
91993     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
Note: See TracBrowser for help on using the repository browser.