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-2013 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

implicit none
integer(i4b) :: i, j
real(dp) :: time, delta_ts, glac_index, z_sl
real(dp) :: time_val, delta_ts_val, glac_index_val, z_sl_val
real(dp) :: h_val(0:jmax,0:imax)

!-------- 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(13,1011) time_val, delta_ts_val, z_sl_val
else if (forcing_flag == 2) then
   write(13,1011) time_val, glac_index_val, z_sl_val
else if (forcing_flag == 3) then
   write(13,1011) time_val, 1.11e+11_dp, z_sl_val
end if

do j=jmax, 0, -1
   write(13,1012) (h_val(j,i), i=0,imax)
end do

 1011 format(/,1pe13.6,2(1pe13.4),/)
 1012 format(imax(f8.1),f8.1)

end subroutine output3
!