Planet Trap on 2D version
This is a research project by @gpichier and @bbitsch , few words of description below, sorry for the italian, Alain knows italian and I guess Bert did a lot of progress ...
La planet trap si basa su Flock et al. 2019 (https://ui.adsabs.harvard.edu/abs/2019A%26A...630A.147F/abstract). L'idea è di usare una transizione della viscosità al bordo interno del disco dove si passa dalla zona attiva in MRI alla dead zone, implementando la formula (3) di Flock et al. 2019. Su Fargo2D1D avevo modificato le funzioni Viscosity.c e Viscosity1D.c, e invece di usare una viscosità costante implementavo quella formula. Localmente funzionava molto bene, c'era una bella planet trap. Il problema era invece al bordo tra le griglie 1D e 2D, motivo per cui vogliamo passare a fargOCA. Inoltre ho l'impressione che con fargOCA sarà più naturale inserire tutti i parametri della viscosity transition nel file *.conf, il che rende il tutto molto più coerente con il resto del codice.