Changeset 3387 for LMDZ6/trunk/libf

Timestamp:

Sep 4, 2018, 7:42:50 PM (6 years ago)

Author:

oboucher

Message:

adding infocfields_init routine and
its call from physiq. This routines reads
which fields need to be transferred between
model components in ESM configuration.
No impact whatsoever in LMDZ mode.

Location:

LMDZ6/trunk/libf/phylmd

Files:

• carbon_cycle_mod.F90 (modified) (8 diffs)
• physiq_mod.F90 (modified) (9 diffs)

LMDZ6/trunk/libf/phylmd/carbon_cycle_mod.F90

@@ -1 +1 @@
 MODULE carbon_cycle_mod
-! Controle module for the carbon CO2 tracers :
+!=======================================================================
+!   Authors: Patricia Cadule and Laurent Fairhead
+!            base sur un travail anterieur mene par Patricia Cadule et Josefine Ghattas
+!
+!  Purpose and description:
+!  -----------------------
+! Control module for the carbon CO2 tracers :
 !   - Identification
 !   - Get concentrations comming from coupled model or read from file to tracers
@@ -6 +12 @@
 !   - Calculate new carbon flux for sending to coupled models (PISCES and ORCHIDEE)
 !
-! Author : Josefine GHATTAS, Patricia CADULE
+! Module permettant de mettre a jour les champs (puits et sources) pour le
+! transport de CO2 en online (IPSL-CM et LMDZOR) et offline (lecture de carte)
+!
+! Le cas online/offline est defini par le flag carbon_cycle_cpl (y/n)
+! Le transport du traceur CO2 est defini par le flag carbon_cycle_tr (y/n)
+! la provenance des champs (termes de puits) est denini par le flag level_coupling_esm
+!
+! level_coupling_esm : level of coupling of the biogeochemical fields between
+! LMDZ, ORCHIDEE and NEMO
+! Definitions of level_coupling_esm in physiq.def
+! level_coupling_esm = 0  ! No field exchange between LMDZ and ORCHIDEE models
+!                         ! No field exchange between LMDZ and NEMO
+! level_coupling_esm = 1  ! Field exchange between LMDZ and ORCHIDEE models
+!                         ! No field exchange between LMDZ and NEMO models
+! level_coupling_esm = 2  ! No field exchange between LMDZ and ORCHIDEE models
+!                         ! Field exchange between LMDZ and NEMO models
+! level_coupling_esm = 3  ! Field exchange between LMDZ and ORCHIDEE models
+!                         ! Field exchange between LMDZ and NEMO models
+!=======================================================================
 
   IMPLICIT NONE
   SAVE
   PRIVATE
-  PUBLIC :: carbon_cycle_init, carbon_cycle
+  PUBLIC :: carbon_cycle_init, carbon_cycle, infocfields_init
 
 ! Variables read from parmeter file physiq.def
@@ -18 +42 @@
   LOGICAL, PUBLIC :: carbon_cycle_cpl       ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) 
 !$OMP THREADPRIVATE(carbon_cycle_cpl)
-  INTEGER, SAVE, PUBLIC :: level_coupling_esm ! Level of coupling for the ESM - 0, 1, 2, 3
+  INTEGER, PUBLIC :: level_coupling_esm ! Level of coupling for the ESM - 0, 1, 2, 3
 !$OMP THREADPRIVATE(level_coupling_esm)
 
@@ -64 +88 @@
 !$OMP THREADPRIVATE(co2_send)
 
+! nbfields : total number of fields 
+  INTEGER, PUBLIC :: nbcf
+!$OMP THREADPRIVATE(nbcf)
+
+! nbcf_in : number of fields IN
+  INTEGER, PUBLIC  :: nbcf_in
+!$OMP THREADPRIVATE(nbcf_in)
+
+! nbcf_in_orc : number of fields IN 
+  INTEGER, PUBLIC  :: nbcf_in_orc
+!$OMP THREADPRIVATE(nbcf_in_orc)
+
+! nbcf_in_inca : number of fields IN (from INCA)
+  INTEGER, PUBLIC  :: nbcf_in_inca
+!$OMP THREADPRIVATE(nbcf_in_inca)
+
+! nbcf_in_nemo : number of fields IN (from nemo)
+  INTEGER, PUBLIC  :: nbcf_in_nemo
+!$OMP THREADPRIVATE(nbcf_in_nemo)
+
+! nbcf_in_ant : number of fields IN (from anthropogenic sources)
+  INTEGER, PUBLIC  :: nbcf_in_ant
+!$OMP THREADPRIVATE(nbcf_in_ant)
+
+! nbcf_out : number of fields OUT
+  INTEGER, PUBLIC :: nbcf_out
+!$OMP THREADPRIVATE(nbcf_out)
+
+! Name of variables
+  CHARACTER(len=25), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfname     ! coupling field short name for restart (?) and diagnostics
+!$OMP THREADPRIVATE(cfname)
+
+  CHARACTER(len=25), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfname_in  ! coupling field short name for restart (?) and diagnostics
+!$OMP THREADPRIVATE(cfname_in)
+
+  CHARACTER(len=25), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfname_out ! coupling field short name for restart (?) and diagnostics
+!$OMP THREADPRIVATE(cfname_out)
+
+  CHARACTER(len=15), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfunits_in  !  coupling field units for diagnostics
+!$OMP THREADPRIVATE(cfunits_in)
+
+  CHARACTER(len=15), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfunits_out !  coupling field units for diagnostics
+!$OMP THREADPRIVATE(cfunits_out)
+
+  CHARACTER(len=120), ALLOCATABLE, DIMENSION(:), PUBLIC :: cftext_in  ! coupling field long name for diagnostics
+!$OMP THREADPRIVATE(cftext_in)
+
+  CHARACTER(len=120), ALLOCATABLE, DIMENSION(:), PUBLIC :: cftext_out ! coupling field long name for diagnostics
+!$OMP THREADPRIVATE(cftext_out)
+
+  CHARACTER(len=5), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfmod1 ! model 1 (rreference) : LMDz
+!$OMP THREADPRIVATE(cfmod1)
+
+  CHARACTER(len=5), ALLOCATABLE, DIMENSION(:), PUBLIC :: cfmod2 ! model 2
+!$OMP THREADPRIVATE(cfmod2)
+
+  REAL, ALLOCATABLE, DIMENSION(:,:), PUBLIC :: zcfields_in
+!$OMP THREADPRIVATE(zcfields_in)
+
+  REAL, ALLOCATABLE, DIMENSION(:,:), PUBLIC :: zcfields_out
+!$OMP THREADPRIVATE(zcfields_out)
 
   TYPE, PUBLIC ::   co2_trac_type
@@ -111 +196 @@
     INTEGER               :: teststop
 
-
-
 ! 1) Read controle parameters from .def input file
 ! ------------------------------------------------
@@ -125 +208 @@
        WRITE(lunout,*) 'carbon_cycle_fos_fuel = ', fos_fuel_s 
     END IF
-
 
     ! Read parmeter for calculation compatible emission
@@ -297 +379 @@
   END SUBROUTINE carbon_cycle_init
 
-
   SUBROUTINE carbon_cycle(nstep, pdtphys, pctsrf, tr_seri, source)
 ! Subroutine for injection of co2 in the tracers
@@ -483 +564 @@
   END SUBROUTINE carbon_cycle
   
+  SUBROUTINE infocfields_init
+
+    USE control_mod, ONLY: planet_type
+    USE phys_cal_mod, ONLY : mth_cur
+    USE mod_synchro_omp
+    USE mod_phys_lmdz_para, ONLY: is_mpi_root, is_omp_root
+    USE mod_phys_lmdz_transfert_para
+    USE mod_phys_lmdz_omp_transfert
+    USE dimphy, ONLY: klon
+
+    IMPLICIT NONE
+
+!=======================================================================
+!
+!   Authors: Patricia Cadule and Laurent Fairhead  
+!   -------
+!
+!  Purpose and description:
+!  -----------------------
+!
+! Infofields
+! this routine enables to define the field exchanges in both directions between
+! the atmospheric circulation model (LMDZ) and ORCHIDEE. In the future this
+! routing might apply to other models (e.g., NEMO, INCA, ...).
+! Therefore, currently with this routine, it is possible to define the coupling
+! fields only between LMDZ and ORCHIDEE.
+! The coupling_fields.def file enables to define the name of the exchanged
+! fields at the coupling interface.
+! field_in_names : the set of names of the exchanged fields in input to ORCHIDEE
+! (LMDZ to ORCHIDEE)
+! field_out_names : the set of names of the exchanged fields in output of
+! ORCHIDEE (ORCHIDEE to LMDZ)
+! n : the number of exchanged fields at th coupling interface
+! nb_fields_in : number of inputs fields to ORCHIDEE (LMDZ to ORCHIDEE)
+! nb_fields_out : number of ouput fields of ORCHIDEE (ORCHIDEE to LMDZ)
+!
+! The syntax for coupling_fields.def is as follows:
+! IMPORTANT: each column entry must be separated from the previous one by 3
+! spaces and only that
+! field name         coupling          model 1         model 2         long_name
+!                    direction
+!   10char  -3spaces-  3char  -3spaces- 4char -3spaces- 4char -3spaces-  30char 
+!
+! n
+! FIELD1 IN LMDZ ORC
+! ....
+! FIELD(j) IN LMDZ ORC
+! FIELD(j+1) OUT LMDZ ORC
+! ...
+! FIELDn OUT LMDZ ORC   
+! 
+!=======================================================================
+!   ... 22/12/2017 ....
+!-----------------------------------------------------------------------
+! Declarations
+
+  INCLUDE "clesphys.h"
+  INCLUDE "dimensions.h"
+  INCLUDE "iniprint.h"
+
+! Local variables
+
+  INTEGER :: iq,  ierr, stat, error
+
+  CHARACTER(LEN=20), ALLOCATABLE, DIMENSION(:), SAVE  :: cfname_root
+  CHARACTER(LEN=120), ALLOCATABLE, DIMENSION(:), SAVE :: cftext_root
+  CHARACTER(LEN=15), ALLOCATABLE, DIMENSION(:), SAVE  :: cfunits_root
+
+  CHARACTER(len=3), ALLOCATABLE, DIMENSION(:) :: cfintent_root
+  CHARACTER(len=5), ALLOCATABLE, DIMENSION(:) :: cfmod1_root
+  CHARACTER(len=5), ALLOCATABLE, DIMENSION(:) :: cfmod2_root
+
+  LOGICAL, ALLOCATABLE, DIMENSION(:), SAVE :: mask_in_root
+  LOGICAL, ALLOCATABLE, DIMENSION(:), SAVE :: mask_out_root
+
+  CHARACTER(len=*),parameter :: modname="infocfields"
+
+!-----------------------------------------------------------------------
+
+nbcf=0
+nbcf_in=0
+nbcf_out=0
+
+IF (planet_type=='earth') THEN
+
+    IF (is_mpi_root .AND. is_omp_root) THEN
+
+       IF (level_coupling_esm.GT.0) THEN
+
+          OPEN(200,file='coupling_fields.def',form='formatted',status='old', iostat=ierr)
+
+          IF (ierr.EQ.0) THEN
+
+             WRITE(lunout,*) trim(modname),': Open coupling_fields.def : ok'
+             READ(200,*) nbcf
+             WRITE(lunout,*) 'infocfields_mod.F90 --- nbcf=',nbcf
+             IF (.NOT. ALLOCATED(cfname_root))   ALLOCATE(cfname_root(nbcf))
+             IF (.NOT. ALLOCATED(cfintent_root)) ALLOCATE(cfintent_root(nbcf))
+             IF (.NOT. ALLOCATED(cfmod1_root))   ALLOCATE(cfmod1_root(nbcf))
+             IF (.NOT. ALLOCATED(cfmod2_root))   ALLOCATE(cfmod2_root(nbcf))
+             IF (.NOT. ALLOCATED(cftext_root))   ALLOCATE(cftext_root(nbcf))
+             IF (.NOT. ALLOCATED(cfunits_root))  ALLOCATE(cfunits_root(nbcf))
+             IF (.NOT. ALLOCATED(mask_in_root))  ALLOCATE(mask_in_root(nbcf))
+             IF (.NOT. ALLOCATED(mask_out_root)) ALLOCATE(mask_out_root(nbcf))
+
+             nbcf_in=0
+             nbcf_out=0
+             DO iq=1,nbcf
+                WRITE(lunout,*) 'infofields : field=',iq
+                READ(200,'(A15,3X,A3,3X,A5,3X,A5,3X,A120,3X,A15)',IOSTAT=ierr) &
+                   cfname_root(iq),cfintent_root(iq),cfmod1_root(iq),cfmod2_root(iq),cftext_root(iq),cfunits_root(iq)
+                cfname_root(iq)=TRIM(cfname_root(iq))
+                cfintent_root(iq)=TRIM(cfintent_root(iq))
+                cfmod1_root(iq)=TRIM(cfmod1_root(iq))
+                cfmod2_root(iq)=TRIM(cfmod2_root(iq))
+                cftext_root(iq)=TRIM(cftext_root(iq))
+                cfunits_root(iq)=TRIM(cfunits_root(iq))
+                WRITE(lunout,*) 'coupling field: ',cfname_root(iq), &  
+                               ', number: ',iq,', INTENT: ',cfintent_root(iq)
+                WRITE(lunout,*) 'coupling field: ',cfname_root(iq), & 
+                               ', number: ',iq,', model 1 (ref): ',cfmod1_root(iq),', model 2: ',cfmod2_root(iq)
+                WRITE(lunout,*) 'coupling field: ',cfname_root(iq), & 
+                               ', number: ',iq,', long name: ',cftext_root(iq),', units ',cfunits_root(iq)
+                IF ((TRIM(cfintent_root(iq)).NE.'OUT') .AND. (nbcf_in.LE.nbcf)) THEN
+                  nbcf_in=nbcf_in+1
+                  mask_in_root(iq)=.TRUE.
+                  mask_out_root(iq)=.FALSE.
+                ELSE IF ((TRIM(cfintent_root(iq)).EQ.'OUT') .AND. (nbcf_out.LE.nbcf)) THEN
+                  nbcf_out=nbcf_out+1
+                  mask_in_root(iq)=.FALSE.
+                  mask_out_root(iq)=.TRUE.
+                ELSE
+                  WRITE(lunout,*) 'abort_gcm --- nbcf: ',nbcf
+                  WRITE(lunout,*) 'abort_gcm --- nbcf_in: ',nbcf_in
+                  WRITE(lunout,*) 'abort_gcm --- nbcf_out: ',nbcf_out
+                  CALL abort_gcm('infocfields_init','Problem in the definition of the coupling fields',1)
+               ENDIF
+             ENDDO !DO iq=1,nbcf
+          ELSE
+             WRITE(lunout,*) trim(modname),': infocfields_mod.F90 --- Problem in opening coupling_fields.def'
+             WRITE(lunout,*) trim(modname),': infocfields_mod.F90 --- WARNING using defaut values'
+          ENDIF ! ierr
+          CLOSE(200)
+
+       ENDIF ! level_coupling_esm
+
+       IF (nbcf_out .EQ. 0) nbcf_out=-1
+       IF (nbcf_in .EQ. 0) nbcf_in=-1
+
+    ENDIF !   (is_mpi_root .AND. is_omp_root)
+!$OMP BARRIER
+
+    CALL bcast(nbcf)
+    CALL bcast(nbcf_in)
+    CALL bcast(nbcf_out)
+
+    WRITE(lunout,*) 'infocfields_mod.F90 --- nbcf=',nbcf
+    WRITE(lunout,*) 'infocfields_mod.F90 --- nbcf_in=',nbcf_in
+    WRITE(lunout,*) 'infocfields_mod.F90 --- nbcf_out=',nbcf_out
+
+    IF ((nbcf.GT.0) .AND. (.NOT. ALLOCATED(cfname))) ALLOCATE(cfname(nbcf))
+    IF ((nbcf_in.GT.0) .AND. (.NOT. ALLOCATED(cfname_in))) ALLOCATE(cfname_in(nbcf_in))
+    IF ((nbcf_in.GT.0) .AND. (.NOT. ALLOCATED(cftext_in))) ALLOCATE(cftext_in(nbcf_in))
+    IF ((nbcf_out.GT.0) .AND. (.NOT. ALLOCATED(cfname_out))) ALLOCATE(cfname_out(nbcf_out))
+    IF ((nbcf_out.GT.0) .AND. (.NOT. ALLOCATED(cftext_out))) ALLOCATE(cftext_out(nbcf_out))
+    IF ((nbcf.GT.0) .AND. (.NOT. ALLOCATED(cfmod1))) ALLOCATE(cfmod1(nbcf))
+    IF ((nbcf.GT.0) .AND. (.NOT. ALLOCATED(cfmod2))) ALLOCATE(cfmod2(nbcf))
+    IF ((nbcf_in.GT.0) .AND. (.NOT. ALLOCATED(cfunits_in))) ALLOCATE(cfunits_in(nbcf_in))
+    IF ((nbcf_out.GT.0) .AND. (.NOT. ALLOCATED(cfunits_out))) ALLOCATE(cfunits_out(nbcf_out))
+       
+    IF (is_mpi_root .AND. is_omp_root) THEN
+
+        IF (nbcf.GT.0) cfname=cfname_root
+        IF (nbcf_in.GT.0) cfname_in=PACK(cfname_root,mask_in_root)
+        IF (nbcf_out.GT.0) cfname_out=PACK(cfname_root,mask_out_root)
+        IF (nbcf_in.GT.0) cftext_in=PACK(cftext_root,mask_in_root)
+        IF (nbcf_out.GT.0) cftext_out=PACK(cftext_root,mask_out_root)
+        IF (nbcf.GT.0) cfmod1=cfmod1_root
+        IF (nbcf.GT.0) cfmod2=cfmod2_root
+        IF (nbcf_in.GT.0) cfunits_in=PACK(cfunits_root,mask_in_root)
+        IF (nbcf_out.GT.0) cfunits_out=PACK(cfunits_root,mask_out_root)
+
+        nbcf_in_orc=0
+        nbcf_in_nemo=0
+        nbcf_in_inca=0
+        nbcf_in_ant=0
+
+        DO iq=1,nbcf
+            IF (cfmod1(iq) == "ORC") nbcf_in_orc=nbcf_in_orc+1  
+            IF (cfmod1(iq) == "NEMO") nbcf_in_nemo=nbcf_in_nemo+1  
+            IF (cfmod1(iq) == "INCA") nbcf_in_inca=nbcf_in_inca+1
+            IF (cfmod1(iq) == "ALL") nbcf_in_orc=nbcf_in_orc+1    ! ALL = ORC/NEMO/INCA
+            IF (cfmod1(iq) == "ALL") nbcf_in_nemo=nbcf_in_nemo+1  ! ALL = ORC/NEMO/INCA 
+            IF (cfmod1(iq) == "ALL") nbcf_in_inca=nbcf_in_inca+1  ! ALL = ORC/NEMO/INCA 
+            IF (cfmod1(iq) == "ANT") nbcf_in_ant=nbcf_in_ant+1  
+        ENDDO
+
+    ENDIF !   (is_mpi_root .AND. is_omp_root)
+!$OMP BARRIER
+
+    CALL bcast(nbcf_in_orc)
+    CALL bcast(nbcf_in_nemo)
+    CALL bcast(nbcf_in_inca)
+    CALL bcast(nbcf_in_ant)
+
+    WRITE(lunout,*) 'nbcf_in_orc ',nbcf_in_orc
+    WRITE(lunout,*) 'nbcf_in_nemo ',nbcf_in_nemo
+    WRITE(lunout,*) 'nbcf_in_inca ',nbcf_in_inca
+    WRITE(lunout,*) 'nbcf_in_ant ',nbcf_in_ant
+
+    IF (nbcf_in.GT.0) THEN 
+        DO iq=1,nbcf_in
+          CALL bcast(cfname_in(iq))
+          CALL bcast(cftext_in(iq))
+          CALL bcast(cfunits_in(iq))
+        ENDDO
+    ENDIF
+
+    IF (nbcf_out.GT.0) THEN
+        DO iq=1,nbcf_out
+          CALL bcast(cfname_out(iq))
+          CALL bcast(cftext_out(iq))
+          CALL bcast(cfunits_out(iq))
+        ENDDO
+    ENDIF
+
+    IF (nbcf.GT.0) THEN
+        DO iq=1,nbcf
+          CALL bcast(cfmod1(iq))
+          CALL bcast(cfmod2(iq))
+        ENDDO
+    ENDIF
+
+    IF (nbcf_in.GT.0) WRITE(lunout,*)'infocfields_mod --- cfname_in: ',cfname_in
+    IF (nbcf_out.GT.0) WRITE(lunout,*)'infocfields_mod --- cfname_out: ',cfname_out
+
+    IF (nbcf_in.GT.0) WRITE(lunout,*)'infocfields_mod --- cftext_in: ',cftext_in
+    IF (nbcf_out.GT.0) WRITE(lunout,*)'infocfields_mod --- cftext_out: ',cftext_out
+
+    IF (nbcf.GT.0) WRITE(lunout,*)'infocfields_mod --- cfmod1: ',cfmod1 
+    IF (nbcf.GT.0) WRITE(lunout,*)'infocfields_mod --- cfmod2: ',cfmod2 
+
+    IF (nbcf_in.GT.0) WRITE(lunout,*)'infocfunits_mod --- cfunits_in: ',cfunits_in
+    IF (nbcf_out.GT.0) WRITE(lunout,*)'infocfunits_mod --- cfunits_out: ',cfunits_out
+
+    IF (nbcf_in.GT.0) WRITE(*,*)'infocfields_init --- number of fields in to LMDZ: ',nbcf_in
+    IF (nbcf_out.GT.0) WRITE(*,*)'infocfields_init --- number of fields out of LMDZ: ',nbcf_out
+
+ ELSE
+ ! Default values for other planets
+    nbcf=0 
+    nbcf_in=0
+    nbcf_out=0
+ ENDIF ! planet_type
+
+ IF ((nbcf_in.GT.0) .AND. (.NOT. ALLOCATED(zcfields_in))) ALLOCATE(zcfields_in(klon,nbcf_in),stat=error)
+        IF (error /= 0)  CALL abort_gcm(modname,'Pb in allocation zcfields_in',1)
+ IF ((nbcf_out.GT.0) .AND. (.NOT. ALLOCATED(zcfields_out))) ALLOCATE(zcfields_out(klon,nbcf_out),stat=error)
+        IF (error /= 0)  CALL abort_gcm(modname,'Pb in allocation zcfields_out',1)
+
+END SUBROUTINE infocfields_init
+
 END MODULE carbon_cycle_mod

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -16 +16 @@
        d_u, d_v, d_t, d_qx, d_ps)
 
-    use assert_m, only: assert
+    USE assert_m, only: assert
     USE ioipsl, only: histbeg, histvert, histdef, histend, histsync, &
          histwrite, ju2ymds, ymds2ju, getin
@@ -168 +168 @@
        omega,  &
        epmax_diag,  &
-
        !    Deep convective variables used in phytrac
        pmflxr, pmflxs,  &
@@ -184 +183 @@
        epmlmMm, eplaMm, &
        sij, &
-
+       !
        cldemi,  &
        cldfra, cldtau, fiwc,  &
@@ -199 +198 @@
        prfl, psfl, fraca, Vprecip,  &
        zw2,  &
-       
+       !
        fluxu, fluxv,  &
        fluxt,  &
-
+       !
        uwriteSTD, vwriteSTD, &                !pour calcul_STDlev.h
        wwriteSTD, phiwriteSTD, &              !pour calcul_STDlev.h
        qwriteSTD, twriteSTD, rhwriteSTD, &    !pour calcul_STDlev.h
-       
+       ! 
        beta_prec,  &
        rneb,  &
@@ -213 +212 @@
     USE phys_state_var_mod ! Variables sauvegardees de la physique
 #ifdef CPP_Dust
-  USE phys_output_write_spl_mod
+    USE phys_output_write_spl_mod
 #else
     USE phys_output_var_mod ! Variables pour les ecritures des sorties
@@ -222 +221 @@
     USE phys_output_mod
     USE phys_output_ctrlout_mod
-    use open_climoz_m, only: open_climoz ! ozone climatology from a file
-    use regr_pr_time_av_m, only: regr_pr_time_av
-    use netcdf95, only: nf95_close
+    USE open_climoz_m, only: open_climoz ! ozone climatology from a file
+    USE regr_pr_time_av_m, only: regr_pr_time_av
+    USE netcdf95, only: nf95_close
     !IM for NMC files
-    !     use netcdf, only: nf90_fill_real
-    use netcdf, only: nf90_fill_real
-    use mod_phys_lmdz_mpi_data, only: is_mpi_root
+    USE netcdf, only: nf90_fill_real
+    USE mod_phys_lmdz_mpi_data, only: is_mpi_root
     USE aero_mod
-    use ozonecm_m, only: ozonecm ! ozone of J.-F. Royer
-    use conf_phys_m, only: conf_phys
-    use radlwsw_m, only: radlwsw
-    use phyaqua_mod, only: zenang_an
+    USE ozonecm_m, only: ozonecm ! ozone of J.-F. Royer
+    USE conf_phys_m, only: conf_phys
+    USE radlwsw_m, only: radlwsw
+    USE phyaqua_mod, only: zenang_an
     USE time_phylmdz_mod, only: day_step_phy, annee_ref, day_ref, itau_phy, &
          start_time, pdtphys, day_ini
@@ -246 +244 @@
     USE indice_sol_mod
     USE phytrac_mod, ONLY : phytrac
+    USE carbon_cycle_mod, ONLY : infocfields_init
 
 #ifdef CPP_RRTM
@@ -263 +262 @@
     !IM stations CFMIP
     USE CFMIP_point_locations
-    use FLOTT_GWD_rando_m, only: FLOTT_GWD_rando
-    use ACAMA_GWD_rando_m, only: ACAMA_GWD_rando
+    USE FLOTT_GWD_rando_m, only: FLOTT_GWD_rando
+    USE ACAMA_GWD_rando_m, only: ACAMA_GWD_rando
     USE VERTICAL_LAYERS_MOD, ONLY: aps,bps
 
-
-    IMPLICIT none
+    IMPLICIT NONE
     !>======================================================================
     !!
@@ -1305 +1303 @@
        nvm_lmdz = 13
        CALL getin_p('NVM',nvm_lmdz)
+
+       !--PC: defining fields to be exchanged between LMDz, ORCHIDEE and NEMO
+       WRITE(lunout,*) 'Call to infocfields from physiq'
+       CALL infocfields_init
+
     ENDIF
 
     IF (prt_level.ge.1) print *,'CONVERGENCE PHYSIQUE THERM 1 '
-
 
     !======================================================================