open TELEMAC-MASCARET The mathematically superior suite of solvers

Thanks for reading my question,
1.When using C-VSM in GAIA，I defined 100 sections,but in VSP.CSV files,there only at most 58 sections？
2.Thanks to Pham at my last question,I tried the three ways to open the vsp_res.slf file,but I found that it can only be displayed using 2D format,I wonder if it can be seen in 3D,especially I really need the vertical profile in C-VSM...if it can't open directly,I wonder how to define to printout the VSP.CVS file of a specific point even it'not a node with a number?
Thank you for your reply.

Hello,

the definition of 100 section is the max. number of possible vertical sections. During simulation new sections are build in case of sedimentation. If 100 sections are reached an algorithm is used to reduce the number of sections <= 100.

Here are 2 papers about cvsm from the telemac conferences:
henry.baw.de/bitstream/handle/20.500.119...quence=1&isAllowed=y
www.uwe-merkel.com/wp-content/uploads/20...5/02_Merkel_2012.pdf

The vsp file is a 3D SERAFIN / SERAFIND file. I can visualise it with 3D in tecplot. But be careful, contrary to the normal 3D plot of water, the 3D "underground" cells are only connected in the vertical and not necessarily in the horizontal. So a normal interpolation will not show a correct picture. The easiest way to visualise is a scatter plot in which no interpolation between points will be made.


More complicated and nicer of course is a visulisation with voronoi-cells. See the last figure on the website of Uwe Merkel: www.uwe-merkel.com/2020/05/07/a-new-cont...lemac-suite-cvsm-ii/
I guess you can do a similar plot with python but I haven't done it yet.

Best regards,
Rebekka

Thank you very much for your reply,It really helped much. And sorry for bothering,I have another question to ask.
There are too many VSP.CVS files in the -tmp file... When I set the time step to 0.5,it generates a chart every 0.5,in the gaia.cas I only see C-VSM FULL PRINTOUT PERIOD,I wonder how to set the interval of output file?
Thank you again.

Hello,

I have checked and think it was probably wanted to have an output every time step for the output per node. If you want to have it only for the CVSM PRINTOUT PERIOD you can do the following:

Insert the 2 lines to cvsp_main_gaia.f in the beginning of the subroutine (e.g. line 75)
CVSM_OUT = .FALSE.
IF(CVSMPPERIOD*(LT/CVSMPPERIOD).EQ.LT) CVSM_OUT = .TRUE.
and add CVSMPPERIOD to the USE command:
USE DECLARATIONS_GAIA, ONLY: CVSMOUTPUT,CVSM_OUT,CVSM_OUT_FULL,
& PRO_D,PRO_MAX,PRO_MAX_MAX,PERCOU,
& HN,LT,DT,MESH,Z,PRO_F,ZR,
& CVSMPPERIOD

Best regards,
Rebekka

Hi,
Thank you for the reply. Did you tried the solution you provided and did it work? I tried to add the content you told me,but it failed,it still genarates a file every time step...

Best wishes,
Xue

Hi Xue,

I tried the solution in the newest version and it worked in my case for C-VSM FULL PRINTOUT PERIOD=0 (then the graphic printout period is used) or with an own printout period (>0).
I have added you the changed subroutine. Maybe I missed something in my description.
Best regards,
Rebekka

Hello again,
I used your program and set the C-VSM FULL PRINTOUT PERIOD=0 in gaia.cas,also GRAPHIC PRINTOUT PERIOD=18 in t2d.cas. But it still genarates a VSP.CSV file every time step in the tmp folder.
Sorry for my little talent and less learing I really don't know where the problem is...
Thanks for your precious time.

Sincerely,
Xue.

Hi,
Ithink I made mistakes,maybe I can't change the output in tmp file...What you told me is about the change in .slf file. I referred to the papers,I may change the calculate total time to save the disk space.

Thank you very much and wish you have a good day.
Xue

Hi Xue,

it should work, I don't see why it is not working... And it is about the files in the tmp-directory. Are you sure that you use the new subroutine?
You can check by using the print command. The output should occur in the listing otherwise you don't use the subroutine: print*,'cvsmout',cvsm_out,cvsp_db_gaia(-1,-1),cvsm_out_full,
& LT,CVSMPPERIOD,CVSMPPERIOD*(LT/CVSMPPERIOD)
cvsm_out should only be true in case of printout period.

The files in the tmp files are written from the commands in line 259-264 and line 277-284

Best regards,
Rebekka

Hi,
Last days I have a ddl so I didn't pay attention to it .Sorry for repling after such a long time..I really used the subroutine and set the C-VSM FULL PRINTOUT PERIOD=0...I trid the print sentence,but it didn't work.Here I put the files on the attachments.Please help me if you have any time.
Thank you,
Xue.