Compare commits

...

2 Commits

Author SHA1 Message Date
sfilippone 926409bf0e Fix matrix generation 2 months ago
sfilippone 6c9337c728 Fix sample matrix generators 2 months ago

@ -467,8 +467,8 @@ contains
! compute gridpoint coordinates ! compute gridpoint coordinates
call idx2ijk(ix,iy,glob_row,idim,idim) call idx2ijk(ix,iy,glob_row,idim,idim)
! x, y coordinates ! x, y coordinates
x = (ix-1)*deltah x = (ix)*deltah
y = (iy-1)*deltah y = (iy)*deltah
zt(k) = f_(x,y) zt(k) = f_(x,y)
! internal point: build discretization ! internal point: build discretization

@ -496,9 +496,9 @@ contains
! compute gridpoint coordinates ! compute gridpoint coordinates
call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim) call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim)
! x, y, z coordinates ! x, y, z coordinates
x = (ix-1)*deltah x = (ix)*deltah
y = (iy-1)*deltah y = (iy)*deltah
z = (iz-1)*deltah z = (iz)*deltah
zt(k) = f_(x,y,z) zt(k) = f_(x,y,z)
! internal point: build discretization ! internal point: build discretization
! !

@ -467,8 +467,8 @@ contains
! compute gridpoint coordinates ! compute gridpoint coordinates
call idx2ijk(ix,iy,glob_row,idim,idim) call idx2ijk(ix,iy,glob_row,idim,idim)
! x, y coordinates ! x, y coordinates
x = (ix-1)*deltah x = (ix)*deltah
y = (iy-1)*deltah y = (iy)*deltah
zt(k) = f_(x,y) zt(k) = f_(x,y)
! internal point: build discretization ! internal point: build discretization

@ -496,9 +496,9 @@ contains
! compute gridpoint coordinates ! compute gridpoint coordinates
call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim) call idx2ijk(ix,iy,iz,glob_row,idim,idim,idim)
! x, y, z coordinates ! x, y, z coordinates
x = (ix-1)*deltah x = (ix)*deltah
y = (iy-1)*deltah y = (iy)*deltah
z = (iz-1)*deltah z = (iz)*deltah
zt(k) = f_(x,y,z) zt(k) = f_(x,y,z)
! internal point: build discretization ! internal point: build discretization
! !

Loading…
Cancel
Save