change cut_off_T_cells to min_dy

README.md

@@ -137,10 +137,10 @@ double-precision topography file.
 Options
   * `--vgrid <vgrid>`  vertical grid (default 'ocean_vgrid.nc')
 
-## cut_off_T_cells
+## min_dy
 
 ```
-usage: topogtools cut_off_T_cells --input <input_file> --output <output_file> --hgrid <grid> --cutoff <cutoff_value>       
+usage: topogtools min_dy --input <input_file> --output <output_file> --hgrid <grid> --cutoff <cutoff_value>       
 ```
 
 Cut off T cells with size smaller than <cutoff_value>. Cut off should be in kilometers

src/topography.f90

@@ -35,7 +35,7 @@ module topography
     procedure :: nonadvective => topography_nonadvective
     procedure :: min_max_depth => topography_min_max_depth
     procedure :: mask => topography_mask
-    procedure :: cut_off_T_cells => topography_cut_off_T_cells
+    procedure :: min_dy => topography_min_dy
   end type topography_t
 
   interface topography_t
@@ -465,7 +465,7 @@ subroutine topography_min_max_depth(this, vgrid_file, vgrid_type, level)
 
   end subroutine topography_min_max_depth
 
- subroutine topography_cut_off_T_cells(this, hgrid, cutoff)
+ subroutine topography_min_dy(this, hgrid, cutoff)
     class(topography_t), intent(inout) :: this
     character(len=*), intent(in) :: hgrid
     real(real64), intent(in) :: cutoff
@@ -502,7 +502,7 @@ subroutine topography_cut_off_T_cells(this, hgrid, cutoff)
         end do
     end do   
 
-end subroutine topography_cut_off_T_cells
+end subroutine topography_min_dy
 
   !-------------------------------------------------------------------------
   subroutine topography_fill_fraction(this, sea_area_fraction)

src/topogtools.f90

@@ -144,7 +144,7 @@ program topogtools
       help_check_nonadvective, version_text)
   case ('mask')
     call set_args('--input:i "unset" --output:o "unset"', help_mask, version_text)
-  case ('cut_off_T_cells')
+  case ('min_dy')
     call set_args('--input:i "unset" --output:o "unset" --hgrid "ocean_hgrid.nc" --cutoff 0.0', help_cutoff, version_text)
   case ('')
     ! Print help in case the user specified the --help flag
@@ -226,7 +226,7 @@ program topogtools
     topog = topography_t(file_in)
     call topog%mask(file_out)
 
-  case ('cut_off_T_cells')
+  case ('min_dy')
     hgrid = sget('hgrid')
     call check_file_exist(hgrid)
     cutoff = rget('cutoff')
@@ -240,7 +240,7 @@ program topogtools
       error stop
     end if  
     topog = topography_t(file_in)
-    call topog%cut_off_T_cells(hgrid, cutoff)
+    call topog%min_dy(hgrid, cutoff)
     call topog%update_history(get_mycommand())
     call topog%write(file_out)