1. Glenn Hammond
  2. pflotran-dev

Commits

Glenn Hammond  committed e2af493

K_H -> K_H_tilde

  • Participants
  • Parent commits b388f3b
  • Branches default

Comments (0)

Files changed (1)

File src/pflotran/general_aux.F90

View file
  • Ignore whitespace
   PetscReal :: den_gp, den_gt, hgp, hgt, dgp, dgt, u
   PetscReal :: krl, visl, dkrl_Se
   PetscReal :: krg, visg, dkrg_Se
-  PetscReal :: K_H, Ps
+  PetscReal :: K_H_tilde, Ps
   PetscReal :: guess, dummy
   PetscInt :: apid, cpid, vpid
   PetscErrorCode :: ierr
 
       call psat(gen_aux_var%temp,Ps,ierr)
       call Henry_air_noderiv(gen_aux_var%pres(lid),gen_aux_var%temp, &
-                             Ps,K_H)
-      gen_aux_var%pres(apid) = K_H*gen_aux_var%xmol(acid,lid)
+                             Ps,K_H_tilde)
+      gen_aux_var%pres(apid) = K_H_tilde*gen_aux_var%xmol(acid,lid)
       gen_aux_var%pres(vpid) = gen_aux_var%pres(lid) - gen_aux_var%pres(apid)
 
     case(GAS_STATE)
       gen_aux_var%pres(lid) = gen_aux_var%pres(gid) - gen_aux_var%pres(cpid)
 
       call Henry_air_noderiv(gen_aux_var%pres(lid),gen_aux_var%temp, &
-                             gen_aux_var%pres(vpid),K_H)
-      gen_aux_var%xmol(acid,lid) = gen_aux_var%pres(apid) / &
-                                  (gen_aux_var%pres(gid)*K_H)
+                             gen_aux_var%pres(vpid),K_H_tilde)
+      gen_aux_var%xmol(acid,lid) = gen_aux_var%pres(apid) / K_H_tilde
       gen_aux_var%xmol(wid,lid) = 1.d0 - gen_aux_var%xmol(acid,lid)
       gen_aux_var%xmol(acid,gid) = gen_aux_var%pres(apid) / gen_aux_var%pres(gid)
       gen_aux_var%xmol(wid,gid) = 1.d0 - gen_aux_var%xmol(acid,gid)