Juan 24/10/2023:ZSOLVER/dyn_sources.f90 , for cce/16.0.0 , add new present_cr...

Juan 24/10/2023:ZSOLVER/dyn_sources.f90 , for cce/16.0.0 , add new present_cr directive to avoid FALSE RECURRENCE on GPU

src/ZSOLVER/dyn_sources.f90

@@ -329,7 +329,7 @@ CALL MYF_DEVICE(ZRVTT,ZTMP2_DEVICE ) !MYF(ZRVT*PVT)
 CALL MXF_DEVICE(ZRUTT,ZTMP3_DEVICE ) !MXF(ZRUT*PUT)
 CALL MXF_DEVICE(PUT,ZTMP4_DEVICE )   !MXF(PUT)
 
-!$acc kernels present_cr(ZTMP1_DEVICE,ZTMP1_DEVICE,ZTMP1_DEVICE,ZTMP1_DEVICE)
+!$acc kernels present_cr(ZTMP1_DEVICE,ZTMP2_DEVICE,ZTMP3_DEVICE,ZTMP4_DEVICE)
 ZTMP1_DEVICE(:,:,:)=ZTMP1_DEVICE(:,:,:)* ZWORK1(:,:,:)   !MYF(PVT) * ZWORK1
 ZTMP2_DEVICE(:,:,:)=ZTMP2_DEVICE(:,:,:)* ZWORK2(:,:,:)   !MYF(ZRVT*PVT) * ZWORK2
 ZTMP3_DEVICE(:,:,:)=ZTMP3_DEVICE(:,:,:)* ZWORK1(:,:,:)   !MXF(ZRUT*PUT) * ZWORK1
@@ -564,7 +564,7 @@ IF( .NOT.L1D ) THEN
 !
 ! stores the specific heat capacity (Cph) in ZWORK1
 
-!$acc kernels present_cr( ZWORK1)
+!$acc kernels present_cr(ZWORK1)
     ZWORK1(:,:,:) = XCPD + XCPV * PRT(:,:,:,1)  ! gas mixing
 !$acc loop seq
     DO JWATER = 2,1+KRRL             !  loop on the liquid components