Changeset 3843 for trunk/LMDZ.MARS/libf

Timestamp:

Jul 9, 2025, 3:22:31 PM (5 days ago)

Author:

jbclement

Message:

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

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (2 diffs)

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