Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 154 trang
THÔNG TIN TÀI LIỆU
Nội dung
I H C QU C GIA THÀNH PH H CHÍ MINH I H C BÁCH KHOA KHOA K THU T XÂY D NG NGUY N TH MINH LÝ N VÙNG B N BI N VI T NAM Chuyên ngành: Xây d ng cơng trình bi n LU TP H Chí Minh, tháng 12 2013 NG C HOÀN THÀNH T I I H C BÁCH KHOA Cán b -HCM ng d n khoa h c : TS NGUY N DANH TH O Cán b ch m nh n xét : C Cán b ch m nh n xét : Lu ngày 15 tháng 01 Thành ph n H i NG c b o v t i Tr 2014 ng ng i h c Bách Khoa, HQG Tp HCM nh giá lu n v n th c s g m: PSG TS LÊ SONG GIANG C NG TS PH M NG C TS NGUY N DANH TH O Xác nh n c a Ch t ch H i ngành sau lu n v n ã CH T CH H NG ng h giá LV c s a ch a (n u có) ng Khoa qu n lý chuyên NG KHOA I H C BÁCH KHOA TP.HCM C NG HÒA XÃ H I CH IH C T NAM c L p - T Do - H nh Phúc Tp.HCM, ngày NHI M V LU H tên h c viên: NGUY N TH MINH LÝ Phái: N sinh: 31 / 07 / 1989 An Giang Chun ngành: Cơng Trình Bi n I- MSHV: 12020445 TÀI: N N VÙNG B BI N VI T NAM II- NHI M V VÀ N I DUNG: - Nghiên c u lý thuy t thành t o lan truy n sóng th n - Áp d ng tính tốn cho 13 k ch b n t - t gây sóng th n c a B Tài nguyên Môi ng cho vùng b bi n Vi t Nam Kh o sát th c tr ng phòng tránh sóng th n Vi t Nam III- NGÀY GIAO NHI M V : 24/06/2013 IV- NGÀY HOÀN THÀNH NHI M V : 22/11/2013 V- CÁN B NG D N: TS NGUY N DANH TH O N ch CÁN B NG D N ng Chuyên Ngành thông qua CH NHI M B MÔN QU N LÝ CHUYÊN NGÀNH TS NGUY N DANH TH O TS TR N THU TÂM L IC Trong su t trình làm lu cs Danh Th o ng d n nhi t nh t kh Em xin chân thành c a th y Nguy n a y Nguy n Danh Th o hồn thành t t lu u ki n thu n l i nh i nh ng kinh nghi m ki n th c quý báo mà th em t cho em s hành trang cho em sau t t nghi p Em r t c n th y Tr n Thu Tâm tồn th th y ngành C ng bi n khoa K thu t Xây d ng ti ng d n em lu ng d em q góp khơng nh giúp em hồn thành t t lu il ic góp ý ki i h c Bách Khoa TP.H Chí Minh, khơng tr c n tình d y b trình h c t p Cơng trình n b ng viên, em hoàn ch nh lu c bi t, xin bày t lòng bi ng su t th M g ng nhi ch nên lu n ba m ng h ng viên, quan c i gian làm lu i ng n, ki n th c h n a em khơng tránh kh i sai sót, mong q th y cô ch d n thêm Cu i cùng, em xin chúc ng l n c a c c, m ng g t hái nhi u thành công, x ng danh m t ng mà t t c h c Em kính chúc ba m , quý th y cô, b n bè l i chúc s c kh e, thành công Em xin chân thành c TP.H Chí Minh, ngày 02 tháng 12 H c viên th c hi n NGUY N TH MINH LÝ LVThS - n n vùng b bi n Vi t Nam M CL C TÓM T i DANH M C B NG BI U ii DANH M C HÌNH ii KÝ HI U S v D NG GI I THI 1.1 1.2 TV TÀI : .1 M M VI NGHIÊN C U: 1.2.1 M u: 1.2.2 Ph m vi nghiên c u: T NG QUAN V 2.1 TÀI T NG QUAN SÓNG TH N: 2.1.1 Khái ni m sóng th n: 2.1.2 Nguyên nhân hình thành sóng th n: 2.1.3 m c a sóng th n: 2.1.4 M t s tr n sóng th n l ch s th gi i: 2.2 CÁC NGU T CÓ KH N TRÊN BI N 2.2.1 2.2.2 11 t gãy vùng rìa Tây Tây b c Bi 2.2.3 2.3 13 14 TÌNH HÌNH NGHIÊN C U SÓNG TH C VÀ TH GI I: 15 2.3.1 M t s nghiên c c v sóng th n: 15 2.3.2 M t s nghiên c u th gi i v sóng th n: 21 M cL c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 LVThS CH n NG LÝ THUY T VÀ PH NG PHÁP TÍNH TỐN 24 3.1 LÝ THUY 3.2 LÝ THUY T SÓNG: 29 3.3 PH T: 24 NG PHÁP SAI PHÂN LEAP-FROG: 32 3.3.1 Công th c sai phân theo Taylor: 32 3.3.2 Công th c sai phân leap 3.3.3 Công th 3.4 CH n vùng b bi n Vi t Nam frog: 33 c thi t l p theo ph U KI ng pháp sai phân leap frog: 39 U KI N BIÊN: 41 3.4.1 u ki 3.4.2 u ki n biên: 42 NG CH u: 41 NG TRÌNH MƠ PH NG S PHÁT SINH VÀ LAN TRUY N SÓNG TH N 44 4.1 S MÔ PH NG SÓNG TH N: 44 4.2 S LI A HÌNH: 45 4.2.1 Ph m vi tính tốn lan truy n sóng th n 4.2.2 X lý s li 4.3 bi t Nam: 45 a hình: 46 MƠ HÌNH D CH CHUY N PHAY: 48 4.3.1 S li 4.3.2 Ch 4.3.3 S 4.3.4 K t qu tính tốn: 50 4.4 u vào: 48 ng trình (subroutine): 48 tính d ch chuy n phay: 49 MÔ PH NG LAN TRUY N SÓNG TH N: 50 4.4.1 S li 4.4.2 Ch 4.4.3 Ch u vào: 50 ng trình (subroutine): 50 ng trình mơ ph ng cho khu v c l n v i b M cL c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 cl i: 51 LVThS CH n NG K T QU TÍNH TỐN .52 5.1 CH K T QU T DO D CH CHUY N PHAY: 52 5.1.1 i hút chìm Manila: 52 5.1.2 i hút chìm RyuKyu: 53 5.1.3 t gãy Nam Trung B : 55 5.1.4 t gãy H i Nam: 55 5.2 CH n vùng b bi n Vi t Nam K T QU LAN TRUY N SÓNG TH N: 56 5.2.1 K ch b n 4: 56 5.2.2 K ch b n 5: 61 5.2.3 K ch b n 10: 69 5.2.4 K ch b n 11: 79 5.2.5 K ch b n 13: 81 5.2.6 Th i gian lan truy n t ngu nb 5.2.7 Chi u cao sóng l n nh ng v i 13 k ch b nb ng v i 13 k ch b NG TH C TR NG PHỊNG CH NG SĨNG TH N t: 87 t: 89 VI T NAM 91 6.1 H TH NG C NH BÁO S M SÓNG TH N: 91 6.2 DI N T P PHÒNG CH NG SÓNG TH N: 96 NG K T LU N VÀ KI N NGH 99 7.1 K T LU N: 99 7.2 KI N NGH : 100 TÀI LI U THAM KH 101 PH L 104 M cL c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 LVThS - n n vùng b bi n Vi t Nam TÓM T T LU Trong nh m h a t nhiên xu t hi nh ng th m h a gây nguy hi gi i M t n khu v c ven b bi n sóng th n N m giáp bi n bi n Vi t Nam có nguy c ti m tàng gánh ch u th m h a sóng th n, v y c n có s quan tâm c sóng th n c M c tiêu c a Lu n b bi n Vi t Nam t k t qu mơ ph ng s hình thành lan truy n sóng th n b bi d ng lý thuy c ta theo k ch b c nông v n d sóng th i gian lan truy t gây sóng th n s y cóc K t qu phân b chi u cao n b bi n cho th y r ng hi m h a sóng th n t i Vi n có thêm nhi u nghiên c u c th ABSTRACT The number of natural disasters around the world has trended to increase in recent years In coast area, tsunami is one among the main natural disasters Facing the East Sea, the Vietnam coast has been recognized to have the potential to generate devastating tsunami, thus it requires a consideration of government This thesis aims to evaluate possibilities of tsunami occurrence at Vietnam coast according to the results of the simulation computation of generation and propagation of tsunami in the East Sea under tsunami source scenarios using numerical model which is based on shallow water equation with leap-frog scheme to compute propagation wave The numerical results show the distribution of the maximum tsunami height and arrival time along the Vietnam coast The results reveal that tsunami hazards in Vietnam are significant and it is necessary to have more careful researches in detail Tóm t t lu HVTH: Nguy n Th Minh Lý - MSHV: 12020445 i LVThS - n n vùng b bi n Vi t Nam DANH M C B NG BI U B ng 1.1 Các thơng s t sóng th n theo k ch b n B ng 2.1 Tóm t t k t qu sóng th n 13 k ch b n 18 B ng 5.1 Th i gian lan truy n sóng th n ven b bi n Vi t Nam theo 13 k ch b n 87 DANH M C HÌNH V Hình 1.1 V t c a 13 k ch b Hình 2.1 B t gây sóng th n cho bi m ng ki n t o c t Nam t Hình 2.2 Các ki u ranh gi i m ng Hình 2.3 Quá trình lan truy n sóng th n t ngu n t i b (minh h a) Hình 2.4 B ki n t o khu v Hình 2.5 Các ngu Hình 2.6 t có kh n Bi i hút chìm RyuKyu Hình 2.7 Phân b Hình 2.8 t Nam 10 11 sâu ngu n Sulu 12 i hút chìm Manila 12 Hình 2.9 Vùng ngu n B c Bi 13 Hình 2.10 Vùng ngu n Tây Bi 14 Hình 2.11 Phân b sâu c i Palawan 14 Hình 2.12 S d ch chuy n phay gây sóng th n (minh h a, ngu n USGS) 21 Hình 2.13 B phân b th m h a sóng th n th gi i 23 Hình 3.1 Hình h c phay h t Hình 3.2 T phay t 25 nh phay 26 Hình 3.3 Sai phân trung tâm 32 Hình 3.4 S s p x Hình 3.5 S s p x m tính tốn ph m cho thông s ng pháp sai phân 33 i l u 35 Danh m c b ng - hình HVTH: Nguy n Th Minh Lý - MSHV: 12020445 ii LVThS - n n vùng b bi n Vi t Nam Hình 4.1 Khu v c tính tốn lan truy n sóng th n .45 Hình 4.2 Hình minh h a v trí phay vùng tính tốn 48 Hình 5.1 Bi n d ng m t bi n phay d ch chuy n ng v i k ch b n t i vùng ngu n Manila 53 Hình 5.2 Bi n d ng m t bi n phay d ch chuy n ng v i k ch b n t i vùng ngu n RyuKyu 54 Hình 5.3 Bi n d ng m t bi n phay d ch chuy n ng v i k ch b n t i vùng ngu n Nam Trung B 55 Hình 5.4 Bi n d ng m t bi n phay d ch chuy n ng v i k ch b n t i vùng ngu n H i Nam 55 Hình 5.5 K t qu th i gian lan truy n sóng th n ng v i k ch b n 58 Hình 5.6 ng m c n c ven b t i m t s a danh theo k ch b n 60 Hình 5.7 K t qu th i gian lan truy n sóng th n ng v i k ch b n 63 Hình 5.8 ng m c n c ven b t i m t s a danh theo k ch b n 65 Hình 5.9 ng m c n c t i Phú Yên ng v i k ch b n 1-5 68 Hình 5.10 K t qu th i gian lan truy n sóng th n ng v i k ch b n 10 71 Hình 5.11 ng m c n c ven b t i m t s a danh theo k ch b n 10 73 Hình 5.12 ng m c n c t i Phú Yên ng v i k ch b n 6-10 75 Hình 5.13 K t qu th i gian lan truy n sóng th n ng v i k ch b n 11 78 Hình 5.14 ng m c n c ven b t i m t s a danh theo k ch b n 11 80 Hình 5.15 K t qu th i gian lan truy n sóng th n ng v i k ch b n 13 83 Hình 5.16 ng m c n c ven b t i m t s a danh theo k ch b n 13 85 Hình 5.17 ng m c n c t i Phú Yên ng v i k ch b n 12-13 86 Hình 5.18 Chi u cao sóng th n l n nh n ven b Hình 6.1 V trí tr m c nh báo sóng th u tiên ng v i 13 k ch b n 89 Vi t Nam t ng ch p t Google Earth 93 Danh m c b ng - hình HVTH: Nguy n Th Minh Lý - MSHV: 12020445 iii LVThS 1.2 C C C C C C n n vùng b bi n Vi t Nam Tính tốn bi n d ng m t bi n s d ch chuy n phay gây ra: ** PROGRAM FOR CALCULATING THE DISPLACEMENT OF GRAND ** ** PROGRAM NAME = TSUNA(SYOKI*) ** ** BY WADA ** In VietNam Sea, revised MinhLy ** *********************************************************** PARAMETER (NI=2880,NJ=2400,ITA=960,JTA=800) C DIMENSION ZDL(NI,NJ),ET(NJ),ZDLA(ITA,JTA) C C C C ! DATA YOMIKOMI ! DATA # XL, YL, XF, YF # / 2592.0, 2160.0, 745.2, 918.0 / DATA # AFL, L, W, ZD, GTHET, SLP, # / 93.0, 646.0, 105.0, 27.0, 15.0, 90.0, UO 0.01749 / IWRIT=30 SLP=SLP*0.017453 US=UO*cos(SLP) UD=UO*sin(SLP) DX=XL/FLOAT(NI) DY=YL/FLOAT(NJ) GTHET=GTHET*0.017453 ALF=ALF*0.017453 FL=L/2 C OPEN(10,FILE='SHOKI_VNS_K5_A.DAT') OPEN(40,FILE='SHOKI_CHART_K5.DAT') C DO 100 I=1,NI XP=DX*(FLOAT(I-1)+0.5) C C C DO 200 J=1,NJ YP=DY*(FLOAT(J-1)+0.5) ! ZAHYO HENKAN ! CALL PTROTE(XP,YP,XF,YF,ALF,X1,Y1) ! HENI NO KEISAN ! CALL EARTH(X1,-Y1,ZDSP,FL,W,ZD,GTHET,US,UD) C ZDL(I,J)=ZDSP ET(J)=ZDSP C 200 CONTINUE C C C ! ZU NI KAKU ! CALL CHART(NJ,ET,12,2.0,-1.0) C 100 CONTINUE C DO 300 I=1,IT DO 300 J=1,JT II=2+(I-1)*3 JJ=2+(J-1)*3 ZDLA(I,J)=ZDL(II,JJ) 300 CONTINUE Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 127 LVThS - n n vùng b bi n Vi t Nam C WRITE(10,600)((ZDLA(I,J),J=1,JTA),I=1,ITA) 600 FORMAT(800F8.3) C STOP END C -C ! ZAHYO HENKAN ! SUBROUTINE PTROTE(XP,YP,XF,YF,ALF,X1,Y1) C X1=(XP-XF)*COS(ALF)+(YP-YF)*SIN(ALF) Y1=-(XP-XF)*SIN(ALF)+(YP-YF)*COS(ALF) C RETURN END C SUBROUTINE EARTH(XX1,XX2,ZDSP,FL,W,ZD,GTHET,US,UD) C DIMENSION STD(4),DPS(4) C D1=ZD/SIN(GTHET) D2=D1+W C CALL STSLIP(XX1,XX2,GTHET,US,UD,D2,FL,STD(1),DPS(1)) CALL STSLIP(XX1,XX2,GTHET,US,UD,D2,-FL,STD(2),DPS(2)) CALL STSLIP(XX1,XX2,GTHET,US,UD,D1,FL,STD(3),DPS(3)) CALL STSLIP(XX1,XX2,GTHET,US,UD,D1,-FL,STD(4),DPS(4)) C ZDDP=(STD(1)-STD(2)-STD(3)+STD(4)) #+(DPS(1)-DPS(2)-DPS(3)+DPS(4)) ZDSP=-1000.*ZDDP C RETURN END C -C ! HENI NO KEISAN ! C SUBROUTINE STSLIP(X1,X2,GTHET,US,UD,GZ,GZ1,US3,UD3) C DATA PAI/3.141593/ C CT=COS(GTHET) ST=SIN(GTHET) TA=TAN(GTHET) R2=X2*ST R3=X2*CT Q2=X2*ST Q3=-X2*CT R=(X1-GZ1)**2+R2**2+(R3-GZ)**2 R=SQRT(R) Q=(X1-GZ1)**2+Q2**2+(Q3+GZ)**2 Q=SQRT(Q) GZ2=GZ*CT GZ3=GZ*ST H=SQRT(Q2**2+(Q3+GZ)**2) C IF(US.EQ.0.0) GO TO 300 C C ! YOKOSUBERI DANSOU ! C Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 128 LVThS - n n vùng b bi n Vi t Nam T1=(ALOG(R+R3-GZ)+(1+3*TA**2)*ALOG(Q+Q3+GZ)-3*TA/CT*ALOG(Q+GZ3)) #*CT T2=2*R2*ST/R+2*ST*(Q2+X2*ST)/Q-2*R2**2*CT/R/(R+R3-GZ) T3=-2*(Q2+X2*ST)*Q3*ST/Q/(Q+Q3+GZ) US3=US*(T1+T2+T3)/PAI/12 GO TO 500 300 US3=0.0 500 CONTINUE C ! TATEZURE DANSOU ! C1=-2*GZ3/R/(R+X1-GZ1)-4*GZ3/Q/(Q+X1-GZ1) A1=(X1-GZ1)*(X2-GZ2) A2=(H+GZ3)*(Q+H) A3=(X1-GZ1)*(R3-GZ) A4=R2*R A5=(X1-GZ1)*(Q3+GZ) A6=Q2*Q C2=-6*ATAN2(A1,A2) C3=3*ATAN2(A3,A4)-6*ATAN2(A5,A6) E1=ST*((X2-GZ2)*C1+C2+C3) C4=ALOG(R+X1-GZ1)-ALOG(Q+X1-GZ1)-2*GZ3**2/R/(R+X1-GZ1) C5=-4*GZ3**2/Q/(Q+X1-GZ1) E2=CT*(C4+C5) UD3=UD*(E1+E2)/PAI/12 C RETURN END C C ! ZU NI KAKU ! SUBROUTINE CHART(N,U,NCT,UMAX,UMIN) C DIMENSION U(N),ST(20) C CHARACTER*1 C(3000),A(16)*1 DATA (A(I),I=1,16)/'1','2','3','4','5','6','7','8','9', #'A','B','C','D','E','F','G'/ C IF(N.GE.3000) RETURN IF(NCT.GE.17) RETURN C NK=NCT+1 DH=(UMAX-UMIN)/FLOAT(NCT) C DO 100 K=1,NK ST(K)=UMIN+DH*FLOAT(K-1) 100 CONTINUE C DO 200 I=1,N C IF(U(I).LT.ST(1)) C(I)='*' IF(U(I).GE.ST(NK)) C(I)='+' C DO 500 JJ=1,NCT IF(U(I).GE.ST(JJ).AND.U(I).LT.ST(JJ+1)) C(I)=A(JJ) 500 CONTINUE C 200 CONTINUE C WRITE(40,300) (C(I),I=1,N) 300 FORMAT(1H ,3000A1) C 20 CONTINUE Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 129 LVThS - n n vùng b bi n Vi t Nam C RETURN END 1.3 Tính tốn lan truy n sóng th n b bi n Vi t Nam: C C C C C C C C C C C C ** TSUNAMINO KEISAN ** ** FILE NAME = TSUNA.FORT(TSUA-H) ** ** MODIFIED BY WADA ** ** ** ******************************************************** ** KR :SUISHIN.DATA ** ** KRI:IB-MAP.DATA ** ** IZR:SYOKI-SUII.DATA ** ** Revised by Minhlies for VietNam Sea ** ******************************************************** CHARACTER*100 fname INTEGER fnumber C REAL MA,NA,MZA C PARAMETER (IT=960,JT=800,NPL=14) C DIMENSION ZA(IT,JT,2),MA(IT,JT,2),NA(IT,JT,2) DIMENSION DZA(IT,JT,2),DMA(IT,JT,2),DNA(IT,JT,2) DIMENSION HZA(IT,JT),HMA(IT,JT),HNA(IT,JT) C DIMENSION IBA(IT,JT),VA(IT,JT,2),MZA(IT,JT) DIMENSION LPI(NPL),LPJ(NPL) C DATA INIZ,KR,KRI,IZR,KV / 2,10,12,16,31 / C OPEN(10,FILE='VietNam_A.DAT') OPEN(12,FILE='IBMAP_A.DAT') OPEN(16,FILE='SHOKI_VNS_K5_A.DAT') OPEN(20,FILE='CH.DAT') OPEN(41,FILE='POINT_Z.DAT') OPEN(42,FILE='POINT_M.DAT') OPEN(43,FILE='POINT_N.DAT') OPEN(44,FILE='POINT_V.DAT') OPEN(51,FILE='MAXZ.DAT') OPEN(61,FILE='MAXV.DAT') C READ(20,100) (LPI(I),LPJ(I),I=1,NPL) 100 FORMAT(14I5) C ITM=IT-1 JTM=JT-1 DXA=2700.0 DTA=3.0 KTA=3 RA=DTA/DXA C fnumber=0 C KE=28800 C C ! SUISHIN-NYURYOKU ! CALL DEPTHO(IT,JT,HZA,IBA,KR,KRI) Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 130 LVThS C C n n vùng b bi n Vi t Nam ! SUISHIN NO KEISAN ! CALL DEPTH(IT,JT,HZA,HMA,HNA) ! SYOKI SUII WO YOMU ! CALL SETZRO(IT,JT,ZA,MZA,MA,NA,VA,DZA,HZA,INIZ,IZR) C DO 10 K=1,KE K1=K WRITE(*,*) K C C C ! RENZOKU-SHIKI ! CALL CONTIN(IT,JT,2,ITM,2,JTM,ZA,MA,NA,DZA,HZA,RA,IBA,K1,KTA) C C C ! UNDOU-SHIKI ! CALL MOTION(IT,JT,2,ITM,2,JTM,ZA,MA,NA,VA,DZA,DMA,DNA,HZA,HMA,HNA, #IBA,RA,0.025,DTA,K1,KTA) C C C C ! KYOUKAI-JYOUKEN ! CALL BOUND(IT,JT,HZA,ZA,MA,NA,RA,K1,KTA) ! OUTPUT ! CALL VISIT(ZA,IT,JT,fnumber,fname,120,KR,KV,K1) CALL OUTPUTA(LPI,LPJ,NPL,IT,JT,ZA,MA,NA,VA,K1) ! TIME STEP WO SUSUMERU ! CALL CHANGE(IT,JT,ZA,MZA,MA,NA,DZA,K1,KTA) C C C 10 CONTINUE ! SAIDAITI WO KAKU ! CALL OUTZMAX(MZA,IT,JT) CALL OUTVMAX(VA,IT,JT) C C C C C C STOP END ** ** ** SUISHIN NO NYURYOKU ** ** ** **************************************************** SUBROUTINE DEPTHO(IG,JG,HZ,IB,KR,KRI) C DIMENSION HZ(IG,JG),IB(IG,JG) C DO 10 I=1,IG DO 10 J=1,JG HZ(I,J)=0.0 10 CONTINUE C IF(KR.EQ.11) GO TO 100 C READ(KR,200) ((HZ(I,J),J=1,JG),I=1,IG) 200 FORMAT(10F7.1) READ(KRI,250) ((IB(I,J),J=1,JG),I=1,IG) 250 FORMAT(800I1) GO TO 20 C 100 READ(KR,300) ((HZ(I,J),J=1,JG),I=1,IG) 300 FORMAT(10F7.1) READ(KRI,350) ((IB(I,J),J=1,JG),I=1,IG) 350 FORMAT(242I1) Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 131 LVThS - n n vùng b bi n Vi t Nam 20 CONTINUE C DO 500 I=1,IG DO 500 J=1,JG IB(I,J)=10*IB(I,J)+1 500 CONTINUE C RETURN END C C C C C ** ** ** M,N,(KEISANTEN) NO SUISHIN WO MOTOMERU ** ** ** **************************************************** SUBROUTINE DEPTH(IG,JG,HZ,HM,HN) C DIMENSION HZ(IG,JG),HM(IG,JG),HN(IG,JG) C DO 10 I=1,IG DO 10 J=1,JG C IF(I.EQ.IG) GO TO 20 HM(I,J)=0.5*(HZ(I,J)+HZ(I+1,J)) GO TO 30 20 HM(I,J)=HZ(I,J) 30 IF(J.EQ.JG) GO TO 40 HN(I,J)=0.5*(HZ(I,J)+HZ(I,J+1)) GO TO 10 40 HN(I,J)=HZ(I,J) C 10 CONTINUE C C C C C C END ** ** ** SYOKI SUII WO YOMU ** ** ** **************************************************** SUBROUTINE SETZRO(IG,JG,Z,M,MZ,N,V,DZ,HZ,INIZ,IZR) C REAL M,N,MZ C DIMENSION Z(IG,JG,2),M(IG,JG,2),N(IG,JG,2),MZ(IG,JG) DIMENSION DZ(IG,JG,2),HZ(IG,JG),V(IG,JG,2) C DO 10 K=1,2 DO 10 I=1,IG DO 10 J=1,JG C Z(I,J,K)=0.0 M(I,J,K)=0.0 N(I,J,K)=0.0 V(I,J,K)=0.0 DZ(I,J,K)=HZ(I,J) C IF(HZ(I,J).LT.1.0E-5) DZ(I,J,K)=0.0 C Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 132 LVThS - n n vùng b bi n Vi t Nam 10 CONTINUE C IF(INIZ.EQ.3) GO TO 16 C READ(IZR,100) ((Z(I,J,1),J=1,JG),I=1,IG) 100 FORMAT(800F8.3) GO TO 18 C 16 READ(IZR,200) ((Z(I,J,1),J=1,JG),I=1,IG) 200 FORMAT(239F8.3) 18 CONTINUE C DO 20 I=1,IG DO 20 J=1,JG C MZ(I,J)=0.0 DZ(I,J,1)=HZ(I,J)+Z(I,J,1) IF(HZ(I,J).LT.1.0E-5) DZ(I,J,1)=0.0 C 20 CONTINUE C C C C C C RETURN END ** ** ** RENZOKU-SHIKI ** ** ** **************************************************** SUBROUTINE CONTIN(IG,JG,IS,IE,JS,JE,Z,M,N,DZ,HZ,R,IB,KK,KT) C REAL M,N C DIMENSION Z(IG,JG,2),M(IG,JG,2),N(IG,JG,2),HZ(IG,JG) DIMENSION DZ(IG,JG,2),IB(IG,JG) C DATA GX/1.0E-5/ C IF(KT.EQ.1) GO TO 200 IF(MOD(KK,KT).NE.1) RETURN 200 CONTINUE C DO 10 I=IS,IE DO 10 J=JS,JE C IF(IB(I,J).GE.40) GO TO 11 C XM=0.0 XN=0.0 C IF(I.NE.1) XM=M(I-1,J,2) IF(J.NE.1) XN=N(I,J-1,2) C ZZ=Z(I,J,1)-R*(M(I,J,2)-XM+N(I,J,2)-XN) IF(ABS(ZZ).LT.1.0E-10) ZZ=0.0 C DD=ZZ+HZ(I,J) IF(DD.LT.GX) GO TO 12 C DZ(I,J,2)=DD Z(I,J,2)=ZZ Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 133 LVThS - n n vùng b bi n Vi t Nam GO TO 10 C 11 DD=0.0 HZ(I,J)=0.0 DZ(I,J,2)=DD Z(I,J,2)=DD GO TO 10 C 12 DD=0.0 DZ(I,J,2)=DD Z(I,J,2)=DD-HZ(I,J) C 10 CONTINUE C C C C C C RETURN END ** ** ** UNDO HOUTEISIKI ** ** ** **************************************************** SUBROUTINE MOTION(IG,JG,IS,IE,JS,JE,Z,M,N,V,DZ,DM,DN,HZ,HM,HN, #IB,R,FM,DT,KK,KT) C REAL M,N C DIMENSION DIMENSION DIMENSION DIMENSION Z(IG,JG,2),M(IG,JG,2),N(IG,JG,2) DZ(IG,JG,2),DM(IG,JG,2),DN(IG,JG,2) HZ(IG,JG),HM(IG,JG),HN(IG,JG) IB(IG,JG),V(IG,JG,2) C DATA GG,GX/9.8,1.0E-7/ C IF(KT.EQ.1) GO TO 200 IF(MOD(KK,KT).NE.KT/2+1) RETURN 200 CONTINUE C DO 10 I=IS,IE DO 10 J=JS,JE C DM1=0.25*(Z(I,J,1)+Z(I+1,J,1)+Z(I,J,2)+Z(I+1,J,2))+HM(I,J) DM2=0.5*(Z(I,J,2)+Z(I+1,J,2))+HM(I,J) DN1=0.25*(Z(I,J,1)+Z(I,J,2)+Z(I,J+1,1)+Z(I,J+1,2))+HN(I,J) C DN2=0.5*(Z(I,J,2)+Z(I,J+1,2))+HN(I,J) C IF(DM1.LT.GX) IF(DM2.LT.GX) IF(DN1.LT.GX) IF(DN2.LT.GX) DM1=0.0 DM2=0.0 DN1=0.0 DN2=0.0 C DM(I,J,1)=DM1 DM(I,J,2)=DM2 DN(I,J,1)=DN1 DN(I,J,2)=DN2 C 10 CONTINUE C FN=0.5*DT*GG*FM**2 C Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 134 LVThS - n n vùng b bi n Vi t Nam DO 20 I=IS,IE DO 20 J=JS,JE C IBB=IB(I,J)/10 IBR=MOD(IB(I,J),10) C IF(IBB.GT.3)GO TO 20 IF(IBB.EQ.1.OR.IBB.EQ.3) GO TO 30 IF(I.EQ.IG)GO TO 30 C IF(IBR.EQ.0)GO TO 33 IF(DZ(I,J,2))31,31,32 C 31 32 34 35 IF(DZ(I+1,J,2))30,30,34 IF(DZ(I+1,J,2))35,35,36 IF(Z(I+1,J,2)+HZ(I,J))30,30,37 IF(Z(I,J,2)+HZ(I+1,J))30,30,38 36 DD=DM(I,J,2) GO TO 39 DD=0.5*(Z(I,J,2)+Z(I+1,J,2))+HZ(I,J) GO TO 39 DD=0.5*(Z(I,J,2)+Z(I+1,J,2))+HZ(I+1,J) GO TO 39 DD=HM(I,J) IF(DD.LT.0.1) GO TO 30 XNN=0.25*(N(I,J,1)+N(I+1,J,1)+N(I,J-1,1)+N(I+1,J-1,1)) C 37 38 33 39 C DF=DD IF(DF.LT.1.0E-2) DF=1.0E-2 FF=FN*SQRT(M(I,J,1)**2+XNN**2)/DF**(7./3.) IF(DD.LT.GX) GO TO 30 XM=(1.0-FF)*M(I,J,1)-GG*R*DD*(Z(I+1,J,2)-Z(I,J,2)) C IF(IBR.EQ.0)GO TO 40 IF(DM(I,J,1).LT.GX) GO TO 40 C 41 IF(M(I,J,1)) 41,42,42 IF(DM(I+1,J,1).LT.GX) GO TO 40 C XM=XM-R*(M(I+1,J,1)**2/DM(I+1,J,1)-M(I,J,1)**2/DM(I,J,1)) GO TO 43 42 IF(DM(I-1,J,1).LT.GX) GO TO 40 C XM=XM-R*(M(I,J,1)**2/DM(I,J,1)-M(I-1,J,1)**2/DM(I-1,J,1)) C 43 IF(XNN)44,45,45 44 XNE=0.25*(N(I,J+1,1)+N(I+1,J+1,1)+N(I,J,1)+N(I+1,J,1)) C IF(DM(I,J+1,1).LT.GX) GO TO 40 XM=XM-R*(M(I,J+1,1)*XNE/DM(I,J+1,1)-M(I,J,1)*XNN/DM(I,J,1)) GO TO 40 45 IF(J.EQ.2)GO TO 46 XNE=0.25*(N(I,J-1,1)+N(I+1,J-1,1)+N(I,J-2,1)+N(I+1,J-2,1)) GO TO 47 46 XNE=0.5*(N(I,J-1,1)+N(I+1,J-1,1)) C 47 IF(DM(I,J-1,1).LT.GX) GO TO 40 XM=XM-R*(M(I,J,1)*XNN/DM(I,J,1)-M(I,J-1,1)*XNE/DM(I,J-1,1)) 40 XM=XM/(1.0+FF) Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 135 LVThS - n n vùng b bi n Vi t Nam C C C IF(ABS(XM).LT.1.0E-10)XM=0.0 M(I,J,2)=XM GO TO 100 30 M(I,J,2)=0.0 100 CONTINUE 100 IF(IBB.GE.2)GO TO 130 IF(J.EQ.JG) GO TO 130 C 131 132 134 135 136 137 138 133 139 IF(IBR.EQ.0) GO TO 133 IF(DZ(I,J,2)) 131,131,132 IF(DZ(I,J+1,2)) 130,130,134 IF(DZ(I,J+1,2)) 135,135,136 IF(Z(I,J+1,2)+HZ(I,J)) 130,130,137 IF(Z(I,J,2)+HZ(I,J+1)) 130,130,138 DD=DN(I,J,2) GO TO 139 DD=0.5*(Z(I,J,2)+Z(I,J+1,2))+HZ(I,J) GO TO 139 DD=0.5*(Z(I,J,2)+Z(I,J+1,2))+HZ(I,J+1) GO TO 139 DD=HN(I,J) IF(DD.LT.0.1) GO TO 130 XMM=0.25*(M(I,J,1)+M(I,J+1,1)+M(I-1,J,1)+M(I-1,J+1,1)) C DF=DD IF(DF.LT.1.0E-2) DF=1.0E-2 FF=FN*SQRT(N(I,J,1)**2+XMM**2)/DF**(7./3.) IF(DD.LT.GX) GO TO 130 XN=(1.-FF)*N(I,J,1)-GG*R*DD*(Z(I,J+1,2)-Z(I,J,2)) IF(IBR.EQ.0) GO TO 140 IF(DN(I,J,1).LT.GX) GO TO 140 C IF(N(I,J,1)) 141,142,142 141 IF(DN(I,J+1,1).LT.GX) GO TO 140 XN=XN-R*(N(I,J+1,1)**2/DN(I,J+1,1)-N(I,J,1)**2/DN(I,J,1)) GO TO 143 142 IF(DN(I,J-1,1).LT.GX) GO TO 140 XN=XN-R*(N(I,J,1)**2/DN(I,J,1)-N(I,J-1,1)**2/DN(I,J-1,1)) C 143 144 IF(XMM)144,145,145 XME=0.25*(M(I+1,J,1)+M(I+1,J+1,1)+M(I,J,1)+M(I,J+1 ,1)) IF(DN(I+1,J,1).LE.GX)GO TO 140 XN=XN-R*(N(I+1,J,1)*XME/DN(I+1,J,1)-N(I,J,1)*XMM/DN(I,J,1)) GO TO 140 145 IF(I.EQ.2)GO TO 146 XME=0.25*(M(I-1,J,1)+M(I-1,J+1,1)+M(I-2,J,1)+M(I-2,J+1,1)) GO TO 147 XME=0.5*(M(I-1,J,1)+M(I-1,J+1,1)) IF(DN(I-1,J,1).LT.GX)GO TO 140 XN=XN-R*(N(I,J,1)*XMM/DN(I,J,1)-N(I-1,J,1)*XME/DN(I-1,J,1)) XN=XN/(1.0+FF) IF(ABS(XN).LT.1.0E-10) XN=0.0 N(I,J,2)=XN GO TO 20 N(I,J,2)=0.0 CONTINUE C 146 147 140 130 20 C Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 136 LVThS - n n vùng b bi n Vi t Nam DO 60 I=1,IG DO 60 J=1,JG C IF(DZ(I,J,2).LE.0.0)GO TO 61 V(I,J,1)=SQRT(M(I,J,2)**2+N(I,J,2)**2)/DZ(I,J,2) GO TO 62 61 V(I,J,1)=-1.0 62 IF(V(I,J,1).GT.V(I,J,2)) V(I,J,2)=V(I,J,1) C 60 C C C C C CONTINUE RETURN END ** ** ** KYOKAI JOKEN ** ** ** **************************************************** SUBROUTINE BOUND(IG,JG,HZ,Z,M,N,R,KK,KT) C REAL M,N C DIMENSION M(IG,JG,2),N(IG,JG,2),Z(IG,JG,2),HZ(IG,JG) C DATA GG/9.80/ C IF(MOD(KK,KT).NE.KT/2+1) RETURN C IT=IG-1 JT=JG-1 C DO 10 J=1,JG C IF(HZ(1,J).LT.0.0) HZ(1,J)=0.0 IF(HZ(IG,J).LT.0.0) HZ(IG,J)=0.0 C IF(M(2,J,2)**2+N(2,J,2)**2.EQ.0.0) GO TO 15 C M(1,J,2)=-SQRT(GG*HZ(1,J))*Z(2,J,2)*ABS(M(2,J,2))/ #SQRT(M(2,J,2)**2+N(2,J,2)**2) N(1,J,2)=SQRT(GG*HZ(1,J))*Z(2,J,2)*ABS(N(2,J,2))/ #SQRT(M(2,J,2)**2+N(2,J,2)**2) C C C CC=SQRT(GG*HZ(1,J)) M(1,J,2)=R*CC*M(2,J,1)+(1-R*CC)*(M(1,J,1)) N(1,J,2)=R*CC*N(2,J,1)+(1-R*CC)*(N(1,J,1)) Z(1,J,2)=R*CC*Z(2,J,1)+(1-R*CC)*(Z(1,J,1)) C GO TO 17 C 15 M(1,J,2)=0.0 Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 137 LVThS - n n vùng b bi n Vi t Nam N(1,J,2)=0.0 C 17 IF(M(IT,J,2)**2+N(IT,J,2)**2.EQ.0.0) GO TO 16 C M(IG,J,2)=SQRT(GG*HZ(IG,J))*Z(IT,J,2)*ABS(M(IT,J,2))/ #SQRT(M(IT,J,2)**2+N(IT,J,2)**2) N(IG,J,2)=SQRT(GG*HZ(IG,J))*Z(IT,J,2)*ABS(N(IT,J,2))/ #SQRT(M(IT,J,2)**2+N(IT,J,2)**2) C C C CC=SQRT(GG*HZ(IG,J)) M(IG,J,2)=R*CC*M(IT,J,1)+(1-R*CC)*(M(IG,J,1)) N(IG,J,2)=R*CC*N(IT,J,1)+(1-R*CC)*(N(IG,J,1)) Z(IG,J,2)=R*CC*Z(IT,J,1)+(1-R*CC)*(Z(IG,J,1)) C GO TO 10 C 16 M(IG,J,2)=0.0 N(IG,J,2)=0.0 C 10 CONTINUE C DO 20 I=1,IG C IF(HZ(I,1).LT.0.0) HZ(I,1)=0.0 IF(HZ(I,JG).LT.0.0) HZ(I,JG)=0.0 C IF(M(I,2,2)**2+N(I,2,2)**2.EQ.0.0) GO TO 25 C M(I,1,2)=SQRT(GG*HZ(I,1))*Z(I,2,2)*ABS(M(I,2,2))/ #SQRT(M(I,2,2)**2+N(I,2,2)**2) N(I,1,2)=-SQRT(GG*HZ(I,1))*Z(I,2,2)*ABS(N(I,2,2))/ #SQRT(M(I,2,2)**2+N(I,2,2)**2) C C C CC=SQRT(GG*HZ(I,1)) M(I,1,2)=R*CC*M(I,2,1)+(1-R*CC)*(M(I,1,1)) N(I,1,2)=R*CC*N(I,2,1)+(1-R*CC)*(N(I,1,1)) Z(I,1,2)=R*CC*Z(I,2,1)+(1-R*CC)*(Z(I,1,1)) C GO TO 27 C 25 M(I,1,2)=0.0 N(I,1,2)=0.0 C 27 IF(M(I,JT,2)**2+N(I,JT,2)**2.EQ.0.0) GO TO 26 C M(I,JG,2)=SQRT(GG*HZ(I,JG))*Z(I,JT,2)*ABS(M(I,JT,2))/ #SQRT(M(I,JT,2)**2+N(I,JT,2)**2) N(I,JG,2)=SQRT(GG*HZ(I,JG))*Z(I,JT,2)*ABS(N(I,JT,2))/ #SQRT(M(I,JT,2)**2+N(I,JT,2)**2) C C C CC=SQRT(GG*HZ(I,JG)) M(I,JG,2)=R*CC*M(I,JT,1)+(1-R*CC)*(M(I,JG,1)) N(I,JG,2)=R*CC*N(I,JT,1)+(1-R*CC)*(N(I,JG,1)) Z(I,JG,2)=R*CC*Z(I,JT,1)+(1-R*CC)*(Z(I,JG,1)) C GO TO 20 C 26 M(I,JG,2)=0.0 N(I,JG,2)=0.0 C Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 138 LVThS - n n vùng b bi n Vi t Nam 20 CONTINUE C C C C C C RETURN END ** ** ** TIME STEP WO SUSUMERU ** ** ** **************************************************** SUBROUTINE CHANGE(IG,JG,Z,MZ,M,N,DZ,KK,KT) C REAL M,N,MZ C DIMENSION Z(IG,JG,2),M(IG,JG,2),N(IG,JG,2),DZ(IG,JG,2),MZ(IG,JG) C LL=MOD(KK,KT) LX=KT/2 IF(LL.NE.0.AND.LL-LX.NE.0) RETURN C DO 10 I=1,IG DO 10 J=1,JG IF(LL.NE.0) GO TO 20 Z(I,J,1)=Z(I,J,2) DZ(I,J,1)=DZ(I,J,2) IF(Z(I,J,2).GT.MZ(I,J)) MZ(I,J)=Z(I,J,2) 20 IF(LL-LX.NE.0) GO TO 10 M(I,J,1)=M(I,J,2) N(I,J,1)=N(I,J,2) 10 CONTINUE C C C C C C C C RETURN END ** ** ** OUTPUT SUBROUTINE ** ** ** **************************************************** ! For VisIt ! SUBROUTINE VISIT(Z,IG,JG,fnumber,fname,OUT,KR,KV,KK) C CHARACTER*100 fname INTEGER fnumber,OUT C DIMENSION Z(IG,JG,2) REAL ZZZ(960,800) C C C Check whether sub folder outputL exists or not LOGICAL EX INQUIRE (DIRECTORY = "outputL", EXIST = EX) IF (.NOT.EX) THEN If there is not a folder \outputL\ then create it CALL system('md ' // "outputL") END IF C IF(MOD(KK,OUT).EQ.1) THEN C fnumber=fnumber+1 C WRITE(fname,"('outputL\grid',i3.3,'.plt')"), fnumber Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 139 LVThS - n n vùng b bi n Vi t Nam C OPEN(KV,file=fname) WRITE(KV,*)'TITLE="Graph - "' WRITE(KV,*)'VARIABLES=x,y,e' WRITE(KV,*)'ZONE T="ZONE1 ",I=', JG,',J=', IG,',F=POINT' C DO 100 I=1,IG DO 100 J=1,JG II=IG+1-I JJ=JG+1-J ZZZ(I,J)=Z(II,JJ,2) 100 CONTINUE C DO 200 I=1,IG DO 200 J=1,JG WRITE(KV,250)I,J,ZZZ(I,J) 250 FORMAT(1H ,2I6,F10.3) 200 CONTINUE C CLOSE(KV) C END IF C RETURN END C -C ! OUTPUTA ! C SUBROUTINE OUTPUTA(LPI,LPJ,NPL,IT,JT,Z,M,N,V,KK) C REAL Z,M,N C DIMENSION Z(IT,JT,2),M(IT,JT,2),N(IT,JT,2),V(IT,JT,2) DIMENSION LPI(NPL),LPJ(NPL) C IF(MOD(KK,6).NE.0) GO TO 405 C WRITE(41,401) KK, #Z(LPI(1),LPJ(1),2),Z(LPI(2),LPJ(2),2), #Z(LPI(3),LPJ(3),2),Z(LPI(4),LPJ(4),2), #Z(LPI(5),LPJ(5),2),Z(LPI(6),LPJ(6),2), #Z(LPI(7),LPJ(7),2),Z(LPI(8),LPJ(8),2), #Z(LPI(9),LPJ(9),2),Z(LPI(10),LPJ(10),2), #Z(LPI(11),LPJ(11),2),Z(LPI(12),LPJ(12),2), #Z(LPI(13),LPJ(13),2),Z(LPI(14),LPJ(14),2) 401 FORMAT(I7,14F12.4) C WRITE(42,402) KK, #M(LPI(1),LPJ(1),2),M(LPI(2),LPJ(2),2), #M(LPI(3),LPJ(3),2),M(LPI(4),LPJ(4),2), #M(LPI(5),LPJ(5),2),M(LPI(6),LPJ(6),2), #M(LPI(7),LPJ(7),2),M(LPI(8),LPJ(8),2), #M(LPI(9),LPJ(9),2),M(LPI(10),LPJ(10),2), #M(LPI(11),LPJ(11),2),M(LPI(12),LPJ(12),2), #M(LPI(13),LPJ(13),2),M(LPI(14),LPJ(14),2) 402 FORMAT(I7,14F12.4) C Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 140 LVThS - n n vùng b bi n Vi t Nam WRITE(43,403) KK, #N(LPI(1),LPJ(1),2),N(LPI(2),LPJ(2),2), #N(LPI(3),LPJ(3),2),N(LPI(4),LPJ(4),2), #N(LPI(5),LPJ(5),2),N(LPI(6),LPJ(6),2), #N(LPI(7),LPJ(7),2),N(LPI(8),LPJ(8),2), #N(LPI(9),LPJ(9),2),N(LPI(10),LPJ(10),2), #N(LPI(11),LPJ(11),2),N(LPI(12),LPJ(12),2), #N(LPI(13),LPJ(13),2),N(LPI(14),LPJ(14),2) 403 FORMAT(I7,14F12.4) C WRITE(44,404) KK, #V(LPI(1),LPJ(1),1),V(LPI(2),LPJ(2),1), #V(LPI(3),LPJ(3),1),V(LPI(4),LPJ(4),1), #V(LPI(5),LPJ(5),1),V(LPI(6),LPJ(6),1), #V(LPI(7),LPJ(7),1),V(LPI(8),LPJ(8),1), #V(LPI(9),LPJ(9),2),V(LPI(10),LPJ(10),2), #V(LPI(11),LPJ(11),2),V(LPI(12),LPJ(12),2), #V(LPI(13),LPJ(13),2),V(LPI(14),LPJ(14),2) 404 FORMAT(I7,14F12.4) 405 CONTINUE C RETURN END C C C ! OUTZMAX ! C SUBROUTINE OUTZMAX(MZA,IT,JT) REAL MZA DIMENSION MZA(IT,JT) C WRITE(51,501) ((MZA(I,J),J=1,JT),I=1,IT) 501 FORMAT(800F8.3) C RETURN END C ! OUTVMAX ! C SUBROUTINE OUTVMAX(VA,IT,JT) DIMENSION VA(IT,JT,2) C WRITE(61,601) ((VA(I,J,2),J=1,JT),I=1,IT) 601 FORMAT(800F8.3) C RETURN END Ph l c HVTH: Nguy n Th Minh Lý - MSHV: 12020445 141 ... Trong vùng ngu t gây sóng th n n vùng b bi n Vi t Nam t Bi n có kh t Nam là: RK Hình 2.5 Các ngu t có kh sóng th n Bi t Nam (ngu n Tr n Th M Thành c ng s , 2011) Chú thích: - Vùng ngu n B c Bi - Vùng. .. 12020445 17 LVThS - n n vùng b bi n Vi t Nam Tóm t t k t qu sóng th n 13 k ch b n K ch b n 1,2 Khu v c Các vùng b Th i gian lan sóng truy n sóng th n bi n Vi t Nam Chi u cao Ghí cao sóng T k ch b n 1-5... lý thuy t sóng n tính v - c nông cho i chia c TUNAMI-N2 : lý thuy n sóng n tính lên b v i chia c nh c nơng sóng leo nh - TUNAMI-N3 : lý thuy t sóng n tính v i - TUNAMI-F1 : lý thuy t sóng n tính