Quentin 25/06/2021: bugfix introduced by commit 09f31 (OCEAN LES). Missing a...

Quentin 25/06/2021: bugfix introduced by commit 09f31 (OCEAN LES). Missing a CALL for the atmospheric case for computation of the pressure

src/MNH/set_mass.f90

@@ -572,6 +572,9 @@ ELSE
 ! ATmos: rho = P/ (R Tv)
  XRHODREFZ(:) = XP00/ (XRD* XTHVREFZ(:))
  ZTHVREF3D(:,:,:)=XTHVREFZ(2)
+   CALL COMPUTE_EXNER_FROM_GROUND(ZTHVREF3D,PZFLUX_MX,&
+          ZEXNSURF2D_MX,ZHEXNFLUX,ZHEXNMASS)
+
  XEXNTOP=SUM_DD_R2_ll(ZHEXNFLUX(IIB:IIE,IJB:IJE,IKE+1))/REAL(NIMAX_ll*NJMAX_ll)
  ZEXNTOP2D=ZHEXNFLUX(:,:,IKE+1)
  CALL COMPUTE_EXNER_FROM_TOP(ZTHVREF3D,XZZ,ZEXNTOP2D,ZHEXNFLUX,ZHEXNMASS)