subroutine sed_rest_in (fname, ids, ide, jds, jde)
#if defined O_sed
!=======================================================================
! input routine for sediment model restarts
! data may be sized differently in x and y from the global fields.
! fields may be written with or without a time dimension. data
! should be defined with the routine defvar and written with putvar.
! if no time dimension, then data is only written once per file.
! make sure the it, iu, ib, and ic arrays and are defining the
! correct dimensions. ln may also need to be recalculated.
! inputs:
! fname = file name
! ids, ide ... = start and end index for data domain
!=======================================================================
implicit none
character(*) :: fname
character(32) :: nstamp
character(3) :: a3
character(120) :: var1, var2
integer i, iou, j, ln, ntrec, ids, ide, jds, jde, l, m, n
integer ils, ile, jls, jle, lls, lle, ib(10), ic(10)
integer nyear, nmonth, nday, nhour, nmin, nsec
logical inqvardef, warning
real tmp
real, allocatable :: tmpij(:,:)
real, allocatable :: tmpip(:)
include "size.h"
include "param.h"
include "pconst.h"
include "stdunits.h"
include "coord.h"
include "grdvar.h"
include "iounit.h"
include "levind.h"
include "tmngr.h"
include "switch.h"
include "csbc.h"
include "sed.h"
integer map_tmp(imt,jmt)
real data(imt,jmt)
!-----------------------------------------------------------------------
! open file
!-----------------------------------------------------------------------
call openfile (fname, iou)
ntrec = 1
if (.not. inqvardef('atsed', iou)) then
print*, "==> Warning: sediment model restart not found"
return
endif
!-----------------------------------------------------------------------
! local domain size (minimum of data domain and global read domain)
!-----------------------------------------------------------------------
lls = 1
lle = ipmax
allocate ( tmpip(lls:lle) )
ils = max(ids,1)
ile = min(ide,imt)
jls = max(jds,1)
jle = min(jde,jmt)
allocate ( tmpij(ils:ile,jls:jle) )
map_tmp(:,:) = map_sed(:,:)
!-----------------------------------------------------------------------
! read 1d data (t)
!-----------------------------------------------------------------------
tmp = itt
call getvars ('itt', iou, ntrec, tmp, c1, c0)
itt = tmp
tmp = irstdy
call getvars ('irstdy', iou, ntrec, tmp, c1, c0)
irstdy = tmp
tmp = msrsdy
call getvars ('msrsdy', iou, ntrec, tmp, c1, c0)
msrsdy = tmp
tmp = year0
call getvars ('year', iou, 1, tmp, c1, c0)
nyear = tmp
tmp = month0
call getvars ('month', iou, 1, tmp, c1, c0)
nmonth = tmp
tmp = day0
call getvars ('day', iou, 1, tmp, c1, c0)
nday = tmp
tmp = hour0
call getvars ('hour', iou, 1, tmp, c1, c0)
nhour = tmp
tmp = min0
call getvars ('minute', iou, 1, tmp, c1, c0)
nmin = tmp
tmp = sec0
call getvars ('second', iou, 1, tmp, c1, c0)
nsec = tmp
call mkstmp (stamp, nyear, nmonth, nday, nhour, nmin, nsec)
tmp = carblith
call getvars ('carblith', iou, 1, tmp, c1, c0)
carblith = tmp
if (init_time_in) then
itt = 0
irstdy = 0
msrsdy = 0
relyr = 0.
call mkstmp (stamp, year0, month0, day0, hour0, min0, sec0)
carblith = 0.
endif
tmp = atsed
call getvars ('atsed', iou, 1, tmp, c1, c0)
atsed = tmp
tmp = weathflx
call getvars ('weathflx', iou, 1, tmp, c1, c0)
weathflx = tmp
tmp = sed_year
call getvars ('sed_year', iou, 1, tmp, c1, c0)
sed_year = tmp
!-----------------------------------------------------------------------
! read 2d data (x,y)
!-----------------------------------------------------------------------
ib(1) = 1
ic(1) = ile-ils+1
ib(2) = 1
ic(2) = jle-jls+1
ln = ic(1)*ic(2)
tmpij(:,:) = kmt(ils:ile,jls:jle)
call getvara ('kmt', iou, ln, ib, ic, tmpij, c1, c0)
map_tmp(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call getvara ('map_sed', iou, ln, ib, ic, tmpij, c1, c0)
map_tmp(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
warning = .false.
do j=jls,jle
do i=ils,ile
if (map_tmp(i,j) .ne. map_sed(i,j)) warning = .true.
enddo
enddo
if (warning) then
print*, ' '
print*, '==> Warning: restart sediment mask is different.'
endif
!-----------------------------------------------------------------------
! read 3d data (x,y,t)
!-----------------------------------------------------------------------
ib(1) = 1
ic(1) = ile-ils+1
ib(2) = 1
ic(2) = jle-jls+1
ib(3) = 1
ic(3) = 1
ln = ic(1)*ic(2)*ic(3)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibtemp)
call getvara ('sbc_btemp', iou, ln, ib, ic, tmpij, c1, c0)
sbc(ils:ile,jls:jle,ibtemp) = tmpij(ils:ile,jls:jle)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibsalt)
call getvara ('sbc_bsalt', iou, ln, ib, ic, tmpij, c1, c0)
sbc(ils:ile,jls:jle,ibsalt) = tmpij(ils:ile,jls:jle)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ircal)
call getvara ('sbc_rcal', iou, ln, ib, ic, tmpij, c1, c0)
sbc(ils:ile,jls:jle,ircal) = tmpij(ils:ile,jls:jle)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,irorg)
call getvara ('sbc_rorg', iou, ln, ib, ic, tmpij, c1, c0)
sbc(ils:ile,jls:jle,irorg) = tmpij(ils:ile,jls:jle)
# if defined O_carbon
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibdic)
call getvara ('sbc_bdic', iou, ln, ib, ic, tmpij, c1, c0)
sbc(ils:ile,jls:jle,ibdic) = tmpij(ils:ile,jls:jle)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibdicfx)
call getvara ('sbc_bdicfx', iou, ln, ib, ic, tmpij, c1, c0)
sbc(ils:ile,jls:jle,ibdicfx) = tmpij(ils:ile,jls:jle)
# endif
# if defined O_npzd_alk
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibalk)
call getvara ('sbc_balk', iou, ln, ib, ic, tmpij, c1, c0)
sbc(ils:ile,jls:jle,ibalk) = tmpij(ils:ile,jls:jle)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibalkfx)
call getvara ('sbc_balkfx', iou, ln, ib, ic, tmpij, c1, c0)
sbc(ils:ile,jls:jle,ibalkfx) = tmpij(ils:ile,jls:jle)
# endif
# if defined O_npzd_o2
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibo2)
call getvara ('sbc_bo2', iou, ln, ib, ic, tmpij, c1, c0)
sbc(ils:ile,jls:jle,ibo2) = tmpij(ils:ile,jls:jle)
# endif
call unloadsed (ipmax, zrct, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_zrct', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
zrct(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, water_z_p, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_water_z_p', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
water_z_p(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, k1, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_k1', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
k1(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, k2, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_k2', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
k2(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, k3, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_k3', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
k3(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, csat, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_csat', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
csat(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, rc, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_rc', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
rc(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, ttrorg, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_ttrorg', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
ttrorg(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, ttrcal, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_ttrcal', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
ttrcal(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, sed_ml_mass, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_ml_mass', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
sed_ml_mass(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, c_advect, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_c_advect', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
c_advect(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, rain_cal_p, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_rain_cal_p', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
rain_cal_p(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, rain_org_p, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_rain_org_p', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
rain_org_p(1:ipmax) = tmpip(1:ipmax)
call unloadsed (ipmax, co3_p, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_co3_p', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
co3_p(1:ipmax) = tmpip(1:ipmax)
!-----------------------------------------------------------------------
! read 4d data (x,y,z,t)
!-----------------------------------------------------------------------
ib(1) = 1
ic(1) = ile-ils+1
ib(2) = 1
ic(2) = jle-jls+1
ib(3) = 1
ic(3) = 1
ib(4) = 1
ic(4) = 1
ln = ic(1)*ic(2)*ic(3)*ic(4)
n = 0
call getdimlen ('slev', iou, n)
n = min(n, nzmax)
do l=1,n
ib(3) = l
do m=1,3
write(a3, '(i1)') m
tmpip(1:ipmax) = carb(l,m,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_carb_'//trim(a3), iou, ln, ib, ic, tmpij
&, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
carb(l,m,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = dcpls(l,m,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_dcpls_'//trim(a3), iou, ln, ib, ic, tmpij
&, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
dcpls(l,m,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = dcmin(l,m,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_dcmin_'//trim(a3), iou, ln, ib, ic, tmpij
&, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
dcmin(l,m,1:ipmax) = tmpip(1:ipmax)
enddo
tmpip(1:ipmax) = pore(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_pore', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
pore(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = form(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_form', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
form(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = o2(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_o2', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
o2(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = orggg(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_orggg', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
orggg(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = orgml(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_orgml', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
orgml(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = calgg(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_calgg', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
calgg(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = calml(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_calml', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
calml(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = dopls(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_dopls', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
dopls(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = domin(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_domin', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
domin(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = dbpls(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_dbpls', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
dbpls(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = dbmin(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_dbmin', iou, ln, ib, ic, tmpij, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
dbmin(l,1:ipmax) = tmpip(1:ipmax)
enddo
n = 0
call getdimlen ('blev', iou, n)
n = min(n, ibmax)
do l=1,n
ib(3) = l
tmpip(1:ipmax) = buried_mass(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_buried_mass', iou, ln, ib, ic, tmpij
&, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
buried_mass(l,1:ipmax) = tmpip(1:ipmax)
tmpip(1:ipmax) = buried_calfrac(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call getvara ('sed_buried_calfrac', iou, ln, ib, ic, tmpij
&, c1, c0)
where (map_tmp(ils:ile,jls:jle) .ne. 0)
& data(ils:ile,jls:jle) = tmpij(ils:ile,jls:jle)
call loadsed (ipmax, tmpip, imt, jmt, map_sed, data)
buried_calfrac(l,1:ipmax) = tmpip(1:ipmax)
enddo
deallocate (tmpip)
call rdstmp (stamp, nyear, nmonth, nday, nhour, nmin, nsec)
nyear = year0 + accel_yr0 + (relyr - accel_yr0)*accel
call mkstmp (nstamp, nyear, nmonth, nday, nhour, nmin, nsec)
print*, '=> Sed restart read from ',trim(fname),' on ', nstamp
deallocate (tmpij)
return
end
subroutine sed_rest_def (fname)
!=======================================================================
! definition routine for sediment model restarts
! inputs:
! fname = file name
!=======================================================================
implicit none
character(*) :: fname
character(3) :: a3
integer iou, ntrec, igs, ige, ig, jgs, jge, jg, m
integer it(10), id_time, id_xt, id_yt, id_xt_e, id_yt_e
integer id_slev, id_slev_e, id_blev, id_blev_e
real c100, c1e3, c1e6, c1e20
include "size.h"
include "param.h"
include "pconst.h"
include "stdunits.h"
include "iounit.h"
include "mw.h"
include "tmngr.h"
include "sed.h"
c100 = 100.
c1e3 = 1.e3
c1e6 = 1.e6
c1e20 = 1.e20
!-----------------------------------------------------------------------
! open file
!-----------------------------------------------------------------------
call openfile (fname, iou)
ntrec = 1
!-----------------------------------------------------------------------
! set global write domain size
!-----------------------------------------------------------------------
igs = 1
ige = imt
ig = ige-igs+1
jgs = 1
jge = jmt
jg = jge-jgs+1
!-----------------------------------------------------------------------
! start definitions
!-----------------------------------------------------------------------
call redef (iou)
!-----------------------------------------------------------------------
! write global atributes
!-----------------------------------------------------------------------
call putatttext (iou, 'global', 'Conventions', 'CF-1.0')
call putatttext (iou, 'global', 'experiment_name', expnam)
call putatttext (iou, 'global', 'run_stamp', runstamp)
!-----------------------------------------------------------------------
! define dimensions
!-----------------------------------------------------------------------
call defdim ('time', iou, 0, id_time)
call defdim ('longitude', iou, ig, id_xt)
call defdim ('latitude', iou, jg, id_yt)
call defdim ('longitude_edges', iou, ig+1, id_xt_e)
call defdim ('latitude_edges', iou, jg+1, id_yt_e)
call defdim ('slev', iou, nzmax, id_slev)
call defdim ('slev_edges', iou, nzmax+1, id_slev_e)
call defdim ('blev', iou, ibmax, id_blev)
call defdim ('blev_edges', iou, ibmax+1, id_blev_e)
!-----------------------------------------------------------------------
! define 1d data (t)
!-----------------------------------------------------------------------
it(1) = id_time
call defvar ('time', iou, 1, it, c0, c0, 'T', 'D'
&, 'time', 'time', 'years since 0-1-1')
call putatttext (iou, 'time', 'calendar', calendar)
call defvar ('itt', iou, 1, it, c0, c0, ' ', 'D'
&, 'itt', ' ',' ')
call defvar ('irstdy', iou, 1, it, c0, c0, ' ', 'D'
&, 'irstdy', ' ',' ')
call defvar ('msrsdy', iou, 1, it, c0, c0, ' ', 'D'
&, 'msrsdy', ' ',' ')
call defvar ('year', iou, 1, it, c0, c0, ' ', 'D'
&, 'year', ' ',' ')
call defvar ('month', iou, 1, it, c0, c0, ' ', 'D'
&, 'month', ' ',' ')
call defvar ('day', iou, 1, it, c0, c0, ' ', 'D'
&, 'day', ' ',' ')
call defvar ('hour', iou, 1, it, c0, c0, ' ', 'D'
&, 'hour', ' ',' ')
call defvar ('minute', iou, 1, it, c0, c0, ' ', 'D'
&, 'minute', ' ',' ')
call defvar ('second', iou, 1, it, c0, c0, ' ', 'D'
&, 'second', ' ',' ')
call defvar ('carblith', iou, 1, it, c0, c0, ' ', 'D'
&, 'carblith', ' ',' ')
call defvar ('atsed', iou, 1, it, c0, c0, ' ', 'D'
&, 'atsed', ' ',' ')
call defvar ('weathflx', iou, 1, it, c0, c0, ' ', 'D'
&, 'weathflx', ' ',' ')
call defvar ('sed_year', iou, 1, it, c0, c0, ' ', 'D'
&, 'sed_year', ' ',' ')
!-----------------------------------------------------------------------
! define 1d data (x, y or z)
!-----------------------------------------------------------------------
it(1) = id_xt
call defvar ('longitude', iou, 1, it, c0, c0, 'X', 'D'
&, 'longitude', 'longitude', 'degrees_east')
it(1) = id_yt
call defvar ('latitude', iou, 1, it, c0, c0, 'Y', 'D'
&, 'latitude', 'latitude', 'degrees_north')
it(1) = id_xt_e
call defvar ('longitude_edges', iou, 1, it, c0, c0, ' ', 'D'
&, 'longitude edges', ' ', 'degrees_east')
it(1) = id_yt_e
call defvar ('latitude_edges', iou, 1, it, c0, c0, ' ', 'D'
&, 'latitude edges', ' ', 'degrees_north')
it(1) = id_slev
call defvar ('slev', iou, 1, it, c0, c0, ' ', 'D'
&, 'sediment level', ' ','1')
it(1) = id_slev_e
call defvar ('slev_edges', iou, 1, it, c0, c0, ' ', 'D'
&, 'sediment level edges', ' ','1')
it(1) = id_blev
call defvar ('blev', iou, 1, it, c0, c0, ' ', 'D'
&, 'burial level', ' ','1')
it(1) = id_blev_e
call defvar ('blev_edges', iou, 1, it, c0, c0, ' ', 'D'
&, 'burial level edges', ' ','1')
!-----------------------------------------------------------------------
! define 2d data (x,y)
!-----------------------------------------------------------------------
it(1) = id_xt
it(2) = id_yt
call defvar ('kmt', iou, 2, it, c0, c1e3, ' ', 'I'
&, 'kmt', ' ' ,' ')
call defvar ('map_sed', iou , 2, it, c0, c1e6, ' ', 'I'
&, 'map_sed', ' ', ' ')
!-----------------------------------------------------------------------
! define 3d data (x,y,t)
!-----------------------------------------------------------------------
it(1) = id_xt
it(2) = id_yt
it(3) = id_time
call defvar ('sbc_btemp', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sbc btemp', ' ', ' ')
call defvar ('sbc_bsalt', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sbc bsalt', ' ', ' ')
call defvar ('sbc_rcal', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sbc rcal', ' ', ' ')
call defvar ('sbc_rorg', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sbc rorg', ' ', ' ')
# if defined O_carbon
call defvar ('sbc_bdic', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sbc bdic', ' ', ' ')
call defvar ('sbc_bdicfx', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sbc bdicfx', ' ', ' ')
# endif
# if defined O_npzd_alk
call defvar ('sbc_balk', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sbc balk', ' ', ' ')
call defvar ('sbc_balkfx', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sbc balkfx', ' ', ' ')
# endif
# if defined O_npzd_o2
call defvar ('sbc_bo2', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sbc bo2', ' ', ' ')
# endif
call defvar ('sed_zrct', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'zrct', ' ', ' ')
call defvar ('sed_water_z_p', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'water_z_p', ' ', ' ')
call defvar ('sed_k1', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'k1', ' ', ' ')
call defvar ('sed_k2', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'k2', ' ', ' ')
call defvar ('sed_k3', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'k3', ' ', ' ')
call defvar ('sed_csat', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'csat', ' ', ' ')
call defvar ('sed_rc', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'rc', ' ', ' ')
call defvar ('sed_ttrorg', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'ttrorg', ' ', ' ')
call defvar ('sed_ttrcal', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'ttrcal', ' ', ' ')
call defvar ('sed_ml_mass', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'sed_ml_mass', ' ', ' ')
call defvar ('sed_c_advect', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'c_advect', ' ', ' ')
call defvar ('sed_rain_cal_p', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'rain_cal_p', ' ', ' ')
call defvar ('sed_rain_org_p', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'rain_org_p', ' ', ' ')
call defvar ('sed_co3_p', iou, 3, it, -c1e20, c1e20, ' '
&, 'D', 'co3_p', ' ', ' ')
!-----------------------------------------------------------------------
! define 4d data (x,y,z,t)
!-----------------------------------------------------------------------
it(1) = id_xt
it(2) = id_yt
it(3) = id_slev
it(4) = id_time
do m=1,3
write(a3, '(i1)') m
call defvar ('sed_carb_'//trim(a3), iou ,4, it, -c1e6, c1e6
&, ' ', 'D', 'sed_carb '//trim(a3), ' ', ' ')
call defvar ('sed_dcpls_'//trim(a3), iou ,4, it, -c1e6, c1e6
&, ' ', 'D', 'sed_dcpls_'//trim(a3), ' ', ' ')
call defvar ('sed_dcmin_'//trim(a3), iou ,4, it, -c1e6, c1e6
&, ' ', 'D', 'sed_dcmin_'//trim(a3), ' ', ' ')
enddo
call defvar ('sed_pore', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_pore', ' ', ' ')
call defvar ('sed_form', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_form', ' ', ' ')
call defvar ('sed_o2', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_o2', ' ', ' ')
call defvar ('sed_orggg', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_orggg', ' ', ' ')
call defvar ('sed_orgml', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_orgml', ' ', ' ')
call defvar ('sed_calgg', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_calgg', ' ', ' ')
call defvar ('sed_calml', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_calml', ' ', ' ')
call defvar ('sed_dopls', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_dopls', ' ', ' ')
call defvar ('sed_domin', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_domin', ' ', ' ')
call defvar ('sed_dbpls', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_dbpls', ' ', ' ')
call defvar ('sed_dbmin', iou, 4, it, -c1e20, c1e20, ' ', 'D'
&, 'sed_dbmin', ' ', ' ')
it(3) = id_blev
call defvar ('sed_buried_mass', iou, 4, it, -c1e20, c1e20
&, ' ', 'D', 'buried_mass', ' ', ' ')
call defvar ('sed_buried_calfrac', iou, 4, it, -c1e20, c1e20
&, ' ', 'D', 'buried_calfrac', ' ', ' ')
!-----------------------------------------------------------------------
! end definitions
!-----------------------------------------------------------------------
call enddef (iou)
return
end
subroutine sed_rest_out (fname, ids, ide, jds, jde)
!=======================================================================
! output routine for sediment model restarts
! data may be sized differently in x and y from the global fields.
! fields may be written with or without a time dimension. data
! should be defined with the routine defvar and written with putvar.
! if no time dimension, then data is only written once per file.
! make sure the it, iu, ib, and ic arrays and are defining the
! correct dimensions. ln may also need to be recalculated.
! inputs:
! fname = file name
! ids, ide ... = start and end index for data domain
!=======================================================================
implicit none
character(*) :: fname
character(3) :: a3
character(32) :: nstamp
character(120) :: var1, var2
integer i, iou, j, ln, l, m, n, ntrec, ids, ide, jds, jde, igs
integer ige, ig, jgs, jge, jg, ils, ile, jls, jle, lls, lle
integer ib(10), ic(10), id_xt, id_yt, nyear, nmonth, nday, nhour
integer nmin, nsec
real tmp, c100, c1e3, c1e20
real, allocatable :: tmpij(:,:), tmpi(:), tmpj(:), tmpie(:)
real, allocatable :: tmpje(:), tmpip(:)
include "size.h"
include "param.h"
include "pconst.h"
include "stdunits.h"
include "coord.h"
include "grdvar.h"
include "iounit.h"
include "levind.h"
include "switch.h"
include "tmngr.h"
include "csbc.h"
include "sed.h"
real xt_e(imt+1), yt_e(jmt+1), data(imt,jmt), slev(nzmax)
real slev_e(nzmax+1), blev(ibmax), blev_e(ibmax+1)
lls = 1
lle = ipmax
allocate ( tmpip(lls:lle) )
c100 = 100.
c1e3 = 1.e3
c1e20 = 1.e20
nstamp = stamp
!-----------------------------------------------------------------------
! open file
!-----------------------------------------------------------------------
call openfile (fname, iou)
ntrec = 1
!-----------------------------------------------------------------------
! set global write domain size
!-----------------------------------------------------------------------
igs = 1
ige = imt
ig = ige-igs+1
jgs = 1
jge = jmt
jg = jge-jgs+1
!-----------------------------------------------------------------------
! local domain size (minimum of data domain and global write domain)
!-----------------------------------------------------------------------
ils = max(ids,igs)
ile = min(ide,ige)
jls = max(jds,jgs)
jle = min(jde,jge)
allocate ( tmpij(ils:ile,jls:jle) )
allocate ( tmpi(igs:ige) )
allocate ( tmpj(jgs:jge) )
allocate ( tmpie(igs:ige+1) )
allocate ( tmpje(jgs:jge+1) )
!-----------------------------------------------------------------------
! write 1d data (t)
!-----------------------------------------------------------------------
if (init_time_out) then
tmp = 0.
call putvars ('time', iou, ntrec, tmp, c1, c0)
tmp = 0.
call putvars ('itt', iou, ntrec, tmp, c1, c0)
tmp = 0.
call putvars ('irstdy', iou, ntrec, tmp, c1, c0)
tmp = 0.
call putvars ('msrsdy', iou, ntrec, tmp, c1, c0)
call mkstmp (nstamp, year0, month0, day0, hour0, min0, sec0)
carblith = 0.
else
tmp = year0 + accel_yr0 + (relyr - accel_yr0)*accel
call putvars ('time', iou, ntrec, tmp, c1, c0)
tmp = itt
call putvars ('itt', iou, ntrec, tmp, c1, c0)
tmp = iday(imodeltime)
call putvars ('irstdy', iou, ntrec, tmp, c1, c0)
tmp = msday(imodeltime)
call putvars ('msrsdy', iou, ntrec, tmp, c1, c0)
endif
call rdstmp (nstamp, nyear, nmonth, nday, nhour, nmin, nsec)
tmp = nyear
call putvars ('year', iou, 1, tmp, c1, c0)
tmp = nmonth
call putvars ('month', iou, 1, tmp, c1, c0)
tmp = nday
call putvars ('day', iou, 1, tmp, c1, c0)
tmp = nhour
call putvars ('hour', iou, 1, tmp, c1, c0)
tmp = nmin
call putvars ('minute', iou, 1, tmp, c1, c0)
tmp = nsec
call putvars ('second', iou, 1, tmp, c1, c0)
tmp = carblith
call putvars ('carblith', iou, 1, tmp, c1, c0)
tmp = atsed
call putvars ('atsed', iou, 1, tmp, c1, c0)
tmp = weathflx
call putvars ('weathflx', iou, 1, tmp, c1, c0)
tmp = sed_year
call putvars ('sed_year', iou, 1, tmp, c1, c0)
!-----------------------------------------------------------------------
! write 1d data (x, y or z)
!-----------------------------------------------------------------------
ib(1) = 1
ic(1) = ig
tmpi(igs:ige) = xt(igs:ige)
call putvara ('longitude', iou, ig, ib, ic, tmpi, c1, c0)
tmpi(igs:ige) = xu(igs:ige)
call putvara ('longitude_V', iou, ig, ib, ic, tmpi, c1, c0)
ic(1) = jg
tmpj(jgs:jge) = yt(jgs:jge)
call putvara ('latitude', iou, jg, ib, ic, tmpj, c1, c0)
tmpj(jgs:jge) = yu(jgs:jge)
call putvara ('latitude_V', iou, jg, ib, ic, tmpj, c1, c0)
ic(1) = ig + 1
call edge_maker (1, xt_e, xt, dxt, xu, dxu, imt)
tmpie(igs:ige+1) = xt_e(igs:ige+1)
call putvara ('longitude_edges', iou, ig+1, ib, ic, tmpie, c1, c0)
ic(1) = jg + 1
call edge_maker (1, yt_e, yt, dyt, yu, dyu, jmt)
tmpje(jgs:jge+1) = yt_e(jgs:jge+1)
call putvara ('latitude_edges', iou, jg+1, ib, ic, tmpje, c1, c0)
ib(1) = 1
do l=1,nzmax
slev(l) = float(l)
slev_e(l) = float(l) - 0.5
enddo
slev_e(nzmax+1) = slev(nzmax) + 0.5
ic(1) = nzmax
call putvara ('slev', iou, nzmax, ib, ic, slev, 1., 0.)
ic(1) = nzmax+1
call putvara ('slev_edges', iou, nzmax+1, ib, ic, slev_e, 1., 0.)
ib(1) = 1
do l=1,ibmax
blev(l) = float(l)
blev_e(l) = float(l) - 0.5
enddo
blev_e(ibmax+1) = blev(ibmax) + 0.5
ic(1) = ibmax
call putvara ('blev', iou, ibmax, ib, ic, blev, 1., 0.)
ic(1) = ibmax+1
call putvara ('blev_edges', iou, ibmax+1, ib, ic, blev_e
&, 1., 0.)
!-----------------------------------------------------------------------
! write 2d data (x,y)
!-----------------------------------------------------------------------
ib(1) = 1
ic(1) = ile-ils+1
ib(2) = 1
ic(2) = jle-jls+1
ln = ic(1)*ic(2)
tmpij(ils:ile,jls:jle) = kmt(ils:ile,jls:jle)
call putvara ('kmt', iou, ln, ib, ic, tmpij, c1, c0)
tmpij(ils:ile,jls:jle) = map_sed(ils:ile,jls:jle)
call putvara ('map_sed', iou, ln, ib, ic, tmpij, c1, c0)
!-----------------------------------------------------------------------
! write 3d data (x,y,t)
!-----------------------------------------------------------------------
ib(1) = 1
ic(1) = ile-ils+1
ib(2) = 1
ic(2) = jle-jls+1
ib(3) = ntrec
ic(3) = 1
ln = ic(1)*ic(2)*ic(3)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibtemp)
call putvara ('sbc_btemp', iou, ln, ib, ic, tmpij, c1, c0)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibsalt)
call putvara ('sbc_bsalt', iou, ln, ib, ic, tmpij, c1, c0)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ircal)
call putvara ('sbc_rcal', iou, ln, ib, ic, tmpij, c1, c0)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,irorg)
call putvara ('sbc_rorg', iou, ln, ib, ic, tmpij, c1, c0)
# if defined O_carbon
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibdic)
call putvara ('sbc_bdic', iou, ln, ib, ic, tmpij, c1, c0)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibdicfx)
call putvara ('sbc_bdicfx', iou, ln, ib, ic, tmpij, c1, c0)
# endif
# if defined O_npzd_alk
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibalk)
call putvara ('sbc_balk', iou, ln, ib, ic, tmpij, c1, c0)
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibalkfx)
call putvara ('sbc_balkfx', iou, ln, ib, ic, tmpij, c1, c0)
# endif
# if defined O_npzd_o2
tmpij(ils:ile,jls:jle) = sbc(ils:ile,jls:jle,ibo2)
call putvara ('sbc_bo2', iou, ln, ib, ic, tmpij, c1, c0)
# endif
tmpip(1:ipmax) = zrct(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_zrct', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = water_z_p(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_water_z_p', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = k1(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_k1', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = k2(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_k2', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = k3(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_k3', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = csat(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_csat', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = rc(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_rc', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = ttrorg(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_ttrorg', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = ttrcal(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_ttrcal', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = sed_ml_mass(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_ml_mass', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = c_advect(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_c_advect', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = rain_cal_p(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_rain_cal_p', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = rain_org_p(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_rain_org_p', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = co3_p(1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_co3_p', iou, ln, ib, ic, tmpij, c1, c0)
!-----------------------------------------------------------------------
! write 4d data (x,y,z,t)
!-----------------------------------------------------------------------
ib(1) = 1
ic(1) = ile-ils+1
ib(2) = 1
ic(2) = jle-jls+1
ib(3) = 1
ic(3) = 1
ib(4) = 1
ic(4) = 1
ln = ic(1)*ic(2)*ic(3)*ic(4)
do l=1,nzmax
ib(3) = l
do m=1,3
write(a3, '(i1)') m
tmpip(1:ipmax) = carb(l,m,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_carb_'//trim(a3), iou, ln, ib, ic, tmpij
&, c1, c0)
tmpip(1:ipmax) = dcpls(l,m,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_dcpls_'//trim(a3), iou, ln, ib, ic, tmpij
&, c1, c0)
tmpip(1:ipmax) = dcmin(l,m,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_dcmin_'//trim(a3), iou, ln, ib, ic, tmpij
&, c1, c0)
enddo
tmpij(ils:ile,jls:jle) = map_sed(ils:ile,jls:jle)
call putvara ('map_sed', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = pore(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_pore', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = form(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_form', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = o2(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_o2', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = orggg(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_orggg', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = orgml(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_orgml', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = calgg(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_calgg', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = calml(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_calml', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = dopls(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_dopls', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = domin(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_domin', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = dbpls(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_dbpls', iou, ln, ib, ic, tmpij, c1, c0)
tmpip(1:ipmax) = dbmin(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_dbmin', iou, ln, ib, ic, tmpij, c1, c0)
enddo
do l=1,ibmax
ib(3) = l
tmpip(1:ipmax) = buried_mass(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_buried_mass', iou, ln, ib, ic, tmpij
&, c1, c0)
tmpip(1:ipmax) = buried_calfrac(l,1:ipmax)
call unloadsed (ipmax, tmpip, imt, jmt, map_sed, data)
tmpij(ils:ile,jls:jle) = data(ils:ile,jls:jle)
call putvara ('sed_buried_calfrac', iou, ln, ib, ic, tmpij
&, c1, c0)
enddo
deallocate (tmpip)
call rdstmp (nstamp, nyear, nmonth, nday, nhour, nmin, nsec)
nyear = year0 + accel_yr0 + (relyr - accel_yr0)*accel
call mkstmp (nstamp, nyear, nmonth, nday, nhour, nmin, nsec)
print*, '=> Sed restart written to ',trim(fname),' on ', nstamp
deallocate ( tmpij )
deallocate ( tmpi )
deallocate ( tmpj )
deallocate ( tmpie )
deallocate ( tmpje )
#endif
return
end