Pierre T. 12/07/2024: bugfix use of JP_DMS even if not present in BASIC for DMS fluxes

src/MNH/ground_paramn.f90

@@ -1072,14 +1072,14 @@ END IF
 IF (LUSECHEM) THEN
    DO JSV=NSV_CHEMBEG,NSV_CHEMEND
       PSFSV(:,:,JSV) = ZSFTS(:,:,JSV) * XMD / ( XAVOGADRO * XRHODREF(:,:,IKB))
-      IF ((LCHEMDIAG).AND.(CPROGRAM == 'DIAG  ')) XCHFLX(:,:,JSV-NSV_CHEMBEG+1) = PSFSV(:,:,JSV)    
-   END DO
-    IF (CPROGRAM == 'MESONH') THEN
+    IF ((LCHEMDIAG).AND.(CPROGRAM == 'DIAG  ')) XCHFLX(:,:,JSV-NSV_CHEMBEG+1) = PSFSV(:,:,JSV)    
+    IF ((CPROGRAM == 'MESONH').AND.(TRIM(CNAMES(JSV-NSV_CHEMBEG+1))=='DMS')) THEN
     ! net flux of so2 in kg/m2/s
-    XFLX_DMS(:,:) = ZSFTS(:,:,NSV_CHEMBEG-1+JP_DMS) *62.13E-3 / XAVOGADRO
+     XFLX_DMS(:,:) = ZSFTS(:,:,JSV) *62.13E-3 / XAVOGADRO
     ! total net flux of so2 in kg/m2
-    XFLXT_DMS(:,:) = XFLXT_DMS(:,:)+ ZSFTS(:,:,NSV_CHEMBEG-1+JP_DMS)*XTSTEP
+     XFLXT_DMS(:,:) = XFLXT_DMS(:,:)+ ZSFTS(:,:,JSV)*XTSTEP
     END IF
+   END DO
 
 ELSE
   PSFSV(:,:,NSV_CHEMBEG:NSV_CHEMEND) = 0.