add i_water to netcdf output and input (#192)

compdyn · Nov 1, 2023 · f871c30 · f871c30
1 parent 9607437
commit f871c30
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/aero_data.F90 b/src/aero_data.F90
@@ -760,6 +760,9 @@ subroutine aero_data_output_netcdf(aero_data, ncid)
     call pmc_nc_write_real_1d(ncid, aero_data%kappa, &
          "aero_kappa", (/ dimid_aero_species /), unit="1", &
          long_name="hygroscopicity parameters (kappas) of aerosol species")
+    call pmc_nc_write_integer(ncid, aero_data%i_water, &
+         "aero_i_water", long_name="Index of aerosol water or " &
+         // "0 if water does not exist.")
     call fractal_output_netcdf(aero_data%fractal, ncid)
 
   end subroutine aero_data_output_netcdf
@@ -840,7 +843,7 @@ subroutine aero_data_input_netcdf(aero_data, ncid)
     end do
     call assert(377166446, aero_source_names == "")
 
-    call aero_data_set_water_index(aero_data)
+    call pmc_nc_read_integer(ncid, aero_data%i_water, "aero_i_water")
 
     call fractal_input_netcdf(aero_data%fractal, ncid)