! Ecriture des champs de forcage radiatif

      ndex2d = 0
      ndex3d = 0

c Champs 2D:

	print*,'itap',m
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1, albsfi,zx_tmp_2d)
      CALL histwrite(nid_day,"albs",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
!
!     Profil de reference
c toa
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1, toa_sw_ref,zx_tmp_2d)
      CALL histwrite(nid_day,"toa_sw_ref",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
C
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,toa_sw_ref0,zx_tmp_2d)
      CALL histwrite(nid_day,"toa_sw_ref0",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,toa_lw_ref,zx_tmp_2d)
      CALL histwrite(nid_day,"toa_lw_ref",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,toa_lw_ref0,zx_tmp_2d)
      CALL histwrite(nid_day,"toa_lw_ref0",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c tps
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1, tps_sw_ref,zx_tmp_2d)
      CALL histwrite(nid_day,"tps_sw_ref",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
C
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,tps_sw_ref0,zx_tmp_2d)
      CALL histwrite(nid_day,"tps_sw_ref0",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,tps_lw_ref,zx_tmp_2d)
      CALL histwrite(nid_day,"tps_lw_ref",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,tps_lw_ref0,zx_tmp_2d)
      CALL histwrite(nid_day,"tps_lw_ref0",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c srf
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1, srf_sw_ref,zx_tmp_2d)
      CALL histwrite(nid_day,"srf_sw_ref",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
C
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,srf_sw_ref0,zx_tmp_2d)
      CALL histwrite(nid_day,"srf_sw_ref0",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,srf_lw_ref,zx_tmp_2d)
      CALL histwrite(nid_day,"srf_lw_ref",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,srf_lw_ref0,zx_tmp_2d)
      CALL histwrite(nid_day,"srf_lw_ref0",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
!
! Flux initiaux après perturbation, sans ajustement strato
!
c toa
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,toa_sw_ini,zx_tmp_2d)
      CALL histwrite(nid_day,"toa_sw_ini",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,toa_sw_ini0,zx_tmp_2d)
      CALL histwrite(nid_day,"toa_sw_ini0",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,toa_lw_ini,zx_tmp_2d)
      CALL histwrite(nid_day,"toa_lw_ini",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,toa_lw_ini0,zx_tmp_2d)
      CALL histwrite(nid_day,"toa_lw_ini0",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c tps
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1, tps_sw_ini,zx_tmp_2d)
      CALL histwrite(nid_day,"tps_sw_ini",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,tps_sw_ini0,zx_tmp_2d)
      CALL histwrite(nid_day,"tps_sw_ini0",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1, tps_lw_ini,zx_tmp_2d)
      CALL histwrite(nid_day,"tps_lw_ini",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,tps_lw_ini0,zx_tmp_2d)
      CALL histwrite(nid_day,"tps_lw_ini0",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
!
! Forcages initiaux après perturbation, sans ajustement strato
!
c toa
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,d_toa_sw_ini,zx_tmp_2d)
      CALL histwrite(nid_day,"d_toa_sw_ini",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,d_toa_sw_ini0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_toa_sw_ini0",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,d_toa_lw_ini,zx_tmp_2d)
      CALL histwrite(nid_day,"d_toa_lw_ini",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,d_toa_lw_ini0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_toa_lw_ini0",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c tps
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,d_tps_sw_ini,zx_tmp_2d)
      CALL histwrite(nid_day,"d_tps_sw_ini",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,d_tps_sw_ini0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_tps_sw_ini0",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,d_tps_lw_ini,zx_tmp_2d)
      CALL histwrite(nid_day,"d_tps_lw_ini",m,zx_tmp_2d,
     s                iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1, ngridmx,iim,jjm+1,d_tps_lw_ini0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_tps_lw_ini0",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
!
! Forcages après ajustement strato
!
c toa
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_toa_sw_adj,zx_tmp_2d)
      CALL histwrite(nid_day,"d_toa_sw_adj",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_toa_sw_adj0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_toa_sw_adj0",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_toa_lw_adj,zx_tmp_2d)
      CALL histwrite(nid_day,"d_toa_lw_adj",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_toa_lw_adj0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_toa_lw_adj0",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c tps
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_tps_sw_adj,zx_tmp_2d)
      CALL histwrite(nid_day,"d_tps_sw_adj",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_tps_sw_adj0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_tps_sw_adj0",m,zx_tmp_2d,
     .              iim*(jjm+1), ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_tps_lw_adj,zx_tmp_2d)
      CALL histwrite(nid_day,"d_tps_lw_adj",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_tps_lw_adj0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_tps_lw_adj0",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c srf
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_srf_sw_adj,zx_tmp_2d)
      CALL histwrite(nid_day,"d_srf_sw_adj",m,zx_tmp_2d,
     s               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_srf_sw_adj0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_srf_sw_adj0",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_srf_lw_adj,zx_tmp_2d)
      CALL histwrite(nid_day,"d_srf_lw_adj",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,d_srf_lw_adj0,zx_tmp_2d)
      CALL histwrite(nid_day,"d_srf_lw_adj0",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
!
! Bilan radiatif de la stratosphere
!
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,bilq_ref,zx_tmp_2d)
      CALL histwrite(nid_day,"bilq_ref",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,bilq_ini,zx_tmp_2d)
      CALL histwrite(nid_day,"bilq_ini",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
      CALL gr_fi_ecrit(1,ngridmx,iim,jjm+1,bilq_adj,zx_tmp_2d)
      CALL histwrite(nid_day,"bilq_adj",m,zx_tmp_2d,
     .               iim*(jjm+1),ndex2d)
c
c Champs 3D:
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,tfi,zx_tmp_3d)
      CALL histwrite(nid_day,"temp",m,zx_tmp_3d,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,tfi_adj,zx_tmp_3d)
      CALL histwrite(nid_day,"temp_adj",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,qfi,zx_tmp_3d)
      CALL histwrite(nid_day,"ovap",m,zx_tmp_3d,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,dHrad_dT,zx_tmp_3d)
      CALL histwrite(nid_day,"dHrad_dT",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
!     Heating and cooling rates
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,heat_ref,zx_tmp_3d)
      CALL histwrite(nid_day,"heat_ref",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,heat0_ref,zx_tmp_3d)
      CALL histwrite(nid_day,"heat0_ref",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,cool_ref,zx_tmp_3d)
      CALL histwrite(nid_day,"cool_ref",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,cool0_ref,zx_tmp_3d)
      CALL histwrite(nid_day,"cool0_ref",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,d_heat_ini,zx_tmp_3d)
      CALL histwrite(nid_day,"d_heat_ini",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,d_heat0_ini,zx_tmp_3d)
      CALL histwrite(nid_day,"d_heat0_ini",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,d_cool_ini,zx_tmp_3d)
      CALL histwrite(nid_day,"d_cool_ini",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,d_cool0_ini,zx_tmp_3d)
      CALL histwrite(nid_day,"d_cool0_ini",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,d_heat_adj,zx_tmp_3d)
      CALL histwrite(nid_day,"d_heat_adj",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,d_heat0_adj,zx_tmp_3d)
      CALL histwrite(nid_day,"d_heat0_adj",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,d_cool_adj,zx_tmp_3d)
      CALL histwrite(nid_day,"d_cool_adj",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      CALL gr_fi_ecrit(llm,ngridmx,iim,jjm+1,d_cool0_adj,zx_tmp_3d)
      CALL histwrite(nid_day,"d_cool0_adj",m,zx_tmp_3d
     $    ,iim*(jjm+1)*llm,ndex3d)
c
      call histsync(nid_day)
	print*,'ok histwrite'