sico_end.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  s i c o _ e n d
!
!> @file
!!
!! Ending of SICOPOLIS.
!!
!! @section Copyright
!!
!! Copyright 2009-2016 Ralf Greve
!!
!! @section License
!!
!! This file is part of SICOPOLIS.
!!
!! SICOPOLIS is free software: you can redistribute it and/or modify
!! it under the terms of the GNU General Public License as published by
!! the Free Software Foundation, either version 3 of the License, or
!! (at your option) any later version.
!!
!! SICOPOLIS is distributed in the hope that it will be useful,
!! but WITHOUT ANY WARRANTY; without even the implied warranty of
!! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!! GNU General Public License for more details.
!!
!! You should have received a copy of the GNU General Public License
!! along with SICOPOLIS.  If not, see <http://www.gnu.org/licenses/>.
!<
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

!-------------------------------------------------------------------------------
!> Ending of SICOPOLIS.
!<------------------------------------------------------------------------------
subroutine sico_end()

use sico_types
use sico_variables
use sico_vars
#if (NETCDF > 1)
use netcdf
use nc_check
#endif

implicit none

integer(i4b) :: ierr

character(len=64), parameter :: thisroutine = 'sico_end'

close(unit=12, status='keep')  ! Close time-series files
#if (OUTSER==2)
close(unit=13, status='keep')
#endif
#if (OUTSER==3)
close(unit=14, status='keep')
deallocate(lambda_core, phi_core, x_core, y_core)
#endif
#if (OUTSER==4)
close(unit=41, status='keep')
close(unit=42, status='keep')
close(unit=43, status='keep')
close(unit=44, status='keep')
close(unit=45, status='keep')
close(unit=46, status='keep')
close(unit=47, status='keep')
close(unit=48, status='keep')
close(unit=49, status='keep')
close(unit=50, status='keep')
deallocate(lambda_surf, phi_surf, x_surf, y_surf)
#endif
#if (defined(XYZ))
#if (defined(HEINO))
close(unit=15, status='keep')
#endif
#endif

#if ((TSURFACE==6) && (ACCSURFACE==6))
call check( nf90_close(ncid_temp_precip), thisroutine )
            ! Closing of NetCDF file containing surface-temperature
            ! and precipitation data as functions of time
#endif

#if (CALCZS==4 || MARGIN==3)
call lis_finalize(ierr)   ! Finalise execution environment of the
                          ! Library of Iterative Solvers Lis, if required
#endif

write(unit=6, fmt='(a)') ' '
write(unit=6, fmt='(a)') ' '
write(unit=6, fmt='(a)') &
'             * * * sicopolis.F90  r e a d y * * *'
write(unit=6, fmt='(a)') ' '
write(unit=6, fmt='(a)') ' '

end subroutine sico_end
!