Changeset 6110 for LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Mar 17, 2026, 4:00:47 PM (3 weeks ago)

Author:

Sebastien Nguyen

Message:

fix bug that can make LMDZISO crash when trying to write Rland_ice before variable allocation and initialisation.

File:

• LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -3215 +3215 @@
     !write(*,*) 'tmp phys_output_write: ntiso=',ntiso
 
+    zx_tmp_fi2d=0. ; zx_tmp_fi3d=0.
     DO ixt = 1, ntiso
         !write(*,*) 'ixt,o_xtovap(ixt)=',ixt,o_xtovap(ixt)
@@ -3230 +3231 @@
         CALL histwrite_phy(o_xtpluc(ixt), zx_tmp_fi2d)
         CALL histwrite_phy(o_xtevap(ixt),   xtevap(ixt,:))
+
+        IF (vars_defined) zx_tmp_fi2d(:) = Rland_ice(ixt,:)
         CALL histwrite_phy(o_Rland_ice(ixt), Rland_ice(ixt,:))
         CALL histwrite_phy(o_xtovap(ixt),  xt_seri(ixt,:,:))
@@ -3236 +3239 @@
         !! runoff land bucket + soil isotopes ratio + river and costal flows PRSN 2025
         CALL histwrite_phy(o_xtrunoff_diag(ixt), xtrunoff_diag(ixt,:))
+          !! LMDZORISO variables set to 0 and not updating in bucket
         CALL histwrite_phy(o_xtriverflow(ixt),   xtriverflow(ixt,:))
         CALL histwrite_phy(o_xtcoastalflow(ixt),   xtcoastalflow(ixt,:))