add meteo tsunami

simple_cases/meteo_tsunami/input.txt

@@ -0,0 +1,90 @@
+!INPUT FILE FOR FUNWAVE_TVD
+  ! NOTE: all input parameter are capital sensitive
+  ! --------------------TITLE-------------------------------------
+  ! title only for log file
+TITLE = MeteoTsunami
+  ! -------------------PARALLEL INFO-----------------------------
+  ! 
+  !    PX,PY - processor numbers in X and Y
+  !    NOTE: make sure consistency with mpirun -np n (px*py)
+  !    
+PX = 4
+PY = 1
+  ! --------------------DEPTH-------------------------------------
+  ! Depth types, DEPTH_TYPE=DATA: from depth file
+  !              DEPTH_TYPE=FLAT: idealized flat, need depth_flat
+  !              DEPTH_TYPE=SLOPE: idealized slope, 
+  !                                 need slope,SLP starting point, Xslp
+  !                                 and depth_flat
+DEPTH_TYPE = FLAT
+DEPTH_FLAT = 40.0
+  ! -------------------PRINT---------------------------------
+  ! PRINT*,
+  ! result folder
+RESULT_FOLDER = output/
+
+  ! ------------------DIMENSION-----------------------------
+  ! global grid dimension
+Mglob = 600
+Nglob = 100
+
+  ! ----------------- TIME----------------------------------
+  ! time: total computational time/ plot time / screen interval 
+  ! all in seconds
+TOTAL_TIME = 30000.0
+PLOT_INTV = 360.0
+PLOT_INTV_STATION = 50000.0
+SCREEN_INTV = 360.0
+  ! -----------------GRID----------------------------------
+  ! if use spherical grid, in decimal degrees
+  ! cartesian grid sizes
+DX = 500.0
+DY = 500.0
+  ! ---------------- SPONGE LAYER ------------------------
+  ! DHI type sponge layer
+  ! need to specify widths of four boundaries and parameters
+  ! set width=0.0 if no sponge
+  ! R_sponge: decay rate
+  ! A_sponge: maximum decay rate
+  ! e.g., sharp: R=0.85
+  !       mild:  R=0.90, A=5.0
+  !       very mild, R=0.95, A=5.0
+DIFFUSION_SPONGE = F
+FRICTION_SPONGE = F
+DIRECT_SPONGE = F
+Csp = 0.10
+CDsponge = 1.0
+Sponge_west_width =  7.5
+Sponge_east_width =  5.0
+Sponge_south_width = 0.0
+Sponge_north_width = 0.0
+R_sponge = 0.85
+A_sponge = 5.0
+  !----------------Friction-----------------------------
+Cd = 0.002
+  ! ----------------NUMERICS----------------------------
+  ! CFL
+CFL = 0.5
+  ! Froude Number Cap (to avoid jumping drop, set 10)
+FroudeCap = 1.0
+
+  ! --------------WET-DRY-------------------------------
+  ! MinDepth for wetting-drying
+MinDepth=0.01
+  ! -----------------OUTPUT-----------------------------
+  ! stations 
+  ! if NumberStations>0, need input i,j in STATION_FILE
+NumberStations = 0
+!STATIONS_FILE = gauges.txt
+  ! output variables, T=.TRUE, F = .FALSE.
+DEPTH_OUT = T
+U = F
+V = F
+ETA = T
+Hmax = T
+Hmin = T
+!   ----- meteo tsunami application
+MeteoGausian = T
+METEO_GAUSIAN_FILE = meteo_data.txt
+OUT_METEO = T
+

simple_cases/meteo_tsunami/meteo_data.txt

@@ -0,0 +1,5 @@
+Meteo data file
+NJ
+time, x, y, dP(mb), SigmaX, SigmaY, Angle(0 = +x direction)
+0.0       0.0 25000.0   5.0 10000.0 100000.0 0.0
+36000.0  720000.0 25000.0  5.0 10000.0 10000.0 0.0 

simple_cases/meteo_tsunami/plot_wave.m

@@ -0,0 +1,44 @@
+clear all
+fdir='/Users/fengyanshi15/tmp1/';
+
+eta=load([fdir 'eta_00001']);
+
+[n,m]=size(eta);
+dx=500.0;
+dy=500.0;
+x=[0:m-1]*dx;
+y=[0:n-1]*dy;
+
+
+nfile=[11 16 21 26 31 36];
+min={'1.0' '1.5' '2.0' '2.5' '3.0' '3.5'};
+
+wid=6;
+len=8;
+set(gcf,'units','inches','paperunits','inches','papersize', [wid len],'position',[1 1 wid len],'paperposition',[0 0 wid len]);
+clf
+
+
+for num=1:length(nfile)
+
+fnum=sprintf('%.5d',nfile(num));
+eta=load([fdir 'eta_' fnum]);
+
+subplot(6,1,num)
+plot(x(:),eta(50,:),'LineWidth',1.5)
+axis([0 300000 -0.5 1])
+grid
+
+title([' Time = ' min{num} ' hr '])
+
+
+ylabel(' eta (m) ')
+
+if(num==length(nfile))
+xlabel(' x (m) ')
+end
+
+
+set(gcf,'Renderer','zbuffer')
+
+end

src/io.F

@@ -3492,6 +3492,11 @@ SUBROUTINE PREVIEW
         call PutFile(TMP_NAME,WindV2D)
     ENDIF
 
+    IF(MeteoGausian)THEN
+        TMP_NAME = TRIM(FDIR)//'Pstorm_'//TRIM(FILE_NAME)
+        call PutFile(TMP_NAME,StormPressureTotal)
+    ENDIF
+
     IF(SlideModel)THEN
         TMP_NAME = TRIM(FDIR)//'Pstorm_'//TRIM(FILE_NAME)
         call PutFile(TMP_NAME,StormPressureTotal)