Minor edits (more fatal info)

schism-dev · Sep 26, 2023 · 1e60664 · 1e60664
1 parent 4ba9298
commit 1e60664
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 8 deletions.
diff --git a/src/Hydro/misc_subs.F90 b/src/Hydro/misc_subs.F90
@@ -475,10 +475,10 @@ subroutine other_hot_init(time)
       enddo !k
 
 !     The following to init th_dt*, th_time* and ath* is only done by
-!     rank 0 and but bcast'ed, b/c in _step we'll continue the reading
-!     by rank0 and only bcast the final
+!     rank 0, not bcast'ed, b/c in _step we'll continue the reading
+!     by rank 0 and only bcast the final
 !     products of eth, trth (since they use global indices) etc, 
-!     and the th_dt[12], th_time[12] and ath[12] are not used further 
+!     and the th_dt[-,12], th_time[-,12] and ath[-,12] are not used further 
       if(myrank==0) then
 !-----------------------------------------------------------------------------
 !...  Init reading t.h. files 

diff --git a/src/Hydro/schism_step.F90 b/src/Hydro/schism_step.F90
@@ -1333,7 +1333,7 @@ subroutine schism_step(it)
 
         rat=(time-th_time(1,1,1))/th_dt(1,1)
         if(rat<-small1.or.rat>1.d0+small1) then
-          write(errmsg,*) 'STEP: rat out in elev.th:',rat,time,th_time(1,1:2,1)
+          write(errmsg,*) 'STEP: rat out in elev.th:',rat,time,th_time(1,1:2,1),th_dt(1,1)
           call parallel_abort(errmsg)
         endif
         icount=0
@@ -1391,7 +1391,7 @@ subroutine schism_step(it)
 
             rat=(time-th_time(m,1,5))/th_dt(m,5)
             if(rat<-small1.or.rat>1.d0+small1) then
-              write(errmsg,*) 'STEP: rat out in htr_.th:',rat,time,th_time(m,1:2,5)
+              write(errmsg,*) 'STEP: rat out in htr_.th:',rat,time,th_time(m,1:2,5),th_dt(m,5)
               call parallel_abort(errmsg)
             endif
             icount=0
@@ -1427,7 +1427,7 @@ subroutine schism_step(it)
 
         rat=(time-th_time2(1,1))/th_dt2(1)
         if(rat<-small1.or.rat>1.d0+small1) then
-          write(errmsg,*) 'STEP: rat out in elev2D.th:',rat,time,th_time2(1:2,1)
+          write(errmsg,*) 'STEP: rat out in elev2D.th:',rat,time,th_time2(1:2,1),th_dt2(1)
           call parallel_abort(errmsg)
         endif
         icount=0
@@ -1478,7 +1478,7 @@ subroutine schism_step(it)
 
         rat=(time-th_time2(1,2))/th_dt2(2)
         if(rat<-small1.or.rat>1.d0+small1) then
-          write(errmsg,*) 'STEP: rat out in uv3D.th:',rat,time,th_time2(1:2,2)
+          write(errmsg,*) 'STEP: rat out in uv3D.th:',rat,time,th_time2(1:2,2),th_dt2(2)
           call parallel_abort(errmsg)
         endif
         icount=0
@@ -1534,7 +1534,7 @@ subroutine schism_step(it)
         if(ntrs(i)>0.and.nnode_tr2(i)>0) then
           rat=(time-th_time2(1,5))/th_dt2(5)
           if(rat<-small1.or.rat>1.d0+small1) then
-            write(errmsg,*) 'STEP: rat out in tr3D.th:',rat,time,th_time2(1:2,5)
+            write(errmsg,*) 'STEP: rat out in tr3D.th:',rat,time,th_time2(1:2,5),th_dt2(5)
             call parallel_abort(errmsg)
           endif
 !          icount=0