Add check in util script for vectors

schism-dev · Aug 15, 2024 · a8cca85 · a8cca85
1 parent 915ef34
commit a8cca85
Show file tree

Hide file tree

Showing 5 changed files with 20 additions and 5 deletions.
diff --git a/src/Utility/Post-Processing-Fortran/read_output10_allnodes.f90 b/src/Utility/Post-Processing-Fortran/read_output10_allnodes.f90
@@ -221,7 +221,10 @@ program read_out
 
       !Find nc file
       file63=varname(1:len_var)//'_'//it_char(1:leng)//'.nc'
-      if(ivs==2) varname2=varname(1:len_var-1)//'Y'
+      if(ivs==2) then
+        if(varname(len_var:len_var).ne."X") stop 'Vector name must end with X'
+        varname2=varname(1:len_var-1)//'Y'
+      endif
       inquire(file=file63,exist=lexist)
       if(lexist) then
         i23d=2 !3D var

diff --git a/src/Utility/Post-Processing-Fortran/read_output10_transect.f90 b/src/Utility/Post-Processing-Fortran/read_output10_transect.f90
@@ -234,7 +234,10 @@ program read_out
 
       !Find nc file
       file63=varname(1:len_var)//'_'//it_char(1:leng)//'.nc'
-      if(ivs==2) varname2=varname(1:len_var-1)//'Y'
+      if(ivs==2) then
+        if(varname(len_var:len_var).ne."X") stop 'Vector name must end with X'
+        varname2=varname(1:len_var-1)//'Y'
+      endif
       inquire(file=file63,exist=lexist)
       if(lexist) then !3D var
         i23d=2 

diff --git a/src/Utility/Post-Processing-Fortran/read_output10_xyt.f90 b/src/Utility/Post-Processing-Fortran/read_output10_xyt.f90
@@ -228,7 +228,10 @@ program read_out
 
           !Find nc file
           file63=varname(1:len_var)//'_'//it_char(1:leng)//'.nc'
-          if(ivs==2) varname2=varname(1:len_var-1)//'Y'
+          if(ivs==2) then
+            if(varname(len_var:len_var).ne."X") stop 'Vector name must end with X'
+            varname2=varname(1:len_var-1)//'Y'
+          endif
           inquire(file=file63,exist=lexist)
           if(lexist) then !3D var
             i23d=2 

diff --git a/src/Utility/Post-Processing-Fortran/read_output10_xyz.f90 b/src/Utility/Post-Processing-Fortran/read_output10_xyz.f90
@@ -208,7 +208,10 @@ program read_out
 
       !Find nc file
       file63=varname(1:len_var)//'_'//it_char(1:leng)//'.nc'
-      if(ivs==2) varname2=varname(1:len_var-1)//'Y'
+      if(ivs==2) then
+        if(varname(len_var:len_var).ne."X") stop 'Vector name must end with X'
+        varname2=varname(1:len_var-1)//'Y'
+      endif
       inquire(file=file63,exist=lexist)
       if(lexist) then !3D var
         i23d=2 

diff --git a/src/Utility/Post-Processing-Fortran/read_output10_xyzt.f90 b/src/Utility/Post-Processing-Fortran/read_output10_xyzt.f90
@@ -210,7 +210,10 @@ program read_out
 
           !Find nc file
           file63=varname(1:len_var)//'_'//it_char(1:leng)//'.nc'
-          if(ivs==2) varname2=varname(1:len_var-1)//'Y'
+          if(ivs==2) then
+            if(varname(len_var:len_var).ne."X") stop 'Vector name must end with X'
+            varname2=varname(1:len_var-1)//'Y'
+          endif
           inquire(file=file63,exist=lexist)
           if(lexist) then !3D var
             i23d=2