Skip to content
Snippets Groups Projects
Commit 722c7f98 authored by Juan Escobar's avatar Juan Escobar
Browse files

Juan 27/09/2022:prandtl.f90, Bug, missing (:,:,:) in array explansion for no reprod side

parent a08cb598
No related branches found
No related tags found
No related merge requests found
......@@ -985,7 +985,7 @@ ELSE IF (L2D) THEN ! 3D case in a 2D model
!$acc kernels
!$mnh_expand_array(JI=1:JIU,JJ=1:JJU,JK=1:JKU)
#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
ZTMP1_DEVICE(:,:,:) = (XG / PTHVREF * PLM * PLEPS / PTKEM)**2
ZTMP1_DEVICE(:,:,:) = (XG / PTHVREF(:,:,:) * PLM(:,:,:) * PLEPS(:,:,:) / PTKEM(:,:,:))**2
#else
ZTMP1_DEVICE(:,:,:) = BR_P2(XG / PTHVREF(:,:,:) * PLM(:,:,:) * PLEPS(:,:,:) / PTKEM(:,:,:))
#endif
......@@ -1066,7 +1066,7 @@ ELSE ! 3D case in a 3D model
!$acc kernels
!$mnh_expand_array(JI=1:JIU,JJ=1:JJU,JK=1:JKU)
#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP)
ZTMP1_DEVICE(:,:,:) = (XG / PTHVREF * PLM * PLEPS / PTKEM)**2
ZTMP1_DEVICE(:,:,:) = (XG / PTHVREF(:,:,:) * PLM(:,:,:) * PLEPS(:,:,:) / PTKEM(:,:,:))**2
#else
ZTMP1_DEVICE(:,:,:) = BR_P2(XG / PTHVREF(:,:,:) * PLM(:,:,:) * PLEPS(:,:,:) / PTKEM(:,:,:))
#endif
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment