Modified inputs for heatpump, added swithc for fan on/off.

Author: madDipanjan

src/heatpump.f90

@@ -10,11 +10,12 @@ module modheatpump
   logical, allocatable :: lhpptsrank(:)  ! Logical array indicating if heat pump points are on this rank
 
   real :: thl_dot_hp  ! Rate of change of temperature due to heat loss from the heat pump in [Km^3/s]
+  real :: w_hp_exhaust ! Exhaust velocity at heat pump points in [m/s]
 
 contains
 
   subroutine init_heatpump
-    use modglobal, only : lheatpump, nhppoints, Q_dot_hp, rhoa, cp, ifinput, cexpnr, ltempeq
+    use modglobal, only : lheatpump, nhppoints, Q_dot_hp, QH_dot_hp, rhoa, cp, ifinput, cexpnr, ltempeq, dxi, dyi
     use modmpi,    only : myid, comm3d, mpierr
     use decomp_2d, only : zstart, zend
     implicit none
@@ -50,12 +51,14 @@ subroutine init_heatpump
       end if
     end do
 
-    thl_dot_hp = Q_dot_hp / (nhppoints*rhoa*cp) ! Calculate temperature change rate from heat loss [Km^3/s]
+    thl_dot_hp = QH_dot_hp / (nhppoints*rhoa*cp) ! Calculate temperature change rate from heat loss [Km^3/s]
 
+    w_hp_exhaust = Q_dot_hp / (nhppoints*dxi*dyi) ! Calculate exhaust velocity at heat pump points [m/s]
+
   end subroutine init_heatpump
 
   subroutine heatpump
-    use modglobal,  only : lheatpump, nhppoints, w_hp_exhaust, dxi, dyi, dzfi, ltempeq
+    use modglobal,  only : lheatpump, lfan_hp, nhppoints, dxi, dyi, dzfi, ltempeq
     use modfields,  only : wm, w0, wp, thlp
     use modmpi,     only : myidx, myidy
     use decomp_2d,  only : zsize
@@ -71,9 +74,12 @@ subroutine heatpump
         j = idhppts_global(n,2) - myidy*zsize(2)
         k = idhppts_global(n,3)
 
-        wm(i,j,k+1) = w_hp_exhaust ! Set exhaust velocity at heat pump point [m/s], at input 'w' cell face k+1
-        w0(i,j,k+1) = w_hp_exhaust
-        !wp(i,j,k+1) = 0.
+        if (lfan_hp) then ! Heat pump fan is on
+          wm(i,j,k+1) = w_hp_exhaust ! Set exhaust velocity at heat pump point [m/s], at input 'w' cell face k+1
+          w0(i,j,k+1) = w_hp_exhaust
+          !wp(i,j,k+1) = 0.
+        end if
+        
         thlp(i,j,k) = thlp(i,j,k) - thl_dot_hp * dxi * dyi * dzfi(k)  ! [K/s], at cell center k
       end if
     end do

src/modglobal.f90

@@ -363,9 +363,10 @@ module modglobal
 
    ! Heat pump
    logical :: lheatpump = .false.   !< switch for heat pump
+   logical :: lfan_hp = .true.      !< switch for heat pump fan on/off
    integer :: nhppoints = 0         !< number of heat pump grid points
-   real    :: Q_dot_hp = 0.         !< rate of heat extracted from the ambient air by the heat pump (W)
-   real    :: w_hp_exhaust = 0.     !< exhaust vertical velocity of the heat pump (m/s)
+   real    :: QH_dot_hp = 0.        !< Total rate of heat extracted from the ambient air by the heat pump (W)
+   real    :: Q_dot_hp = 0.         !< Total volume flux going out from the heat pump in vertical direction(m^3/s)
 
    ! Poisson solver
    integer, parameter :: POISS_FFT2D = 0, &

src/modstartup.f90

@@ -77,7 +77,7 @@ subroutine readnamelists
                                     lrandomize, prandtlturb, fkar, lwritefac, dtfac, tfac, tnextfac, &
                                     ltrees,ntrees,Qstar,dQdt,lad,lsize,r_s,cd,dec,ud,ltreedump, &
                                     lpurif,npurif,Qpu,epu, &
-                                    lheatpump, nhppoints, Q_dot_hp, w_hp_exhaust
+                                    lheatpump,lfan_hp,nhppoints,Q_dot_hp,QH_dot_hp
       use modsurfdata,       only : z0, z0h,  wtsurf, wttop, wqtop, wqsurf, wsvsurf, wsvtop, wsvsurfdum, wsvtopdum, ps, thvs, thls, thl_top, qt_top, qts
       use modfields,         only : initfields, dpdx, ncname
       use modpois,           only : initpois
@@ -170,7 +170,7 @@ subroutine readnamelists
       namelist/PURIFS/ &
          lpurif, npurif, Qpu, epu
       namelist/HEATPUMP/ &
-         lheatpump, nhppoints, Q_dot_hp, w_hp_exhaust
+         lheatpump, lfan_hp, nhppoints, Q_dot_hp, QH_dot_hp
 
       if (myid == 0) then
          if (command_argument_count() >= 1) then
@@ -603,9 +603,10 @@ subroutine readnamelists
       call MPI_BCAST(Qpu, 1, MY_REAL, 0, comm3d, mpierr)
       call MPI_BCAST(epu, 1, MY_REAL, 0, comm3d, mpierr)
       call MPI_BCAST(lheatpump, 1, MPI_LOGICAL, 0, comm3d, mpierr)
+      call MPI_BCAST(lfan_hp, 1, MPI_LOGICAL, 0, comm3d, mpierr)
       call MPI_BCAST(nhppoints, 1, MPI_INTEGER, 0, comm3d, mpierr)
       call MPI_BCAST(Q_dot_hp, 1, MY_REAL, 0, comm3d, mpierr)
-      call MPI_BCAST(w_hp_exhaust, 1, MY_REAL, 0, comm3d, mpierr)
+      call MPI_BCAST(QH_dot_hp, 1, MY_REAL, 0, comm3d, mpierr)
 
       ! ! Allocate and initialize core modules
       ! call initglobal