Should be like this. Will check after if some matrix is rotated between Matlab and Python code. l’Idea dietro a questa operazione e’ la seguente:

- lo jacobiano e’ la ∂R/∂x  (dove R e’ la radianza e x e uno dei livelli del vapor d’acqua);

- il forwad model calcola lo jacobiano usando per x le unita’ di Kg/Kg (ovvero q il mixing ratio), mentre nello state vector usato nell’inversione il vapor d’acqua e’ espresso il log(q).

- per trasformare lo jacobiano in q, nello jacobiano in log(q) si usa  ∂R/∂(ln(x))= x*∂R/∂x;

quindi lo jacobiano in output dal forward model, che ha per ogni canale 61 livelli, deve essere moltiplicato per il profilo stesso del vapor d’acqua (per ogni canale).