Philippe 10/01/2019: IO: use NEWUNIT argument of OPEN for chemistry files + bug: close all files

6b6d917a · WAUTELET Philippe · 9a39855d · 6b6d917a · 6b6d917a · 6b6d917a
Commit 6b6d917a authored 6 years ago by WAUTELET Philippe
--- a/src/MNH/ch_f77.fx90
+++ b/src/MNH/ch_f77.fx90
--- a/src/MNH/ch_init_diagnostics.f90
+++ b/src/MNH/ch_init_diagnostics.f90
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 chimie 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !!    ##############################
      SUBROUTINE CH_INIT_DIAGNOSTICS
 !!    ##############################
@@ -38,6 +33,7 @@
 !!    Original 14/05/95
 !!    27/07/96 (K. Suhre) restructured
 !!    01/12/03  (Gazen)   change Chemical scheme interface
+!!    Philippe Wautelet: 10/01/2019: use newunit argument to open files
 !!
 !!    EXTERNAL
 !!    --------
@@ -82,11 +78,11 @@ ALLOCATE(IINDEX(2,NNONZEROTERMS))
 !     ---------------------
 !
 ! open the file
+OPEN(NEWUNIT =  NDIAGIO,    &
+     FILE    =  CDIAGFILE,  &
+     FORM    = "FORMATTED", &
+     STATUS  = "UNKNOWN"    )
 PRINT *, "CH_INIT_DIAGNOSTICS: opening unit ", NDIAGIO, " for file ", CDIAGFILE
-OPEN(UNIT   =  NDIAGIO,    &
-     FILE   =  CDIAGFILE,  &
-     FORM   = "FORMATTED", &
-     STATUS = "UNKNOWN"    )
 !
 ! write the headder
 CALL DATE_AND_TIME(YDATE, YTIME)

--- a/src/MNH/ch_init_model0d.f90
+++ b/src/MNH/ch_init_model0d.f90
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 chimie 2006/12/21 15:21:50
-!-----------------------------------------------------------------
 !!    ###########################
      MODULE MODI_CH_INIT_MODEL0D
 !!    ###########################
@@ -50,6 +45,7 @@ END MODULE MODI_CH_INIT_MODEL0D
 !!    19/03/96 (K. Suhre) insert namelists directly for compiler reasons
 !!    31/08/96 (K. SUHRE) read MODD_BLANK namelist
 !!    01/08/01 (C. Mari)  change NAM_CH_SOLVER to $n 
+!!    Philippe Wautelet: 10/01/2019: use newunit argument to open files
 !!
 !!    EXTERNAL
 !!    --------
@@ -74,6 +70,9 @@ USE MODN_CH_SOLVER_n
 IMPLICIT NONE

 CHARACTER*(*), INTENT(IN) :: HNAMELISTFILE ! name of namelist input file
+
+INTEGER :: ILU ! unit number for IO
+
 NAMELIST /NAM_CH_MODEL0D/ XTBEGIN, XTEND, XDTACT,          &
                       XDTOUT, XDTDIAG,                 &
                       CRUNID,                          &
@@ -107,10 +106,10 @@ CALL INIT_NAM_CH_SOLVERn


 PRINT *, 'CH_INIT_MODEL0D: opening namelist file: ', HNAMELISTFILE
-OPEN(UNIT=NFILEIO,FILE=HNAMELISTFILE,STATUS="OLD",FORM="FORMATTED")
-READ(NFILEIO,NAM_CH_MODEL0D)
-READ(NFILEIO,NAM_CH_SOLVERn)
-CLOSE(NFILEIO)
+OPEN(NEWUNIT=ILU,FILE=HNAMELISTFILE,STATUS="OLD",FORM="FORMATTED")
+READ(UNIT=ILU,NML=NAM_CH_MODEL0D)
+READ(UNIT=ILU,NML=NAM_CH_SOLVERn)
+CLOSE(UNIT=ILU)

 ! quick fix for use of the surface scheme in the box model
 CCHEM_INPUT_FILE = 'CHCONTROL1.nam'

--- a/src/MNH/ch_init_output.f90
+++ b/src/MNH/ch_init_output.f90
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 chimie 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !!    ##########################
      MODULE MODI_CH_INIT_OUTPUT
 !!    ##########################
@@ -54,6 +49,7 @@ END MODULE MODI_CH_INIT_OUTPUT
 !!    Original 25/04/95
 !!    27/07/96 (K. Suhre) restructured
 !!    01/12/03 (D. Gazen) change Chemical scheme interface
+!!    Philippe Wautelet: 10/01/2019: use newunit argument to open files
 !!
 !!    EXTERNAL
 !!    --------
@@ -88,11 +84,11 @@ INTEGER                :: NAERO
 !*       1.   OPEN OUTPUT FILE
 !        ----------------------
 !
+OPEN(NEWUNIT = NRESULTIO,   &
+     FILE    = CRESULTFILE, &
+     FORM    = "FORMATTED", &
+     STATUS  = "UNKNOWN"    )
 PRINT *, "CH_INIT_OUTPUT: opening unit ", NRESULTIO, " for file ", CRESULTFILE
-OPEN(UNIT   = NRESULTIO,   &
-     FILE   = CRESULTFILE, &
-     FORM   = "FORMATTED", &
-     STATUS = "UNKNOWN"    )
 !
 !*       2.   WRITE HEADDER
 !        ------------------

--- a/src/MNH/ch_make_lookup.f90
+++ b/src/MNH/ch_make_lookup.f90
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1999-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 chimie 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !      ######################
       PROGRAM CH_MAKE_LOOKUP
 !      ######################
@@ -57,6 +52,7 @@
 !!    MODIFICATIONS
 !!    -------------
 !!    Original 01/03/99
+!!    Philippe Wautelet: 10/01/2019: use newunit argument to open files
 !!
 !!    EXTERNAL
 !!    --------
@@ -92,6 +88,7 @@ CHARACTER*40, DIMENSION(NJOUT)        :: JLABELOUT
 !
 CHARACTER*120 :: HEADDER
 REAL          :: UT
+INTEGER       :: ILU ! unit number for IO
 INTEGER       :: I, J, K, NJIO
 CHARACTER*40  :: YFMT = '(2F11.2,5E11.4/99(7E11.4/))'
 !
@@ -122,11 +119,11 @@ DOBNEW = -1.
 !
 !      read the namelist
 WRITE(*,*) "TRYING TO OPEN FILE LOOKUP1.nam ..."
-OPEN(UNIT=42,FILE="LOOKUP1.nam",STATUS="OLD",FORM="FORMATTED", ERR=100)
-READ(42,NAM_TUV)
+OPEN(NEWUNIT=ILU,FILE="LOOKUP1.nam",STATUS="OLD",FORM="FORMATTED", ERR=100)
+READ(UNIT=ILU,NML=NAM_TUV)
 WRITE(*,*) "NAMELIST NAM_TUV INITIALIZED TO:"
 WRITE(*,NAM_TUV)
-CLOSE(42)
+CLOSE(UNIT=ILU)
 GOTO 200
 !
 100    CONTINUE
@@ -171,21 +168,21 @@ ENDDO
 !*       3.   WRITE LOOKUP TABLE FILE
 !        -----------------
 !
-OPEN(66,FILE="PHOTO.TUV39",STATUS="UNKNOWN",FORM="FORMATTED")
-WRITE(66,'(A)') HEADDER
-WRITE(66,'(I6,A)') NJIO,   "  NUMBER OF PHOTOLYSIS REACTIONS"
-WRITE(66,'(I6, F10.0, A)') NLEVEL, DZ, &
+OPEN(NEWUNIT=ILU,FILE="PHOTO.TUV39",STATUS="UNKNOWN",FORM="FORMATTED")
+WRITE(ILU,'(A)') HEADDER
+WRITE(UNIT=ILU, FMT='(I6,A)') NJIO,   "  NUMBER OF PHOTOLYSIS REACTIONS"
+WRITE(UNIT=ILU, FMT='(I6, F10.0, A)') NLEVEL, DZ, &
      "  NUMBER OF LEVELS, VERTICAL INCREMENT (M)"
-WRITE(66,'(I6, F10.4, A)') NTIME,  DT, &
+WRITE(UNIT=ILU, FMT='(I6, F10.4, A)') NTIME,  DT, &
      "  NUMBER OF TEMPORAL RECS, TEMPORAL INCREMENT (H)"
-WRITE(66,'(A)') (JLABELOUT(K), K=1, NJIO)
-WRITE(66,'(A)') YFMT
+WRITE(UNIT=ILU, FMT='(A)') (JLABELOUT(K), K=1, NJIO)
+WRITE(UNIT=ILU, FMT='(A)') YFMT
 DO J = 1, NLEVEL
  DO I = 1, NTIME
-    WRITE(66,YFMT) AZ(J), ATIME(I), (JDATA(J,K,I), K=1, NJIO)
+    WRITE(UNIT=ILU, FMT=YFMT) AZ(J), ATIME(I), (JDATA(J,K,I), K=1, NJIO)
  ENDDO
 ENDDO
-CLOSE(66)
+CLOSE(UNIT=ILU)
 !
 PRINT *, 'Lookup table file PHOTO.TUV39 has been generated.'
 PRINT *, 'CH_MAKE_LOOKUP ended correctly.'

--- a/src/MNH/ch_read_chem.f90
+++ b/src/MNH/ch_read_chem.f90
-!MNH_LIC Copyright 1994-2018 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
 !MNH_LIC for details. version 1.
@@ -56,6 +56,7 @@ END MODULE MODI_CH_READ_CHEM
 !!    01/12/03 (D. Gazen)   change Chemical scheme interface
 !!    M.Leriche 2015 : masse molaire Black carbon à 12 g/mol
 !!    Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
+!!    Philippe Wautelet: 10/01/2019: use newunit argument to open files
 !!
 !!    EXTERNAL
 !!    --------
@@ -66,7 +67,7 @@ USE MODE_FM,    ONLY: IO_FILE_CLOSE_ll
 !!
 !!    IMPLICIT ARGUMENTS
 !!    ------------------
-USE MODD_CH_MODEL0D, ONLY: NFILEIO, NVERB
+USE MODD_CH_MODEL0D, ONLY: NVERB
 USE MODD_CH_M9_n, ONLY:      NEQ, CNAMES
 USE MODD_CH_AEROSOL
 !!
@@ -83,6 +84,7 @@ CHARACTER(LEN=*), INTENT(IN)      :: HFILE ! name of the file to be read from
 !!    ------------------------------
 CHARACTER(LEN=32) :: YVARNAME
 CHARACTER(LEN=80) :: YINPUT
+INTEGER :: ILU ! unit number for IO
 INTEGER :: JI, JJ, IIN
 REAL :: ZMD
 REAL, DIMENSION(NSP+NCARB+NSOA) :: ZMI ! aerosol molecular mass in g/mol
@@ -114,15 +116,15 @@ ELSE
 ! open file 
 !
  PRINT *, 'CH_READ_CHEM: opening file ', HFILE
-  OPEN(UNIT = NFILEIO, &
-       FILE = HFILE, &
-       FORM = 'FORMATTED', &
-       STATUS = 'OLD')
+  OPEN(NEWUNIT = ILU,         &
+       FILE    = HFILE,       &
+       FORM    = 'FORMATTED', &
+       STATUS  = 'OLD')
 !
 ! read line by line and check variable names
 !
  outer_loop : DO JI = 1, NEQ
-    READ(NFILEIO,*,END=999) YVARNAME, PCONC(JI)
+    READ(UNIT=ILU,FMT=*,END=999) YVARNAME, PCONC(JI)
    check_loop : DO JJ = 1, 32
      IF (YVARNAME(JJ:JJ).NE.CNAMES(JI)(JJ:JJ)) THEN
        PRINT *, 'CH_READ_CHEM: Error: variable names do not match:'
@@ -139,7 +141,7 @@ ELSE
 PCONC(:) =  PCONC(:) * 1E-9
 IF (LORILAM) THEN
  outer_loop2 : DO JI = 1, SIZE(PAERO,1)
-    READ(NFILEIO,*,END=997) YVARNAME, PAERO(JI)
+    READ(UNIT=ILU,FMT=*,END=997) YVARNAME, PAERO(JI)
    check_loop2 : DO JJ = 1, 32
      IF (YVARNAME(JJ:JJ).NE.CAERONAMES(JI)(JJ:JJ)) THEN
        PRINT *, 'CH_READ_CHEM: Error: variable names do not match:'
@@ -260,7 +262,7 @@ END IF
 END IF
 !
 ! close file
-  CLOSE(NFILEIO)
+  CLOSE(UNIT=ILU)
 !

 END IF

--- a/src/MNH/ch_surface0d.f90
+++ b/src/MNH/ch_surface0d.f90
-!MNH_LIC Copyright 1994-2018 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1999-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
 !     ####################################################
@@ -39,6 +39,7 @@
 !!    -------------
 !!    Original 03/03/99
 !!  Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
+!!  Philippe Wautelet: 10/01/2019: use newunit argument to open files
 !!
 !!    EXTERNAL
 !!    --------
@@ -51,7 +52,7 @@ USE MODI_CH_EMISSION_FLUX0D
 !!
 !!    IMPLICIT ARGUMENTS
 !!    ------------------
-USE MODD_CH_MODEL0D, ONLY: NFILEIO, XDTOUT, XTBEGIN, NVERB
+USE MODD_CH_MODEL0D, ONLY: XDTOUT, XTBEGIN, NVERB
 USE MODD_CH_M9_n,    ONLY: CNAMES
 !!
 !!    EXPLICIT ARGUMENTS
@@ -64,6 +65,7 @@ REAL, DIMENSION(KEQ), INTENT(INOUT) :: PCONC ! concentration vector
 !
 !!    DECLARATION OF LOCAL VARIABLES
 !!    ------------------------------
+INTEGER, SAVE :: ILU = -1 ! unit number for IO
 INTEGER :: JI
 REAL, DIMENSION(KEQ) :: ZDEPO  ! deposition velocity
 REAL, DIMENSION(KEQ) :: ZEMIS  ! emission flux after multiplication 
@@ -174,49 +176,48 @@ PCONC(:) = (PCONC(:) + PDTACT*ZEMIS(:)/(100.*ZHEIGHT)) &
 !
 IF (LSFIRST) THEN
  ZSTNEXTOUT = XTBEGIN
-  OPEN(UNIT   = 86,            &
-       FILE   = "SURFACE.out", &
-       FORM   = "FORMATTED",   &
-       STATUS = "UNKNOWN"      )
-  WRITE(86,'(A)') 'parameters from the Wesley scheme'
-  WRITE(86,'(A)') 'quick & dirty coding :-)'
-  WRITE(86,'(2I5,A)') 1+2*KEQ+21, KEQ, ' 0 0'
-  WRITE(86,'(A)') 'XTSIMUL'
+  OPEN(NEWUNIT = ILU,            &
+       FILE    = "SURFACE.out", &
+       FORM    = "FORMATTED",   &
+       STATUS  = "UNKNOWN"      )
+  WRITE(UNIT=ILU, FMT='(A)') 'parameters from the Wesley scheme'
+  WRITE(UNIT=ILU, FMT='(A)') 'quick & dirty coding :-)'
+  WRITE(UNIT=ILU, FMT='(2I5,A)') 1+2*KEQ+21, KEQ, ' 0 0'
+  WRITE(UNIT=ILU, FMT='(A)') 'XTSIMUL'
  DO JI = 1, KEQ
-    WRITE (86,'(2A)') "VDEP-", CNAMES(JI)
+    WRITE (UNIT=ILU, FMT='(2A)') "VDEP-", CNAMES(JI)
  END DO
  DO JI = 1, KEQ
-    WRITE (86,'(2A)') "EMIS-", CNAMES(JI)
+    WRITE (UNIT=ILU, FMT='(2A)') "EMIS-", CNAMES(JI)
  END DO
-  WRITE(86,'(A)') 'VMOD'
-  WRITE(86,'(A)') 'RHODREF'
-  WRITE(86,'(A)') 'RS'
-  WRITE(86,'(A)') 'RESA'
-  WRITE(86,'(A)') 'THT2D'
-  WRITE(86,'(A)') 'PABST2D'
-  WRITE(86,'(A)') 'CD'
-  WRITE(86,'(A)') 'CLAY'
-  WRITE(86,'(A)') 'SAND'
-  WRITE(86,'(A)') 'LAND'
-  WRITE(86,'(A)') 'VEG'
-  WRITE(86,'(A)') 'LAI'
-  WRITE(86,'(A)') 'SST'
-  WRITE(86,'(A)') 'TS'
-  WRITE(86,'(A)') 'HU'
-  WRITE(86,'(A)') 'PSN'
-  WRITE(86,'(A)') 'PSNG'
-  WRITE(86,'(A)') 'PSNV'
-  WRITE(86,'(A)') 'Z0VEG'
-  WRITE(86,'(A)') 'HEIGHT'
-  WRITE(86,'(A)') 'EMISFACTOR'
-print *, "209"
-  WRITE(86,'(A)') '(5E16.8)'
+  WRITE(UNIT=ILU, FMT='(A)') 'VMOD'
+  WRITE(UNIT=ILU, FMT='(A)') 'RHODREF'
+  WRITE(UNIT=ILU, FMT='(A)') 'RS'
+  WRITE(UNIT=ILU, FMT='(A)') 'RESA'
+  WRITE(UNIT=ILU, FMT='(A)') 'THT2D'
+  WRITE(UNIT=ILU, FMT='(A)') 'PABST2D'
+  WRITE(UNIT=ILU, FMT='(A)') 'CD'
+  WRITE(UNIT=ILU, FMT='(A)') 'CLAY'
+  WRITE(UNIT=ILU, FMT='(A)') 'SAND'
+  WRITE(UNIT=ILU, FMT='(A)') 'LAND'
+  WRITE(UNIT=ILU, FMT='(A)') 'VEG'
+  WRITE(UNIT=ILU, FMT='(A)') 'LAI'
+  WRITE(UNIT=ILU, FMT='(A)') 'SST'
+  WRITE(UNIT=ILU, FMT='(A)') 'TS'
+  WRITE(UNIT=ILU, FMT='(A)') 'HU'
+  WRITE(UNIT=ILU, FMT='(A)') 'PSN'
+  WRITE(UNIT=ILU, FMT='(A)') 'PSNG'
+  WRITE(UNIT=ILU, FMT='(A)') 'PSNV'
+  WRITE(UNIT=ILU, FMT='(A)') 'Z0VEG'
+  WRITE(UNIT=ILU, FMT='(A)') 'HEIGHT'
+  WRITE(UNIT=ILU, FMT='(A)') 'EMISFACTOR'
+  WRITE(UNIT=ILU, FMT='(A)') '(5E16.8)'
  LSFIRST = .FALSE.
 END IF

 IF (PTSIMUL >= ZSTNEXTOUT) THEN 
  ZSTNEXTOUT = ZSTNEXTOUT + XDTOUT
-  WRITE(86, '(5E16.8)') PTSIMUL, &
+  WRITE(UNIT=ILU, FMT= '(5E16.8)') PTSIMUL, &
        (ZDEPO(JI),JI=1,KEQ), &
        (ZEMIS(JI),JI=1,KEQ), &
        ZVMOD(1,1), &   
@@ -240,6 +241,7 @@ IF (PTSIMUL >= ZSTNEXTOUT) THEN
        ZZ0VEG(1,1), &   
        ZHEIGHT, &
        ZEMISFACTOR
+  CALL FLUSH(UNIT=ILU)
 END IF

 RETURN

--- a/src/MNH/ch_write_chem.f90
+++ b/src/MNH/ch_write_chem.f90
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 chimie 2006/06/30 11:47:01
-!-----------------------------------------------------------------
 !!    ######################### 
      MODULE MODI_CH_WRITE_CHEM
 !!    ######################### 
@@ -55,6 +50,7 @@ END MODULE MODI_CH_WRITE_CHEM
 !!    Original 21/04/95
 !!    01/12/03 (D. Gazen)   change Chemical scheme interface
 !!    M.Leriche 2015 : masse molaire Black carbon à 12 g/mol
+!!    Philippe Wautelet: 10/01/2019: use newunit argument to open files
 !!
 !!    EXTERNAL
 !!    --------
@@ -62,7 +58,6 @@ END MODULE MODI_CH_WRITE_CHEM
 !!
 !!    IMPLICIT ARGUMENTS
 !!    ------------------
-USE MODD_CH_MODEL0D, ONLY: NFILEIO
 USE MODD_CH_M9_n, ONLY:      NEQ, CNAMES
 USE MODD_CH_AEROSOL
 !!
@@ -76,7 +71,8 @@ CHARACTER(LEN=*), INTENT(IN)     :: HFILE ! HFILE: name of the file for write
 !
 !!    DECLARATION OF LOCAL VARIABLES
 !!    ------------------------------
-INTEGER :: JI, NAERO
+INTEGER :: ILU ! unit number for IO
+INTEGER :: JI, IAER
 REAL :: ZMD
 REAL, DIMENSION(NSP+NCARB+NSOA) :: ZMI ! aerosol molecular mass in g/mol

@@ -84,20 +80,20 @@ REAL, DIMENSION(NSP+NCARB+NSOA) :: ZMI ! aerosol molecular mass in g/mol
 !*    EXECUTABLE STATEMENTS
 !     ---------------------
 IF (LORILAM) THEN
-NAERO=SIZE(PAERO)
+IAER=SIZE(PAERO)
 ELSE
-NAERO=0
+IAER=0
 END IF
 !
 ! open file 
 PRINT *, 'CH_WRITE_CHEM: opening file ', HFILE
-OPEN(UNIT   =  NFILEIO,    &
-     FILE   =  HFILE,      &
-     FORM   = 'FORMATTED', &
-     STATUS = 'UNKNOWN'    )
+OPEN(NEWUNIT =  ILU,        &
+     FILE    =  HFILE,      &
+     FORM    = 'FORMATTED', &
+     STATUS  = 'UNKNOWN'    )
 !
 DO JI = 1, NEQ
-  WRITE(NFILEIO,'(3A,E20.8)') "'", CNAMES(JI), "' ", PCONC(JI)*1.E9 ! convert ppp to ppb
+  WRITE(UNIT=ILU,FMT='(3A,E20.8)') "'", CNAMES(JI), "' ", PCONC(JI)*1.E9 ! convert ppp to ppb
 ENDDO
 IF (LORILAM) THEN
 !Conversion  ppp to microgram/m3
@@ -206,12 +202,12 @@ IF (NSOA .EQ. 10) THEN
 END IF
 END IF

-DO JI = 1, NAERO
-  WRITE(NFILEIO,'(3A,E20.8)') "'", CAERONAMES(JI), "' ", PAERO(JI) 
+DO JI = 1, IAER
+  WRITE(UNIT=ILU,FMT='(3A,E20.8)') "'", CAERONAMES(JI), "' ", PAERO(JI)
 ENDDO

 !
 ! close file
-CLOSE(NFILEIO)
+CLOSE(UNIT=ILU)
 !
 END SUBROUTINE CH_WRITE_CHEM
--- a/src/MNH/modd_ch_model0d.f90
+++ b/src/MNH/modd_ch_model0d.f90
-!MNH_LIC Copyright 1995-2018 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
 !!    ###################### 
@@ -32,6 +32,7 @@
 !!              CCH_TUV_CLOUDS  : method for the impact of clouds on radiation
 !!    18/02/99  LCH_SURFACE0D   : apply surface emission / deposition fluxes
 !!  Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
+!!  Philippe Wautelet: 10/01/2019: use newunit argument to open files
 !!
 !!    IMPLICIT ARGUMENTS
 !!    ------------------
@@ -74,9 +75,8 @@ CHARACTER*80  :: CRUNID        = "no runid specified" ! runid for output file
 CHARACTER*40  :: CRESULTFORMAT = "(5E16.8)" ! Format for results
 CHARACTER*40  :: CDIAGFORMAT   = "(5E16.8)" ! Format for diagnostics
 !
-INTEGER :: NFILEIO   = 42 ! channel to be used for all intermediate file I/O
-INTEGER :: NRESULTIO = 44 ! channel to be used for all regular result file I/O
-INTEGER :: NDIAGIO   = 45 ! channel to be used for all diagnostics file I/O
+INTEGER :: NRESULTIO = -1 ! channel to be used for all regular result file I/O (set in CH_INIT_OUTPUT)
+INTEGER :: NDIAGIO   = -1 ! channel to be used for all diagnostics file I/O (set in CH_INIT_DIAGNOSTICS)
 !
 !*       0.3  verbosity level
 !