only write deltas if not doing transport opacities

evanoconnor · Oct 11, 2017 · 5f8a477 · 5f8a477
1 parent 0bb6087
commit 5f8a477
Showing 1 changed file with 9 additions and 7 deletions.
diff --git a/src/make_table_example.F90 b/src/make_table_example.F90
@@ -936,13 +936,15 @@ subroutine write_h5(filename,timestamp)
     call h5sclose_f(dspace_id, error)  
     cerror = cerror + error   
 
-    call h5screate_simple_f(rank, dims5, dspace_id, error)
-    call h5dcreate_f(file_id, "scattering_delta", H5T_NATIVE_DOUBLE, &
-         dspace_id, dset_id, error)
-    call h5dwrite_f(dset_id, H5T_NATIVE_DOUBLE,table_delta, dims5, error)
-    call h5dclose_f(dset_id, error)
-    call h5sclose_f(dspace_id, error)  
-    cerror = cerror + error   
+    if(.not. do_transport_opacities) then
+       call h5screate_simple_f(rank, dims5, dspace_id, error)
+       call h5dcreate_f(file_id, "scattering_delta", H5T_NATIVE_DOUBLE, &
+            dspace_id, dset_id, error)
+       call h5dwrite_f(dset_id, H5T_NATIVE_DOUBLE,table_delta, dims5, error)
+       call h5dclose_f(dset_id, error)
+       call h5sclose_f(dspace_id, error)
+       cerror = cerror + error
+    endif
 
     if (doing_inelastic.or.doing_epannihil) then