Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
xe=`awk {'print $1'} tmp/domain_input.gdat`
ye=`awk {'print $2'} tmp/domain_input.gdat`
south=`awk {'print $3'} tmp/domain_input.gdat`
est=`awk {'print $4'} tmp/domain_input.gdat`
north=`awk {'print $5'} tmp/domain_input.gdat`
west=`awk {'print $6'} tmp/domain_input.gdat`
reso=`awk {'print 0.01*$7'} tmp/domain_input.gdat`
reso_dim=`awk {'print $7'} tmp/domain_input.gdat`
xmin=`awk {'print $8'} tmp/domain_input.gdat`
xmax=`awk {'print $9'} tmp/domain_input.gdat`
ymin=`awk {'print $10'} tmp/domain_input.gdat`
ymax=`awk {'print $11'} tmp/domain_input.gdat`
rm tmp/domain_input.gdat
echo "x_source = "$xe
echo "y_source = "$ye
echo "resolution = "$reso_dim
echo "xmin = "$xmin
echo "xmax = "$xmax
echo "ymin = "$ymin
echo "ymax = "$ymax
awk '{if (($1 >= '$xmin')&&($1<='$xmax')&&($2 >= '$ymin')&&($2<='$ymax')) {printf "%0.6f %0.6f %0.6f \n",0.01*($1-'$xmin'),0.01*($2-'$ymin'),0.01*$3}}' $1 > DEMadim.xyz
nxny=`wc -l DEMadim.xyz | awk '{print $1}'`
echo $nxny
x0=`head -n 1 DEMadim.xyz | awk '{print $2}'`
nx=`awk '{if ($2 != '$x0') {nl_val=NR; exit}} END{print nl_val-1;}' DEMadim.xyz`
ny=$((($nxny)/$nx))
echo $nx
echo $ny
xmin=`head -n 1 DEMadim.xyz | awk '{print $1}'`
xmax=`tail -n 1 DEMadim.xyz | awk '{print $1}'`
ymin=`tail -n 1 DEMadim.xyz | awk '{print $2}'`
ymax=`head -n 1 DEMadim.xyz | awk '{print $2}'`
$2/mkgeo_grid_2d_var -t $(($nx-1)) $(($ny-1)) -a $xmin -b $xmax -c $ymin -d $ymax -order-ji -v4 > DEM.geo
geo -upgrade DEM.geo > tmp.geo
mv tmp.geo DEM.geo
> DEM.field
echo "field" >> DEM.field
echo "1 $nxny" >> DEM.field
echo "DEM.geo" >> DEM.field
echo "P1" >> DEM.field
awk '{print $3}' DEMadim.xyz > z_tmp
cat DEM.field z_tmp > DEM_tmp
mv DEM_tmp DEM.field
rm z_tmp
rm DEMadim.xyz
> grille.dmn
echo "EdgeDomainNames" >> grille.dmn
echo "4" >> grille.dmn
echo "bottom" >> grille.dmn
echo "right" >> grille.dmn
echo "top" >> grille.dmn
echo "left" >> grille.dmn
> grille.bamgcad
echo "MeshVersionFormatted" >> grille.bamgcad
echo " 0" >> grille.bamgcad
echo " Dimension" >> grille.bamgcad
echo " 2" >> grille.bamgcad
echo " Vertices" >> grille.bamgcad
echo " 4" >> grille.bamgcad
echo " $xmin $ymin 1" >> grille.bamgcad
echo " $xmax $ymin 2" >> grille.bamgcad
echo " $xmax $ymax 3" >> grille.bamgcad
echo " $xmin $ymax 4" >> grille.bamgcad
echo " Edges" >> grille.bamgcad
echo " 4" >> grille.bamgcad
echo " 1 2 101" >> grille.bamgcad
echo " 2 3 102" >> grille.bamgcad
echo " 3 4 103" >> grille.bamgcad
echo " 4 1 104" >> grille.bamgcad
echo " hVertices" >> grille.bamgcad
echo " $reso $reso $reso $reso" >> grille.bamgcad
bamg -g grille.bamgcad -o grille.bamg -nbv 5000000
bamg2geo grille.bamg grille.dmn > grille.geo