Changeset 3843

Timestamp:

Jul 9, 2025, 3:22:31 PM (26 hours ago)

Author:

jbclement

Message:

Mars PCM 1D:
Better error handling for reading "start1D.txt" and "profile_prescribed_h2o_vap" files.
JBC

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4901 +4901 @@
 == 07/07/2025 == JBC
 In "display_netcdf.py", addition of the possibility to display altitude as function of time for 1D variables + showing value at cursor for 2D heatmaps.
+
+== 09/07/2025 == JBC
+Better error handling for reading "start1D.txt" and "profile_prescribed_h2o_vap" files.

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -518 +518 @@
 else
     do iq = 1,nq
-        read(3,*) header, (qsurf(1,iq,j), j = 1,size(qsurf,3)), (q(1,ilayer,iq), ilayer = 1,nlayer)
+        read(3,*,iostat = ierr) header, (qsurf(1,iq,j), j = 1,size(qsurf,3)), (q(1,ilayer,iq), ilayer = 1,nlayer)
+        if (ierr /= 0) error stop 'Not enough atmospheric layers defined in the file "'//trim(start1Dname)//'" for the tracer "'//trim(header)//'"!'
         if (trim(tname(iq)) /= trim(header)) then
             write(*,*) 'Tracer names between "traceur.def" and "'//trim(start1Dname)//'" do not match!'
@@ -793 +794 @@
         if (ierr == 0) then
             do ilayer = 1,nlayer
-                read(10,*) q_prescribed_h2o_vap(ilayer)
+                read(10,*,iostat=ierr) q_prescribed_h2o_vap(ilayer)
+                if (ierr /= 0) error stop 'Not enough atmospheric layers defined in the file "profile_prescribed_h2o_vap"!'
             enddo
         else