Commits

Anonymous committed ae71c88

Refactored function/subroutine headers

Comments (0)

Files changed (189)

src/pflotran/auxiliary.F90

 contains
 
 ! ************************************************************************** !
-!
-! AuxInit: Nullifies pointers in auxiliary object
-! author: Glenn Hammond
-! date: 04/09/08
-!
-! ************************************************************************** !
+
 subroutine AuxInit(aux)
+  ! 
+  ! Nullifies pointers in auxiliary object
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 04/09/08
+  ! 
 
   implicit none
   
 end subroutine AuxInit
 
 ! ************************************************************************** !
-!
-! AuxDestroy: Deallocates any allocated pointers in auxiliary object
-! author: Glenn Hammond
-! date: 04/09/08
-!
-! ************************************************************************** !
+
 subroutine AuxDestroy(aux)
+  ! 
+  ! Deallocates any allocated pointers in auxiliary object
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 04/09/08
+  ! 
 
   implicit none
   

src/pflotran/block_solve.F90

             bl3dsolb
       
 contains
-      
+
+! ************************************************************************** !
+
 subroutine bl3dfac(n, k, E, D, F, pivot)
 
 !	This version of bl3dfac dated March 15, 2000.
   
 end subroutine bl3dfac
 
-!************************ END OF bl3dfac ********************************
+! ************************************************************************** !
 
 subroutine bl3dsol(n, k, E, D, F, pivot, nrhs, rhs)
 
       
 end subroutine bl3dsol
 
-!************************ END OF bl3dsol ********************************
-
+! ************************************************************************** !
 
 subroutine bl3dsolf(n, k, E, D, F, pivot, nrhs, rhs)
 
 
 end subroutine bl3dsolf
 
+! ************************************************************************** !
+
 subroutine bl3dsolb(n, k, E, D, F, pivot, nrhs, rhs)
 
 !       This version of bl3dsol dated April 20, 2000.

src/pflotran/block_tridiag.F90

             solbtb
      
 contains
-      
+
+! ************************************************************************** !
+
 subroutine decbt (m, n, ndim, a, b, c, ip, ier)
       
   use PFLOTRAN_Constants_module
      return
 
 end subroutine decbt
-!-----------------------  end of subroutine decbt  ----------
+
+! ************************************************************************** !
 
 subroutine solbt (m, n, ndim, a, b, c, ip, y)
 
   return
 end subroutine solbt
 
-!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
+! ************************************************************************** !
 
 subroutine solbtf (m, n, ndim, a, b, c, ip, y)
 
 
   return
 end subroutine solbtf
-      
+
+! ************************************************************************** !
+
 subroutine solbtb (m, n, ndim, a, b, c, ip, y)
 
   implicit none
   return
 end subroutine solbtb
 
-!-----------------------  end of subroutine solbt  -----
+! ************************************************************************** !
 
 subroutine dec (n, ndim, a, ip, ier)
 
 
 end subroutine dec
 
-!----------------------- end of subroutine dec -----------
+! ************************************************************************** !
 
 subroutine sol (n, ndim, a, b, ip)
 

src/pflotran/checkpoint.F90

 
 #ifndef PROCESS_MODEL
 
+! ************************************************************************** !
+
 subroutine Checkpoint(realization, &
                       flow_time_steps,flow_cumulative_newton_iterations, &
                       flow_cumulative_time_step_cuts, &
 
 end subroutine Checkpoint
 
+! ************************************************************************** !
+
 subroutine Restart(realization, &
                    flow_time_steps, &
                    flow_cumulative_newton_iterations, &
 end subroutine Restart
 
 ! ************************************************************************** !
-!
-! CheckpointRegisterBagHeader: Registers entities within the PETSc bag to
-!                              header
-! author: Glenn Hammond
-! date: 01/12/12
-!
-! ************************************************************************** !
+
 subroutine CheckpointRegisterBagHeader(bag,header)
+  ! 
+  ! Registers entities within the PETSc bag to
+  ! header
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 01/12/12
+  ! 
 
   implicit none
   
 #endif
 
 ! ************************************************************************** !
-!
-! OpenCheckpointFile: Opens checkpoint file; sets format
-! author: Glenn Hammond
-! date: 07/26/13
-!
-! ************************************************************************** !
+
 subroutine OpenCheckpointFile(viewer,id,option,id_stamp)
+  ! 
+  ! Opens checkpoint file; sets format
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 07/26/13
+  ! 
 
   use Option_module
   use String_module, only : StringNull
 end subroutine OpenCheckpointFile
 
 ! ************************************************************************** !
-!
-! CloseCheckpointFile: Closes checkpoint file
-! author: Glenn Hammond
-! date: 07/26/13
-!
-! ************************************************************************** !
+
 subroutine CloseCheckpointFile(viewer)
+  ! 
+  ! Closes checkpoint file
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 07/26/13
+  ! 
 
   use Option_module
   
 end subroutine CloseCheckpointFile
 
 ! ************************************************************************** !
-!
-! CheckpointFlowProcessModel: Checkpoints flow process model vectors
-! author: Glenn Hammond
-! date: 07/26/13
-!
-! ************************************************************************** !
+
 subroutine CheckpointFlowProcessModel(viewer,realization)
+  ! 
+  ! Checkpoints flow process model vectors
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 07/26/13
+  ! 
 
   use Option_module
   use Realization_class
 end subroutine CheckpointFlowProcessModel
 
 ! ************************************************************************** !
-!
-! RestartFlowProcessModel: Restarts flow process model
-! author: Glenn Hammond
-! date: 07/26/13
-!
-! ************************************************************************** !
+
 subroutine RestartFlowProcessModel(viewer,realization)
+  ! 
+  ! Restarts flow process model
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 07/26/13
+  ! 
 
   use Option_module
   use Realization_class

src/pflotran/co2_span_wagner.F90

       private
 
       contains
-    
+
+! ************************************************************************** !
+
 subroutine initialize_span_wagner(itable,myrank,option)
 
       use Option_module
   endif
 
 end subroutine initialize_span_wagner
-      
+
+! ************************************************************************** !
+
 subroutine co2_span_wagner(pl,tl,rho,dddt,dddp,fg,dfgdp,dfgdt, &
       eng,ent,dhdt,dhdp,visc,dvdt,dvdp,iflag,itable)
       
       
 end subroutine co2_span_wagner
 
-
+! ************************************************************************** !
 
 subroutine guess(lguess,uguess)
 
 
 end subroutine guess
 
+! ************************************************************************** !
+
 subroutine co2den(den,f,df)
      
     IMPLICIT NONE
       return
 end subroutine co2den
 
+! ************************************************************************** !
+
       PetscReal function psi(i,del2,tau2)
       implicit none
 !     PetscReal :: psi
 
       end function psi
 
+! ************************************************************************** !
+
 subroutine dphiodtau(dr,del2,tau2)
       implicit none
       PetscInt :: i
 
 end subroutine dphiodtau
 
+! ************************************************************************** !
+
 subroutine dphiodtautau(dr,del2,tau2)
 
       implicit none
 
 end subroutine dphiodtautau
 
+! ************************************************************************** !
+
 subroutine phir(r_helm,del2,tau2)
 
 !     residual helmholtz energy: Span & Wagner (1996), p. 1544, eq. (6.5)
       return
 end subroutine phir
 
+! ************************************************************************** !
+
 subroutine dphirddel(dr,del2,tau2)
 
 !     Span & Wagner (1996) Table 32
       
 end subroutine dphirddel
 
+! ************************************************************************** !
+
 subroutine dphirdddel(dpdd,del2,tau2)
       implicit none
       PetscInt :: i
       dpdd = derdr_helm
       
 end subroutine dphirdddel
-      
+
+! ************************************************************************** !
+
 subroutine dphirdtau(dpdtau,del2,tau2)
 
       implicit none
       
 end subroutine dphirdtau
 
+! ************************************************************************** !
+
 subroutine dphirdtautau(dpdtt,del2,tau2)
 
       implicit none
       
 end subroutine dphirdtautau
 
+! ************************************************************************** !
+
 subroutine dphirddeldtau(dpddt,del2,tau2)
 
       implicit none
       
 end subroutine dphirddeldtau
 
+! ************************************************************************** !
+
 function dpsiddel(i,del2,tau2)
       implicit none
       PetscReal :: dpsiddel
 
       end function dpsiddel
 
+! ************************************************************************** !
+
       function d2psiddel2(i,del2,tau2)
       implicit none
       PetscReal :: d2psiddel2
 
       end function d2psiddel2
 
+! ************************************************************************** !
+
       function dpsidtau(i,del2,tau2)
       implicit none
       PetscReal :: dpsidtau
 
       end function dpsidtau
 
+! ************************************************************************** !
+
       function d2psidtau2(i,del2,tau2)
       implicit none
       PetscReal :: d2psidtau2
 
       end function d2psidtau2
 
+! ************************************************************************** !
+
       function d2psiddeltau(i,del2,tau2)
       implicit none
       PetscReal :: d2psiddeltau
 
       end function d2psiddeltau
 
+! ************************************************************************** !
+
       function theta(i,del2,tau2)
       implicit none
       PetscReal :: theta
 
       end function theta
 
+! ************************************************************************** !
+
       function capdel(i,del2,tau2)
       implicit none
       PetscReal :: capdel
       capdel=theta1*theta1+capb(i)*(((del2-1.d0)**2)**aco2(i))
       end function capdel
 
+! ************************************************************************** !
+
       function dcapdelddel(i,del2,tau2)
       implicit none
       PetscReal :: dcapdelddel
 
       end function dcapdelddel
 
+! ************************************************************************** !
+
       function d2capdelddel2(i,del2,tau2) ! d^2 Delta/d delta^2
       implicit none
       PetscReal :: d2capdelddel2
 
       end function d2capdelddel2
 
+! ************************************************************************** !
+
       function ddelbiddel(i,del2,tau2)
       implicit none
       PetscReal :: ddelbiddel
 
       end function ddelbiddel
 
+! ************************************************************************** !
+
       function d2delbiddel2(i,del2,tau2)
       implicit none
       PetscReal :: d2delbiddel2
 
       end function d2delbiddel2
 
+! ************************************************************************** !
+
       function ddelbidtau(i,del2,tau2)
       implicit none
       PetscReal :: ddelbidtau
 
       end function ddelbidtau
 
+! ************************************************************************** !
+
       function d2delbidtau2(i,del2,tau2)
       implicit none
       PetscReal :: d2delbidtau2
 
       end function d2delbidtau2
 
+! ************************************************************************** !
+
       function d2delbiddeltau(i,del2,tau2)
       implicit none
       PetscReal :: d2delbiddeltau
 
 end function d2delbiddeltau
 
+! ************************************************************************** !
+
 subroutine vappr(tm,ps,dertp,derpt,ifl1)
 
       implicit none
       endif
       
 end subroutine vappr
-      
+
+! ************************************************************************** !
+
 subroutine viscosity(p,t,rho,drhodp,drhodt,mu,dmudt,dmudp)
 
 ! Fenghour, A., W. A. Wakeham, and V. Vesovic,
 
 end subroutine viscosity
 
+! ************************************************************************** !
+
 subroutine dissco2(p,t,mco2,fg,mol)
 
       implicit none

src/pflotran/co2_span_wagner_orig.F90

       private
 
       contains
-    
+
+! ************************************************************************** !
+
 subroutine initialize_span_wagner(itable,myrank)
 
       implicit none
   endif
 
 end subroutine initialize_span_wagner      
-      
+
+! ************************************************************************** !
+
 subroutine co2_span_wagner(pl,tl,rho,dddt,dddp,fg,dfgdp,dfgdt, &
       eng,ent,dhdt,dhdp,visc,dvdt,dvdp,itable)
       
       
 end subroutine co2_span_wagner
 
-
+! ************************************************************************** !
 
 subroutine guess(lguess,uguess)
 
 
 end subroutine guess
 
+! ************************************************************************** !
+
 subroutine co2den(den,f,df)
      
 	  IMPLICIT NONE
       return
 end subroutine co2den
 
+! ************************************************************************** !
+
 double precision function psi(i,del2,tau2)
       implicit none
 !     real*8 :: psi
 
       end function psi
 
+! ************************************************************************** !
+
       subroutine dphiodtau(dr,del2,tau2)
       implicit none
       integer :: i
 
 end subroutine dphiodtau
 
+! ************************************************************************** !
+
 subroutine dphiodtautau(dr,del2,tau2)
 
       implicit none
 
 end subroutine dphiodtautau
 
+! ************************************************************************** !
+
 subroutine phir(r_helm,del2,tau2)
 
 !     residual helmholtz energy: Span & Wagner (1996), p. 1544, eq. (6.5)
       return
 end subroutine phir
 
+! ************************************************************************** !
+
 subroutine dphirddel(dr,del2,tau2)
 
 !     Span & Wagner (1996) Table 32
       
 end subroutine dphirddel
 
+! ************************************************************************** !
+
 subroutine dphirdddel(dpdd,del2,tau2)
       implicit none
       integer :: i
       dpdd = derdr_helm
       
 end subroutine dphirdddel
-      
+
+! ************************************************************************** !
+
 subroutine dphirdtau(dpdtau,del2,tau2)
 
       implicit none
       
 end subroutine dphirdtau
 
+! ************************************************************************** !
+
 subroutine dphirdtautau(dpdtt,del2,tau2)
 
       implicit none
       
 end subroutine dphirdtautau
 
+! ************************************************************************** !
+
 subroutine dphirddeldtau(dpddt,del2,tau2)
 
       implicit none
       
 end subroutine dphirddeldtau
 
+! ************************************************************************** !
+
 function dpsiddel(i,del2,tau2)
       implicit none
       real*8 :: dpsiddel
 
       end function dpsiddel
 
+! ************************************************************************** !
+
       function d2psiddel2(i,del2,tau2)
       implicit none
       real*8 :: d2psiddel2
 
       end function d2psiddel2
 
+! ************************************************************************** !
+
       function dpsidtau(i,del2,tau2)
       implicit none
       real*8 :: dpsidtau
 
       end function dpsidtau
 
+! ************************************************************************** !
+
       function d2psidtau2(i,del2,tau2)
       implicit none
       real*8 :: d2psidtau2
 
       end function d2psidtau2
 
+! ************************************************************************** !
+
       function d2psiddeltau(i,del2,tau2)
       implicit none
       real*8 :: d2psiddeltau
 
       end function d2psiddeltau
 
+! ************************************************************************** !
+
       function theta(i,del2,tau2)
       implicit none
       real*8 :: theta
 
       end function theta
 
+! ************************************************************************** !
+
       function capdel(i,del2,tau2)
       implicit none
       real*8 :: capdel
       capdel=theta1*theta1+capb(i)*(((del2-1.d0)**2)**aco2(i))
       end function capdel
 
+! ************************************************************************** !
+
       function dcapdelddel(i,del2,tau2)
       implicit none
       real*8 :: dcapdelddel
 
       end function dcapdelddel
 
+! ************************************************************************** !
+
       function d2capdelddel2(i,del2,tau2)
       implicit none
       real*8 :: d2capdelddel2
 
       end function d2capdelddel2
 
+! ************************************************************************** !
+
       function ddelbiddel(i,del2,tau2)
       implicit none
       real*8 :: ddelbiddel
 
       end function ddelbiddel
 
+! ************************************************************************** !
+
       function d2delbiddel2(i,del2,tau2)
       implicit none
       real*8 :: d2delbiddel2
 
       end function d2delbiddel2
 
+! ************************************************************************** !
+
       function ddelbidtau(i,del2,tau2)
       implicit none
       real*8 :: ddelbidtau
 
       end function ddelbidtau
 
+! ************************************************************************** !
+
       function d2delbidtau2(i,del2,tau2)
       implicit none
       real*8 :: d2delbidtau2
 
       end function d2delbidtau2
 
+! ************************************************************************** !
+
       function d2delbiddeltau(i,del2,tau2)
       implicit none
       real*8 :: d2delbiddeltau
 
 end function d2delbiddeltau
 
+! ************************************************************************** !
+
 subroutine vappr(tm,ps,dertp,derpt,ifl1)
 
       implicit none
       endif
       
 end subroutine vappr
-      
+
+! ************************************************************************** !
+
 subroutine viscosity(p,t,rho,drhodp,drhodt,mu,dmudt,dmudp)
 
 ! Fenghour, A., W. A. Wakeham, and V. Vesovic, 
 
 end subroutine viscosity
 
+! ************************************************************************** !
+
 subroutine dissco2(p,t,mco2,fg,mol)
 
       implicit none

src/pflotran/co2_span_wagner_spline.F90

 
   return
 end subroutine sw_spline_read
-  
-  
+
+! ************************************************************************** !
+
 subroutine sw_prop(tx,px,rho,h,u,fg)
       
        use spline_module

src/pflotran/co2_sw.F90

 
       public initialize_sw_interp, co2_sw_interp      
  contains
- 
- 
-! prepare data table for interpolation 
+
+! ************************************************************************** !
+
 subroutine initialize_sw_interp(itable,myrank)
+  ! 
+  ! prepare data table for interpolation
+  ! 
       use co2_span_wagner_module, only: co2_span_wagner, vappr
         
         implicit none
      end select
       
 end subroutine initialize_sw_interp 
-     
 
- 
-     
+! ************************************************************************** !
+
 PetscReal function co2_prop_spwag(ip,it,iv)
      implicit none 
  !    PetscReal co2_prop_spwag
       endif 
      
   end function co2_prop_spwag
-           
-                     
- ! 2-d function interpolation     
+
+! ************************************************************************** !
+
 subroutine interp(x1,x2,y)
+  ! 
+  ! 2-d function interpolation
+  ! 
 
       use co2_span_wagner_module, only: vappr, co2_span_wagner
       implicit none 
       endif  
 end subroutine 
 
+! ************************************************************************** !
 
 subroutine co2_sw_interp(p,tc,rho,dddt,dddp,fg,dfgdp,dfgdt, &
         eng,ent,dhdt,dhdp,visc,dvdt,dvdp,itable)

src/pflotran/co2_sw_rtsafe.F90

 
   contains
 
+! ************************************************************************** !
+
   FUNCTION rtsafe(funcd,x1,x2,xacc)
 
   IMPLICIT NONE
   end do
   print *,'rtsafe: exceeded maximum iterations'
   END FUNCTION rtsafe
-  
+
+! ************************************************************************** !
+
 subroutine bracket(func,x1,x2)
   
   implicit none

src/pflotran/co2eos.F90

 
 contains
 
-!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+! ************************************************************************** !
 
-! VERSION/REVISION HISTORY
- 
-! $Id: co2eos.F90,v 1.1.1.1 2004/07/30 21:49:42 lichtner Exp $
-! $Log: co2eos.F90,v $
-! Revision 1.1.1.1  2004/07/30 21:49:42  lichtner
-! initial import
-!
-! Revision 1.2  2004/01/10 18:32:06  lichtner
-! Began work on 2 phase capability.
-!
-! Revision 1.1.1.1  2003/11/23 20:12:46  lichtner
-! initial entry
-!
-! Revision 1.2  2003/05/09 15:22:41  lichtner
-! commented out icall statements
-!
-! Revision 1.1.1.1  2003/03/03 01:33:27  lichtner
-! PFLOTRAN initial implementation
-!
-! Revision 1.6  2002/09/28 17:25:49  lichtner
-! Improved fit of dissolved CO2.
-!
-! Revision 1.5  2002/05/19 18:53:01  lichtner
-! Added documentation of CO2 EOS
-!
-! Revision 1.4  2002/05/19 00:21:46  lichtner
-! Modified Crovetto (1991) fit to Henry's law to be consistent with 
-! Duan et al. CO2 EOS.
-!
-! Revision 1.3  2002/05/07 03:14:39  lichtner
-! Modified Henry's law subroutine to output the Poynting term.
-!
-! Revision 1.2  2002/05/04 18:28:34  lichtner
-! Added Duan and Weare CO2 eos and new Henry's law.
-!
-! Revision 1.1  2002/04/12 19:03:10  lichtner
-! Initial entry
-!
-
-!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-!********1*********2*********3*********4*********5*********6*********7**
 subroutine CO2(TX,PCX,DC,FC,PHI,HC)
+  ! 
+  ! VERSION/REVISION HISTORY
+  ! $Id: co2eos.F90,v 1.1.1.1 2004/07/30 21:49:42 lichtner Exp $
+  ! $Log: co2eos.F90,v $
+  ! Revision 1.1.1.1  2004/07/30 21:49:42  lichtner
+  ! initial import
+  ! Revision 1.2  2004/01/10 18:32:06  lichtner
+  ! Began work on 2 phase capability.
+  ! Revision 1.1.1.1  2003/11/23 20:12:46  lichtner
+  ! initial entry
+  ! Revision 1.2  2003/05/09 15:22:41  lichtner
+  ! commented out icall statements
+  ! Revision 1.1.1.1  2003/03/03 01:33:27  lichtner
+  ! PFLOTRAN initial implementation
+  ! Revision 1.6  2002/09/28 17:25:49  lichtner
+  ! Improved fit of dissolved CO2.
+  ! Revision 1.5  2002/05/19 18:53:01  lichtner
+  ! Added documentation of CO2 EOS
+  ! Revision 1.4  2002/05/19 00:21:46  lichtner
+  ! Modified Crovetto (1991) fit to Henry's law to be consistent with
+  ! Duan et al. CO2 EOS.
+  ! Revision 1.3  2002/05/07 03:14:39  lichtner
+  ! Modified Henry's law subroutine to output the Poynting term.
+  ! Revision 1.2  2002/05/04 18:28:34  lichtner
+  ! Added Duan and Weare CO2 eos and new Henry's law.
+  ! Revision 1.1  2002/04/12 19:03:10  lichtner
+  ! Initial entry
+  ! 
       
   use PFLOTRAN_Constants_module
 
       
       RETURN
 end subroutine CO2
-      
-!********1*********2*********3*********4*********5*********6*********7**
+
+! ************************************************************************** !
+
 subroutine MRK(Y,T,PCX,V,DV,AT)
       
       implicit none
            -(AT/(DSQRT(T)*V*(V+B))))
       RETURN
 end subroutine MRK
-      
-!********1*********2*********3*********4*********5*********6*********7**
+
+! ************************************************************************** !
+
 subroutine FUGACITY(Y,T,V,Z,PHI)
       
       implicit none
 
       RETURN
 end subroutine FUGACITY
-      
-!********1*********2*********3*********4*********5*********6*********7**
+
+! ************************************************************************** !
+
 subroutine ENTHALPY(T,V,Z,H)
 
       implicit none
       H = (URES+Z-(TREF/T)) * R * T + 8.1858447D5*.044D0
       RETURN
 end subroutine ENTHALPY
-      
-!********1*********2*********3*********4*********5*********6*********7**
-!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-! Subroutine: duanco2.f
 
-! Input: tt   [C]        temperature
-!        p    [Pa]       CO2 partial pressure
-
-! Output: fc   [bars]    CO2 fugacity
-!         phi  [-]       CO2 fugacity coefficient
-!         dc   [g/cm^3]  CO2 density
-
-! Duan, Z., Moller, N., and Weare, J.H. (1992) An equation of state for 
-! the CH4-CO2-H2O system: I. Pure systems from 0 to 1000 oC and 0 to   
-! 8000 bar. Geochimica Cosmochimica Acta, 56, 2605-2617.
-!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+! ************************************************************************** !
 
 subroutine duanco2 (tt,p,dc,fc,phi)
+  ! 
+  ! Subroutine: duanco2.f
+  ! Input: tt   [C]        temperature
+  ! p    [Pa]       CO2 partial pressure
+  ! Output: fc   [bars]    CO2 fugacity
+  ! phi  [-]       CO2 fugacity coefficient
+  ! dc   [g/cm^3]  CO2 density
+  ! Duan, Z., Moller, N., and Weare, J.H. (1992) An equation of state for
+  ! the CH4-CO2-H2O system: I. Pure systems from 0 to 1000 oC and 0 to
+  ! 8000 bar. Geochimica Cosmochimica Acta, 56, 2605-2617.
+  ! 
 
       implicit none
       
       return
 
 end subroutine duanco2
-!----------------------------------------------
+
+! ************************************************************************** !
+
 subroutine Henry_duan_sun_0NaCl (p,tc,henry)
 
   implicit none
 
   return
 end subroutine Henry_duan_sun_0NaCl
-  
+
+! ************************************************************************** !
 
 subroutine Henry_duan_sun(tc,p,keqco2,phico2,lngamco2,mc,ma,psat,co2_aq_actcoef)
   
   !print *, 'keqco2: ', mu0,lngamco2,phico2,psat,yco2,t,p
   return
 end subroutine Henry_duan_sun
-  
+
+! ************************************************************************** !
+
 subroutine duan_sun_param(t,p,c,par)
   
   implicit none
   return
 end subroutine duan_sun_param
 
-
-!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-! Subroutine: henry.f
-
-! Input: tx   [C]       temperature
-!        pco2 [Pa]      CO2 partial pressure
-!        psys [Pa]      total system pressure ***not implemented***
-!        phi  [-]       CO2 fugacity coefficient
-
-! Output: xmole [-]     mole fraction CO2 in aqueous phase
-!         x1m   [-]     mass fraction CO2
-!         rkh   [Pa]  Henry constant
-!         poyn  [-]     Poynting factor
-
-! Crovetto, R. (1991) Evaluation of solubility data of the system CO2-H2O 
-! from 273�Z K to the critical point of water. Journal of Physical and 
-! Chemical Reference Data, 20(3), 575-589.
-!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+! ************************************************************************** !
 
 subroutine HENRY_co2_noderiv(xmole,x1m,tx,pcx,xphi,rkh,poyn)
+  ! 
+  ! Subroutine: henry.f
+  ! Input: tx   [C]       temperature
+  ! pco2 [Pa]      CO2 partial pressure
+  ! psys [Pa]      total system pressure ***not implemented***
+  ! phi  [-]       CO2 fugacity coefficient
+  ! Output: xmole [-]     mole fraction CO2 in aqueous phase
+  ! x1m   [-]     mass fraction CO2
+  ! rkh   [Pa]  Henry constant
+  ! poyn  [-]     Poynting factor
+  ! Crovetto, R. (1991) Evaluation of solubility data of the system CO2-H2O
+  ! from 273�Z K to the critical point of water. Journal of Physical and
+  ! Chemical Reference Data, 20(3), 575-589.
+  ! 
 
 !     input:
 !     tx   [C]  temperature
 
       RETURN
 end subroutine HENRY_CO2_NODERIV
-      
+
+! ************************************************************************** !
+
 subroutine HENRY_sullivan (TX,PCX,PS,FC,X1M,XCO2,HP)
 
       implicit none
 
       RETURN
 end subroutine HENRY_sullivan
-      
-!********1*********2*********3*********4*********5*********6*********7**
+
+! ************************************************************************** !
+
 subroutine SOLUT(PCX,TX,HSOL)
       implicit none
       
       HSOL = HSOL * 1.D6
       RETURN
 end subroutine SOLUT
-      
-!********1*********2*********3*********4*********5*********6*********7**
+
+! ************************************************************************** !
+
 subroutine DENMIX(TX,DW,X1M,D1M)
       implicit none
       
 
       RETURN
 end subroutine DENMIX
-      
-!********1*********2*********3*********4*********5*********6*********7**
+
+! ************************************************************************** !
+
 subroutine VISCO2(TX,DC,VC)
       implicit none
 !
       VC = VC * 1.0D-06
       RETURN
 end subroutine VISCO2
-                  
-!********1*********2*********3*********4*********5*********6*********7**
+
+! ************************************************************************** !
+
 subroutine SAT(T,P)
 !--------- Fast SAT M.J.O'Sullivan - 17 SEPT 1990 ---------
 !
     1 FORMAT(A1,'TEMPERATURE = ',E12.6,'  OUT OF RANGE IN SAT ')
       RETURN
 end subroutine SAT
-      
+
+! ************************************************************************** !
+
 subroutine COWAT0(TF,PP,D,U)
 !--------- Fast COWAT M.J.O'Sullivan - 17 SEPT 1990 ---------
 
   102 FORMAT(1H ,5X,A6,5X,I2,2X,E20.10)
       RETURN
 end subroutine COWAT0
-      
+
+! ************************************************************************** !
+
 subroutine SUPST(T,P,D,U)
 !--------- Fast SUPST M.J.O'Sullivan - 17 SEPT 1990 ---------
 ! SUPST    1.0 S     1 February  1991
       RETURN
 end subroutine SUPST
 
+! ************************************************************************** !
+
 subroutine TSAT(PX,TX00,TS)
       implicit none
       
       goto 1
       
 end subroutine TSAT
-      
+
+! ************************************************************************** !
+
 subroutine SIGMA(T,ST)
       implicit none
 !
       ST=0.
       RETURN
 end subroutine SIGMA
-      
+
+! ************************************************************************** !
+
 subroutine VIS(T,P,D,VW,VS,PS)
       implicit none
       
       IF(T.GT.350.) VS=1.E-7*(V1+.353*D+676.5E-6*D**2+102.1E-9*D**3)
       RETURN
 end subroutine VIS
-      
+
+! ************************************************************************** !
+
 subroutine VISW0(T,P,PS,VW)
       implicit none
 
 
       RETURN
 end subroutine VISW0
-      
+
+! ************************************************************************** !
+
 subroutine VISS(T,P,D,VS)
       implicit none
 
       RETURN
 end subroutine VISS
 
+! ************************************************************************** !
+
 subroutine THERC(T,P,D,CONW,CONS,PS)
       implicit none
 

src/pflotran/communicator_base.F90

 contains
 
 ! ************************************************************************** !
-!
-! CommCreateProcessorGroups: Splits MPI_COMM_WORLD into N separate
-!                            processor groups
-! author: Glenn Hammond
-! date: 08/11/09
-!
-! ************************************************************************** !
+
 subroutine CommCreateProcessorGroups(option,num_groups)
+  ! 
+  ! Splits MPI_COMM_WORLD into N separate
+  ! processor groups
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 08/11/09
+  ! 
 
   use Option_module
   

src/pflotran/condition.F90

 contains
 
 ! ************************************************************************** !
-!
-! FlowConditionCreate: Creates a condition
-! author: Glenn Hammond
-! date: 10/23/07
-!
-! ************************************************************************** !
+
 function FlowConditionCreate(option)
+  ! 
+  ! Creates a condition
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/23/07
+  ! 
 
   use Option_module
   
 end function FlowConditionCreate
 
 ! ************************************************************************** !
-!
-! TranConditionCreate: Creates a transport condition
-! author: Glenn Hammond
-! date: 10/23/07
-!
-! ************************************************************************** !
+
 function TranConditionCreate(option)
+  ! 
+  ! Creates a transport condition
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/23/07
+  ! 
 
   use Option_module
   
 end function TranConditionCreate
 
 ! ************************************************************************** !
-!
-! FlowGeneralConditionCreate: Creates a condition for general mode
-! author: Glenn Hammond
-! date: 05/26/11
-!
-! ************************************************************************** !
+
 function FlowGeneralConditionCreate(option)
+  ! 
+  ! Creates a condition for general mode
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 05/26/11
+  ! 
 
   use Option_module
   
 end function FlowGeneralConditionCreate
 
 ! ************************************************************************** !
-!
-! FlowGeneralSubConditionPtr: Returns a pointer to a subcondition, creating
-!                             them if necessary
-! author: Glenn Hammond
-! date: 06/09/11
-!
-! ************************************************************************** !
+
 function FlowGeneralSubConditionPtr(sub_condition_name,general, &
                                     option)
+  ! 
+  ! Returns a pointer to a subcondition, creating
+  ! them if necessary
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 06/09/11
+  ! 
 
   use Option_module
 
 end function FlowGeneralSubConditionPtr
 
 ! ************************************************************************** !
-!
-! FlowSubConditionCreate: Creates a sub_condition
-! author: Glenn Hammond
-! date: 02/04/08
-!
-! ************************************************************************** !
+
 function FlowSubConditionCreate(ndof)
+  ! 
+  ! Creates a sub_condition
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 02/04/08
+  ! 
 
   use Dataset_Ascii_class
   use Option_module
 end function FlowSubConditionCreate
 
 ! ************************************************************************** !
-!
-! GetFlowSubCondFromArrayByName: returns a pointer to a subcondition with
-!                                 matching name
-! author: Glenn Hammond
-! date: 06/02/08
-!
-! ************************************************************************** !
+
 function GetFlowSubCondFromArrayByName(sub_condition_ptr_list,name)
+  ! 
+  ! returns a pointer to a subcondition with
+  ! matching name
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 06/02/08
+  ! 
 
   use Input_Aux_module
   use String_module
 end function GetFlowSubCondFromArrayByName
 
 ! ************************************************************************** !
-!
-! FlowSubConditionVerify: Verifies the data in a subcondition
-! author: Glenn Hammond
-! date: 02/04/08
-!
-! ************************************************************************** !
+
 subroutine FlowSubConditionVerify(option, condition, sub_condition_name, &
                                   sub_condition, default_time_storage, &
                                   destroy_if_null)
+  ! 
+  ! Verifies the data in a subcondition
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 02/04/08
+  ! 
   use Time_Storage_module
   use Option_module
   use Dataset_module
 end subroutine FlowSubConditionVerify
 
 ! ************************************************************************** !
-!
-! FlowConditionRead: Reads a condition from the input file
-! author: Glenn Hammond
-! date: 10/31/07
-!
-! ************************************************************************** !
+
 subroutine FlowConditionRead(condition,input,option)
+  ! 
+  ! Reads a condition from the input file
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/31/07
+  ! 
 
   use Option_module
   use Input_Aux_module
 end subroutine FlowConditionRead
 
 ! ************************************************************************** !
-!
-! FlowConditionGeneralRead: Reads a condition from the input file for
-!                           general mode
-! author: Glenn Hammond
-! date: 09/14/11
-!
-! ************************************************************************** !
+
 subroutine FlowConditionGeneralRead(condition,input,option)
+  ! 
+  ! Reads a condition from the input file for
+  ! general mode
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 09/14/11
+  ! 
 
   use Option_module
   use Input_Aux_module
 end subroutine FlowConditionGeneralRead
 
 ! ************************************************************************** !
-!
-! TranConditionRead: Reads a transport condition from the input file
-! author: Glenn Hammond
-! date: 10/14/08
-!
-! ************************************************************************** !
+
 subroutine TranConditionRead(condition,constraint_list,reaction,input,option)
+  ! 
+  ! Reads a transport condition from the input file
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/14/08
+  ! 
 
   use Option_module
   use Input_Aux_module
 end subroutine TranConditionRead
 
 ! ************************************************************************** !
-!
-! ConditionReadValues: Read the value(s) of a condition variable
-! author: Glenn Hammond
-! date: 10/31/07
-!
-! ************************************************************************** !
+
 subroutine ConditionReadValues(input,option,keyword,string,dataset_base,units)
+  ! 
+  ! Read the value(s) of a condition variable
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/31/07
+  ! 
 
   use Input_Aux_module
   use String_module
 end subroutine ConditionReadValues
 
 ! ************************************************************************** !
-!
-! FlowConditionPrint: Prints flow condition info
-! author: Glenn Hammond
-! date: 12/04/08
-!
-! ************************************************************************** !
+
 subroutine FlowConditionPrint(condition,option)
+  ! 
+  ! Prints flow condition info
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 12/04/08
+  ! 
 
   use Option_module
   use Dataset_module
 end subroutine FlowConditionPrint
 
 ! ************************************************************************** !
-!
-! FlowConditionPrintSubCondition: Prints flow subcondition info
-! author: Glenn Hammond
-! date: 12/04/08
-!
-! ************************************************************************** !
+
 subroutine FlowConditionPrintSubCondition(subcondition,option)
+  ! 
+  ! Prints flow subcondition info
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 12/04/08
+  ! 
 
   use Option_module
   use Dataset_module
 end subroutine FlowConditionPrintSubCondition
 
 ! ************************************************************************** !
-!
-! SubConditionName: Return name of subcondition
-! author: Gautam Bisht
-! date: 10/16/13
-!
-! ************************************************************************** !
+
 function GetSubConditionName(subcon_itype)
+  ! 
+  ! SubConditionName: Return name of subcondition
+  ! 
+  ! Author: Gautam Bisht
+  ! Date: 10/16/13
+  ! 
 
   implicit none
 
 end function GetSubConditionName
 
 ! ************************************************************************** !
-!
-! FlowConditionUpdate: Updates a transient condition
-! author: Glenn Hammond
-! date: 11/02/07
-!
-! ************************************************************************** !
+
 subroutine FlowConditionUpdate(condition_list,option,time)
+  ! 
+  ! Updates a transient condition
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 11/02/07
+  ! 
 
   use Option_module
   use Dataset_module
 end subroutine FlowConditionUpdate
 
 ! ************************************************************************** !
-!
-! TranConditionUpdate: Updates a transient transport condition
-! author: Glenn Hammond
-! date: 11/02/07
-!
-! ************************************************************************** !
+
 subroutine TranConditionUpdate(condition_list,option,time)
+  ! 
+  ! Updates a transient transport condition
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 11/02/07
+  ! 
 
   use Option_module
   
 end subroutine TranConditionUpdate
 
 ! ************************************************************************** !
-!
-! FlowConditionInitList: Initializes a condition list
-! author: Glenn Hammond
-! date: 11/01/07
-!
-! ************************************************************************** !
+
 subroutine FlowConditionInitList(list)
+  ! 
+  ! Initializes a condition list
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 11/01/07
+  ! 
 
   implicit none
 
 end subroutine FlowConditionInitList
 
 ! ************************************************************************** !
-!
-! FlowConditionAddToList: Adds a new condition to a condition list
-! author: Glenn Hammond
-! date: 11/01/07
-!
-! ************************************************************************** !
+
 subroutine FlowConditionAddToList(new_condition,list)
+  ! 
+  ! Adds a new condition to a condition list
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 11/01/07
+  ! 
 
   implicit none
   
 end subroutine FlowConditionAddToList
 
 ! ************************************************************************** !
-!
-! FlowConditionGetPtrFromList: Returns a pointer to the condition matching &
-!                          condition_name
-! author: Glenn Hammond
-! date: 11/01/07
-!
-! ************************************************************************** !
+
 function FlowConditionGetPtrFromList(condition_name,condition_list)
+  ! 
+  ! Returns a pointer to the condition matching &
+  ! condition_name
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 11/01/07
+  ! 
 
   use String_module
   
 end function FlowConditionGetPtrFromList
 
 ! ************************************************************************** !
-!
-! TranConditionInitList: Initializes a transport condition list
-! author: Glenn Hammond
-! date: 10/13/08
-!
-! ************************************************************************** !
+
 subroutine TranConditionInitList(list)
+  ! 
+  ! Initializes a transport condition list
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/13/08
+  ! 
 
   implicit none
 
 end subroutine TranConditionInitList
 
 ! ************************************************************************** !
-!
-! TranConditionAddToList: Adds a new condition to a transport condition list
-! author: Glenn Hammond
-! date: 10/13/08
-!
-! ************************************************************************** !
+
 subroutine TranConditionAddToList(new_condition,list)
+  ! 
+  ! Adds a new condition to a transport condition list
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/13/08
+  ! 
 
   implicit none
   
 end subroutine TranConditionAddToList
 
 ! ************************************************************************** !
-!
-! TranConditionGetPtrFromList: Returns a pointer to the condition matching
-!                              condition_name
-! author: Glenn Hammond
-! date: 10/13/08
-!
-! ************************************************************************** !
+
 function TranConditionGetPtrFromList(condition_name,condition_list)
+  ! 
+  ! Returns a pointer to the condition matching
+  ! condition_name
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/13/08
+  ! 
 
   use String_module
 
 end function TranConditionGetPtrFromList
 
 ! ************************************************************************** !
-!
-! FlowConditionIsTransient: Returns PETSC_TRUE
-! author: Glenn Hammond
-! date: 01/12/11
-!
-! ************************************************************************** !
+
 function FlowConditionIsTransient(condition)
+  ! 
+  ! Returns PETSC_TRUE
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 01/12/11
+  ! 
 
   use Dataset_module
 
 end function FlowConditionIsTransient
 
 ! ************************************************************************** !
-!
-! FlowSubConditionIsTransient: Returns PETSC_TRUE
-! author: Glenn Hammond
-! date: 10/26/11
-!
-! ************************************************************************** !
+
 function FlowSubConditionIsTransient(sub_condition)
+  ! 
+  ! Returns PETSC_TRUE
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/26/11
+  ! 
 
   use Dataset_module
 
 end function FlowSubConditionIsTransient
 
 ! ************************************************************************** !
-!
-! FlowConditionDestroyList: Deallocates a list of conditions
-! author: Glenn Hammond
-! date: 11/01/07
-!
-! ************************************************************************** !
+
 subroutine FlowConditionDestroyList(condition_list)
+  ! 
+  ! Deallocates a list of conditions
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 11/01/07
+  ! 
 
   implicit none
   
 end subroutine FlowConditionDestroyList
 
 ! ************************************************************************** !
-!
-! FlowConditionDestroy: Deallocates a condition
-! author: Glenn Hammond
-! date: 10/23/07
-!
-! ************************************************************************** !
+
 subroutine FlowConditionDestroy(condition)
+  ! 
+  ! Deallocates a condition
+  ! 
+  ! Author: Glenn Hammond
+  ! Date: 10/23/07
+  ! 
 
   use Dataset_module
   use Dataset_Ascii_class
 end subroutine FlowConditionDestroy
 
 ! ************************************************************************** !
-!
-! FlowGeneralConditionDestroy:Destroys a general mode condition
-! author: Glenn Hammond
-! date: 05/26/11
-!
-! ************************************************************************** !
+
 subroutine FlowGeneralConditionDestroy(general_condition)
+  ! 
+  ! Destroys a general mode condition