open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

We are implementing a new two-class population balance model for flocculation in Telemac-3D. Now the most part is done and we are trying to put additional variables, such as floc size and floc desity in the result files.

We checked the telemac3d.dico. Inside, we found there are "private variables (P1, P2, P3 and P4)" available for the VARIABLES FOR 3D GRAPHIC PRINTOUTS. Can these be used as customized variables for output?

We also found that the way of writing the variables to the result files:
In telemac3d.f, it calls BIEF_DESIMP --> WRITE_DATA --> ADD_DATA --> ADD_DATA_SRF

It seems that the values are stored in VARSO3(3D) and VARSOR(2D) and then they are passed into BIEF_DESIMP and so on. But We couldn't find the way of giving values to VARSO3 and VARSOR in telemac3d.f. If we have TA1 (tracer values) in the keyword VARIABLES FOR 3D GRAPHIC PRINTOUTS, in which subroutine they are passed to VARSO3? How does Telemac write them in the result files?

Thanks in advance!

Qilong

Hello Qilong,

VARSO3 and VARSOR are blocks which are filled in subroutine point_telemac3d.f by successively calling subroutine ADDBLO. This is where you can add a new variable, preferably after all the others, otherwise it will shift all others and it can be a nightmare to deal with. See also NOMVAR_TELEMAC3D and NOMVAR_2D_IN_3D where the same order is taken for names and wre MNEMOs are chose for the VARIABLES FOR GRAOHIC PRINTOUTS. Whether these variables should be put in the results file is decided in lecdon_telemac3d.f (see the 2 calls to subroutine SORTIE (French word for exit...))

With best regards,

Jean-Michel Hervouet

Thank you very much again, Jean-Michel!

I followed your instruction and now I understand better. I'm able to generate the variables I want by using the work array PRIVE defined in the keyword VARIABLES FOR 3D GRAPHIC PRINTOUTS.

I simply give P1, P2... as additional variable for the VARIABLES FOR 3D GRAPHIC PRINTOUTS. Then in the code, I give values of other variables to PRIVE%ADR(1)%P, PRIVE%ADR(2)%P, and so on. And it works!

Best regards,
Qilong

Hello Qilong,

Good idea, we could put P1,... P4 instead of N, O, R Z which I never remember. Note that P* can exit all variables P1 up to P9...

With best regards,

Jean-Michel Hervouet