README 2.21 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
aim: read output from flexpart default and parameter sweep with fortran

history: 
2018-08-16
cloned from /Users/ignacio/todo/doing/flexpart/2018-06-03_Sun22/PREPARE_FOOTPRINT

modified read header and read grid from f77 to f90

added

readheader10.f
 from readheader10.f modified to be more verbose
readheader_stdout.f
 wrapper to call the above function and display the read variables

readheader1090.f90
readheader_stdout90.f90
 translated into free f90 format from above



examples from before
sum_avg
 compilation:
	make sum_avg	
 run:
	# sum_avg needs dir to be in file
	./sum_avg conc

flexout_sum_max
 compilation:
	make flexout_sum_max	
 run:
	./flexout_sum_max /Users/ignacio/repos/flex_cases/output/ conc


read and display header

f77
make readheader_stdout
./readheader_stdout 
f90
make readheader_stdout90
./readheader_stdout90

read and display grid



read and display particles

read and verify a fingerprint

particles: some central measure or other moments


2018-08-22
can read (header and grids) group 1 of standard use cases

08.22 22.58 ignacio@Kjeller:~/repos/flex_read_fortran$./flexout_sum_max /Users/ignacio/repos/flex_gen_input/output_2/ pptv
08.22 22.58 ignacio@Kjeller:~/repos/flex_read_fortran$./flexout_sum_max /Users/ignacio/repos/flex_gen_input/output_3/ conc
08.22 22.58 ignacio@Kjeller:~/repos/flex_read_fortran$./flexout_sum_max /Users/ignacio/repos/flex_gen_input/output_3/ pptv
08.22 22.58 ignacio@Kjeller:~/repos/flex_read_fortran$./flexout_sum_max /Users/ignacio/repos/flex_gen_input/output_5/ conc
 

2018-08-23

compare with matlab

case 1A _2

./flexout_sum_max /Users/ignacio/repos/flex_cases_reference/output_2/ pptv
 sum_sum_grid_out=   90265168.0    
 max_max_grid_out=   28875172.0   

with matlab:
% diagnostics
% sumgrid=90265170.9066
% maxgrid=28875172

CAse 1A _3 

./flexout_sum_max /Users/ignacio/repos/flex_cases_reference/output_3/ pptv
 sum_sum_grid_out=   90265168.0    
 max_max_grid_out=   28875172.0   

./flexout_sum_max /Users/ignacio/repos/flex_cases_reference/output_3/ conc
 sum_sum_grid_out=   101040392.    
 max_max_grid_out=   32318520.0    


./flexout_sum_max /Users/ignacio/repos/flex_cases_reference/output_5/ conc
 sum_sum_grid_out=   101040392.    
 max_max_grid_out=   32318520.0