! defaultmesh , routine for old VSV grid compatibility
! Copyright (C) 1999,2010 Gauthier Delerce
!
! This program 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 2, or (at your option)
! any later version.
!
! This program 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 this program; if not, write to the Free Software Foundation,
! Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
!
! Gauthier Delerce <gauthier@delerce.fr>
subroutine defaultmesh
use modvsv
use modgrille
implicit none
integer::num,ix,iy
!real xp(40000),yp(40000)
print*,ibsx,ibsy,isrx,isry,new_default,dxshift,dyshift
if (new_default) then
num=0
do iy=(isry+1)/2,(ny-(isry-1)/2),dyshift
do ix=(isrx+1)/2,(nx-(isrx-1)/2),dxshift
num=num+1
! print*,num
! xp(num)=float(ix)
! yp(num)=float(iy)
end do
end do
allocate(grille(num,7))
mg=num
ng=1
num=0
do iy=(isry+1)/2,(ny-(isry-1)/2),dyshift
do ix=(isrx+1)/2,(nx-(isrx-1)/2),dxshift
num=num+1
! print*,num
grille(num,1)=float(ix)
grille(num,2)=float(iy)
end do
end do
!goto 124
!print*,'writing grid: ',trim(fgrille)
!open(unit=16,file=fgrille,form='formatted',status='replace',err=10)
!write(16,*) 1,num
!do ix=1,num
! write(16,200) grille(ix,1),grille(ix,2)
!end do
!close(16)
print*,'default grid generated '
else
open(unit=16,file=trim(adjustl(fgrille)),form='formatted',status='old',ACTION='READ' )
Print*,'open user grid',fgrille
read(16,*) mg,ng
! print*,mg,ng
num=mg*ng
allocate(grille(num,7))
do ix=1,num
read(16,*) grille(ix,1),grille(ix,2)
end do
close(16)
endif
goto 124
!10 print*,'can not open ',fgrille
!* formats
200 format(f7.2,2x,f7.2)
124 return
end