From 4a64f87d770ec176c68636f2da706aadb44ed521 Mon Sep 17 00:00:00 2001 From: Lionel GUEZ <guez@lmd.ens.fr> Date: Wed, 7 Jun 2023 23:07:22 +0200 Subject: [PATCH] Do not write `speed.csv` Since we now have speed in the shapefile `cont_list`. --- Inst_eddies/set_max_speed.F90 | 14 ++------------ 1 file changed, 2 insertions(+), 12 deletions(-) diff --git a/Inst_eddies/set_max_speed.F90 b/Inst_eddies/set_max_speed.F90 index 322d7d62..05abace1 100644 --- a/Inst_eddies/set_max_speed.F90 +++ b/Inst_eddies/set_max_speed.F90 @@ -19,7 +19,7 @@ contains IEEE_QUIET_NAN ! Libraries: - use jumble, only: new_unit, rad_to_deg + use jumble, only: rad_to_deg use shapelib, only: shpt_polygon use shapelib_03, only: shp_append_object_03, dbf_write_attribute_03 @@ -70,7 +70,7 @@ contains ! Local: real, allocatable:: speed(:) ! (n_cont) speed on the contour - integer i, i_outer, unit, ishape + integer i, i_outer, ishape integer sort_ssh(size(cont_list)) ! (n_max_cont) ! Sorts cont_list%ssh in ascending order. Defined only for @@ -102,16 +102,6 @@ contains end do #ifdef DEBUG - call new_unit(unit) - open(unit, file = "speed.csv", status = "replace", action = "write") - write(unit, fmt = *) "speed(sort_ssh) cont_list%ssh(sort_ssh)" - - do i = 1, n_cont - write(unit, fmt = *) speed(sort_ssh(i)), cont_list(sort_ssh(i))%ssh - end do - - close(unit) - if (write_cont_list) then do i = 1, n_cont call shp_append_object_03(ishape, hshp_cont_list, shpt_polygon, & -- GitLab