source: CONFIG/publications/ICOLMDZORINCA_CO2_Transport_GMD_2023/INCA/src/INCA_SRC/chem_write_restart.F90 @ 6610

!Subroutine permettant d'écrire dans le fichier netcdf chemrestart.nc la valeur de la variable source à la fin du premier jour

#include <inca_define.h>

SUBROUTINE chem_write_restart

  USE XIOS_INCA
  USE SPECIES_NAMES
  USE INCA_DIM
  USE SFLX, ONLY : flux_source,eflux
  USE AIRPLANE_SRC, ONLY : itrop
  USE CHEM_MODS, ONLY : invariants, hrates_cv, h2oc
  USE RATE_INDEX_MOD
  USE AEROSOL_DIAG, ONLY : tau_allaer, piz_allaer, cg_allaer, mrfname, naero_grp, ccm_aer
  IMPLICIT NONE

  INTEGER :: i

  CALL xios_inca_change_context("inca")
  
  DO i=1,PCNST
     CALL xios_inca_send_field("source_"//tracnam(i),flux_source(:,i))
  ENDDO
  

  DO i=1, NFS
     CALL xios_inca_send_field("inv_"//trim(invname(i)), invariants(:,:,i))
  ENDDO

#ifdef AER
  
  DO i=1, naero_grp
     CALL xios_inca_send_field("tau_inca1_"//trim(mrfname(i)), tau_allaer(:,:,i,1) )
     CALL xios_inca_send_field("tau_inca2_"//trim(mrfname(i)), tau_allaer(:,:,i,2) )
                                         
     CALL xios_inca_send_field("piz_inca1_"//trim(mrfname(i)), piz_allaer(:,:,i,1) )
     CALL xios_inca_send_field("piz_inca2_"//trim(mrfname(i)), piz_allaer(:,:,i,2) )

     CALL xios_inca_send_field("cg_inca1_"//trim(mrfname(i)), cg_allaer(:,:,i,1) )
     CALL xios_inca_send_field("cg_inca2_"//trim(mrfname(i)), cg_allaer(:,:,i,2) )
  ENDDO

  CALL xios_inca_send_field("ccm1", ccm_aer(:,:,1) )
  CALL xios_inca_send_field("ccm2", ccm_aer(:,:,2) ) 
#endif

#ifdef STRAT
  call xios_inca_send_field("h2oc", h2oc) 
#endif

  CALL xios_inca_change_context("LMDZ")

END SUBROUTINE chem_write_restart