Page 1 of 1

Almost zero mass

Posted: September 20th, 2020, 8:10 pm
by rosatrancoso
Hi

Just wondering if there is a way to increase the threshold of "zero mass".

This is because I have a lot of simulation runtime with masses as low as 1e-30 it's not considered "zero" because for fortran real*8 is below 1e-38 (with option -funsafe-math-optimizations in gfortan otherwise the threshold is as low as 1e-45).

In the code it seems particles are only deleted when rmass <0.0, i.e. a hard zero.

If no hysplit option, what Fortran compilation flags do you recommend for intel compiler?

Thanks in advance
Cheers
Rosa

Re: Almost zero mass

Posted: September 23rd, 2020, 1:07 pm
by alicec
You can use the
FRMR
namelist variable.
https://ready.arl.noaa.gov/hysplitusersguide/S630.htm

FRMR=0.01
would remove particles that had mass less than 0.01*total_mass.