Reverted to intended 'rhod' reading

GEOS-ESM · Feb 26, 2025 · e3b264f · e3b264f · jrjoshi1 · Feb 26, 2025
1 parent 7cf6f60
commit e3b264f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/Process_Library/GOCART2G_MieMod.F90 b/Process_Library/GOCART2G_MieMod.F90
@@ -65,7 +65,7 @@ module GOCART2G_MieMod
       real, pointer  :: pmom(:,:,:,:,:) => Null() ! (r,c,b,m,p) moments of phase function
       real, pointer  :: gf(:,:) => Null()         ! (r,b) hygroscopic growth factor
       real, pointer  :: rhop(:,:) => Null()       ! (r,b) wet particle density [kg m-3]
-      real, pointer  :: rhod(:,:) => Null()       ! (r,b) wet particle density [kg m-3]
+      real, pointer  :: rhod(:,:) => Null()       ! (r,b) dry particle density [kg m-3]
       real, pointer  :: vol(:,:) => Null()        ! (r,b) wet particle volume [m3 kg-1]
       real, pointer  :: area(:,:) => Null()       ! (r,b) wet particle cross section [m2 kg-1]
       real, pointer  :: refr(:,:,:) => Null()     ! (r,c,b) real part of refractive index
@@ -291,7 +291,7 @@ type(GOCART2G_Mie) function GOCART2G_MieCreate ( MieFile, wavelengths, nmom, rc
       endif
 
 !     Dry particle density (will be pulled from wet particle radius)
-      rc = nf90_inq_varid(ncid,'rhod',ivarid)
+      rc = nf90_inq_varid(ncid,'rhop',ivarid)
       if(rc .ne. NF90_NOERR) then   ! not in table, fill in dummy variable
         rhod_table = -999.
       else