From 84175ef48d3f5d986d889ae3c4c761235e71a6dc Mon Sep 17 00:00:00 2001
From: Thomas Gastine <gastine@mps.mpg.de>
Date: Fri, 13 Sep 2024 15:27:41 +0200
Subject: [PATCH] remove trailing white spaces

---
 src/LMmapping.f90          |   8 +--
 src/Makefile               |  16 +++---
 src/README.md              |   4 +-
 src/dirk_schemes.f90       |   2 +-
 src/dt_fieldsLast.f90      |   2 +-
 src/multistep_schemes.f90  |   6 +-
 src/nonlinear_bcs.f90      |   2 +-
 src/parallel_solver.f90    |   6 +-
 src/radial_data.f90        |  10 ++--
 src/radial_derivatives.f90 | 114 ++++++++++++++++++-------------------
 src/storeCheckPoints.f90   |   6 +-
 src/truncation.f90         |  14 ++---
 12 files changed, 95 insertions(+), 95 deletions(-)

diff --git a/src/LMmapping.f90 b/src/LMmapping.f90
index 0599b9cc..26b5d598 100644
--- a/src/LMmapping.f90
+++ b/src/LMmapping.f90
@@ -6,15 +6,15 @@ module LMmapping
    use parallel_mod, only: load
 
    implicit none
- 
+
    private
- 
+
    type, public :: mappings
       integer :: l_max, m_max, lm_max ,m_min
       integer, allocatable :: lm2(:,:),lm2l(:),lm2m(:)
       integer, allocatable :: lm2lmS(:), lm2lmA(:)
    end type mappings
- 
+
    type, public :: subblocks_mappings
       integer :: nLMBs,l_max,m_max,sizeLMB2max,m_min
       integer, allocatable :: nLMBs2(:)
@@ -23,7 +23,7 @@ module LMmapping
       integer, allocatable :: lm22l(:,:,:)
       integer, allocatable :: lm22m(:,:,:)
    end type subblocks_mappings
- 
+
    public :: allocate_mappings, deallocate_mappings, &
    &         allocate_subblocks_mappings, deallocate_subblocks_mappings
 
diff --git a/src/Makefile b/src/Makefile
index 06fe08a8..12aaa485 100644
--- a/src/Makefile
+++ b/src/Makefile
@@ -76,12 +76,12 @@ ifeq ($(COMPILER),intel)
 	else
 		COMP_MPFC = mpiifort
 		COMP_MPCC = mpiicc
-	endif	
+	endif
 
 	FFLAG_MKL_LIBS = -L$(MKLROOT)/lib/intel64 -lmkl_intel_lp64 \
                          -lmkl_lapack95_lp64\
 			 -lmkl_sequential -lmkl_core -lpthread -lm \
-			 -Wl,-rpath,$(MKLROOT)/lib/intel64 
+			 -Wl,-rpath,$(MKLROOT)/lib/intel64
 
 endif
 ifeq ($(COMPILER),gnu)
@@ -94,7 +94,7 @@ ifeq ($(COMPILER),gnu)
 	FFLAG_PREPROCESSOR = -cpp
 	FFLAG_WARN_ALL = -Wall -Wno-maybe-uninitialized
 	FFLAG_CHECK_ALL = -fcheck=all
-	FFLAG_FPE = -ffpe-trap=invalid,zero,overflow 
+	FFLAG_FPE = -ffpe-trap=invalid,zero,overflow
 	COMP_FC = gfortran
 	COMP_CC = gcc
 	COMP_MPFC = mpif90
@@ -212,7 +212,7 @@ ifeq ($(PRODRUN),yes)
 	DEBUG=no
 	OPT += -O3 $(FFLAG_ARCH_OPT)
 else
-	OPT += -O 
+	OPT += -O
 endif
 ifeq ($(USE_OMP),yes)
 	FFLAGS += $(FFLAG_OPENMP)
@@ -238,7 +238,7 @@ endif
 
 INCPATHS += -I.
 
-FFLAGS += $(FFLAG_CONVERT) 
+FFLAGS += $(FFLAG_CONVERT)
 FFLAGS  += $(FFLAG_PREPROCESSOR) $(OPT)
 LDFLAGS =
 
@@ -296,7 +296,7 @@ endif
 OBJS   := $(addsuffix .o, $(basename $(RED_SOURCES)))
 OBJS += truncation.o
 
-.SUFFIXES: 
+.SUFFIXES:
 
 #ifeq ($(COMPILER),gnu)
 #ifeq ($(USE_MPI),yes)
@@ -513,7 +513,7 @@ preCalculations.o     : truncation.o phys_param.o num_param.o constants.o\
                         useful.o integration.o parallel.o special.o
 RMS_helpers.o         : truncation.o blocking.o radial.o constants.o\
                         horizontal.o  useful.o integration.o\
-                        LMmapping.o 
+                        LMmapping.o
 kinetic_energy.o      : truncation.o blocking.o num_param.o radial.o\
                         horizontal.o output_data.o logic.o constants.o\
                         phys_param.o useful.o integration.o\
@@ -539,7 +539,7 @@ RMS.o	              : truncation.o blocking.o mem_alloc.o
 init_fields.o         : truncation.o blocking.o radial.o horizontal.o\
                         constants.o logic.o $(FFT_OBJS) $(SHT_OBJS)\
                         useful.o phys_param.o mpi_transpose.o\
-                        $(DCT_OBJS) mem_alloc.o parallel.o 
+                        $(DCT_OBJS) mem_alloc.o parallel.o
 mean_sd.o             : precision_mod.o mem_alloc.o matrices.o communications.o
 movie.o               : truncation.o parallel.o radial_data.o\
                         output_data.o logic.o radial.o mem_alloc.o\
diff --git a/src/README.md b/src/README.md
index 493866c7..9e6b04ff 100644
--- a/src/README.md
+++ b/src/README.md
@@ -5,7 +5,7 @@
 Set a suitable compiler in the first line of the Makefile: `COMPILER = <compiler_phrase>`. The options are `intel`, `gnu` or `portland` - depending on your available compilers.
 
 **Table** : List of default compilers
- 
+
 | Compiler Option |    Normal     |     With MPI     |
 |:---------------:|:-------------:|:----------------:|
 |      intel      |   ifort, icc  | mpiifort, mpiicc |
@@ -17,7 +17,7 @@ Set a suitable compiler in the first line of the Makefile: `COMPILER = <compiler
 * `PRODRUN` Set it to `yes` for production run, `no` for debugging.
 * `PRECISION` Set it to 'dble' for double-precision calculations or to 'sngl' for single-precision calculations
 * `OUT_PREC` Set it to 'dble' for double-precision in binary outputs or to 'sngl' for single precision
-* `DEBUG`   Set to `yes` to run in debugging mode. *While running in debugging mode, set* `PRODRUN` *to* `no`. The debug mode with intel compilers uses `marmot90`. 
+* `DEBUG`   Set to `yes` to run in debugging mode. *While running in debugging mode, set* `PRODRUN` *to* `no`. The debug mode with intel compilers uses `marmot90`.
 * `USE_MPI` Set to `yes` to use MPI
 * `USE_OMP`  Set to `yes` to use openmp (cannot work without MPI)
 * `USE_PRECOND` Set to `yes` to perform some pre-conditioning  of the matrices
diff --git a/src/dirk_schemes.f90 b/src/dirk_schemes.f90
index 6035dcaf..e37c63c6 100644
--- a/src/dirk_schemes.f90
+++ b/src/dirk_schemes.f90
@@ -784,7 +784,7 @@ subroutine set_weights(this, lMatNext)
 
       this%butcher_ass_imp(:) = this%dt(1)*this%butcher_ass_imp(:)
       this%butcher_ass_exp(:) = this%dt(1)*this%butcher_ass_exp(:)
-         
+
    end subroutine set_weights
 !------------------------------------------------------------------------------
    subroutine set_dt_array(this, dt_new, dt_min, time, n_log_file,  &
diff --git a/src/dt_fieldsLast.f90 b/src/dt_fieldsLast.f90
index 2d5085f7..68266011 100644
--- a/src/dt_fieldsLast.f90
+++ b/src/dt_fieldsLast.f90
@@ -151,7 +151,7 @@ subroutine initialize_fieldsLast(nold, nexp, nimp)
          dVSrLM_Rloc(:,:) =zero
          dVxBhLM_Rloc(:,:)=zero
          bytes_allocated = bytes_allocated+                               &
-         &                 6*lm_max*(nRstop-nRstart+1)*SIZEOF_DEF_COMPLEX+& 
+         &                 6*lm_max*(nRstop-nRstart+1)*SIZEOF_DEF_COMPLEX+&
          &                 3*lm_maxMag*(nRstopMag-nRstartMag+1)*SIZEOF_DEF_COMPLEX
       else
          if ( l_double_curl ) then
diff --git a/src/multistep_schemes.f90 b/src/multistep_schemes.f90
index 520f3a43..4495206d 100644
--- a/src/multistep_schemes.f90
+++ b/src/multistep_schemes.f90
@@ -190,13 +190,13 @@ subroutine set_weights(this, lMatNext)
       !-- Local variables
       real(cp) :: delta, delta_n, delta_n_1, delta_n_2
       real(cp) :: a0, a1, a2, a3, a4, b0, b1, b2, b3, c0, c1, c2, c3
-      real(cp) :: gam, theta, c 
+      real(cp) :: gam, theta, c
       real(cp) :: wimp_old
 
       wimp_old = this%wimp_lin(1)
 
       select case ( this%time_scheme )
-         case ('CNAB2') 
+         case ('CNAB2')
             this%wimp(1)    =one
             this%wimp_lin(1)=alpha*this%dt(1)
             this%wimp_lin(2)=(1-alpha)*this%dt(1)
@@ -213,7 +213,7 @@ subroutine set_weights(this, lMatNext)
 
             this%wexp(1)=(one+delta)*this%dt(1)
             this%wexp(2)=0.0_cp
-         case ('MODCNAB') 
+         case ('MODCNAB')
             delta = this%dt(1)/this%dt(2)
             this%wimp(1)    =one
             this%wimp(2)    =0.0_cp
diff --git a/src/nonlinear_bcs.f90 b/src/nonlinear_bcs.f90
index be843926..598714a2 100644
--- a/src/nonlinear_bcs.f90
+++ b/src/nonlinear_bcs.f90
@@ -90,7 +90,7 @@ subroutine get_b_nl_bcs(bc,br_vt_lm,br_vp_lm,b_nl_bc,aj_nl_bc)
 
       if ( bc == 'CMB' ) then
 
-         b_nl_bc(1) =zero                                             
+         b_nl_bc(1) =zero
          aj_nl_bc(1)=zero
          fac=conductance_ma*prmag
          !$omp parallel do default(shared)
diff --git a/src/parallel_solver.f90 b/src/parallel_solver.f90
index 77ca5085..318634e9 100644
--- a/src/parallel_solver.f90
+++ b/src/parallel_solver.f90
@@ -17,10 +17,10 @@ module parallel_solvers
 
    private
 
-   type, public :: type_tri_par   
+   type, public :: type_tri_par
       integer :: nRMin
       integer :: nRMax
-      integer :: lMin 
+      integer :: lMin
       integer :: lMax
       real(cp), allocatable :: low(:,:)
       real(cp), allocatable :: diag(:,:)
@@ -578,7 +578,7 @@ subroutine solver_finish_3(this, x, lms_block, nlm_block, nRstart, nRstop, tag,
       integer,     intent(inout) :: array_req(:)
       complex(cp), intent(inout) :: x(1:lm_max, nRstart-1:nRstop+1)
       integer,     intent(inout) :: req
-      
+
       !-- Local variables
       integer :: lmb, lmu
 
diff --git a/src/radial_data.f90 b/src/radial_data.f90
index 1752bd62..9f302e1a 100644
--- a/src/radial_data.f90
+++ b/src/radial_data.f90
@@ -6,16 +6,16 @@ module radial_data
    use iso_fortran_env, only: output_unit
    use parallel_mod, only: rank, n_procs, nR_per_rank, load, getBlocks
    use logic, only: l_mag, lVerbose, l_finite_diff
- 
+
    implicit none
- 
+
    private
 
    type(load), public, allocatable :: radial_balance(:)
- 
+
    integer, protected, public :: nRstart,nRstop,nRstartMag,nRstopMag
    integer, protected, public :: n_r_cmb,n_r_icb
- 
+
    public :: initialize_radial_data, finalize_radial_data
 
 contains
@@ -31,7 +31,7 @@ subroutine initialize_radial_data(n_r_max)
       n_r_icb=n_r_max
 
       allocate(radial_balance(0:n_procs-1))
-      call getBlocks(radial_balance, n_r_max, n_procs)   
+      call getBlocks(radial_balance, n_r_max, n_procs)
 
       nRstart = radial_balance(rank)%nStart
       nRstop = radial_balance(rank)%nStop
diff --git a/src/radial_derivatives.f90 b/src/radial_derivatives.f90
index 9d70d26a..9e1df886 100644
--- a/src/radial_derivatives.f90
+++ b/src/radial_derivatives.f90
@@ -66,9 +66,9 @@ end subroutine finalize_der_arrays
 !------------------------------------------------------------------------------
    subroutine get_dcheb_complex(f,df,n_f_max,n_f_start,n_f_stop,n_r_max,n_cheb_max)
       !
-      !  Returns Chebyshev coeffitients of first derivative df and second  
-      !  derivative ddf for a function whose cheb-coeff. are given as     
-      !  columns in array f(n_f_max,n_r_max).                             
+      !  Returns Chebyshev coeffitients of first derivative df and second
+      !  derivative ddf for a function whose cheb-coeff. are given as
+      !  columns in array f(n_f_max,n_r_max).
       !
 
       !-- Input variables:
@@ -153,11 +153,11 @@ end subroutine get_dcheb_real_1d
 !------------------------------------------------------------------------------
    subroutine get_ddcheb(f,df,ddf,n_f_max,n_f_start,n_f_stop,n_r_max,n_cheb_max)
       !
-      !  Returns Chebyshev coefficients of first derivative df and second  
-      !  derivative ddf for a function whose cheb-coeff. are given as     
-      !  columns in array f(n_c_tot,n_r_max).                             
+      !  Returns Chebyshev coefficients of first derivative df and second
+      !  derivative ddf for a function whose cheb-coeff. are given as
+      !  columns in array f(n_c_tot,n_r_max).
       !
-    
+
       !-- Input variables:
       integer,     intent(in) :: n_f_start  ! No of column to start with
       integer,     intent(in) :: n_f_stop   ! No of column to stop with
@@ -165,15 +165,15 @@ subroutine get_ddcheb(f,df,ddf,n_f_max,n_f_start,n_f_stop,n_r_max,n_cheb_max)
       integer,     intent(in) :: n_r_max    ! second dimension of f,df,ddf
       integer,     intent(in) :: n_cheb_max ! Number of cheb modes
       complex(cp), intent(in) :: f(n_f_max,n_r_max)
-    
+
       !-- Output variables:
       complex(cp), intent(out) ::  df(n_f_max,n_r_max)
       complex(cp), intent(out) ::  ddf(n_f_max,n_r_max)
-    
+
       !-- local variables:
       integer :: n_f,n_cheb
       real(cp) :: fac_cheb
-    
+
       !----- initialize derivatives:
       do n_cheb=n_cheb_max,n_r_max
          do n_f=n_f_start,n_f_stop
@@ -193,7 +193,7 @@ subroutine get_ddcheb(f,df,ddf,n_f_max,n_f_start,n_f_stop,n_r_max,n_cheb_max)
          df(n_f,n_cheb) =fac_cheb*f(n_f,n_cheb+1)
          ddf(n_f,n_cheb)=zero
       end do
-    
+
       !----- recursion
       do n_cheb=n_cheb_max-2,1,-1
          fac_cheb=real(2*n_cheb,kind=cp)
@@ -208,11 +208,11 @@ end subroutine get_ddcheb
    subroutine get_dddcheb(f,df,ddf,dddf,n_f_max,n_f_start,n_f_stop,n_r_max, &
               &           n_cheb_max)
       !
-      !  Returns chebychev coeffitiens of first derivative df and second  
-      !  derivative ddf for a function whose cheb-coeff. are given as     
-      !  columns in array f(n_c_tot,n_r_max).                             
+      !  Returns chebychev coeffitiens of first derivative df and second
+      !  derivative ddf for a function whose cheb-coeff. are given as
+      !  columns in array f(n_c_tot,n_r_max).
       !
-         
+
       !-- Input variables:
       integer,     intent(in) :: n_f_start  ! No of column to start with
       integer,     intent(in) :: n_f_stop   ! No of column to stop with
@@ -265,34 +265,34 @@ subroutine get_dddcheb(f,df,ddf,dddf,n_f_max,n_f_start,n_f_stop,n_r_max, &
    end subroutine get_dddcheb
 !---------------------------------------------------------------------------
    subroutine get_dr_real_1d(f,df,n_r_max,r_scheme)
- 
+
       !-- Input variables:
       integer,             intent(in) :: n_r_max    ! number of radial grid points
       real(cp),            intent(in) :: f(n_r_max)
       class(type_rscheme), intent(in) :: r_scheme
-    
+
       !-- Output variables:
       real(cp), intent(out) :: df(n_r_max)   ! first derivative of f
-    
+
       !-- Local:
       integer :: n_r, od
-    
+
       select type (r_scheme)
 
       type is(type_cheb_odd)
 
          !-- Copy input functions:
          work_1d_real(1:n_r_max)=f(:)
-    
+
          !-- Transform f to cheb space:
          call r_scheme%costf1(work_1d_real)
-    
+
          !-- Get derivatives:
          call get_dcheb(work_1d_real,df,n_r_max,r_scheme%n_max)
-    
+
          !-- Transform back:
          call r_scheme%costf1(df)
-    
+
          !-- New map:
          df(:)=r_scheme%drx(:)*df(:)
 
@@ -324,13 +324,13 @@ end subroutine get_dr_real_1d
    subroutine get_dr_complex(f,df,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme, &
               &              nocopy,l_dct_in)
       !
-      !  Returns first radial derivative df of the input function f.      
-      !  Array f(n_f_max,*) may contain several functions numbered by     
-      !  the first index. The subroutine calculates the derivaties of     
-      !  the functions f(n_f_start,*) to f(n_f_stop) by transforming      
-      !  to a Chebychev representation using n_r_max radial grid points . 
+      !  Returns first radial derivative df of the input function f.
+      !  Array f(n_f_max,*) may contain several functions numbered by
+      !  the first index. The subroutine calculates the derivaties of
+      !  the functions f(n_f_start,*) to f(n_f_stop) by transforming
+      !  to a Chebychev representation using n_r_max radial grid points .
       !
-    
+
       !-- Input variables:
       integer,             intent(in) :: n_r_max          ! number of radial grid points
       integer,             intent(in) :: n_f_max          ! first dim of f
@@ -340,14 +340,14 @@ subroutine get_dr_complex(f,df,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme, &
       class(type_rscheme), intent(in) :: r_scheme
       logical, optional,   intent(in) :: nocopy
       logical, optional,   intent(in) :: l_dct_in
-    
+
       !-- Output variables:
       complex(cp), intent(out) :: df(n_f_max,n_r_max)   ! first derivative of f
-    
+
       !-- Local:
       integer :: n_r,n_f,od
       logical :: copy_array, l_dct_in_loc
-    
+
       select type (r_scheme)
 
       type is(type_cheb_odd)
@@ -376,16 +376,16 @@ subroutine get_dr_complex(f,df,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme, &
                   work(n_f,n_r)=f(n_f,n_r)
                end do
             end do
-       
+
             !-- Transform f to cheb space:
             if ( l_dct_in_loc ) then
                call r_scheme%costf1(work,n_f_max,n_f_start,n_f_stop)
             end if
-          
+
             !-- Get derivatives:
             call get_dcheb(work,df,n_f_max,n_f_start,n_f_stop,n_r_max, &
                  &         r_scheme%n_max)
-          
+
             !-- Transform back:
             call r_scheme%costf1(df,n_f_max,n_f_start,n_f_stop)
 
@@ -393,18 +393,18 @@ subroutine get_dr_complex(f,df,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme, &
 
             !-- Transform f to cheb space:
             if ( l_dct_in_loc ) call r_scheme%costf1(f,n_f_max,n_f_start,n_f_stop)
-          
+
             !-- Get derivatives:
             call get_dcheb(f,df,n_f_max,n_f_start,n_f_stop,n_r_max, &
                  &         r_scheme%n_max)
-          
+
             !-- Transform back:
             if ( l_dct_in_loc ) call r_scheme%costf1(f,n_f_max,n_f_start,n_f_stop)
             call r_scheme%costf1(df,n_f_max,n_f_start,n_f_stop)
 
          end if
 #endif
-       
+
          !-- New map:
          do n_r=1,n_r_max
             do n_f=n_f_start,n_f_stop
@@ -448,14 +448,14 @@ end subroutine get_dr_complex
    subroutine get_ddr(f,df,ddf,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme, &
               &       l_dct_in)
       !
-      !  Returns first radial derivative df and second radial             
-      !  derivative ddf of the input function f.                          
-      !  Array f(n_f_max,*) may contain several functions numbered by     
-      !  the first index. The subroutine calculates the derivatives of    
-      !  the functions f(n_f_start,*) to f(n_f_stop) by transforming      
-      !  to a Chebychev representation using n_r_max radial grid points.  
+      !  Returns first radial derivative df and second radial
+      !  derivative ddf of the input function f.
+      !  Array f(n_f_max,*) may contain several functions numbered by
+      !  the first index. The subroutine calculates the derivatives of
+      !  the functions f(n_f_start,*) to f(n_f_stop) by transforming
+      !  to a Chebychev representation using n_r_max radial grid points.
       !
-    
+
       !-- Input variables:
       integer,             intent(in) :: n_r_max       ! number of radial grid points
       integer,             intent(in) :: n_f_max       ! first dim of f
@@ -464,11 +464,11 @@ subroutine get_ddr(f,df,ddf,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme, &
       integer,             intent(in) :: n_f_stop      ! last function to be treated
       class(type_rscheme), intent(in) :: r_scheme
       logical, optional,   intent(in) :: l_dct_in
-    
+
       !-- Output variables:
       complex(cp), intent(out) :: df(n_f_max,n_r_max)   ! first derivative of f
       complex(cp), intent(out) :: ddf(n_f_max,n_r_max)  ! second derivative of f
-    
+
       !-- Local variables:
       integer :: n_r,n_f,od
       logical :: l_dct_in_loc
@@ -482,7 +482,7 @@ subroutine get_ddr(f,df,ddf,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme, &
          else
             l_dct_in_loc=.true.
          end if
-    
+
 #ifdef USE_DCT_FFT
          call r_scheme%chebt_oc%get_ddr_fft(f,df,ddf,r_scheme%x_cheb,n_f_max, &
               &                             n_f_start,n_f_stop,r_scheme%n_max,&
@@ -499,16 +499,16 @@ subroutine get_ddr(f,df,ddf,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme, &
          if ( l_dct_in_loc ) then
             call r_scheme%costf1(work,n_f_max,n_f_start,n_f_stop)
          end if
-    
+
          !-- Get derivatives:
          call get_ddcheb(work,df,ddf,n_f_max,n_f_start,n_f_stop,n_r_max, &
               &          r_scheme%n_max)
-    
+
          !-- Transform back:
          call r_scheme%costf1(df,n_f_max,n_f_start,n_f_stop)
          call r_scheme%costf1(ddf,n_f_max,n_f_start,n_f_stop)
 #endif
-    
+
          !-- New map:
          do n_r=1,n_r_max
             do n_f=n_f_start,n_f_stop
@@ -568,11 +568,11 @@ subroutine get_dddr(f,df,ddf,dddf,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme,&
               &        l_dct_in)
       !
       !  Returns first radial derivative df, the second radial deriv. ddf,
-      !  and the third radial derivative dddf of the input function f.    
-      !  Array f(n_f_max,*) may contain several functions numbered by     
-      !  the first index. The subroutine calculates the derivaties of     
-      !  the functions f(n_f_start,*) to f(n_f_stop) by transforming      
-      !  to a Chebychev representation using n_r_max radial grid points.  
+      !  and the third radial derivative dddf of the input function f.
+      !  Array f(n_f_max,*) may contain several functions numbered by
+      !  the first index. The subroutine calculates the derivaties of
+      !  the functions f(n_f_start,*) to f(n_f_stop) by transforming
+      !  to a Chebychev representation using n_r_max radial grid points.
       !
 
       !-- Input variables:
@@ -583,7 +583,7 @@ subroutine get_dddr(f,df,ddf,dddf,n_f_max,n_f_start,n_f_stop,n_r_max,r_scheme,&
       integer,             intent(in) :: n_f_stop        ! last function to be treated
       class(type_rscheme), intent(in) :: r_scheme
       logical, optional,   intent(in) :: l_dct_in
-    
+
       !-- Output variables:
       complex(cp), intent(out) :: df(n_f_max,n_r_max)    ! first derivative of f
       complex(cp), intent(out) :: ddf(n_f_max,n_r_max)   ! second derivative of f
diff --git a/src/storeCheckPoints.f90 b/src/storeCheckPoints.f90
index 4704ae23..03e98305 100644
--- a/src/storeCheckPoints.f90
+++ b/src/storeCheckPoints.f90
@@ -87,7 +87,7 @@ subroutine store(time,tscheme,n_time_step,l_stop_time,l_new_rst_file,
       character(len=72) :: string,rst_file
 
       version = 5
-      l_press_store = ( .not. l_double_curl ) 
+      l_press_store = ( .not. l_double_curl )
 
       if ( l_ave_file ) then
          rst_file="checkpoint_ave."//tag
@@ -678,7 +678,7 @@ subroutine store_mpi(time,tscheme,n_time_step,l_stop_time,l_new_rst_file,  &
 
       end if
 
-      !-- Close file 
+      !-- Close file
       call MPI_Info_free(info, ierr)
       call MPI_File_close(fh, ierr)
 
@@ -726,7 +726,7 @@ subroutine write_one_field_mpi(fh, info, datatype, tscheme, w, dwdt, &
 
       !-- Output variables
       integer(lip),        intent(inout) :: disp
-      complex(cp),         intent(inout) :: work(llm:ulm, n_r_max) 
+      complex(cp),         intent(inout) :: work(llm:ulm, n_r_max)
 
       !-- Local variables
       integer :: n_o
diff --git a/src/truncation.f90 b/src/truncation.f90
index 517bef46..6b2f755e 100644
--- a/src/truncation.f90
+++ b/src/truncation.f90
@@ -1,6 +1,6 @@
 module truncation
    !
-   ! This module defines the grid points and the truncation 
+   ! This module defines the grid points and the truncation
    !
 
    use precision_mod, only: cp
@@ -14,7 +14,7 @@ module truncation
    integer :: n_phi_tot     ! number of longitude grid points
    integer :: n_r_ic_max    ! number of grid points in inner core
    integer :: n_cheb_ic_max ! number of chebs in inner core
-   integer :: minc          ! basic wavenumber, longitude symmetry  
+   integer :: minc          ! basic wavenumber, longitude symmetry
    integer :: nalias        ! controls dealiasing in latitude
    logical :: l_axi         ! logical for axisymmetric calculations
    character(len=72) :: radial_scheme ! radial scheme (either Cheybev of FD)
@@ -23,7 +23,7 @@ module truncation
    real(cp) :: rcut_l        ! Cut-off radius when degree varies with radius
    integer :: fd_order       ! Finite difference order (for now only 2 and 4 are safe)
    integer :: fd_order_bound ! Finite difference order on the  boundaries
- 
+
    !-- Derived quantities:
    integer :: n_phi_max   ! absolute number of phi grid-points
    integer :: n_theta_max ! number of theta grid-points
@@ -35,7 +35,7 @@ module truncation
    integer :: n_m_max     ! max number of ms (different oders)
    integer :: lm_max      ! number of l/m combinations
    integer :: n_r_tot     ! total number of radial grid points
- 
+
    !--- Now quantities for magnetic fields:
    !    Set lMag=0 if you want to save this memory (see c_fields)!
    integer :: lMagMem       ! Memory for magnetic field calculation
@@ -44,7 +44,7 @@ module truncation
    integer :: n_r_totMag    ! n_r_maxMag + n_r_ic_maxMag
    integer :: l_maxMag      ! Max. degree for magnetic field calculation
    integer :: lm_maxMag     ! Max. number of l/m combinations for magnetic field calculation
- 
+
 contains
 
    subroutine initialize_truncation
@@ -125,8 +125,8 @@ end subroutine initialize_truncation
 !--------------------------------------------------------------------------------
    subroutine checkTruncation
       !  This function checks truncations and writes it
-      !  into STDOUT and the log-file.                                 
-      !  MPI: called only by the processor responsible for output !  
+      !  into STDOUT and the log-file.
+      !  MPI: called only by the processor responsible for output !
 
       if ( minc < 1 ) then
          call abortRun('! Wave number minc should be > 0!')