how obtaining curves in SWE-SPHysics ???

Issues and examples using this model (available in the Models Repository)

how obtaining curves in SWE-SPHysics ???

Postby rachedelwafi on Sat Nov 23, 2013 3:38 am

Hi everybody,
I have a question about "how obtaining curves (water depth = fct(distane "X") at a section "Y") ?"

I work with the SWE-SPHysics, i have introduced the geometry of my case in SWE-SPHysics and i runed the case with FTN95 compilator. Then i have done the visualisation with PARAVIEW.

My problem is to now how obtaining curves for example (Water depth at section "Y") or for example (velocity at a section "Y" and a time "t").

Thank you very much.
rachedelwafi
 
Posts: 1
Joined: Fri Nov 22, 2013 9:23 am

Re: how obtaining curves in SWE-SPHysics ???

Postby rvacondio on Mon Feb 17, 2014 10:14 am

Besides data on particles, SWE-SPHysics plots water depth and velocity maps on a regular grid (files GRD_d_[timestep_number], GRD_u_[timestep_number], GRD_v_[timestep_number]) those files are
loaded by the SPHYSICS_SWE_Plot_2D.m matlab script and they can be used to evaluate the water depth and/or the velocities at any point.

For example, if you want to know the water depth at point with coordinate (x_point, y_point) you have to include the following command in the SPHYSICS_SWE_Plot_2D.m matlab script:

h_point = interp2(XX,YY,h_int,x_point,y_point)

the same for the velocities

I hope this might help you.

Best Regards
Renato
rvacondio
 
Posts: 1
Joined: Mon Feb 17, 2014 9:41 am

Re: how obtaining curves in SWE-SPHysics ???

Postby panda.wang.1989 on Tue Jun 03, 2014 7:48 pm

Hello, I am Wei Wang, form Kyushu University, Japan.
I recognised you are a SWE_SPHysics user form your post, and i am now using SWE_SPHysics.
While i cant't compile the code successfully using FTN95 as yours.
Could you leave me an e-mail so that i can ask your help.
My e-mail is panda.wang.1989@gmail.com
Thank you very much.
Best regards.
panda.wang.1989
 
Posts: 2
Joined: Tue Jun 03, 2014 10:36 am

Re: how obtaining curves in SWE-SPHysics ???

Postby henri8388 on Sun Nov 16, 2014 4:09 am

Hi, I have the same question as same as Wei above mentioned, I also can not get through SWE-SPHysics by the compiler of FTN95. And it showed as follows(1D-Case1):

Salford MK/Win32 Version 3.89, Copyright (c) Silverfrost Ltd. 1996-2006]
ftn95 global_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
PROCESSING MODULE [<GLOBAL_1D> FTN95/Win32 v7.10.0]
NO ERRORS [<GLOBAL_1D> FTN95/Win32 v7.10.0]
ftn95 ini_divide_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<INI_DIVIDE> FTN95/Win32 v7.10.0]
ftn95 SPHYSICS_SWE_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0049) write(*,*) 'the FreeBSD License;
*** Unterminated character constant
1 ERROR [<SPH> FTN95/Win32 v7.10.0]
*** Compilation failed
ftn95 getdata_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<GETDATA> FTN95/Win32 v7.10.0]
ftn95 check_limits_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<CHECK_LIMITS_1D> FTN95/Win32 v7.10.0]
NO ERRORS [<NEW_BCP> FTN95/Win32 v7.10.0]
0171) integer iflagg, ifl, ipos, i
WARNING - The argument IPOS has not been used
NO ERRORS, 1 WARNING [<NEW_OB_FLP> FTN95/Win32 v7.10.0]
NO ERRORS [<NEW_FL_OB> FTN95/Win32 v7.10.0]
ftn95 poute_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<POUTE> FTN95/Win32 v7.10.0]
ftn95 divide_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<DIVIDE> FTN95/Win32 v7.10.0]
ftn95 variable_time_step_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0030) integer i, i_fa
WARNING - Variable I_FA has been given a value but never used
NO ERRORS, 1 WARNING [<VARIABLE_TIME_STEP> FTN95/Win32 v7.10.0]
ftn95 ac_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<AC> FTN95/Win32 v7.10.0]
ftn95 ac_dw_var_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0035) real Res, hmax_def, adh, dxx, ncou, pippo
WARNING - Variable PIPPO has been declared but not used
NO ERRORS, 1 WARNING [<AC_DW> FTN95/Win32 v7.10.0]
ftn95 celij_dw_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0031) integer i, j1, j2, kind_p2, j, jj
WARNING - The argument J1 has not been used
0032) real drx, rr2, dux
WARNING - Variable DUX has been given a value but never used
NO ERRORS, 2 WARNINGS [<CELIJ_DW> FTN95/Win32 v7.10.0]
ftn95 self_dw_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<SELF_DW> FTN95/Win32 v7.10.0]
ftn95 kernel_cubic_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0032) real h, const_ker(2), one_over_h, pippo
WARNING - Variable PIPPO has been declared but not used
0039) if (const_ker(1).ne.1./6..or.const_ker(2).ne.2./3.) then
WARNING - Comparing floating point quantities for inequality may give
misleading results
NO ERRORS, 2 WARNINGS [<KERNEL> FTN95/Win32 v7.10.0]
ftn95 step_leap_frog_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<STEP> FTN95/Win32 v7.10.0]
ftn95 ghost_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<GHOST> FTN95/Win32 v7.10.0]
ftn95 ac_alpha_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<AC_ALPHA> FTN95/Win32 v7.10.0]
ftn95 celij_alpha_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<CELIJ_ALPHA> FTN95/Win32 v7.10.0]
ftn95 self_alpha_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<SELF_ALPHA> FTN95/Win32 v7.10.0]
ftn95 grid_h_var_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<GRID_H_VAR> FTN95/Win32 v7.10.0]
ftn95 celij_visc_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<CELIJ> FTN95/Win32 v7.10.0]
ftn95 self_visc_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0032) integer j1, j2, ii, jj, i, j, jj_start
WARNING - Variable J2 has been given a value but never used
NO ERRORS, 1 WARNING [<SELF> FTN95/Win32 v7.10.0]
ftn95 celij_ob_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0036) integer i, j1, j2, j, jj
WARNING - The argument J1 has not been used
NO ERRORS, 1 WARNING [<CELIJ_OB> FTN95/Win32 v7.10.0]
ftn95 viscosity_LF_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0027) For more information see http://opensource.org/licenses
*** Invalid characters in label field of FIXED format source
*** Statement not recognised
0033) subroutine viscosity(pi_visc,dot,drx,dux,kind,rr2,
0034) + cbar,rhobar,one_over_rhobar,hbar)
*** SUBROUTINE cannot be declared inside PROGRAM block (perhaps missing
CONTAINS or END statement?)

3 ERRORS [viscosity_LF_1D.F] - Compilation failed.
ftn95 divide_b_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<DIVIDE_B> FTN95/Win32 v7.10.0]
ftn95 bottom_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<BOTTOM> FTN95/Win32 v7.10.0]
ftn95 celij_hb_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0032) integer i, j1, j2, j, jj
WARNING - The argument J1 has not been used
NO ERRORS, 1 WARNING [<CELIJ_HB> FTN95/Win32 v7.10.0]
ftn95 source_slope_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<SOURCE_SLOPE> FTN95/Win32 v7.10.0]
ftn95 ac_b_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<AC_B> FTN95/Win32 v7.10.0]
ftn95 celij_b_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0031) integer i, j1, j2, jj, j
WARNING - The argument J1 has not been used
NO ERRORS, 1 WARNING [<CELIJ_B> FTN95/Win32 v7.10.0]
ftn95 limiter_noMUSCL.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0035) real r, fi
WARNING - The argument R has not been used
NO ERRORS, 1 WARNING [<LIMITER> FTN95/Win32 v7.10.0]
ftn95 open_bc_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<OPEN_BC> FTN95/Win32 v7.10.0]
0176) integer i, j1, j2, jj, j
WARNING - The argument J1 has not been used
NO ERRORS, 1 WARNING [<CELIJ_OPENBC> FTN95/Win32 v7.10.0]
ftn95 recount_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0034) real xp_n(npar),up_n(npar),uo_n(npar),pm_n(npar),h_var_n(npar),
0035) + rhop_n(npar), dw_n(npar), alphap_n(npar), cs_n(npar),
0036) + alphapo_n(npar), dwo_n(npar), h_var0_n(npar), dw0_n(npar),
0037) + ar_n(npar), eta0_n(npar), pm_eta_n(npar)
WARNING - Variable ALPHAPO_N has been declared but not used
WARNING - Variable DWO_N has been declared but not used
WARNING - Variable AR_N has been given a value but never used
WARNING - Variable ETA0_N has been declared but not used
WARNING - Variable PM_ETA_N has been declared but not used
NO ERRORS, 5 WARNINGS [<RECOUNT> FTN95/Win32 v7.10.0]
ftn95 celij_dw_ob_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0031) integer i, j1, j2, j, jj
WARNING - The argument J1 has not been used
NO ERRORS, 1 WARNING [<CELIJ_DW_OB> FTN95/Win32 v7.10.0]
ftn95 celij_alpha_ob_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0031) integer i, j1, j2, j, jj
WARNING - The argument J1 has not been used
NO ERRORS, 1 WARNING [<CELIJ_ALPHA_OB> FTN95/Win32 v7.10.0]
ftn95 divide_ob_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<DIVIDE_OB> FTN95/Win32 v7.10.0]
ftn95 open_bc_pos_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
NO ERRORS [<OPEN_BC_POS> FTN95/Win32 v7.10.0]
ftn95 SPHYSICS_SWE_1D.f /OPTIMISE
[FTN95/Win32 Ver. 7.10.0 Copyright (c) Silverfrost Ltd 1993-2014]
0049) write(*,*) 'the FreeBSD License;
*** Unterminated character constant
1 ERROR [<SPH> FTN95/Win32 v7.10.0]
*** Compilation failed
slink global_1D.obj ini_divide_1D.obj SPHYSICS_SWE_1D.obj getdata_
1D.obj check_limits_1D.obj poute_1D.obj divide_1D.obj variable_time_st
ep_1D.obj ac_1D.obj ac_dw_var_1D.obj celij_dw_1D.obj
self_dw_1D.obj kernel_cubic_1D.obj step_leap_frog_1D.obj ghost_1D.obj
ac_alpha_1D.obj celij_alpha_1D.obj self_alpha_1D.obj grid_h_v
ar_1D.obj celij_visc_1D.obj self_visc_1D.obj celij_ob_1D.obj
viscosity_LF_1D.obj divide_b_1D.obj bottom_1D.obj celij_hb
_1D.obj source_slope_1D.obj ac_b_1D.obj celij_b_1D.obj limiter_
noMUSCL.obj open_bc_1D.obj recount_1D.obj celij_dw_ob_1D.obj celij_al
pha_ob_1D.obj divide_ob_1D.obj open_bc_pos_1D.obj -FILE:SPHYSICS_SWE_1D
.exe
*** Could not open: sphysics_swe_1d.obj
^C
SPHYSICScompilation FAILED
Check you have specified the correct compiler in Case file
henri8388
 
Posts: 8
Joined: Tue Sep 30, 2014 9:14 am


Return to SPHysics User Discussions

Who is online

Users browsing this forum: No registered users and 2 guests

cron