Index: LMDZ6/branches/Portage_acc/libf/phylmd/physiq_mod.F90
===================================================================
--- LMDZ6/branches/Portage_acc/libf/phylmd/physiq_mod.F90	(revision 4441)
+++ LMDZ6/branches/Portage_acc/libf/phylmd/physiq_mod.F90	(revision 4442)
@@ -1271,7 +1271,8 @@
     !$acc & copyout(d_u,d_v,d_t,d_qx,d_ps) &
     !$acc &    create(t_w, q_w, limbas, d_t_adjwk, d_q_adjwk) &
-    !$acc &    create(t_x, q_x) 
+    !$acc &    create(d_deltaq_ajs_cv, d_deltat_ajs_cv) &
+    !$acc &    create(t_x, q_x) &
 !    !$acc &    present(wake_s, wake_deltat, wake_deltaq)  &
-!    !$acc &    present(t_seri,q_seri)
+    !$acc &    present(t_seri,q_seri, ql_seri, v_seri, u_seri, qs_seri, rneb_seri)
 
     
@@ -1803,6 +1804,8 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
        CALL wake_ini(rg,rd,rv,prt_level)
+   CALL iophys_ini(pdtphys) ! replay automatic include
+   open(81,file='dump_param.bin',form='unformatted')  ! replay automatic include
        CALL thermcell_ini(iflag_thermals,prt_level,tau_thermals,lunout, &
-   &    RG,RD,RCPD,RKAPPA,RLVTT,RETV)
+   &    RG,RD,RCPD,RKAPPA,RLVTT,RETV, RTT)
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -3010,5 +3013,5 @@
              !$acc update self(d_t_adjwk, d_q_adjwk)
              !
-                 !$acc kernels default(none) async
+             !$acc kernels default(none) async
                 DO k=1,klev
                    DO i=1,klon
Index: LMDZ6/branches/Portage_acc/libf/phylmd/thermcell_ini_mod.F90
===================================================================
--- LMDZ6/branches/Portage_acc/libf/phylmd/thermcell_ini_mod.F90	(revision 4441)
+++ LMDZ6/branches/Portage_acc/libf/phylmd/thermcell_ini_mod.F90	(revision 4442)
@@ -4,5 +4,5 @@
 save
    integer :: dvdq=1,dqimpl=-1,prt_level=0,lunout
-   real RG,RD,RCPD,RKAPPA,RLVTT,RLvCp,RETV
+   real RG,RD,RCPD,RKAPPA,RLVTT,RLvCp,RETV, RTT
    real           :: r_aspect_thermals,tau_thermals,fact_thermals_ed_dz
    integer        :: iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure
@@ -10,5 +10,5 @@
 
 !$OMP THREADPRIVATE(dvdq,dqimpl,prt_level,lunout)
-!$OMP THREADPRIVATE(RG,RD,RCPD,RKAPPA,RLVTT,RLvCp)
+!$OMP THREADPRIVATE(RG,RD,RCPD,RKAPPA,RLVTT,RLvCp, RTT)
 !$OMP THREADPRIVATE(r_aspect_thermals,tau_thermals,fact_thermals_ed_dz)
 !$OMP THREADPRIVATE(iflag_thermals_ed,iflag_thermals_optflux,iflag_thermals_closure)
@@ -33,10 +33,10 @@
 
 SUBROUTINE thermcell_ini(iflag_thermals,prt_level_in,tau_thermals_in,lunout_in, &
-   &    RG_in,RD_in,RCPD_in,RKAPPA_in,RLVTT_in,RETV_in)
+   &    RG_in,RD_in,RCPD_in,RKAPPA_in,RLVTT_in,RETV_in, RTT_in)
 
    USE ioipsl_getin_p_mod, ONLY : getin_p
 
 integer, intent(in) :: iflag_thermals,prt_level_in,lunout_in
-real, intent(in) :: RG_in,RD_in,RCPD_in,RKAPPA_in,RLVTT_in,RETV_in,tau_thermals_in
+real, intent(in) :: RG_in,RD_in,RCPD_in,RKAPPA_in,RLVTT_in,RETV_in,tau_thermals_in, RTT_in
 
    if (iflag_thermals==15.or.iflag_thermals==16) then
@@ -55,4 +55,5 @@
    RLvCp = RLVTT/RCPD
    RETV=RETV_in
+   RTT=RTT_in
    tau_thermals=tau_thermals_in
    lunout=lunout_in
Index: LMDZ6/branches/Portage_acc/libf/phylmd/thermcell_main.F90
===================================================================
--- LMDZ6/branches/Portage_acc/libf/phylmd/thermcell_main.F90	(revision 4441)
+++ LMDZ6/branches/Portage_acc/libf/phylmd/thermcell_main.F90	(revision 4442)
@@ -130,4 +130,5 @@
 #endif
 
+      
 !
 
@@ -165,7 +166,15 @@
 !   --------------------------------------------------------------------
 !
+
+   !$acc data &
+   !$acc &   copyin (po, pt, pu, pv, pplay, pplev) &
+   !$acc &   copyout(zo, zl, zh, ztv, zthl)        &
+   !$acc &   copyout (zpspsk, zu, zv, zqsat)       &
+   !$acc &
+       
       CALL thermcell_env(ngrid,nlay,po,pt,pu,pv,pplay,  &
      &           pplev,zo,zh,zl,ztv,zthl,zu,zv,zpspsk,zqsat,lev_out)
        
+   !$acc end data      
       if (prt_level.ge.1) print*,'thermcell_main apres thermcell_env'
 
@@ -687,4 +696,5 @@
 
       if (prt_level.ge.1) print*,'thermcell_main FIN  OK'
+      
 
  RETURN
