Initialisations manquantes de mass_solu_aero() et mass_solu_aero_pi() ?

[Ehouarn Millour]

Bug remonté des tests avec dynamico-aquaplanète de Yann. A vérifier, mais il semble qu'il manque dans physiq_mod.F90 les initialisations de mass_solu_aero() et mass_solu_aero_pi()...

Proposition de la modification:

Ligne 3301+ :

ELSE

tausum_aero(:,:,:) = 0.

+ mass_solu_aero(:,:)=0 + mass_solu_aero_pi(:,:)=0

IF (iflag_rrtm .EQ. 0) THEN !--old radiation

tau_aero(:,:,:,:) = 1.e-15 piz_aero(:,:,:,:) = 1. cg_aero(:,:,:,:) = 0.

ELSE

tau_aero_sw_rrtm(:,:,:,:) = 1.e-15 tau_aero_lw_rrtm(:,:,:,:) = 1.e-15 piz_aero_sw_rrtm(:,:,:,:) = 1.0 cg_aero_sw_rrtm(:,:,:,:) = 0.0

ENDIF

ENDIF

[Ehouarn Millour]

Pb de mise en forme ci-dessus, pour la partie code lire

Ligne 3301+ :
        ELSE
           tausum_aero(:,:,:) = 0.
+         mass_solu_aero(:,:)=0
+         mass_solu_aero_pi(:,:)=0
           IF (iflag_rrtm .EQ. 0) THEN !--old radiation
              tau_aero(:,:,:,:) = 1.e-15
              piz_aero(:,:,:,:) = 1.
              cg_aero(:,:,:,:)  = 0.
           ELSE
              tau_aero_sw_rrtm(:,:,:,:) = 1.e-15
              tau_aero_lw_rrtm(:,:,:,:) = 1.e-15
              piz_aero_sw_rrtm(:,:,:,:) = 1.0
              cg_aero_sw_rrtm(:,:,:,:)  = 0.0
           ENDIF
        ENDIF

[Ehouarn Millour]

Solved with r2623.

Note that in practice it was not a bug; mass_solu_aero() and mass_solu_aero() were either not initialized but not used, or initialized (elsewhere) when used.