output3.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  o u t p u t 3
!
!> @file
!!
!! Writing of time-series data of the ice thickness field on file
!! in ASCII format.
!!
!! @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/>.
!<
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

!-------------------------------------------------------------------------------
!> Writing of time-series data of the ice thickness field on file
!! in ASCII format.
!<------------------------------------------------------------------------------

subroutine output3(time, delta_ts, glac_index, z_sl)

use sico_types
use sico_variables
use sico_vars

implicit none

real(dp), intent(in) :: time, delta_ts, glac_index, z_sl

integer(i4b)       :: i, j
real(dp)           :: time_val, delta_ts_val, glac_index_val, z_sl_val
real(dp)           :: h_val(0:jmax,0:imax)

character(len=  8) :: ch_imax
character(len=128) :: fmt5, fmt6

fmt5 = '(/,1pe13.6,2(1pe13.4),/)'

write(ch_imax, fmt='(i8)') imax
write(fmt6,    fmt='(a)')  '('//trim(adjustl(ch_imax))//'(f8.1),f8.1)'

!-------- Conversion --------

#if ( !defined(OUT_TIMES) || OUT_TIMES==1 )
time_val = time /year_sec   ! s -> a
#elif OUT_TIMES == 2
time_val = (time+year_zero) /year_sec   ! s -> a
#else
stop ' output3: OUT_TIMES must be either 1 or 2!'
#endif

delta_ts_val   = delta_ts         ! in deg C
glac_index_val = glac_index
z_sl_val       = z_sl             ! in m

h_val          = h_c + h_t        ! in m

!-------- Write data on time-series file --------

if (forcing_flag == 1) then
   write(unit=13, fmt=trim(fmt5)) time_val, delta_ts_val, z_sl_val
else if (forcing_flag == 2) then
   write(unit=13, fmt=trim(fmt5)) time_val, glac_index_val, z_sl_val
else if (forcing_flag == 3) then
   write(unit=13, fmt=trim(fmt5)) time_val, 1.11e+11_dp, z_sl_val
end if

do j=jmax, 0, -1
   write(unit=13, fmt=trim(fmt6)) (h_val(j,i), i=0,imax)
end do

end subroutine output3
!