Changeset 2494 for LMDZ5/trunk/libf/phylmd

Timestamp:

Apr 12, 2016, 10:00:36 AM (9 years ago)

Author:

Ehouarn Millour

Message:

Fix to handle cases when q_ancien, t_ancien, etc. contain
dummy (constant) fields, as is the case when startphy.nc is
generated by ce0l or iniaqua: even though the fields are
present, they should not be used to compute corresponding
dynamical tendencies (done in physiq) for the first physics step.
EM

File:

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

LMDZ5/trunk/libf/phylmd/phyetat0.F90

@@ -327 +327 @@
   ancien_ok=ancien_ok.AND.phyetat0_get(klev,v_ancien,"VANCIEN","VANCIEN",0.)
 
+  ! Ehouarn: addtional tests to check if t_ancien, q_ancien contain
+  !          dummy values (as is the case when generated by ce0l,
+  !          or by iniaqua)
+  if ( (maxval(q_ancien).eq.minval(q_ancien)) .or. &
+       (maxval(t_ancien).eq.minval(t_ancien)) ) then
+    ancien_ok=.false.
+  endif
+
   found=phyetat0_get(klev,clwcon,"CLWCON","CLWCON",0.)
   found=phyetat0_get(klev,rnebcon,"RNEBCON","RNEBCON",0.)