Skip to content
Snippets Groups Projects
Commit e38b77a2 authored by Lionel GUEZ's avatar Lionel GUEZ
Browse files

Rename `i_shpc` to `i_slice`

parent 4530d358
No related branches found
No related tags found
No related merge requests found
......@@ -23,7 +23,7 @@ contains
! Local:
integer ishape, ishape_first
integer:: i_shpc = huge(0)
integer:: i_slice = huge(0)
integer date_read
integer eddy_i ! eddy index
type(eddy) e
......@@ -31,23 +31,23 @@ contains
!---------------------------------------------------------------------
if (size(hshpc) == 1) then
i_shpc = 1
i_slice = 1
else
call hunt([ssm%d0, huge(0)], k, i_shpc)
call hunt([ssm%d0, huge(0)], k, i_slice)
end if
if (k == ssm(i_shpc)%d0) then
if (k == ssm(i_slice)%d0) then
ishape_first = 0
else
! {k > ssm(i_shpc)%d0}
ishape_first = ssm(i_shpc)%ishape_last(k - 1) + 1
! {k > ssm(i_slice)%d0}
ishape_first = ssm(i_slice)%ishape_last(k - 1) + 1
end if
s%number_extr = ssm(i_shpc)%ishape_last(k) - ishape_first + 1
s%number_extr = ssm(i_slice)%ishape_last(k) - ishape_first + 1
allocate(s%list(s%number_extr))
do ishape = ishape_first, ssm(i_shpc)%ishape_last(k)
call read_eddy(e, date_read, eddy_i, hshpc(i_shpc), ishape)
do ishape = ishape_first, ssm(i_slice)%ishape_last(k)
call read_eddy(e, date_read, eddy_i, hshpc(i_slice), ishape)
! Check that all the eddies have the same date index:
if (date_read /= k) then
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment