Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
flexpart
flex_read_fortran
Commits
0af3b1c0
Commit
0af3b1c0
authored
Dec 17, 2018
by
Ignacio Pisso
Browse files
README is in markdown
parent
151589be
Changes
2
Hide whitespace changes
Inline
Side-by-side
README
deleted
100644 → 0
View file @
151589be
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
README.md
0 → 100644
View file @
0af3b1c0
#Read output from flexpart with fortran functions
## compilation
make (all)
or
make printheader
make printgrid
make flex_read_compare2
## test
make -f makefile_dist test output=$HOME/repos/flexpart/output/
## index
-
juldate.f
-
caldate.f
-
readheader.f
-
readgrid.f
-
readheader_stdout90.f90
-
readheader1090.f90
-
printgrid.f90
-
flex
\_
read_compare2.f90
-
makefile
## Print header to stdout
Compilation
make printheader
Usage
./printheader /path/to/output/
## Print grid stats to stdout
Compilation
make printgrid
Usage
./printgrid /path/to/output/ <unit>
unit= conc/pptv/time
## Compare gridded output
Compilation
make flex_read_compare2
Usage
./flex_read_compare2 /path/to/output/reference /path/to/output/new <unit>
unit= conc/pptv/time
## history
./prepare_footprint_modif
reads dirlist. Only bwd.
./sum_avg
reads dirlist. bwd OK. fwd?
./flexout_sum_max
reads dirlist. bwd OK. fwd?
./flexout_sum_max90
OK with dirlist.
./flexout_compare_default
Segmentation fault
./flexout_compare_default /home/ignacio/repos/flexpart/tests/examples/output_2-bwd/ time
STOP
./flex_read_compare /home/ignacio/repos/flexpart/tests/examples/output_2-bwd/ time
OK
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment