remove trailing whitespace from all files and tabs from Fortran files

Author: bensonr

CHANGELOG.md

@@ -20,7 +20,7 @@ sequential patch number (starting from `01`).
 
 ### Changed
 - The `diag_manager` IO is handled by fms2_io instead of mpp_io. Default behavior assumes that the mpp_io namelist variable is set to *cf_compliant = .true.*
-- The user must specify the diagnostic attributes that they want to write to the output files.  Example: If there is no *units* attribute, then the variable metadata will not include *units*, and it will not automatically add *units = “none”*. 
+- The user must specify the diagnostic attributes that they want to write to the output files.  Example: If there is no *units* attribute, then the variable metadata will not include *units*, and it will not automatically add *units = “none”*.
 - Calls to `register_diag_axis` for an X or Y axis that is shifted from the *CENTER* position need to include the optional argument *domain_position* and should be equal to *EAST* or *NORTH* based on the position relative to the domain. EAST and *NORTH* are exposed through `diag_manager_mod`.
 - Changed the handling of *average_T* and *time_bnds* variables so that they are set to values that are sent in and are not manipulated as was the case in mpp_io.
 - `interpolator`, `xgrid`, `data_override`, and `amip_interp` now call  `fms2_io` routines

amip_interp/amip_interp.F90

@@ -577,8 +577,8 @@ subroutine get_amip_sst (Time, Interp, sst, err_msg, lon_model, lat_model)
 
         if(.not. open_file(fileobj, trim(ncfilename), 'read')) &
              call error_mesg ('get_amip_sst', 'Error in opening file '//trim(ncfilename), FATAL)
- 
-        call get_dimension_size(fileobj, 'TIME', nrecords) 
+
+        call get_dimension_size(fileobj, 'TIME', nrecords)
         if (nrecords < 1) call mpp_error('amip_interp_mod', &
                            'Invalid number of SST records in daily SST data file: '//trim(ncfilename), FATAL)
         allocate(timeval(nrecords), ryr(nrecords), rmo(nrecords), rdy(nrecords))
@@ -1434,7 +1434,7 @@ subroutine read_record (type, Date, Adate, dat)
        Curr_date = Adate
        if (verbose > 2 .and. mpp_pe() == 0)  &
              print *, '....... checking   ', Adate
-       if (Date == Adate) ierr = 0 
+       if (Date == Adate) ierr = 0
        if (yr == 0 .and. mo == Date%month) ierr = 0
        if (ierr == 0) exit
      enddo

axis_utils/axis_utils2.F90

@@ -32,7 +32,7 @@ module axis_utils2_mod
   !<DESCRIPTION>
   !
   ! subroutine get_axis_cart(axis,cart) : Returns X,Y,Z or T cartesian attribute
-  ! subroutine get_axis_bounds(axis,axis_bound,axes) : Return axis_bound either from an array of 
+  ! subroutine get_axis_bounds(axis,axis_bound,axes) : Return axis_bound either from an array of
   !                                                    available axes, or defined based on axis mid-points
   ! function get_axis_modulo : Returns true if axis has the modulo attribute
   ! function get_axis_fold   : Returns is axis is folded at a boundary (non-standard meta-data)
@@ -59,7 +59,7 @@ module axis_utils2_mod
   integer, parameter :: maxatts = 100
   real, parameter    :: epsln= 1.e-10
   real, parameter    :: fp5 = 0.5, f360 = 360.0
-  
+
 ! Include variable "version" to be written to log file.
 #include<file_version.h>
 
@@ -340,13 +340,13 @@ subroutine tranlon(lon, lon_start, istrt)
     istrt=0
     do i=1,len-1
        if (lon(i+1) < lon(i)) then
-          istrt=i+1 
+          istrt=i+1
           exit
        endif
     enddo
 
     if (istrt>1) then ! grid is not monotonic
-       if (abs(lon(len)-lon(1)) < epsln) then 
+       if (abs(lon(len)-lon(1)) < epsln) then
           tmp = cshift(lon(1:len-1),istrt-1)
           lon(1:len-1) = tmp
           lon(len) = lon(1)
@@ -405,15 +405,15 @@ function frac_index (value, array)
     real :: value, frac_index
     real, dimension(:) :: array
     logical keep_going
-    
+
     ia = size(array(:))
 
     do i=2,ia
        if (array(i) < array(i-1)) then
-          unit = stdout() 
+          unit = stdout()
           write (unit,*) '=> Error: "frac_index" array must be monotonically increasing when searching for nearest value to ',&
                               value
-          write (unit,*) '          array(i) < array(i-1) for i=',i 
+          write (unit,*) '          array(i) < array(i-1) for i=',i
           write (unit,*) '          array(i) for i=1..ia follows:'
           do ii=1,ia
              write (unit,*) 'i=',ii, ' array(i)=',array(ii)
@@ -431,7 +431,7 @@ function frac_index (value, array)
        do while (i <= ia .and. keep_going)
           i = i+1
           if (value <= array(i)) then
-             frac_index = float(i-1) + (value-array(i-1))/(array(i)-array(i-1)) 
+             frac_index = float(i-1) + (value-array(i-1))/(array(i)-array(i-1))
              keep_going = .false.
           endif
        enddo
@@ -489,7 +489,7 @@ function nearest_index (value, array)
           unit = stdout()
           write (unit,*) '=> Error: "nearest_index" array must be monotonically increasing &
                          &when searching for nearest value to ',value
-          write (unit,*) '          array(i) < array(i-1) for i=',i 
+          write (unit,*) '          array(i) < array(i-1) for i=',i
           write (unit,*) '          array(i) for i=1..ia follows:'
           do ii=1,ia
              write (unit,*) 'i=',ii, ' array(i)=',array(ii)
@@ -516,7 +516,7 @@ end function nearest_index
 
   !#############################################################################
 
-  subroutine interp_1d_linear(grid1,grid2,data1,data2)  
+  subroutine interp_1d_linear(grid1,grid2,data1,data2)
 
     real, dimension(:),    intent(in) :: grid1, data1, grid2
     real, dimension(:), intent(inout) :: data2
@@ -550,8 +550,8 @@ subroutine interp_1d_linear(grid1,grid2,data1,data2)
              data2(i) = data1(n)
           else
              w = (grid2(i)-grid1(n-1))/(grid1(n)-grid1(n-1))
-             data2(i) = (1.-w)*data1(n-1) + w*data1(n)   
-          endif     
+             data2(i) = (1.-w)*data1(n-1) + w*data1(n)
+          endif
        endif
     enddo
 
@@ -561,7 +561,7 @@ subroutine interp_1d_linear(grid1,grid2,data1,data2)
   end subroutine interp_1d_linear
 
   !###################################################################
-  subroutine interp_1d_cubic_spline(grid1, grid2, data1, data2, yp1, ypn)  
+  subroutine interp_1d_cubic_spline(grid1, grid2, data1, data2, yp1, ypn)
 
     real, dimension(:),    intent(in) :: grid1, grid2, data1
     real, dimension(:), intent(inout) :: data2
@@ -572,7 +572,7 @@ subroutine interp_1d_cubic_spline(grid1, grid2, data1, data2, yp1, ypn)
     integer                           :: n, m, i, k, klo, khi
 
     n = size(grid1(:))
-    m = size(grid2(:))    
+    m = size(grid2(:))
 
     do i=2,n
        if (grid1(i) <= grid1(i-1)) call mpp_error(FATAL, 'grid1 not monotonic')
@@ -622,7 +622,7 @@ subroutine interp_1d_cubic_spline(grid1, grid2, data1, data2, yp1, ypn)
        else
           if(n==1) then
             klo = n
-          else 
+          else
             klo = n -1
           endif
        endif
@@ -637,15 +637,15 @@ end subroutine interp_1d_cubic_spline
 
   !###################################################################
 
-  subroutine interp_1d_1d(grid1,grid2,data1,data2, method, yp1, yp2)  
+  subroutine interp_1d_1d(grid1,grid2,data1,data2, method, yp1, yp2)
 
     real, dimension(:),      intent(in)    :: grid1, data1, grid2
     real, dimension(:),      intent(inout) :: data2
     character(len=*), optional, intent(in) :: method
     real,             optional, intent(in) :: yp1, yp2
 
     real              :: y1, y2
-    character(len=32) :: interp_method    
+    character(len=32) :: interp_method
     integer           :: k2, ks, ke
 
     k2 = size(grid2(:))
@@ -673,7 +673,7 @@ end subroutine interp_1d_1d
   !###################################################################
 
 
-  subroutine interp_1d_2d(grid1,grid2,data1,data2)  
+  subroutine interp_1d_2d(grid1,grid2,data1,data2)
 
     real, dimension(:,:),    intent(in) :: grid1, data1, grid2
     real, dimension(:,:), intent(inout) :: data2
@@ -698,7 +698,7 @@ end subroutine interp_1d_2d
 
   !###################################################################
 
-  subroutine interp_1d_3d(grid1,grid2,data1,data2, method, yp1, yp2)  
+  subroutine interp_1d_3d(grid1,grid2,data1,data2, method, yp1, yp2)
 
     real, dimension(:,:,:),  intent(in)    :: grid1, data1, grid2
     real, dimension(:,:,:),  intent(inout) :: data2
@@ -797,7 +797,7 @@ program test
 
 integer :: n_src = 0
 integer :: n_dst = 0
-real, dimension(MAXSIZE) :: grid_src = 0 
+real, dimension(MAXSIZE) :: grid_src = 0
 real, dimension(MAXSIZE) :: grid_dst = 0
 real, dimension(MAXSIZE) :: data_src = 0
 
@@ -839,7 +839,7 @@ program test
                          371.054289820675, 395.098187506342, 446.150726850039 /)
 
 
-  !---reading namelist 
+  !---reading namelist
 #ifdef INTERNAL_FILE_NML
       read (input_nml_file, test_axis_utils_nml, iostat=io)
       ierr = check_nml_error(io,'test_axis_utils_nml')

column_diagnostics/column_diagnostics.F90

@@ -238,7 +238,7 @@ subroutine initialize_diagnostic_columns     &
       character(len=8)   ::  char
       character(len=32)  ::  filename
       logical            ::  allow_ij_input
-      logical            ::  open_file        
+      logical            ::  open_file
       integer            ::  io
 !--------------------------------------------------------------------
 !    local variables:

diag_manager/diag_manager.F90

@@ -3927,7 +3927,7 @@ SUBROUTINE diag_manager_init(diag_model_subset, time_init, err_msg)
     ALLOCATE(fnum_for_domain(max_files))
     ALLOCATE(pelist(mpp_npes()))
     !> Initialize fnum_for_domain with "dn" which stands for done
-     fnum_for_domain(:) = "dn" 
+     fnum_for_domain(:) = "dn"
     CALL mpp_get_current_pelist(pelist, pelist_name)
 
     ! set the diag_init_time if time_init present.  Otherwise, set it to base_time

diag_manager/diag_output.F90

@@ -169,9 +169,9 @@ SUBROUTINE diag_output_init(file_name, FORMAT, file_title, file_unit,&
        else
           fname_no_tile = trim(file_name)
        endif
-    elseif (lowercase(file_name(len_file_name-4:len_file_name-1)) .eq. "tile") then 
+    elseif (lowercase(file_name(len_file_name-4:len_file_name-1)) .eq. "tile") then
        fname_no_tile = file_name(1:len_file_name-6)
-    elseif (len_file_name < 9) then                        
+    elseif (len_file_name < 9) then
        fname_no_tile = trim(file_name)
     elseif (lowercase(file_name(len_file_name-7:len_file_name-4)) .eq. "tile") then
        fname_no_tile = file_name(1:len_file_name-9)
@@ -190,7 +190,7 @@ SUBROUTINE diag_output_init(file_name, FORMAT, file_title, file_unit,&
 !        trim(fname_no_tile)(len_file_name-3:len_file_name) == ".nc") write (6,*)trim(fname_no_tile)
 !        trim(fname_no_tile)(len(trim(fname_no_tile))-3:len(trim(fname_no_tile))) == ".nc") &
 !        trim(fname_no_tile)(len(trim(fname_no_tile))-3:len(trim(fname_no_tile)))  = "   "
-    
+
 !> Check to make sure that only domain2D or domainUG is used.  If both are not null, then FATAL
     if (domain .NE. NULL_DOMAIN2D .AND. domainU .NE. NULL_DOMAINUG)&
           & CALL error_mesg('diag_output_init', "Domain2D and DomainUG can not be used at the same time in "//&
@@ -380,7 +380,7 @@ SUBROUTINE write_axis_meta_data(file_unit, axes, fileob, time_ops, time_axis_reg
                                    call register_variable_attribute(fptr, axis_name, "positive", "up")
                               case (-1)
                                    call register_variable_attribute(fptr, axis_name, "positive", "down")
-                         end select 
+                         end select
                          call write_data(fptr, axis_name, axis_data(istart:iend) )
                       endif
                     type is (FmsNetcdfFile_t) !< For regional X and Y axes, treat as any other axis
@@ -404,7 +404,7 @@ SUBROUTINE write_axis_meta_data(file_unit, axes, fileob, time_ops, time_axis_reg
                                 "FmsNetcdfFile_t for a X or Y axis, ", FATAL)
                   end select
              endif
-             
+
           ELSE
                select type (fptr)
                     type is (FmsNetcdfDomainFile_t)
@@ -436,7 +436,7 @@ SUBROUTINE write_axis_meta_data(file_unit, axes, fileob, time_ops, time_axis_reg
                               call register_variable_attribute(fptr, axis_name, "positive", "up")
                          case (-1)
                               call register_variable_attribute(fptr, axis_name, "positive", "down")
-                    end select 
+                    end select
                     call write_data(fileob, axis_name, axis_data(istart:iend) )
           END IF
        ELSE
@@ -518,7 +518,7 @@ SUBROUTINE write_axis_meta_data(file_unit, axes, fileob, time_ops, time_axis_reg
                                   call register_variable_attribute(fptr, axis_name, "positive", "up")
                              case (-1)
                                   call register_variable_attribute(fptr, axis_name, "positive", "down")
-                        end select 
+                        end select
                         call write_data(fptr, axis_name, axis_data)
                     endif
                    type is (FmsNetcdfFile_t)
@@ -552,7 +552,7 @@ SUBROUTINE write_axis_meta_data(file_unit, axes, fileob, time_ops, time_axis_reg
 
                         call register_variable_attribute(fptr, axis_name, "long_name", axis_long_name)
                         if(trim(axis_cart_name).ne."N") call register_variable_attribute(fptr, axis_name, "axis",trim(axis_cart_name))
-                        is_time_axis_registered = .true. 
+                        is_time_axis_registered = .true.
                         if (present(time_axis_registered)) time_axis_registered = is_time_axis_registered
                    type is (FmsNetcdfUnstructuredDomainFile_t)
                         call register_axis(fptr, axis_name, size(axis_data) )
@@ -943,7 +943,7 @@ FUNCTION write_field_meta_data ( file_unit, name, axes, units, long_name, range,
                & fill=CMOR_MISSING_VALUE,&
                & pack=ipack, time_method=time_method)
        END IF
-    END IF 
+    END IF
   if (.not. variable_exists(fileob,name)) then
   ! ipack Valid values:
   !        1 = 64bit </LI>
@@ -1144,7 +1144,7 @@ subroutine diag_field_write_field (field, buffer, static, fileob, file_num, file
     TYPE(diag_fieldtype), INTENT(inout) :: Field
     REAL , INTENT(inout) :: buffer(:,:,:,:)
     logical, intent(in), optional :: static
-    class(FmsNetcdfFile_t), optional, intent(inout),target :: fileob 
+    class(FmsNetcdfFile_t), optional, intent(inout),target :: fileob
     class(FmsNetcdfFile_t), pointer :: fptr => null()
     integer, intent(in), optional  :: file_num
     type(FmsNetcdfUnstructuredDomainFile_t),intent(inout), optional :: fileobjU(:)
@@ -1294,10 +1294,10 @@ subroutine diag_write_time (fileob,rtime_value,time_index,time_name)
      endif
 !> Write the time data
      call write_data (fileob, trim(name_time), rtime_value, unlim_dim_level=time_index)
-!> Cleanup     
+!> Cleanup
      if (allocated(name_time)) deallocate(name_time)
      if (associated(fptr)) nullify(fptr)
-  end subroutine diag_write_time 
+  end subroutine diag_write_time
   ! </SUBROUTINE>
 
   ! <FUNCTION NAME="get_axis_index">

diag_manager/diag_util.F90

@@ -2495,7 +2495,7 @@ SUBROUTINE diag_data_out(file, field, dat, time, final_call_in, static_write_in)
     IF ( PRESENT(final_call_in) ) final_call = final_call_in
     static_write = .FALSE.
     IF ( PRESENT(static_write_in) ) static_write = static_write_in
-!> dif is the time as a real that is evaluated 
+!> dif is the time as a real that is evaluated
     dif = get_date_dif(time, base_time, files(file)%time_units)
 
     ! get file_unit, open new file and close curent file if necessary

doxyfile

@@ -820,7 +820,7 @@ RECURSIVE              = YES
 # Note that relative paths are relative to the directory from which doxygen is
 # run.
 
-EXCLUDE                = 
+EXCLUDE                =
 # The EXCLUDE_SYMLINKS tag can be used to select whether or not files or
 # directories that are symbolic links (a Unix file system feature) are excluded
 # from the input.
@@ -852,7 +852,7 @@ EXCLUDE_SYMBOLS        =
 # that contain example code fragments that are included (see the \include
 # command).
 
-EXAMPLE_PATH           = 
+EXAMPLE_PATH           =
 
 # If the value of the EXAMPLE_PATH tag contains directories, you can use the
 # EXAMPLE_PATTERNS tag to specify one or more wildcard pattern (like *.cpp and