output5.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  o u t p u t 5
!
!> @file
!!
!! Writing of time-series data for all defined surface points on file
!! in ASCII format. Modification of Tolly's output7 by Thorben Dunse.
!!
!! @section Copyright
!!
!! Copyright 2009-2016 Thorben Dunse, 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 for all defined surface points on file
!! in ASCII format. Modification of Tolly's output7 by Thorben Dunse.
!<------------------------------------------------------------------------------
subroutine output5(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) :: n, k
real(dp) :: time_val
real(dp), dimension(0:JMAX,0:IMAX) :: field
real(dp), dimension(:), allocatable :: zl_surf, zs_surf, &
                  accum_surf, as_perp_surf, &
                  snowfall_surf, rainfall_surf, runoff_surf, &
                  vx_surf, vy_surf, vz_surf, &
                  vx_base, vy_base, vz_base, &
                  temp_base_pmp

allocate(zl_surf(n_surf), zs_surf(n_surf), &
        accum_surf(n_surf), &
        as_perp_surf(n_surf), snowfall_surf(n_surf), &
        rainfall_surf(n_surf), runoff_surf(n_surf), &
        vx_surf(n_surf), vy_surf(n_surf), vz_surf(n_surf), &
        vx_base(n_surf), vy_base(n_surf), vz_base(n_surf), &
        temp_base_pmp(n_surf))

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

do n=1, n_surf

!  ------ Bedrock elevation

   field = zl
   call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', zl_surf(n))

!  ------ Surface elevation

   field = zs
   call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', zs_surf(n))


!  ------ Accumulation

   field = accum
   call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', accum_surf(n))

!  ------ Surface mass balance

   field = as_perp
   call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', as_perp_surf(n))

!  ------ Snowfall

   field = snowfall
   call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', snowfall_surf(n))

!  ------ Rainfall

   field = rainfall
   call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', rainfall_surf(n))

!  ------ Runoff

   field = runoff
   call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', runoff_surf(n))

! ------ Surface velocities

      field = vx_s_g
      call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', vx_surf(n))

      field = vy_s_g
      call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', vy_surf(n))

      field = vz_s
      call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', vz_surf(n))

! ------ Basal velocities

      field = vx_b_g
      call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', vx_base(n))

      field = vy_b_g
      call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', vy_base(n))

      field = vz_b
      call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', vz_base(n))

! ------ Basal temperature relative to pressure melting point

      field = temph_b
      call borehole(field, x_surf(n), y_surf(n), dxi, deta, &
        'grid', temp_base_pmp(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 ' output5: OUT_TIMES must be either 1 or 2!'
#endif

do n=1, n_surf
   accum_surf(n)    = accum_surf(n)    *year_sec   ! m/s -> m/a
   as_perp_surf(n)  = as_perp_surf(n)  *year_sec   ! m/s -> m/a
   snowfall_surf(n) = snowfall_surf(n) *year_sec   ! m/s -> m/a
   rainfall_surf(n) = rainfall_surf(n) *year_sec   ! m/s -> m/a
   runoff_surf(n)   = runoff_surf(n)   *year_sec   ! m/s -> m/a
   vx_surf(n)       = vx_surf(n)       *year_sec   ! m/s -> m/a
   vy_surf(n)       = vy_surf(n)       *year_sec   ! m/s -> m/a
   vz_surf(n)       = vz_surf(n)       *year_sec   ! m/s -> m/a
   vx_base(n)       = vx_base(n)       *year_sec   ! m/s -> m/a
   vy_base(n)       = vy_base(n)       *year_sec   ! m/s -> m/a
   vz_base(n)       = vz_base(n)       *year_sec   ! m/s -> m/a
end do

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

if (forcing_flag == 1) then
   write(41,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(42,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(43,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(44,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(45,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(46,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(47,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(48,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(49,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(50,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(51,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(52,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(53,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
   write(54,'(1pe13.6,2(1pe13.4))') time_val, delta_ts, z_sl
else if (forcing_flag == 2) then
   write(41,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(42,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(43,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(44,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(45,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(46,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(47,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(48,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(49,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(50,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(51,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(52,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(53,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
   write(54,'(1pe13.6,2(1pe13.4))') time_val, glac_index, z_sl
end if

do n=1, n_surf-1
   write(41,'(1pe13.4)',advance='no') zl_surf(n)
   write(42,'(1pe13.4)',advance='no') zs_surf(n)
   write(43,'(1pe13.4)',advance='no') accum_surf(n)
   write(44,'(1pe13.4)',advance='no') as_perp_surf(n)
   write(45,'(1pe13.4)',advance='no') snowfall_surf(n)
   write(46,'(1pe13.4)',advance='no') rainfall_surf(n)
   write(47,'(1pe13.4)',advance='no') runoff_surf(n)
   write(48,'(1pe13.4)',advance='no') vx_surf(n)
   write(49,'(1pe13.4)',advance='no') vy_surf(n)
   write(50,'(1pe13.4)',advance='no') vz_surf(n)
   write(51,'(1pe13.4)',advance='no') vx_base(n)
   write(52,'(1pe13.4)',advance='no') vy_base(n)
   write(53,'(1pe13.4)',advance='no') vz_base(n)
   write(54,'(1pe13.4)',advance='no') temp_base_pmp(n) 
end do

n=n_surf
   write(41,'(1pe13.4)') zl_surf(n)
   write(42,'(1pe13.4)') zs_surf(n)
   write(43,'(1pe13.4)') accum_surf(n)
   write(44,'(1pe13.4)') as_perp_surf(n)
   write(45,'(1pe13.4)') snowfall_surf(n)
   write(46,'(1pe13.4)') rainfall_surf(n)
   write(47,'(1pe13.4)') runoff_surf(n)
   write(48,'(1pe13.4)') vx_surf(n)
   write(49,'(1pe13.4)') vy_surf(n)
   write(50,'(1pe13.4)') vz_surf(n)
   write(51,'(1pe13.4)') vx_base(n)
   write(52,'(1pe13.4)') vy_base(n)
   write(53,'(1pe13.4)') vz_base(n)
   write(54,'(1pe13.4)') temp_base_pmp(n) 

deallocate(zl_surf, zs_surf, accum_surf, as_perp_surf, &
           snowfall_surf, rainfall_surf, runoff_surf, &
           vx_surf, vy_surf, vz_surf, &
           vx_base, vy_base, vz_base,temp_base_pmp)

end subroutine output5
!