add subroutine changes

Author: sbryngelson

src/post_process/m_data_output.fpp

@@ -727,8 +727,8 @@ contains
 
             if (p > 0) then
                 err = DBMKOPTLIST(2, optlist)
-                err = DBADDIOPT(optlist, DBOPT_LO_OFFSET, lo_offset)
-                err = DBADDIOPT(optlist, DBOPT_HI_OFFSET, hi_offset)
+                err = DBADDIAOPT(optlist, DBOPT_LO_OFFSET, 3, lo_offset)
+                err = DBADDIAOPT(optlist, DBOPT_HI_OFFSET, 3, hi_offset)
                 if (grid_geometry == 3) then
                     err = DBPUTQM(dbfile, 'rectilinear_grid', 16, &
                                   'x', 1, 'y', 1, 'z', 1, &
@@ -745,8 +745,8 @@ contains
                 err = DBFREEOPTLIST(optlist)
             elseif (n > 0) then
                 err = DBMKOPTLIST(2, optlist)
-                err = DBADDIOPT(optlist, DBOPT_LO_OFFSET, lo_offset)
-                err = DBADDIOPT(optlist, DBOPT_HI_OFFSET, hi_offset)
+                err = DBADDIAOPT(optlist, DBOPT_LO_OFFSET, 2, lo_offset)
+                err = DBADDIAOPT(optlist, DBOPT_HI_OFFSET, 2, hi_offset)
                 err = DBPUTQM(dbfile, 'rectilinear_grid', 16, &
                               'x', 1, 'y', 1, 'z', 1, &
                               x_cb, y_cb, DB_F77NULL, dims, 2, &
@@ -755,7 +755,7 @@ contains
                 err = DBFREEOPTLIST(optlist)
             else
                 err = DBMKOPTLIST(2, optlist)
-                err = DBADDIOPT(optlist, DBOPT_LO_OFFSET, lo_offset)
+                err = DBADDIAOPT(optlist, DBOPT_LO_OFFSET, 1, lo_offset)
                 err = DBADDIOPT(optlist, DBOPT_HI_OFFSET, hi_offset)
                 err = DBPUTQM(dbfile, 'rectilinear_grid', 16, &
                               'x', 1, 'y', 1, 'z', 1, &

src/simulation/m_data_output.fpp

@@ -1144,6 +1144,8 @@ contains
         real(wp) :: G_local
         real(wp) :: dyn_p, T
         real(wp) :: damage_state
+        character(LEN=15) :: FMT_glb
+        character(len=30) :: FMT
 
         integer :: i, j, k, l, s, d !< Generic loop iterator
 
@@ -1165,6 +1167,13 @@ contains
 
         T = dflt_T_guess
 
+        ! Set format string based on precision
+        if (precision == 1) then
+            FMT_glb = 'F28.7'
+        else
+            FMT_glb = 'F28.16'
+        end if
+
         ! Non-dimensional time calculation
         if (time_stepper == 23) then
             nondim_time = mytime