how to specify the sediment source? (Solved)

edit

Hello everybody

Thanks to Andres and the team of CROCO developers. Whit your help i did a configuration successfull.

suggesting: in current version of CROCO (v1.1) it is not possible to read sediment variable from croco_runoff.nc. for this reason I suggest that you do not use this option. I tried but had many problems.

Solution:
1.- define the following in your cppdefs.h

# define SEDIMENT
#define BEDLOAD
#define PSOURCE

2.- in croco.in modify this section

psource: Nsrc Isrc Jsrc Dsrc Qbar [m3 / s] Lsrc Tsrc
1
53 51 0 -10000. T T T T 20. 2. 9.9 9.8

where: Nsrc = Number of sources or rivers (1), Isrc and Jsrc = source position in xi and eta (53 51), Dsrc = source direction (0 = zonal), Qbar = flow (-10000, negative value indicate that flow is north / south), Lsrc = tracers (Temp, Salt, Sand_1, Sand_2) and Tsrc = tracers value (20. 2. 9.9 9.8).

3- in sediment.in modify value of CSED

Sd (1-NST), CSED, SRHO, WSED, ERATE, TAU_CE, TAU_CD, BED_FRAC (1: NLAY)
0.125 0.0 2650. 9.4 25.0e-5 0.05 0.14 0.4 0.4
0.050 0.0 2650. 1.6 4.0e-5 0.01 0.14 0.6 0.6

Change both value to zero. with this change i avoided a lot problems in my model.
Clarifying that this configuration it is for when the sediment comes from a river.

Hi Andrés

1.- this is my croco.im file for this configuration model: I defined SEDIMENT, BBL and PSOURCE C:\fakepath\croco.in

2.- on my second try I actived SEDIMENT, BBL, PSOURCE_NCFILES and PSOURCE_NCFILE_TS

this is my croco.in C:\fakepath\croco.in and this is my cppdefs.h C:\fakepath\cppdefs.h and this is mi out of model C:\fakepath\16135786999468359.txt

I specified one river whit variable caudal (~500 m3/s) in a Netcdf file. my grid size is of 54 x 75.

Thaks

Hi Andrés and the whole community

I did what you recommended but i had a few problems.

I defined SEDIMENT, BBL , PSOURCE on my simulation and i modified croco.in specifying tracers. I have this problem: " READ_INP ERROR while reading block with keyword 'psource'." I tryed many things but nothing worked. I attach the files C:\fakepath\out_Psorce_Sediment.txt C:\fakepath\cppdefs_Psorce_Sediment.h

I did can not resolve this problem then i decided use make_runoff.m to create croco_runoff.nc. I defined SEDIMENT, BBL , PSOURCE, PSOURCE_NCFILE_TS on my simulation and i modified croco.in specifying psource_ncfile. I have this problem: " STEP2D: ABNORMAL JOB END BLOW UP " . I attach the files. C:\fakepath\out_ncPsource_Sediment.txt C:\fakepath\cppdefs_ncPsource_Sediment.h

Finally I defined only SEDIMENT and undefine BBL. with this configuration the simulation did successful but it is not what i expected. The river is not sediment source, apparently only the bottom sediment is moving. I attach the files. C:\fakepath\out_ncPsorce_Sediment_NOBBL.txt C:\fakepath\cppdefs_ncPsorce_Sediment_NOBBL.h. C:\fakepath\Screenshot from 2021-02-17 10-38-33.png C:\fakepath\Screenshot from 2021-02-17 10-39-18.png

Now my questions are: How do I add a new tracer in make_runoff.m? (actually the script only add T,S or/and NO3) What name should the sediment tracer variable have?

Note: in croco.in line 234 define that sediment.in is a input file (i not modified nothing)

I attach general files used in my simulation.C:\fakepath\croco.in C:\fakepath\sediment.in

Thank you

Thank you so much, I will try