Problem with aed2.nml running

[Léo Faivre]

Hello everybody, 

I'm currently facing a problem with AED2 Library. Infact hen I try to launch my simulation of GLM thanks to glmatools in R (through R studio) the following error shows and the rest of AED modules are not taken in count. Do you have any idea where does it come from and what can I do ? I also give you my parameter list attached to organic matter. 

The error is also occuring with phytoplankton parameters and zooplankton parameters when i desactivate module one per one to test.

I have to admit that I'm new to R and to GLM model. 

Thanks in advance for your help ! 

[matt.hipsey]

Hi Leo

I must admit I have not seen the IEEE_underflow warning, so I think there is something about the numbers or characters that are being misread, or sometimes there are hidden "gremlin" characters. International keyboards can also use slightly different characters in some cases. Do you have a working example of an add nil file you can swap in for this block? From a quick look it seems OK, but you can cross check your entries with this list int he code to make sure:

https://github.com/AquaticEcoDynamics/libaed2/blob/45e9d95cfc9d128f9c76ec96fccab8e096f22c1f/src/aed2_organic_matter.F90#L259

Also, please note AED2 is now getting a bit outdated as GLM-AED is the current supported version for the past couple of years.

https://github.com/AquaticEcoDynamics/glm-aed

If you are not too far into your setup I would recommend using this latest version, and note also there are some examples in there. This should also work with glmR.

Let us know how you go and happy to help get you running

Best wishes

Matt

[Léo Faivre]

Hi Matt, 

Thanks a lot for your help, finally it seems that indeed my aed2.nml file had too many digits in few parameters, it's no longer the case and my simulation is working fine. 
Thank you for the information about glm-aed latest version but I had some trouble at first when I tried to launch GLM with aed and notre aed2 so I finally preferred to stay on aed2. Since indeed my calibration is not began yet, i'll give it another try and find out if it's possible to use the latest AED version. 

Thanks again for your help and I will tell you if I face any problem with the latest AED version. 

Best wishes,

Leo.

--
You received this message because you are subscribed to a topic in the Google Groups "Aquatic Ecosystem MOdelling Network (AEMON) discussion forum" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/aquaticmodelling/MiDjEIZGKn4/unsubscribe.
To unsubscribe from this group and all its topics, send an email to aquaticmodelli...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/aquaticmodelling/3d1c6792-6ba3-4a7c-bd86-9c1cd0bf1a5fn%40googlegroups.com.

[Léo Faivre]

Hi Matt,

As you advised me last time, I converted my AED2 configuration files to AED and used the latest version of AED. However, for Cyanobacteria, I found in the AED manual and other examples of the AED.nml file that the settling regime is set to 3 (instead of 1 as for other phyto groups). I noticed that this configuration used the PHY rho (density of the group), c1 (rate coefficient for density increase), and c3 (minimum rate of density increase) for each group configured with settling on 3.

However, I was wondering if there is an equation that links all these variables together for a better understanding of the utilization of such parameters. Unfortunately, I couldn't find it in the manual or the GitHub files. Do you have any idea where I could find the equation for settling options for each phytoplankton group where c1, c3, and w_p parameters are used?

Thank you in advance for your precious help.
Best wishes,

Leo.