Debug ppoly0 values

adcroft · Mar 27, 2024 · a6b16f0 · a6b16f0
1 parent fb5fc62
commit a6b16f0
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/ALE/MOM_remapping.F90 b/src/ALE/MOM_remapping.F90
@@ -1792,6 +1792,8 @@ logical function remapping_unit_tests(verbose)
   ! u_sub =     |1|   2     |3|  4    | 5.5 |  6  |
   call PLM_reconstruction(3, (/2.,2.,1./), (/2.,4.,5.5/), ppoly0_E, ppoly0_coefs, h_neglect )
   call PLM_boundary_extrapolation(3, (/2.,2.,1./),(/2.,4.,5.5/), ppoly0_E, ppoly0_coefs, h_neglect)
+  call test%real_arr(3, ppoly0_E(:,1), (/1.,3.,5./), 'ppoly0_E(:,1)')
+  call test%real_arr(3, ppoly0_E(:,2), (/3.,5.,6./), 'ppoly0_E(:,1)')
   allocate(u_sub(6), uh_sub(6))
   call remap_src_to_sub_grid(3, (/2.,2.,1./), (/2.,4.,5.5/), ppoly0_E, ppoly0_coefs, 3, h_sub, &
                              h0_eff, isrc_start, isrc_end, isrc_max, isub_src, &