diff --git a/Inst_eddies/config.f90 b/Inst_eddies/config.f90 index 216360ac8c4ec671644a4446253c0836225d9cab..1bf4432a0604ddfb183ad41e2c2e0679b0da6274 100644 --- a/Inst_eddies/config.f90 +++ b/Inst_eddies/config.f90 @@ -30,13 +30,14 @@ contains !---------------------------------------------------------------------- call new_unit(unit) - + if (exist) then open(unit, file = "SHPC_cyclo/used_nml.txt", status = "old", & action = "read", position = "rewind") read(unit, nml = config_nml) close(unit) else + print *, "config: Reading config_nml.txt..." open(unit, file = "config_nml.txt", status = "old", action = "read", & position = "rewind") read(unit, nml = config_nml) diff --git a/Inst_eddies/input_ssh.f90 b/Inst_eddies/input_ssh.f90 index 7c2234bc6ed8bc5f391ca00cff13da2c264403fc..9a5177e2c5aafeba65210a2bb256559829cd019d 100644 --- a/Inst_eddies/input_ssh.f90 +++ b/Inst_eddies/input_ssh.f90 @@ -76,9 +76,9 @@ contains ! As we are requiring the grid spacing to be uniform, the value of ! "periodic" may be deduced from the values of step(1) and nlon: periodic = nint(twopi / step(1)) == nlon - print *, "inst_eddies: periodic = ", periodic + print *, "input_ssh: periodic = ", periodic - call assert(2 * max_radius_deg(1) < 180., "inst_eddies max_radius_deg") + call assert(2 * max_radius_deg(1) < 180., "input_ssh: max_radius_deg") ! (Let us require this even if not periodic. This is clearer.) max_radius = nint(max_radius_deg / step_deg) diff --git a/Inst_eddies/inst_eddies.py.in b/Inst_eddies/inst_eddies.py.in index 8cc6cc9061e204e571c4147d509f5ea82f06e368..0383e1c2c64c0b8c02c4fb3797668d4589ed5d6f 100644 --- a/Inst_eddies/inst_eddies.py.in +++ b/Inst_eddies/inst_eddies.py.in @@ -28,9 +28,9 @@ def loop_inst_eddies(files, bbox, d): sys.exit(inst_eddies_exe + " not found or not executable") if os.access("config_nml.txt", os.R_OK): - print("Will use config_nml.txt.") + print("inst_eddies.loop_inst_eddies: Will use config_nml.txt.") else: - sys.exit('inst_eddies.py: "config_nml.txt" not found') + sys.exit('"config_nml.txt" not found') nco_instance = nco.Nco() @@ -45,8 +45,9 @@ def loop_inst_eddies(files, bbox, d): writer.writerow(["date", "elapsed_NCO", "elapsed_Fortran"]) for nc_file in files: - print("inst_eddies.py: Date:", d) - print("Input file:", nc_file) + print("inst_eddies.loop_inst_eddies: Date:", d) + print("inst_eddies.loop_inst_eddies: Input file:", nc_file, + flush = True) if os.access(nc_file, os.F_OK): my_pc = time.perf_counter() @@ -73,7 +74,7 @@ def loop_inst_eddies(files, bbox, d): os.remove("uv.nc") writer.writerow([d, elapsed_NCO, elapsed_Fortran]) else: - print("inst_eddies.py: Missing file:", nc_file) + print("inst_eddies.loop_inst_eddies: Missing file:", nc_file) for orient in ["cyclo", "anti"]: with open(f"SHPC_{orient}/ishape_last.txt", "a") \ @@ -98,6 +99,6 @@ if __name__ == "__main__": argparser.add_argument("-b", "--bbox", nargs=4, type = float, metavar = ("xmin", "xmax", "ymin", "ymax")) args = argparser.parse_args() - d = input("Enter first date index (integer value):") + d = input("Enter first date index (integer value): ") d = int(d) loop_inst_eddies(args.file, args.bbox, d) diff --git a/Overlap/Tests/save_snapshot.f90 b/Overlap/Tests/save_snapshot.f90 index 6144cb2dbe60552bdcdb29278a81184ebcc829c2..e3d206b621b4447d7248eba69a43514645e0bf12 100644 --- a/Overlap/Tests/save_snapshot.f90 +++ b/Overlap/Tests/save_snapshot.f90 @@ -16,7 +16,7 @@ contains use write_extr_map_m, only: write_extr_map type(snapshot), intent(in):: s - + real, intent(in):: corner(:) ! (2) ! longitude and latitude of the corner of the whole grid in input ! NetCDF, in degrees diff --git a/Overlap/Tests/test_overlap.f90 b/Overlap/Tests/test_overlap.f90 index 09eae3c16f7dd658e3912675ca8579714817d712..7fef4fbd5be22ec8c2a1da533fa145c4744c9e74 100644 --- a/Overlap/Tests/test_overlap.f90 +++ b/Overlap/Tests/test_overlap.f90 @@ -3,8 +3,8 @@ program test_overlap use, intrinsic:: ISO_FORTRAN_ENV ! Libraries: - use jumble, only: get_command_arg_dyn, read_opcol, new_unit, ediff1d - use jumble, only: deg_to_rad, assert + use jumble, only: get_command_arg_dyn, read_opcol, new_unit, ediff1d, & + deg_to_rad, assert use shapelib_03, only: dbf_read_attribute_03 use derived_types, only: snapshot, shpc_slice_handler, shpc_slice_meta