Changeset 4185 for LMDZ6/branches

Timestamp:

Jun 23, 2022, 6:36:15 PM (2 years ago)

Author:

Laurent Fairhead

Message:

Ported to GPU using the replay_mod method: the phys.nc file that is output show now difference
between a run using GPU and a run without GPU

Location:

LMDZ6/branches/Portage_acc/libf/phylmd

Files:

• thermcell_env.F90 (modified) (5 diffs)
• thermcell_ini_mod.F90 (modified) (1 diff)

LMDZ6/branches/Portage_acc/libf/phylmd/thermcell_env.F90

@@ -25 +25 @@
    logical mask(ngrid,nlay)
 
+   !$acc data create (mask) &
+   !$acc &    copyin (po, pt, pu, pv, pplay, pplev) &
+   !$acc &    copyout (zo, zl, zh, ztv, zthl)     &
+   !$acc &    copyout (zpspsk, zu, zv, pqsat)       &
+   !$acc &
+
 
 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
 ! Initialisations :
 !------------------
+!!! temporary for GPU port   write(81) ngrid,nlay,lev_out,po,pt,pu,pv,pplay,pplev
 
+
+  !$acc kernels default(none) async
    mask(:,:)=.true.
-
+   !$acc end kernels
+   
 !
 ! calcul des caracteristiques de l environnement
+  !$acc kernels default(none) async
    DO  ll=1,nlay
      DO ig=1,ngrid
@@ -41 +52 @@
      enddo
    enddo
-
+  !$acc end kernels
+   
 ! Condensation :
 !---------------
@@ -47 +59 @@
 
    call thermcell_qsat(ngrid*nlay,mask,pplev,pt,po,pqsat)
+  !$acc kernels default(none) async
    do ll=1,nlay
       do ig=1,ngrid
@@ -54 +67 @@
       enddo
    enddo
-
+  !$acc end kernels
 !-----------------------------------------------------------------------
    if (prt_level.ge.1) print*,'0 OK convect8'
 
+  !$acc kernels default(none) async
    do ll=1,nlay
       do ig=1,ngrid
@@ -74 +88 @@
       enddo
    enddo
+  !$acc end kernels
  
- RETURN
+!!! temporary for GPU port !!!$acc update self(zo, zl, zh, ztv, zthl,zpspsk, zu, zv, pqsat)
+!!! temporary for GPU port include "dump_param.h" ! replay automatic include 
+  !$acc end data
+   RETURN
    END

LMDZ6/branches/Portage_acc/libf/phylmd/thermcell_ini_mod.F90

@@ -40 +40 @@
 real, intent(in) :: RG_in,RD_in,RCPD_in,RKAPPA_in,RLVTT_in,RETV_in,tau_thermals_in
 
-print*,'thermcell_ini'
-      if (iflag_thermals==15.or.iflag_thermals==16) then
-         dvdq=0
-         dqimpl=-1
-      else
-         dvdq=1
-         dqimpl=1
-      endif
+   if (iflag_thermals==15.or.iflag_thermals==16) then
+      dvdq=0
+      dqimpl=-1
+   else
+      dvdq=1
+      dqimpl=1
+   endif
    prt_level=prt_level_in
    RG=RG_in