Minor change in aspect ratio script

Author: josephzhang8

src/Readme.beta_notes

@@ -9,7 +9,7 @@ This beta notes are divided into 5 part:
 
 Github versions:
 
-- `8e0b1477`: removed/renamed vclose_surf_frac as vclose_surf_frac0. If vclose_surf_frac0<0, needs vclose.gr3 (depth in [0,1]);
+- `8e0b1477` (30 Sept 2025): removed/renamed vclose_surf_frac as vclose_surf_frac0. If vclose_surf_frac0<0, needs vclose.gr3 (depth in [0,1]);
 
 Tag stofs3d-atl.v3.1.0 (from v5.14): Sept 17, 2025
 - `b3644238` (21 Aug 2025): added istemp as on/off and removed stemp_thick; restored soil_thick.gr3;
@@ -640,6 +640,7 @@ git versions:
 159. `1e2b181209` (21 Aug 2025): merged with shallow_T, quad_init (fixed a bug in quad_int() c/o Jose O.)
 160. `6c2ef8f3` (25 Aug 2025): updated WWM nc outputs to use OR(NF90_NETCDF4,NF90_CLOBBER) (as in hydro)
 161. `b2c4d490`: changed soil conductivity in a cell from min to mean;
+162. `8f61b12e` (10 Oct 2025): fixed a bug (uninit'ed var) in fix_bad_quads.f90
 
 ================================================================================================
 (D) Auto-test history:

src/Utility/Pre-Processing/eq_diameter_aspect_ratio.f90

@@ -1,6 +1,6 @@
 !     Compute equivalent diameter/radius of a tri-quad grid, and dx/dz (aspect ratio; dx is eq. diameter)
 !     Input: hgrid.gr3 (projection)
-!     Output: fort.12 (*.prop format, eq. radius) and asrat.gr3 (aspect ratio for h>5m; -9999 otherwise)
+!     Output: fort.12 (*.prop format, eq. radius) and asrat.gr3 (aspect ratio for h>5m; 9999 otherwise. A.R.<1 indicates problem)
 !     ifx -Bstatic -O3 -o eq_diameter_aspect_ratio eq_diameter_aspect_ratio.f90
 !
 !     Below is matlab code for plotting histogram for radii
@@ -78,7 +78,7 @@
       write(13,*); write(13,*)ne,np
       do i=1,np
         if(dp(i)<=5.) then
-          ar=-9999.
+          ar=9999.
         else
           ar=rad_node(i)*2/dp(i)
         endif