Changeset 2048

Timestamp:

May 20, 2014, 1:38:23 PM (11 years ago)

Author:

lguez

Message:

Added a check. nstra can be 0 for small values of llm. If nstra
computed by sugwd_strato is 0 then there is an array index out of
bounds later in gwprofil_strato. The check does not cost much because
sugwd_strato is called only once per run.

File:

• LMDZ5/trunk/libf/phylmd/orografi_strato.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/orografi_strato.F90

@@ -1873 +1873 @@
   PRINT *, ' DANS SUGWD nktopg=', nktopg
   PRINT *, ' DANS SUGWD nstra=', nstra
+  if (nstra == 0) call abort_gcm("sugwd_strato", "no level in stratosphere", 1)
 
   gsigcr = 0.80