diff options
Diffstat (limited to 'src/EHFinder_Integrate.F90')
| -rw-r--r-- | src/EHFinder_Integrate.F90 | 20 |
1 files changed, 10 insertions, 10 deletions
diff --git a/src/EHFinder_Integrate.F90 b/src/EHFinder_Integrate.F90 index 9e5257d..b8b8fcb 100644 --- a/src/EHFinder_Integrate.F90 +++ b/src/EHFinder_Integrate.F90 @@ -92,9 +92,9 @@ subroutine EHFinder_FindSurfaceElement(CCTK_ARGUMENTS) if ( interp_handle .lt. 0 ) then warn_message = 'Cannot get handle for interpolation. ' - warn_message = warn_message//'Forgot to activate an implementation ' - warn_message = warn_message//'providing interpolation operators?' - call CCTK_WARN( 0, warn_message ) + warn_message = trim(warn_message)//'Forgot to activate an implementation ' + warn_message = trim(warn_message)//'providing interpolation operators?' + call CCTK_WARN( 0, trim(warn_message) ) end if ! Write the interpolation order parameter into a string @@ -110,9 +110,9 @@ subroutine EHFinder_FindSurfaceElement(CCTK_ARGUMENTS) call CCTK_CoordSystemHandle ( coord_system_handle, 'cart3d' ) if ( coord_system_handle .lt. 0) then warn_message = 'Cannot get handle for cart3d coordinate system. ' - warn_message = warn_message//'Forgot to activate an implementation ' - warn_message = warn_message//'providing coordinates?' - call CCTK_WARN( 0, warn_message ) + warn_message = trim(warn_message)//'Forgot to activate an implementation ' + warn_message = trim(warn_message)//'providing coordinates?' + call CCTK_WARN( 0, trim(warn_message) ) endif ! Get the pointers to the interpolation points. @@ -271,7 +271,7 @@ subroutine EHFinder_IntegrateArea(CCTK_ARGUMENTS) ! Write an info message reporting the area found. write(info_message,'(a15,f10.5)') 'Horizon area = ',eh_area_tmp - call CCTK_INFO(info_message) + call CCTK_INFO( trim(info_message) ) ! Figure out the bounds of the area storage grid array. call CCTK_GrouplbndGN ( ierr, cctkGH, 1, lbnd, 'ehfinder::eh_area2' ) @@ -394,7 +394,7 @@ subroutine EHFinder_IntegrateCentroid(CCTK_ARGUMENTS) eh_centroid2_x(sn,l), & eh_centroid2_y(sn,l), & eh_centroid2_z(sn,l) - call CCTK_INFO(info_message) + call CCTK_INFO( trim(info_message) ) end subroutine EHFinder_IntegrateCentroid @@ -509,7 +509,7 @@ subroutine EHFinder_IntegrateCircumference(CCTK_ARGUMENTS) ! Write an info message with the calculated circunference. write(info_message,'(a35,f10.5)') 'Horizon equatorial circumference = ', & eq_circ - call CCTK_INFO(info_message) + call CCTK_INFO( trim(info_message) ) ! Now do the same for the polar circomference. call CCTK_INFO ( 'Integrating polar circumference' ) @@ -543,7 +543,7 @@ subroutine EHFinder_IntegrateCircumference(CCTK_ARGUMENTS) ! Write an info message with the calculated circumference. write(info_message,'(a30,f10.5)') 'Horizon polar circumference = ',pol_circ - call CCTK_INFO(info_message) + call CCTK_INFO( trim(info_message) ) end subroutine EHFinder_IntegrateCircumference |