Changeset 2724 for trunk/LMDZ.GENERIC

Timestamp:

Jun 21, 2022, 11:05:46 AM (3 years ago)

Author:

aslmd

Message:

rneb_generic and RH_generic added

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• condensation_generic_mod.F90 (modified) (4 diffs)
• physiq_mod.F90 (modified) (6 diffs)

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

@@ -5 +5 @@
     
     subroutine condensation_generic(ngrid,nlayer,nq,ptimestep, pplev, pplay,    &
-                pt, pq, pdt, pdq, pdtlsc, pdqvaplsc, pdqliqlsc)
+                pt, pq, pdt, pdq, pdtlsc, pdqvaplsc, pdqliqlsc, rneb)
         use ioipsl_getin_p_mod, only: getin_p !-> to get the metallicity 
         use generic_cloud_common_h
@@ -43 +43 @@
         REAL, intent(out) :: pdqvaplsc(ngrid,nlayer,nq) ! incrementation de la vapeur du traceur
         REAL, intent(out) :: pdqliqlsc(ngrid,nlayer,nq) ! incrementation du traceur liquide
+        REAL, intent(out) :: rneb(ngrid,nlayer,nq)    ! fraction nuageuse
 
 !       Options : 
@@ -95 +96 @@
         qevap(1:ngrid,1:nlayer,1:nq)=0.0
         tevap(1:ngrid,1:nlayer)=0.0
+        rneb(1:ngrid,1:nlayer,1:nq) = 0.0
         ! Let's loop on tracers 
         do iq=1,nq
@@ -138 +140 @@
                                         End do ! niter
                                         zcond(i)=MAX(zcond(i),-(pq(i,k,igcm_generic_ice)+pdq(i,k,igcm_generic_ice)*ptimestep))
+
+                                        if (zcond(i) .gt. 0.) then
+                                                rneb(i,k,iq)=1
+                                        else
+                                                rneb(i,k,iq)=0.
+                                        endif
+
                                         zcond(i) = zcond(i)/ptimestep
                                 ENDDO ! i=1,ngrid

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

@@ -17 +17 @@
       use radinc_h, only : L_NSPECTI,L_NSPECTV,naerkind
       use watercommon_h, only : RLVTT, Psat_water,epsi,su_watercycle, RV, T_h2o_ice_liq
+      use generic_cloud_common_h, only : Psat_generic
       use thermcell_mod, only: init_thermcell_mod
       use gases_h, only: gnom, gfrac
@@ -415 +416 @@
       real RH(ngrid,nlayer)   ! Relative humidity.
       real psat_tmp
+
+      real qsat_generic(ngrid,nlayer,nq) ! generic condensable tracers (GCS) specific concentration at saturation (kg/kg_of_air).
+      real RH_generic(ngrid,nlayer,nq)   ! generic condensable tracers (GCS) Relative humidity.
+      real rneb_generic(ngrid,nlayer,nq) ! GCS cloud fraction (generic condensation).
+      real psat_tmp_generic
+      real metallicity
       
       logical clearsky ! For double radiative transfer call. By BC
@@ -1559 +1566 @@
             call condensation_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,   &
                                           pt,pq,pdt,pdq,dtlscale_generic, &
-                                          dqvaplscale_generic,dqcldlscale_generic)
+                                          dqvaplscale_generic,dqcldlscale_generic,rneb_generic)
             pdt(1:ngrid,1:nlayer) = pdt(1:ngrid,1:nlayer)+dtlscale_generic(1:ngrid,1:nlayer)
             pdq(1:ngrid,1:nlayer,1:nq) = pdq(1:ngrid,1:nlayer,1:nq)+dqvaplscale_generic(1:ngrid,1:nlayer,1:nq)
@@ -1573 +1580 @@
                if (is_master) print*,'In generic condensation atmospheric energy change =',dEtot,' W m-2'
             end if
+
+            if (.not. water) then ! Compute GCS cloud fraction. For now we can not have both water cloud fraction and GCS cloud fraction 
+
+               ! Let's loop on tracers
+               do iq=1,nq
+
+                  call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,call_ice_vap_generic)
+                  
+                  if (call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
+                     do l = 1, nlayer
+                        do ig=1,ngrid
+                           cloudfrac(ig,l)=MAX(rneb_generic(ig,l,iq),rneb_man(ig,l)) 
+                           ! To be improved : for now, cloudfrac is only the cloudfrac of the last generic tracer
+                           ! Because it is rewritten every tracer
+                           ! One should create a cloudfrac_generic(ngrid,nlayer,nq)
+                        enddo
+                     enddo
+                  endif
+            end do ! do iq=1,nq loop on tracers 
+            endif ! .not. water
+
          endif !generic_condensation
 
@@ -2083 +2111 @@
       endif ! end of 'water'
 
+      ! Calculate RH_generic (Generic Relative Humidity) for diagnostic.
+      if(generic_condensation)then
+         do iq=1,nq
+
+            call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,call_ice_vap_generic)
+            
+            if (call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
+      
+               do l = 1, nlayer
+                  do ig=1,ngrid
+                     metallicity = 0.
+                     call Psat_generic(zt(ig,l),pplay(ig,l),metallicity,psat_tmp_generic,qsat_generic(ig,l,iq))
+                     RH_generic(ig,l,iq) = zq(ig,l,igcm_generic_vap) / qsat_generic(ig,l,iq)
+                  enddo
+               enddo
+
+            end if
+         
+         end do ! iq=1,nq
+
+      endif ! end of 'generic_condensation'
+
 
       if (is_master) print*,'--> Ls =',zls*180./pi
@@ -2423 +2473 @@
                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
+
+            if(generic_condensation)then
+               call writediagfi(ngrid,"rneb_generic","GCS cloud fraction (generic condensation)"," ",3,rneb_generic)
+               call writediagfi(ngrid,"CLF","GCS cloud fraction"," ",3,cloudfrac)
+               call writediagfi(ngrid,"RH_generic","GCS relative humidity"," ",3,RH_generic)
             endif