38 subroutine output3(time, delta_ts, glac_index, z_sl)
45 real(dp),
intent(in) :: time, delta_ts, glac_index, z_sl
48 real(dp) :: time_val, delta_ts_val, glac_index_val, z_sl_val
49 real(dp) :: h_val(0:jmax,0:imax)
51 character(len= 8) :: ch_imax
52 character(len=128) :: fmt5, fmt6
54 fmt5 =
'(/,1pe13.6,2(1pe13.4),/)'
56 write(ch_imax, fmt=
'(i8)') imax
57 write(fmt6, fmt=
'(a)')
'('//trim(adjustl(ch_imax))//
'(f8.1),f8.1)'
61 #if ( !defined(OUT_TIMES) || OUT_TIMES==1 )
62 time_val = time /year_sec
64 time_val = (time+year_zero) /year_sec
66 stop
' output3: OUT_TIMES must be either 1 or 2!'
69 delta_ts_val = delta_ts
70 glac_index_val = glac_index
77 if (forcing_flag == 1)
then
78 write(unit=13, fmt=trim(fmt5)) time_val, delta_ts_val, z_sl_val
79 else if (forcing_flag == 2)
then
80 write(unit=13, fmt=trim(fmt5)) time_val, glac_index_val, z_sl_val
81 else if (forcing_flag == 3)
then
82 write(unit=13, fmt=trim(fmt5)) time_val, 1.11e+11_dp, z_sl_val
86 write(unit=13, fmt=trim(fmt6)) (h_val(j,i), i=0,imax)