Custom Query (210 matches)

I have built Flexpart from source on a couple of Linux machines and I always get a compile error due to a problem in erf.f90. The error message is:

gfortran -c -O2 -m64 -mcmodel=medium -fconvert=little-endian -frecord-marker=4 -I/usr/lib64/gfortran/modules erf.f90 erf.f90:106.26:

real, external

gammln 1

Error: Return type mismatch of function 'gammln' at (1) (REAL(4)/REAL(8)) erf.f90:111.13:

gln=gammln(a)

1

Warning: Type mismatch in argument 'xx' at (1); passed REAL(4) to REAL(8)

The same error happens at line 140. The following patch allows Flexpart to compile although it may not be the best way of achieving this:

aurora.hcp>diff baderf.f90 erf.f90 106c106 < real, external :: gammln ---

double precision, external

gammln

111c111 < gln=gammln(a) ---

gln=gammln(DBLE(a))

140c140 < real, external :: gammln ---

double precision, external

gammln

145c145 < gln=gammln(a) ---

gln=gammln(DBLE(a))

The version of gfortran used does not appear to be important; the same thing happens with gfortran 4.8.5 on Scientific Linux 7 and with gfortran 6.3.0 on Debian 9.2.

Dear all,

Anyone interested in implementing the MQUASILAG in backward mode?

Preferably in FP v9.02, which I currently use. But if there is a reason to rather do it already in v10 its also fine, I can wait for it.

I read that some people think that the option in forward mode is not properly implemented. So I think this can be dealt with at the same time.

(Im not sure yet if there is some money for this from the project, but at least I can offer co-authorship.)

Thanks!

I need to get particle positions after I release the particles in my measurement site in order to have the number of trajectory segment endpoints in cell ijth in T time interval. I need particle positions in order to obtain Residence Time Analysis and subsequently, Concentration Field Analysis. My question is how can I get the particle positions? Is CONC output parameter (in netCDF format) with s. m^{3 /kg unit representing the residence time of each grid divided by concentration at the measurement site? If not, what is this parameter representing? Many Thanks in advance for your kind help.}