new feature: add multi faults tsunamis

Feature introduced in paper from TUC 2022. It required:

1. A txt file attached here that contains:
   • in the header the number of subfaults.
   • for each row the okada information for each subfault + dynamic information. Details below in the okada modified file:
   • was introduced as FORMATTED DATA FILE 1 = 'earthquakes/tohoku_fuji_50km_dyn10s.txt' see full CAS file 03_tokohu_2011_fuji_10s.cas
2. A modified version of okada.f with 12 coefs instead of the existing 10, in the same order:
   • TI = COEFS(1) ! TIME INIT - **NEW** DYNAMIC // TSR
   • TR = COEFS(2) ! RISING TIME - **NEW** DYNAMIC // TSR
   • HH = COEFS(3) ! FOCAL DEPTH
   • L = COEFS(4) ! FAULT LENGTH
   • W = COEFS(5) ! FAULT WIDTH
   • D = COEFS(6) ! DISLOCATION
   • TH = COEFS(7) * DEG2RAD ! STRIKE DIRECTION
   • DL = COEFS(8) * DEG2RAD ! DIP ANGLE
   • RD = COEFS(9) * DEG2RAD ! SLIP ANGLE
   • Y0 = COEFS(10) * DEG2RAD ! EPICENTRE LATITUDE
   • X0 = COEFS(11) * DEG2RAD ! EPICENTRE LONGITUDE
   • C0 = COEFS(12) ! SIZE OF THE ELLIPSE
3. A user input file to read the multi-faults info:
   • here is the dynamic implementation, i.e. modifying bottom over time: user_utimp_telemac2d.f
   • here is the hydrostatic implementation, i.e. translating okada fault deformation into free-surface anomaly: user_condin_h.f.
   • the best way forward would be to implement static when the first 2 parameters are equal to 0 (i.e. TIME INIT = 0 and RISING TIME = 0)

It is the first time I implement a new feature, that requires also new keywords: (we should add a new keyword like INPUT FILE FOR TSUNAMI CHARACTERISTICS or TSUNAMI INPUT FILE?)

@chi-tuan.pham, please have a look at the files and let me know what you think + how to proceed:

1. to add new keyword (+ suppress the old one?)
2. static vs dynamic implementation (mentioned above)
3. add a new test case (or modify the existing okada test case)

also pinging @sebastien.bourban, @michael.turnbull and @alexander.breugem for more thoughts on this