Merge pull request #12 from jswhit/bugfix/enkf-readobs

GitHub Issue #15. Fix a bug in EnKF that can cause segfault in a diagnostic print.
NOAA-EMC · Jun 12, 2020 · 93cd14b · 93cd14b
2 parents e2a33cc + afb30bb
commit 93cd14b
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/enkf/mpi_readobs.f90 b/src/enkf/mpi_readobs.f90
@@ -317,8 +317,8 @@ subroutine mpi_getobs(obspath, datestring, nobs_conv, nobs_oz, nobs_sat, nobs_to
                                      maxval(sprd_ob(nobs_conv+nobs_oz+1:nobs_tot))
        do nob =nobs_conv+nobs_oz+1 , nobs_tot
           if (sprd_ob(nob) > 1000.) then 
-             print *, nob, ' sat spread: ', sprd_ob(nob), ', ensmean_ob: ', ensmean_obbc(nob), &
-                           ', anal_ob: ', anal_ob(:,nob), ', mem_ob: ', mem_ob(nob)
+             print *, nob, trim(obtype(nob)),ob(nob),' sat spread: ', sprd_ob(nob), ', ensmean_ob: ', ensmean_obbc(nob), &
+                           ', anal_ob: ', anal_ob(:,nob)
           endif
        enddo
     endif