Adding idealized radiative tendency forcing

Author: natgeo-wong

forcing.f90

@@ -322,6 +322,26 @@ subroutine forcing
 
 endif
 
+!----------------------------------------------------------------------------
+! Temperature Tendency Forcing:
+! Simple Radiative Tendencies taken from Pauluis & Garner [2006]
+
+if(doradtendency.and.time.gt.timelargescale) then
+
+  do k = 1,nzm
+    do j=1,ny
+      do i=1,nx
+        if (tabs(i,j,k)>207.5) then
+            t(i,j,k) = t(i,j,k) - dtn * troptend / 86400
+        else
+            t(i,j,k) = t(i,j,k) + dtn * (200 - tabs(i,j,k)) / (5*86400)
+        end if
+      end do
+    end do
+  end do
+
+endif
+
 !-------------------------------------------------------------------------------
 
 if(.not.dosfcforcing.and.dodynamicocean) call sst_evolve()

params.f90

@@ -134,7 +134,7 @@ module params
   real:: bubble_dtemp = 0.
   real:: bubble_dq = 0.
 
-! The following are added by Kuang Lab at Harvard
+!!!!!  The following are added by Kuang Lab at Harvard !!!!!
 
 ! Options
 logical:: dompiensemble = .false. ! Subdomains defined in domains.f90 are run separately
@@ -146,5 +146,9 @@ module params
 real:: tperturbA = 1.     ! Default perturbation 1 time positive
 real:: qperturbA = 1.
 
+! Radiative tendencies as per Pauluis & Garner [2006]
+logical :: doradtendency = .false. 
+real :: troptend = 1.5 ! Convective tendency in Pauluis & Garner [2006]
+
 
 end module params

setparm.f90

@@ -47,7 +47,8 @@ subroutine setparm
 
 ! Parameters added by Kuang Lab at Harvard
 NAMELIST /KUANG_PARAMS/ dompiensemble, &
-                dolayerperturb, tperturbi, qperturbi, tperturbA, qperturbA ! linear response perturbation: layer by layer (Song Qiyu, 2022)
+                dolayerperturb, tperturbi, qperturbi, tperturbA, qperturbA, &
+                doradtendency, troptend 
 
 !bloss: Create dummy namelist, so that we can figure out error code
 !       for a mising namelist.  This lets us differentiate between