GIZMO ¤ÇÍ·¤ó¤Ç¤ß¤ë
|
|
|
- かずき とりこし
- 5 years ago
- Views:
Transcription
1 GIZMO February 21, 2019 GIZMO February 21, / 17
2 GIZMO Users Guide URL phopkins/site/gizmo_files/gizmo_documentation.html /home/hydro00/gizmo_hydro2018.tar.gz GIZMO_hydro2018/practice User Guide (scripts ) scripts/test_problems/ practice GIZMO February 21, / 17
3 (1) Intel XC50 module list PrgEnv-intel/6.0.4 gsl/2.4o intel-18.0 cray-hdf5/ fftw/ mdule load (or add) GIZMO February 21, / 17
4 (2) gizmo Config.sh cp Template-Config.sh Config.sh editor (vim ) Makefile.systype SYSTYPE= XC-intel GIZMO February 21, / 17
5 1: Config.sh (10 K) n = 0.5 cm erg M 1000 Users Guide n = cm 3, E = J. BOX PERIODIC SELFGRAVITY OFF EOS GAMMA=(5.0/3.0). MFM (HYDRO MESHLESS FINITE MASS) make Config.sh make GIZMO February 21, / 17
6 1: /work/hydro??/sedov GIZMO sedov.params sedov.param sedov.params InitCondFile /path_to_ic_dir/sedov_ics OutputDir. current directory GIZMO February 21, / 17
7 1: run.sh 1 MPI (40 ) GIZMO_hydro2018 batch_gizmo.sh run.sh #PBS -N sedov #PBS -l nodes=1 #PBS -q test-bp cd ${PBS O WORKDIR} aprun -n 40./GIZMO sedov.params > mpi output.txt qsub run.sh OutputDir snapshto???.hdf5 ( ) GIZMO February 21, / 17
8 1: yt /work/hydro??/ ipython ipython (or projection.py) import yt ds = yt.load( /xc-work/okamtotk/hydro/tests/sedov/snapshot_010.hdf5 ) snapshot plot = yt.projectionplot(ds, z, fields=[( gas, density )]) plot.save() projection plot display snapshot_011_projection_z_density.png plot = yt.projectionplot(ds, z, fields=[( gas, temperature )]), weight field=[( gas, density )] projection map. GIZMO February 21, / 17
9 1: yt python ipython (or sedov density.py) import h5py import numpy as np ds = h5py.file( /xc-work/okamtotk/hydro/tests/sedov/snapshot_010.hdf5, r ) # list(ds.kyes()) # Header PartType0 PartType0 (gas) list(ds[ PartType0 ].keys()) #. Header list(ds[ Header ].attrs) # simulation box boxlen = ds[ Header ].attrs[ BoxSize ] GIZMO February 21, / 17
10 1: matplotlib plot ipython (or sedov density.py) import matplotlib.pyplot as plt c = np.full(3, 0.5*boxlen) # pos = np.array(ds[ PartType0/Coordinates ]) r = np.sqrt( (pos[:,0] - c[0])**2 + (pos[:,1] - c[1])**2 + (pos[:,2] - c[2])**2) # plt.plot(r, ds[ PartType0/Density ],,, rasterized=true) plt.savefig( density.png ) exit ipython GIZMO February 21, / 17
11 2: SPH PSPH (Saito & Makino DISPH ) Config.sh HYDRO PRESSURE SPH Cullen & Dehnen 2010 make GIZMO cp sedov.params sedov psph.params sedov psph.params SnapshotFileBase snapshot snapshot psph run.sh sedov.params sedov psph.params qsub run.sh GIZMO February 21, / 17
12 2: sedov.txt sedov.py /home/okamtotk/hydro sedov.py python sedov.py sedov.png GIZMO February 21, / 17
13 : UnitLength in cm UnitMass in g UnitVelocity in cm per s e+21, 1.989e+43, kpc, M, km/s code units GIZMO February 21, / 17
14 3: Config.sh HYDRO PRESSURE SPH The Noh (Spherical Implosion) Problem (noh.params) BOX PERIODIC SELFGRAVITY OFF EOS GAMMA=(5.0/3.0) Kelvin Helmholtz Instabilities (kh mcnally 2d.params) BOX PERIODIC BOX SPATIAL DIMENSION=2 PREVENT PARTICLE MERGE SPLIT SELFGRAVITY OFF KERNEL FUNCTION=3 EOS GAMMA=(5.0/3.0) GIZMO February 21, / 17
15 3: The Blob Test (blob.params) BOX PERIODIC BOX LONG X=1 BOX LONG Y=1 BOX LONG Z=3 SELFGRAVITY OFF EOS GAMMA=(5.0/3.0) The Evrard (Spherical Collapse) Test (evrard.params) EOS GAMMA=(5.0/3.0) ADAPTIVE GRAVSOFT FORGAS bolob.params SofteningGas SofteningGasMaxPhys SPH GIZMO February 21, / 17
16 4: ( m i V tot = ρ i V tot BoxSize 3 SPH MFM, MFV i GIZMO February 21, / 17
17 5: SPH 3 SPH m j W(r i j, h i ) ρ j j 1 i W(r i j, h i ) { } 8 πh (u 1.0)u 2 for u = r 3 i j /h i 0.5 i { W(r i j, h i ) = } 8 πh 2.0(1.0 u) 3 for 0.5 u i 0 for u > 1.0 h i SmoothingLength GIZMO February 21, / 17
Python (Anaconda ) Anaconda 2 3 Python Python IDLE Python NumPy 6
Python (Anaconda ) 2017. 05. 30. 1 1 2 Anaconda 2 3 Python 3 3.1 Python.......................... 3 3.2 IDLE Python....................... 5 4 NumPy 6 5 matplotlib 7 5.1..................................
1 matplotlib matplotlib Python matplotlib numpy matplotlib Installing A 2 pyplot matplotlib 1 matplotlib.pyplot matplotlib.pyplot plt import import nu
Python Matplotlib 2016 ver.0.06 matplotlib python 2 3 (ffmpeg ) Excel matplotlib matplotlib doc PDF 2,800 python matplotlib matplotlib matplotlib Gallery Matplotlib Examples 1 matplotlib 2 2 pyplot 2 2.1
Python ( ) Anaconda 2 3 Python Python IDLE Python NumPy 6 5 matpl
Python ( ) 2017. 11. 21. 1 1 2 Anaconda 2 3 Python 3 3.1 Python.......................... 3 3.2 IDLE Python....................... 5 4 NumPy 6 5 matplotlib 7 5.1.................................. 7 5.2..................................
2
01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 2 3 01 02 03 4 04 05 06 5 07 08 09 6 10 11 12 7 13 14 15 8 16 17 18 9 19 20 21 10 22 23 24 11 FIELD MAP 12 13 http://www.pref.ishikawa.jp/shinrin/zei/index.html
Visual Python, Numpy, Matplotlib
Visual Python, Numpy, Matplotlib 1 / 38 Contents 1 2 Visual Python 3 Numpy Scipy 4 Scipy 5 Matplotlib 2 / 38 Contents 1 2 Visual Python 3 Numpy Scipy 4 Scipy 5 Matplotlib 3 / 38 3 Visual Python: 3D Numpy,
90 27 1 29 ( ) 27 4 26 ( ) 27 4 29 ( ) 4 27 3 30 ( ) 1 30 4 27 4 8 ( ) 16 ( ) 9 P2 3 301 27 4 14 ( ) 15 ( ) 9 5 27 4 21 ( ) 4 27 4 21 ( ) 8 30 5 27 4 23 ( ) 5 10 11 12 13 4 3 4 14 27 4 23 ( 27 4 25 ( )
格子スピン模型の計算科学2018＿実習
e-mail: t-okubo@phys.s.u-tokyo.ac.jp ALPS ALPS mc-09-snapshot Paraview Mac Paraview Mac ITC-LMS (https://itc-lms.ecc.u-tokyo.ac.jp/portal/login) 0618.pdf python python2.7 pyenv shell anaconda-4.0.0 ALPS
Visual Python, Numpy, Matplotlib
Visual Python, Numpy, Matplotlib 1 / 57 Contents 1 2 Visual Python 3 Numpy Scipy 4 Scipy 5 Matplotlib 2 / 57 Contents 1 2 Visual Python 3 Numpy Scipy 4 Scipy 5 Matplotlib 3 / 57 3 Visual Python: 3D Numpy,
1. ( ) SPH 1. 1: 2: 2.
SPH 2014/08/06 2014 1. ( ) 2. 3. 4. SPH 1. 1: 2: 2. 1. ( ) 2. 3. 4. MPI SIMD : 1. MPI : 2. ( ) MPI : ( ) : DRY (Don t Repeat Yourself) ( ) : DRY (Don t Repeat Yourself) ( 2014 7 ) (SPH MPS MLS ) ( : Tree,
appli_HPhi_install
2018/3/7 HΦ version 3.0.0 インストール手順書 (Linux 64 ビット版 ) 目次 1. アプリケーション概要...- 1-2. システム環境...- 1-3. 必要なツール ライブラリのインストール...- 1-1 cmake...- 2-2 numpy...- 3-4. アプリケーションのインストール...- 4-5. 動作確認の実施...- 5 - 本手順書は HΦ
170403_1_FABEX_2017会場案内マップ
OFFICIAL GUIDE MAP 07 4 4 34 000700 3 0 7 J a p a n M e a t I n d u s t r y F a i r 4 t h 4/3 4/4 4/ 4/34 4/3 4/4 4/4 4/4 4/4 4/4 4/4 4/3 4/4 4/3 4/4 4/ 4/4 4/4 4/4 P-0 4 3 4 4 :0 3:0 3:40 4:30 4:50 5:40
200708_LesHouches_02.ppt
Numerical Methods for Geodynamo Simulation Akira Kageyama Earth Simulator Center, JAMSTEC, Japan Part 2 Geodynamo Simulations in a Sphere or a Spherical Shell Outline 1. Various numerical methods used
In [168]: soup.find_all("label")
![In [168]: soup.find_all(label)](/thumbs/72/68060162.jpg "In [168]: soup.find_all(label)")
step 1 kakaku.com/bicycle/bicycle-battery/ web web Chrome cntl + Source Code and Copy cd 'kakaku_com_bicycle_bicycle-battery.html' In [166]: from bs4 import BeautifulSoup In [167]: html = open('kakaku_com_bicycle_bicycle-battery_2017.html')
nakao
Fortran+Python 4 Fortran, 2018 12 12 !2 Python!3 Python 2018 IEEE spectrum https://spectrum.ieee.org/static/interactive-the-top-programming-languages-2018!4 Python print("hello World!") if x == 10: print
Netcommunity SYSTEM X7000 IPコードレス電話機 取扱説明書
4 5 6 7 8 9 . 4 DS 0 4 5 4 4 4 5 5 6 7 8 9 0 4 5 6 7 8 9 4 5 6 4 0 4 4 4 4 5 6 7 8 9 40 4 4 4 4 44 45 4 6 7 5 46 47 4 5 6 48 49 50 5 4 5 4 5 6 5 5 6 4 54 4 5 6 7 55 5 6 4 56 4 5 6 57 4 5 6 7 58 4
.A. D.S
1999-1- .A. D.S 1996 2001 1999-2- -3- 1 p.16 17 18 19 2-4- 1-5- 1~2 1~2 2 5 1 34 2 10 3 2.6 2.85 3.05 2.9 2.9 3.16 4 7 9 9 17 9 25 10 3 10 8 10 17 10 18 10 22 11 29-6- 1 p.1-7- p.5-8- p.9 10 12 13-9- 2
Microsoft Word - Documento2
? 3 2 ... 19 20 20 2 1 1 20 40 20 2 9 20 20 2 40 : http://www.hirahidenobu.com/scw/entry 平 秀 信 : 250 400 2 4 58CD4 CD1: 2 2 CD2 URL CD3: CD4: CD CD4... 250 200 250 http://www.hirahidenobu.com/scw/entry
PowerPoint プレゼンテーション
斎藤貴之東京工業大学地球生命研究所 1 Contents 研究背景 ASURA: 並列 N 体 /SPHコード Density Independent SPH まとめ 2 3 内海洋輔(広島大学) NASA/JPL-Caltech/ESO/R. Hurt 4 Credit: S. Beckwith & the HUDF Working Group (STScI), HST, ESA, NASA 我々の宇宙
宿泊産業活性化のための実証実験
121 32 10 12 12 19 2 15 59 40 33 34 35 36 37 38 3637 20 39 12 19 OFF 2008/12/19 2008/12/25 3 1 1 72,000 2008/12/19 2008/12/26 2 1 1 36,000 2008/12/28 2009/1/5 2 1 1 24,000 2009/1/6 2009/1/16 3 1 1 25,200
901 902 2 40 5 786 30 2 2 100 10100200 903 904 2 3 2 12 905 6765 30 3 61016 1 10162532 253240 2 2 1 2 100 24 45 545 1 2 2 510 1515 1010 50 300 0 10 2942 560 2 1 1 2 24 15 2565 2 10 2942 560 3 3 56 03 18
1037 1038 2 40 5 876 30 2 2 100 10100200 1039 1040 2 3 2 12 1041 6765 30 1 1 2 2 1 2 100 24 45 545 1 2 2 510 1515 1010 50 300 0 10 2942 560 2 1 3 2 10 2942 560 3 61016 1 10162532 253240 1 2 24 15 2565
2015 I ( TA)
2015 I ( TA) Schrödinger PDE Python u(t, x) x t 2 u(x, t) = k u(t, x) t x2 k k = i h 2m Schrödinger h m 1 ψ(x, t) i h ( 1 ψ(x, t) = i h ) 2 ψ(x, t) t 2m x Cauchy x : x Fourier x x Fourier 2 u(x, t) = k
endo.PDF
MAP 18 19 20 21 3 1173 MAP 22 700800 106 3000 23 24 59 1984 358 358 399 25 12 8 1996 3 39 24 20 10 1998 9,000 1,400 5,200 250 12 26 4 1996 156 1.3 1990 27 28 29 8 606 290 250 30 11 24 8 1779 31 22 42 9
1
1 2 3 4 5 0% 20% 40% 60% 80% 100% 6 7 8 0% 20% 40% 60% 80% 100% 9 0% 20% 40% 60% 80% 100% 10 100% 90% 80% 70% 60% 50% 40% 30% 20% 10% 0% 2529 (n=17) 3034 35 (n=21) (n=17) 2529 (n=19) 3034 35 (n=34) (n=64)
1 6/13 2 6/20 3 6/27 4 7/4 5 7/11 6 7/18 N 7 7/25 Warshall-Floyd, Bellman-Ford, Dijkstra TSP DP, 8/1 2 / 36
3 2016 6 27 1 / 36 1 6/13 2 6/20 3 6/27 4 7/4 5 7/11 6 7/18 N 7 7/25 Warshall-Floyd, Bellman-Ford, Dijkstra TSP DP, 8/1 2 / 36 1 2 3 3 / 36 4 / 36 os.urandom(n) n >>> import os >>> r = os.urandom(4) #
PowerPoint Presentation
2010 KEK (Japan) (Japan) (Japan) Cheoun, Myun -ki Soongsil (Korea) Ryu,, Chung-Yoe Soongsil (Korea) 1. S.Reddy, M.Prakash and J.M. Lattimer, P.R.D58 #013009 (1998) Magnetar : ~ 10 15 G ~ 10 17 19 G (?)
RL_tutorial
)! " = $ % & ' "(& &*+ = ' " + %' "(- + %. ' "(. + γ γ=0! " = $ " γ=0.9! " = $ " + 0.9$ " + 0.81$ "+, + ! " #, % #! " #, % # + (( + #,- +. max 2 3! " #,-, % 4! " #, % # ) α ! " #, % ' ( )(#, %)!
untitled
186 17 100160250 1 10.1 55 2 18.5 6.9 100 38 17 3.2 17 8.4 45 3.9 53 1.6 22 7.3 100 2.3 31 3.4 47 OR OR 3 1.20.76 63.4 2.16 4 38,937101,118 17 17 17 5 1,765 1,424 854 794 108 839 628 173 389 339 57 6 18613
untitled
1. 3 14 2. 1 12 9 7.1 3. 5 10 17 8 5500 4. 6 11 5. 1 12 101977 1 21 45.31982.9.4 79.71996 / 1997 89.21983 41.01902 6. 7 5 10 2004 30 16.8 37.5 3.3 2004 10.0 7.5 37.0 2004 8. 2 7 9. 6 11 46 37 25 55 10.
81 34 6 5 6 6 5 5 6 1 5 2 6 176 414
5 6 5 6 5 6 5 4 4 6 5 5 5 413 175 81 34 6 5 6 6 5 5 6 1 5 2 6 176 414 1 5 5 12 3 4 5 1 1 2 3 2 1 415 177 81 34 2 3 4 5 3 1 2 4 5 1 1 1 4 178 416 65,000 35,000 6,000 2 4 417 179 81 34 4 10,000 5,000 750
... 6 1) 2) No. 01 02 03 04 05 06 07 08 09 10 11 12 No. 1 2 2 3 3cm 4 
untitled
21 14 487 2,322 2 7 48 4 15 ( 27) 14 3(1867) 3 () 1 2 3 ( 901923 ) 5 (1536) 3 4 5 6 7 8 ( ) () () 9 10 21 11 12 13 14 16 17 18 20 1 19 20 21 22 23 21 22 24 25 26 27 28 22 5 29 30cm 7.5m 1865 3 1820 5
1948 1907 4024 1925 14 19281929 30 111931 4 3 15 4 16 3 15 4 161933 813 1935 12 17 11 17 1938 1945 2010 14 221 1945 10 1946 11 1947 1048 1947 1949 24
15 4 16 1988 63 28 19314 29 3 15 4 16 19283 15294 16 1930 113132 3 15 4 16 33 13 35 12 3 15 4 16 1945 10 10 10 10 40 1948 1907 4024 1925 14 19281929 30 111931 4 3 15 4 16 3 15 4 161933 813 1935 12 17 11
<91E6825289F1938C966B95FA8ECB90FC88E397C38B5A8F708A778F7091E589EF8EC08D7388CF88F5837D836A83858341838B566572312E696E6464>
Microsoft Word - 入居のしおり.doc
1 1 2 2 2 3 2 4 3 5 3 6 3 7 3 8 4 1 7 2 7 3 7 4 8 5 9 6 9 7 10 8 10 9 11 10 11 11 11 12 12 13 13 1 14 2 17 3 18 4 19 5 20 6 22 (1) 24 (2) 24 (3) 24 (4) 24 (5) 24 (6) 25 (7) 25 (8) 25 (9) 25 1 29 (1) 29
( )
( ) () () 3 cm cm cm cm cm cm 1000 1500 50 500 1000 1000 1500 1000 10 50 300 1000 2000 1000 1500 50 10 1000 2000 300 50 1000 2000 1000 1500 50 10 1000 2000 300 30 10 300 1000 2000 1000 1500 1000
-26-
-25- -26- -27- -28- -29- -30- -31- -32- -33- -34- -35- -36- -37- -38- cm -39- -40- 1 2 3 4 4 3 2 1 5 5-41- -42- -43- -44- -45- -46- -47- -48- -49- -50- cm -51- -52- -53- -54- -55- -56- -57- -58- -59- -60-
-1- 4 1 2 4-2- -- 2 2 cm 0 80cm 2 80cm 80cm 80cm 50cm 80cm 50cm 6 80cm 100 50 50 cm 10 6 4 50cm 4 4 50cm -4- -5- cm 50cm 4 4 4 50cm 50cm 4 80cm 50cm 80cm 50cm 6 cm -6- 20 250cm 1 2 1 4 0cm 60cm cm cm 1
河川砂防技術基準・基本計画編.PDF
4 1 1 1 1 1 2 1 2.1 1 2.2 2 2.3 2 2.4 2 3 2 4 3 2 4 1 4 1.1 4 1.2 4 2 4 2.1 4 2.2 4 2.3 5 2.4 5 2.5 5 2.5.1 5 2.5.2 5 2.6 5 2.6.1 5 2.6.2 5 2.6.3 5 2.6.4 5 2.6.5 6 2.7 6 2.7.1 6 2.7.2 6 2.7.3 6 2.7.4
4 100g
100g 10 20 30 40 50 60 70 80 4 5 7 9 12 15 19 24 60 100 10 80 100 20 10 5 20 195 20-1- 60 60 15 100 60 100 15 15 15 100 15 15 60 100 10 60 100 100 15 10 10 60 100 15 10 15 10 5-2- 80 80 24 100 80 100 24
180 30 30 180 180 181 (3)(4) (3)(4)(2) 60 180 (1) (2) 20 (3)
12 12 72 (1) (2) (3) 12 (1) (2) (3) (1) (2) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (1) (2) 180 30 30 180 180 181 (3)(4) (3)(4)(2) 60 180 (1) (2) 20 (3) 30 16 (1) 31 (2) 31 (3) (1) (2) (3) (4) 30
untitled
1....1 2....2 2.1...2 2.2...2 3....14 3.1...14 3.2...14 4....15 4.1...15 4.2...18 4.3...21 4.4...23 4.5...26 5....27 5.1...27 5.2...35 5.3...54 5.4...64 5.5...75 6....79 6.1...79 6.2...85 6.3...94 6.4...
113 120cm 1120cm 3 10cm 900 500+240 10 1 2 3 5 4 5 3 8 6 3 8 6 7 6 8 4 4 4 4 23 23 5 5 7
裁定審議会における裁定の概要 (平成23年度)
23 23 23 4 24 3 10 11 12 13 14 () 1 23 7 21 23 12 14 (19 ) 30 1.876% 60 8 24 19 78 27 1 (10) 37 (3) 2 22 9 21 23 5 9 21 12 1 22 2 27 89 10 11 6 A B 3 21 12 1 12 10 10 12 5 1 9 1 2 61 ( 21 10 1 11 30 )
和県監査H15港湾.PDF
...1...1...1...1...1...1...1...1...2...2...2...3...3...3...5...5...10...11...12...13...13...13...14...14...14...14...14...14...15...15...15...15...15 ...16...17 14...17...18...18...19...21...23 2...25...27...27...28...28...28
2002 (1) (2) (3) (4) (5) (1) (2) (3) (4) (5) (1) (2) (3) (4) (5) (6) (7) (8) (1) (2) (3) (4) (1) (2) (3) (4) (5) (6) (7) (8) No 2,500 3 200 200 200 200 200 50 200 No, 3 1 2 00 No 2,500 200 7 2,000 7
untitled
() () () () () ( ) () ( ) () ( ) () 2 () () 2 () () ( ) () () () 2 () () 2 3 ( ) () ( ) 2 3 4 () () 2 3 4 () () ( )( ) ( ) 2 ( ) 3 () () 2 3 () () 2 3 () () () () () () () () (( ) ( ) (( ))( )( ) ) 2 3
総研大恒星進化概要.dvi
The Structure and Evolution of Stars I. Basic Equations. M r r =4πr2 ρ () P r = GM rρ. r 2 (2) r: M r : P and ρ: G: M r Lagrange r = M r 4πr 2 rho ( ) P = GM r M r 4πr. 4 (2 ) s(ρ, P ) s(ρ, P ) r L r T
研究紀要 第５号
3 4 5 6 7 8 a s d f a 9 10 s d a 11 12 s d f g 13 h j a d s 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 a 35 s 36 a 37 s 38 a 39 s 40 a 41 s d 42 f 43 44 46 47 48 49 50 a s d as d 51
17. (1) 18. (1) 19. (1) 20. (1) 21. (1) (3) 22. (1) (3) 23. (1) (3) (1) (3) 25. (1) (3) 26. (1) 27. (1) (3) 28. (1) 29. (1) 2
1. (1) 2. 2 (1) 4. (1) 5. (1) 6. (1) 7. (1) 8. (1) 9. (1) 10. (1) 11. (1) 12. (1) 13. (1) 14. (1) 15. (1) (3) 16. (1) 1 17. (1) 18. (1) 19. (1) 20. (1) 21. (1) (3) 22. (1) (3) 23. (1) (3) 24. 1 (1) (3)
untitled
0120-888-089 0120-919-498 160-83381-26-1TEL. 03-3349-3111 URL http://www.sjnk.co.jp/ 160-83381-26-1 TEL. 03-3349-3111 URL http://www.sompo-japan.co.jp/ 100-89653-7-3 TEL. 0120-919-498 URL http://www.nipponkoa.co.jp/
表紙.PDF
2 3 4 2001 2010 33 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 H18 21 22 23 24 25 26 27 28 29 30 31 32 33 300 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 410-8601 055-934-4747 Fax 055-933-1412
