Changeset 411 for codes/icosagcm/trunk/src

Timestamp:

06/09/16 02:20:07 (8 years ago)

Author:

ymipsl

Message:

adapt dcmip 2016 kessler physic
YM

File:

• codes/icosagcm/trunk/src/dcmip2016_kessler_physic.f90 (modified) (6 diffs)

codes/icosagcm/trunk/src/dcmip2016_kessler_physic.f90

@@ -39 +39 @@
 !     z      - heights of thermodynamic levels in the grid column (m)
 !     nz     - number of thermodynamic levels in the column
-!     rainnc - accumulated precip beneath the grid column (mm)
+!     precl  - Precipitation rate (m_water/s)
 !
 ! Output variables:
@@ -66 +66 @@
 !=======================================================================
 
-SUBROUTINE KESSLER(theta, qv, qc, qr, rho, pk, dt, z, nz, rainnc) &
-  BIND(c, name = "kessler")
+SUBROUTINE KESSLER(theta, qv, qc, qr, rho, pk, dt, z, nz, precl)
 
   IMPLICIT NONE
@@ -79 +78 @@
             qv      ,     & ! Water vapor mixing ratio (gm/gm)
             qc      ,     & ! Cloud water mixing ratio (gm/gm)
-            qr      ,     & ! Rain  water mixing ratio (gm/gm)
+            qr              ! Rain  water mixing ratio (gm/gm)
+
+  REAL(8), DIMENSION(nz), INTENT(IN) :: &
             rho             ! Dry air density (not mean state as in KW) (kg/m^3)
 
-  REAL(8), INTENT(INOUT) :: &
-            rainnc          ! Accumulated precip beneath the grid column (mm)
+  REAL(8), INTENT(OUT) :: &
+            precl          ! Precipitation rate (m_water / s)
 
   REAL(8), DIMENSION(nz), INTENT(IN) :: &
@@ -123 +124 @@
 
   ! Maximum time step size in accordance with CFL condition
+  if (dt .le. 0.d0) then
+    write(*,*) 'kessler.f90 called with nonpositive dt'
+    stop
+  end if
+
   dt_max = dt
   do k=1,nz-1
@@ -135 +141 @@
 
   ! Subcycle through rain process
+  precl = 0.d0
+
   do nt=1,rainsplit
 
-    ! Precipitation accumulated beneath the column (mm rain)
-    rainnc = rainnc + 1000.d0*rho(1)*qr(1)*velqr(1)*dt0/rhoqr
+    ! Precipitation rate (m/s)
+    precl = precl + rho(1) * qr(1) * velqr(1) / rhoqr
 
     ! Sedimentation term using upstream differencing
@@ -180 +188 @@
   end do
 
+  precl = precl / dble(rainsplit)
+
 END SUBROUTINE KESSLER
 
 !=======================================================================
 
+
+
 END MODULE dcmip2016_kessler_physic_mod