Fixed NumPy compatibility issue in PyNGA, small improvements and typo fixes to other codes

Author: fabiolsilva

gmsvtoolkit/metrics/calc_gmpe.py

@@ -2,7 +2,7 @@
 """
 BSD 3-Clause License
 
-Copyright (c) 2023, University of Southern California
+Copyright (c) 2026, University of Southern California
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without
@@ -141,9 +141,6 @@ def parse_arguments(self):
                             help="output directory")
         parser.add_argument("--gmpe-group", dest="gmpe_group", required=True,
                             help="GMPE group %s" % (self.gmpe_groups))
-        parser.add_argument("-o", "--output", "--output-file",
-                            dest="output_file",
-                            help="output rd100 file")
         parser.add_argument("--station-list", "-s", dest="station_list", required=True,
                             help="station list for batch processing")
         parser.add_argument("--src-file", "--src", dest="src_file", required=True,

gmsvtoolkit/models/pynga/AS08.py

@@ -710,7 +710,7 @@ def calc_sigma_tau(self):
                       2. * alpha * tauB_T * tauB_PGA * self.rho[indT])
         sigmaTot = np.sqrt(sigma**2 + tau**2)
 
-        return (sigma, tau, sigmaTot)
+        return (np.float64(sigma[0]), np.float64(tau[0]), np.float64(sigmaTot[0]))
 
 def AS08nga_test(T, CoefTerms):
     """

gmsvtoolkit/models/pynga/BA08.py

@@ -434,9 +434,13 @@ def compute_std(self):
         try:
             ind = (np.array(self.periods) == self.T).nonzero()[0]
             if FT == 'U':
-                return (self.sigma_TU[ind], self.tau_U[ind], self.sigma0[ind])
+                return (np.float64(self.sigma_TU[ind][0]),
+                        np.float64(self.tau_U[ind][0]),
+                        np.float64(self.sigma0[ind][0]))
             else:
-                return (self.sigma_TM[ind], self.tau_M[ind], self.sigma0[ind])
+                return (np.float64(self.sigma_TM[ind][0]),
+                        np.float64(self.tau_M[ind][0]),
+                        np.float64(self.sigma0[ind][0]))
         except:
             print('inputed T not found in the available periods list, try to do interpolation')
             raise ValueError

gmsvtoolkit/models/pynga/CB08.py

@@ -526,7 +526,10 @@ def sd_calc(self):
         sigmaArb = np.sqrt(sigmaT**2 + self.sigma_C[indT]**2)
 
         # standard deviations are in logarithm scale !!!
-        return (sigma, tau, sigmaT, sigmaArb)
+        return (np.float64(sigma[0]),
+                np.float64(tau[0]),
+                np.float64(sigmaT[0]),
+                np.float64(sigmaArb[0]))
 
 def CB08nga_test(T, CoefTerms):
     """

gmsvtoolkit/models/pynga/CY08.py

@@ -525,7 +525,9 @@ def calc_sigma_tau(self):
         sigmaT = np.sqrt(sigma**2 + tauNL**2)
 
         #return (sigma, tau, sigmaT)
-        return (sigma, tauNL, sigmaT)
+        return (np.float64(sigma[0]),
+                np.float64(tauNL[0]),
+                np.float64(sigmaT[0]))
 
 def CY08nga_test(T, CoefTerms):
     """

gmsvtoolkit/plots/plot_gmpe.py

@@ -1,8 +1,8 @@
-#!/bin/env python
+#!/bin/env python3
 """
 BSD 3-Clause License
 
-Copyright (c) 2023, University of Southern California
+Copyright (c) 2026, University of Southern California
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without

gmsvtoolkit/plots/plot_gmpe_gof.py

@@ -1,8 +1,8 @@
-#!/usr/bin/env python
+#!/usr/bin/env python3
 """
 BSD 3-Clause License
 
-Copyright (c) 2023, University of Southern California
+Copyright (c) 2026, University of Southern California
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without
@@ -318,9 +318,13 @@ def run():
     output_dir = os.path.abspath(output_dir)
     station_file = os.path.abspath(args.station_list)
 
-    # Plot GMPE GOF plot
-    plot_title = "Comparison between GMPEs and %s" % (args.comp_label)
+    # Set up plot title
+    if args.plot_title is not None:
+        plot_title = args.plot_title
+    else:
+        plot_title = "Comparison between GMPEs and %s" % (args.comp_label)
 
+    # Plot GMPE GOF plot
     plot_gmpe_gof(station_file, args.gmpe_group, args.comp_label,
                   plot_title, input_dir, output_dir, args.run_prefix)
 

gmsvtoolkit/stats/fas_gof.py

@@ -1,8 +1,8 @@
-#!/usr/bin/env python
+#!/usr/bin/env python3
 """
 BSD 3-Clause License
 
-Copyright (c) 2023, University of Southern California
+Copyright (c) 2026, University of Southern California
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without

gmsvtoolkit/stats/gmpe_gof.py

@@ -1,8 +1,8 @@
-#!/usr/bin/env python
+#!/usr/bin/env python3
 """
 BSD 3-Clause License
 
-Copyright (c) 2023, University of Southern California
+Copyright (c) 2026, University of Southern California
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without