Commits

Richard Mills committed dc3a38f

Added a few comments on where I need to figure out some things related to freezing in surface flows.

Comments (0)

Files changed (1)

src/pflotran/surface_th.F90

 
 end subroutine SurfaceTHUpdateSurfBC
 
+! RTM: TODO: Figure out how this needs to be modified for surface freezing.
 ! ************************************************************************** !
 !> This routine updates source/sink term for the subsurface model
 !!
         local_id = cur_connection_set%id_dn(iconn)
         ghosted_id = surf_grid%nL2G(local_id)
 
+        ! RTM: TODO: I think this is incorrect in the case of surface freezing.
         ! Compute densities:
         call density(surf_global_aux_vars(ghosted_id)%temp(1), &
                      option%reference_pressure,den_surf_kg)
   ff_p = 0.d0
   Res  = 0.d0
 
+  ! RTM: Does this computed density get used anywhere?
   call density(option%reference_temperature,option%reference_pressure,rho)
 
   xc => surf_realization%discretization%grid%x