Changeset 2724 for trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90

Timestamp:

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

Author:

aslmd

Message:

rneb_generic and RH_generic added

File:

• trunk/LMDZ.GENERIC/libf/phystd/condensation_generic_mod.F90 (modified) (4 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