Juan 30/05/2024:modd_cst.f90/compute_frac_ice.func.h/mode_tiwmx.f90 , For GPU...

Juan 30/05/2024:modd_cst.f90/compute_frac_ice.func.h/mode_tiwmx.f90 , For GPU routine to work , add PURE & "!$acc routine seq" & "!$acc declare VAR"

src/MNH/modd_cst.f90

@@ -141,6 +141,7 @@ REAL, POINTER :: XCPD=>NULL(), XCPV=>NULL()
 REAL, POINTER :: XRHOLW=>NULL()
 REAL, POINTER :: XCL=>NULL(), XCI=>NULL()
 REAL, POINTER :: XTT=>NULL()
+!$acc declare create(XTT)
 REAL, POINTER :: XLVTT=>NULL()
 REAL, POINTER :: XLSTT=>NULL()
 REAL, POINTER :: XLMTT=>NULL()

src/PHYEX/micro/compute_frac_ice.func.h

@@ -3,7 +3,7 @@
 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
       ELEMENTAL SUBROUTINE COMPUTE_FRAC_ICE(HFRAC_ICE,NEB,PFRAC_ICE,PT,KERR)
-
+      !$acc routine seq
 ! ******* TO BE INCLUDED IN THE *CONTAINS* OF A SUBROUTINE, IN ORDER TO EASE AUTOMATIC INLINING ******
 ! => Don't use drHook !!!
 !

src/PHYEX/micro/mode_tiwmx.f90

@@ -37,6 +37,7 @@ REAL ::  ESTABW(NSTART:NSTOP)
 REAL :: DESTABW(NSTART:NSTOP)
 REAL ::  ESTABI(NSTART:NSTOP)
 REAL :: DESTABI(NSTART:NSTOP)
+!$acc declare create(ESTABW,ESTABI)
 
 ! Ice crystal- or water droplet tables
 REAL ::   A2TAB(NSTART:NSTOP)
@@ -50,6 +51,7 @@ REAL :: REDINTAB(NSTART:NSTOP)
 CONTAINS
 
   REAL FUNCTION ESATW(TT)
+    !$acc routine seq
     REAL,INTENT(IN) :: TT 
     ESATW = ESTABW(NINT(XNDEGR*TT))
   END FUNCTION ESATW
@@ -60,6 +62,7 @@ CONTAINS
   END FUNCTION
 
   REAL FUNCTION ESATI(TT)
+    !$acc routine seq
     REAL,INTENT(IN) :: TT 
     ESATI = ESTABI(NINT(XNDEGR*TT))
   END FUNCTION