output4.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  o u t p u t 4
!
!> @file
!!
!! Writing of time-series data of the deep ice cores 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 deep ice cores on file in ASCII format.
!<------------------------------------------------------------------------------
subroutine output4(time, dxi, deta, delta_ts, glac_index, z_sl)

use sico_types
use sico_variables
use sico_vars

implicit none
real(dp), intent(in) :: time, dxi, deta, delta_ts, glac_index, z_sl

integer(i4b) :: i, j, n
real(dp) :: time_val
real(dp), dimension(0:JMAX,0:IMAX) :: field
real(dp), dimension(:), allocatable :: h_core, temp_core, &
                                       vx_core, vy_core, v_core, &
                                       rbx_core, rby_core, rb_core

allocate(h_core(n_core), temp_core(n_core), &
         vx_core(n_core), vy_core(n_core), v_core(n_core), &
         rbx_core(n_core), rby_core(n_core), rb_core(n_core))

!-------- Determination of ice-core data --------

do n=1, n_core

!  ------ Ice thickness

   field = h_c + h_t
   call borehole(field, x_core(n), y_core(n), dxi, deta, 'grid', h_core(n))

!  ------ Surface velocity

   do i=0, imax
   do j=0, jmax
      field(j,i) = vx_c(kcmax,j,i)
   end do
   end do

   call borehole(field, x_core(n), y_core(n), dxi, deta, 'sg_x', vx_core(n))

   do i=0, imax
   do j=0, jmax
      field(j,i) = vy_c(kcmax,j,i)
   end do
   end do

   call borehole(field, x_core(n), y_core(n), dxi, deta, 'sg_y', vy_core(n))

   v_core(n) = sqrt(vx_core(n)**2+vy_core(n)**2)

!  ------ Basal temperature

   do i=0, imax
   do j=0, jmax
      field(j,i) = temp_r(krmax,j,i)
   end do
   end do

   call borehole(field, x_core(n), y_core(n), dxi, deta, 'grid', temp_core(n))

!  ------ Basal frictional heating

   do i=0, imax
   do j=0, jmax
      field(j,i) = vx_t(0,j,i)*txz_t(0,j,i)
   end do
   end do

   call borehole(field, x_core(n), y_core(n), dxi, deta, 'sg_x', rbx_core(n))

   do i=0, imax
   do j=0, jmax
      field(j,i) = vy_t(0,j,i)*tyz_t(0,j,i)
   end do
   end do

   call borehole(field, x_core(n), y_core(n), dxi, deta, 'sg_y', rby_core(n))

   rb_core(n) = rbx_core(n) + rby_core(n)

end do

!  ------ 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 ' output4: OUT_TIMES must be either 1 or 2!'
#endif

do n=1, n_core
   v_core(n) = v_core(n) *year_sec   ! m/s -> m/a
end do

!-------- Writing of data on file --------

if (forcing_flag == 1) then
   write(unit=14, fmt='(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
else if (forcing_flag == 2) then
   write(unit=14, fmt='(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
else if (forcing_flag == 3) then
   write(unit=14, fmt='(1pe13.6,2(1pe13.4))') time_val, 1.11e+11_dp, z_sl
end if

n=1
write(unit=14, fmt='(13x,1pe13.4)', advance='no') h_core(n)
do n=2, n_core-1
   write(unit=14, fmt='(1pe13.4)', advance='no') h_core(n)
end do
n=n_core
write(unit=14, fmt='(1pe13.4)') h_core(n)

n=1
write(unit=14, fmt='(13x,1pe13.4)', advance='no') v_core(n)
do n=2, n_core-1
   write(unit=14, fmt='(1pe13.4)', advance='no') v_core(n)
end do
n=n_core
write(unit=14, fmt='(1pe13.4)') v_core(n)

n=1
write(unit=14, fmt='(13x,1pe13.4)', advance='no') temp_core(n)
do n=2, n_core-1
   write(unit=14, fmt='(1pe13.4)', advance='no') temp_core(n)
end do
n=n_core
write(unit=14, fmt='(1pe13.4)') temp_core(n)

n=1
write(unit=14, fmt='(13x,1pe13.4)', advance='no') rb_core(n)
do n=2, n_core-1
   write(unit=14, fmt='(1pe13.4)', advance='no') rb_core(n)
end do
n=n_core
write(unit=14, fmt='(1pe13.4,/)') rb_core(n)

deallocate(h_core, vx_core, vy_core, v_core, temp_core, &
           rbx_core, rby_core, rb_core)

end subroutine output4
!