From 85fb11d2f7738ad9444677d2f20fb98f30eff3de Mon Sep 17 00:00:00 2001 From: Jan Mandel Date: Fri, 24 Sep 2021 23:16:30 -0600 Subject: [PATCH] debugging --- femwind/fortran/module_utils.f90 | 2 ++ femwind/fortran/module_wrfout.f90 | 1 + 2 files changed, 3 insertions(+) diff --git a/femwind/fortran/module_utils.f90 b/femwind/fortran/module_utils.f90 index e8f3334..e6c7490 100644 --- a/femwind/fortran/module_utils.f90 +++ b/femwind/fortran/module_utils.f90 @@ -239,6 +239,7 @@ equivalence(rel,iel) !*** executable lb = lbound(a) ub = ubound(a) +print *,'shape ',shape(a),' bounds',lbound(a),ubound(a) lsum=0 do j=lb(2),ub(2) do i=lb(1),ub(1) @@ -247,6 +248,7 @@ lsum=0 !print *,i,j,rel,iel,lsum enddo enddo +print *,'lsum=',lsum get_chsum_2d = lsum diff --git a/femwind/fortran/module_wrfout.f90 b/femwind/fortran/module_wrfout.f90 index 216901c..dfce268 100644 --- a/femwind/fortran/module_wrfout.f90 +++ b/femwind/fortran/module_wrfout.f90 @@ -165,6 +165,7 @@ call ncopen(filename,nf90_write,ncid) call netcdf_write_2d(ncid,uf,"UF",istep) call netcdf_write_2d(ncid,vf,"VF",istep) chsum=ior(get_chsum_2d(uf),get_chsum_2d(vf)) +print *,'chsum=',chsum call netcdf_write_int(ncid,chsum,"CHSUM_FMW") call netcdf_write_int(ncid,frame0_fmw,"FRAME_FMW") call ncclose(ncid) -- 2.11.4.GIT