Changeset 1859 for trunk/LMDZ.GENERIC

Timestamp:

Dec 16, 2017, 3:15:53 PM (7 years ago)

Author:

mturbet

Message:

add diagnostic of reevaporation of precipitation

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• physiq_mod.F90 (modified) (3 diffs)
• rain.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -275 +275 @@
       real zdqsrain(ngrid), zdqssnow(ngrid) ! Rain routine.
       real dqs_hyd(ngrid,nq)                ! Hydrol routine.
+      real reevap_precip(ngrid)             ! re-evaporation flux of precipitation (integrated over the atmospheric column)
                   
       ! For Tracers : (kg/kg_of_air/s)
@@ -1284 +1285 @@
 
                call rain(ngrid,nlayer,nq,ptimestep,pplev,pplay,pt,pdt,pq,pdq,            &
-                         zdtrain,zdqrain,zdqsrain,zdqssnow,cloudfrac)
+                         zdtrain,zdqrain,zdqsrain,zdqssnow,reevap_precip,cloudfrac)
 
                pdq(1:ngrid,1:nlayer,igcm_h2o_vap) = pdq(1:ngrid,1:nlayer,igcm_h2o_vap) &
@@ -2111 +2112 @@
                call writediagfi(ngrid,"rain","rainfall","kg m-2 s-1",2,zdqsrain)
                call writediagfi(ngrid,"snow","snowfall","kg m-2 s-1",2,zdqssnow)
+               call writediagfi(ngrid,"reevap","reevaporation of precipitation","kg m-2 s-1",2,reevap_precip)
             endif
 

trunk/LMDZ.GENERIC/libf/phystd/rain.F90

@@ -1 @@
-subroutine rain(ngrid,nlayer,nq,ptimestep,pplev,pplay,t,pdt,pq,pdq,d_t,dqrain,dqsrain,dqssnow,rneb)
+subroutine rain(ngrid,nlayer,nq,ptimestep,pplev,pplay,t,pdt,pq,pdq,d_t,dqrain,dqsrain,dqssnow,reevap_precip,rneb)
 
 
@@ -38 +38 @@
       real,intent(out) :: dqsrain(ngrid)  ! rain flux at the surface (kg.m-2.s-1)
       real,intent(out) :: dqssnow(ngrid)  ! snow flux at the surface (kg.m-2.s-1)
+      real,intent(out) :: reevap_precip(ngrid)  ! re-evaporation flux of precipitation integrated over the atmospheric column (kg.m-2.s-1)
       real,intent(in) :: rneb(ngrid,nlayer) ! cloud fraction
 
@@ -379 +380 @@
         dqrain(1:ngrid,1:nlayer,i_vap)=d_q(1:ngrid,1:nlayer)/ptimestep
         d_t(1:ngrid,1:nlayer)=d_t(1:ngrid,1:nlayer)/ptimestep
+        do i=1,ngrid
+           reevap_precip(i)=0.
+           do k=1,nlayer
+              reevap_precip(i)=reevap_precip(i)+dqrain(i,k,i_vap)*l2c(i,k)
+           enddo
+        enddo
       else
         dqrain(1:ngrid,1:nlayer,i_vap)=0.0