Commit d79e1653 authored by Ignacio Pisso's avatar Ignacio Pisso
Browse files

prepare for flexpart 10 paper

parent 9c87b65d
......@@ -6,7 +6,7 @@ FC = gfortran
INCPATH = /xnilu_wrk/flex_wrk/bin64/grib_api/include
LIBPATH1 = /xnilu_wrk/flex_wrk/bin64/grib_api/lib
LIBPATH2 = /flex_wrk/flexpart/lib64/gfortran/lib/
FFLAGS = -O2 -m64 -mcmodel=medium -fconvert=little-endian -frecord-marker=4 -I$(INCPATH)
FFLAGS = -g -fcheck=all -Wall -fbacktrace -O2 -m64 -mcmodel=medium -fconvert=little-endian -frecord-marker=4 -I$(INCPATH)
# LDFLAGS = $(FFLAGS) -L$(LIBPATH2) -L$(LIBPATH1) -lgrib_api_f90 -lgrib_api -lm -ljasper
# linker flags not necessary: no need for grib files
......
program flexpartplot_latlon
parameter(nxmax=360,nymax=200,nzmax=3,maxtime=40)
parameter(maxpoint=5000,maxspec=88,maxageclass=1)
parameter(nxmax=720,nymax=360,nzmax=40,maxtime=10)
parameter(maxpoint=1,maxspec=1,maxageclass=1)
double precision jul,julstart,juldate
real outheight(nzmax)
......@@ -83,10 +83,10 @@ C Loop about all times, given in file dates
if (mother_or_nest.eq.1) then
filegrid=dirname(1:len)//'/grid_time_'//aday//atime//'_001'
fileout=dirname(1:len)//'footprint_'//aday//atime
fileout=dirname(1:len)//'footprint_'//aday//atime//'_001'
else
filegrid=dirname(1:len)//'grid_time_nest_'//aday//atime//'_001'
fileout=dirname(1:len)//'footprint_nest_'//aday//atime
fileout=dirname(1:len)//'footprint_nest_'//aday//atime//'_001'
endif
write(*,*) filegrid
......@@ -165,9 +165,9 @@ C Dump time-integrated footprint
********************************
if (mother_or_nest.eq.1) then
fileout=dirname(1:len)//'footprint_total'
fileout=dirname(1:len)//'footprint_total'//'_001'
else
fileout=dirname(1:len)//'footprint_total_nest'
fileout=dirname(1:len)//'footprint_total_nest'//'_001'
endif
open(97,file=fileout,form='unformatted')
write(97) itime
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment