Skip to content
Snippets Groups Projects
Commit 5d7780a9 authored by RODIER Quentin's avatar RODIER Quentin
Browse files

Quentin 16/08/2022: bugfix bit-repro packing turb_ver_dyn_flux

parent 888772b2
No related branches found
No related tags found
No related merge requests found
...@@ -1060,7 +1060,7 @@ IF(HTURBDIM=='3DIM') THEN ...@@ -1060,7 +1060,7 @@ IF(HTURBDIM=='3DIM') THEN
ZWORK31D(IIJB:IIJE ) = - ZFLXZ(IIJB:IIJE,IKB+D%NKL) & ZWORK31D(IIJB:IIJE ) = - ZFLXZ(IIJB:IIJE,IKB+D%NKL) &
* ( ZWORK1(IIJB:IIJE,IKB+D%NKL) - ZWORK51D(IIJB:IIJE ) & * ( ZWORK1(IIJB:IIJE,IKB+D%NKL) - ZWORK51D(IIJB:IIJE ) &
* PDZY(IIJB:IIJE,IKB+D%NKL) ) & * PDZY(IIJB:IIJE,IKB+D%NKL) ) &
/ (0.5*(PDYY(IIJB:IIJE,IKB+D%NKL)+PDYY(IIJB:IIJE,IKE))) / (0.5*(PDYY(IIJB:IIJE,IKB+D%NKL)+PDYY(IIJB:IIJE,IKB)))
!$mnh_end_expand_array(JIJ=IIJB:IIJE) !$mnh_end_expand_array(JIJ=IIJB:IIJE)
CALL MYF2D_PHY(D,ZWORK31D,ZWORK41D) CALL MYF2D_PHY(D,ZWORK31D,ZWORK41D)
ZA(IIJB:IIJE,IKB) = ZWORK41D(IIJB:IIJE) ZA(IIJB:IIJE,IKB) = ZWORK41D(IIJB:IIJE)
......
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