Hgardfou reports NaNs in temperature field - Revision history

Maureenjcohen at 07:43, 18 October 2024

2024-10-18T07:43:18Z

Maureenjcohen: /* Age of air tracer included in chemistry array */

2024-10-10T12:22:08Z

‎Age of air tracer included in chemistry array

Maureenjcohen: /* Wrong value of nbapp_chem */

2024-10-10T10:10:54Z

‎Wrong value of nbapp_chem

Maureenjcohen at 10:09, 10 October 2024

2024-10-10T10:09:50Z

Maureenjcohen at 10:08, 10 October 2024

2024-10-10T10:08:19Z

Maureenjcohen: Created page with "== The symptoms == The log gcm.out reads:

 i,k,temperature =               1                  1               NaN i,k,temperature =               2                  1..."

2024-10-10T10:06:07Z

Created page with "== The symptoms == The log gcm.out reads: <pre> i,k,temperature = 1 1 NaN i,k,temperature = 2 1..."

New page

== The symptoms ==
The log gcm.out reads:
<pre>
i,k,temperature = 1 1 NaN
i,k,temperature = 2 1 NaN
i,k,temperature = 3 1 NaN
etc.
</pre>
At the end of a long report of NaNs, there is a message like:
<pre>
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP 1
</pre>
and
<pre>Stopping in hgardfou
Reason = hgardfou stops debutphy
Houston, we have a problem, ierr = 1
</pre>

== Known causes ==
The message indicates numerical instability in the temperature field, which can be caused by different things. It may happen quickly when the simulation starts or more slowly after significant runtime. Note: even though the reported NaNs started at gridbox i,k 1,1, this doesn't mean the instability began there - the hgardfou subroutine checks the temperature field starting from level 1, so it will report bad values in this order.

=== Wrong value of callthermos ===
(October 2024): In physiq.def, the flag callthermos must be:
<pre>
callthermos = y
</pre>
when running a chemistry simulation with 78 levels. A value of **n** causes an immediate crash.

=== Wrong value of nbapp_chem ===
(October 2024): In physiq.def, the flag nbapp_chem must be:
<pre>
nbapp_chem = 24000
</pre>
Other values cause a crash after 0.1-0.2 Venus days.

=== Age of air tracer included in chemistry array ===
(October 2024): The age of air cannot be included in chemistry runs. The flags in physiq.def should be:
<pre>
ok_chem = y
ok_aoa = n
</pre>
or
<pre>
ok_chem = n
ok_aoa = n
</pre>
Or both can be set to **n**. This also causes a crash after 0.1-0.2 Venus days. There is a check in conf_phys which should warn the user if the settings are incoherent and stop the simulation.

@@ Line 47: / Line 47: @@
 </pre>
 Or both can be set to '''n'''. This also causes a crash after 0.1-0.2 Venus days. There is a check in conf_phys which should warn the user if the settings are incoherent and stop the simulation.
 [[Category:Errors]]
 [[Category:Venus-Model]]
 [[Category:Venus-LMDZ]]

@@ Line 44: / Line 44: @@
 <pre>
 ok_chem = n
-ok_aoa = n
+ok_aoa = y
 </pre>
 Or both can be set to '''n'''. This also causes a crash after 0.1-0.2 Venus days. There is a check in conf_phys which should warn the user if the settings are incoherent and stop the simulation.

← Older revision		Revision as of 10:09, 10 October 2024
Line 47:		Line 47:
	</pre>		</pre>
	Or both can be set to '''n'''. This also causes a crash after 0.1-0.2 Venus days. There is a check in conf_phys which should warn the user if the settings are incoherent and stop the simulation.		Or both can be set to '''n'''. This also causes a crash after 0.1-0.2 Venus days. There is a check in conf_phys which should warn the user if the settings are incoherent and stop the simulation.
		+
		+	[[Category:Errors]]
		+	[[Category:Venus-Model]]
		+	[[Category:Venus-LMDZ]]

@@ Line 33: / Line 33: @@
 nbapp_chem = 24000
 </pre>
-Other values cause a crash after 0.1-0.2 Venus days.
+Other values can cause a slow crash after 0.1-0.2 Venus days.
 === Age of air tracer included in chemistry array ===

@@ Line 26: / Line 26: @@
 callthermos = y
 </pre>
-when running a chemistry simulation with 78 levels. A value of **n** causes an immediate crash.
+when running a chemistry simulation with 78 levels. A value of '''n''' causes an immediate crash.
 === Wrong value of nbapp_chem ===
@@ Line 46: / Line 46: @@
 ok_aoa = n
 </pre>
-Or both can be set to **n**. This also causes a crash after 0.1-0.2 Venus days. There is a check in conf_phys which should warn the user if the settings are incoherent and stop the simulation.
+Or both can be set to '''n'''. This also causes a crash after 0.1-0.2 Venus days. There is a check in conf_phys which should warn the user if the settings are incoherent and stop the simulation.