Changeset 2958 for LMDZ5/trunk/libf/phylmd/cloudth_mod.F90

Timestamp:

Jul 18, 2017, 2:21:01 PM (7 years ago)

Author:

jbmadeleine

Message:

Added new outputs to cloudth_vert (level 10):
s_th, Thermal plume saturation deficit, kg/kg
s_env, Environment saturation deficit, kg/kg
sigma_th, Thermal plume gauss variance, kg/kg
sigma_env, Environment gauss variance, kg/kg

File:

• LMDZ5/trunk/libf/phylmd/cloudth_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/cloudth_mod.F90

@@ -819 +819 @@
 
       USE ioipsl_getin_p_mod, ONLY : getin_p
+      USE phys_output_var_mod, ONLY : cloudth_sth,cloudth_senv, &
+     &                                cloudth_sigmath,cloudth_sigmaenv
 
       IMPLICIT NONE
@@ -1145 +1147 @@
       alenv=(0.622*Lv*zqsatenv(ind1,ind2))/(rdd*zthl(ind1,ind2)**2)  
       aenv=1./(1.+(alenv*Lv/cppd))
-      senv=aenv*(po(ind1)-zqsatenv(ind1,ind2)) 
-      
+      senv=aenv*(po(ind1)-zqsatenv(ind1,ind2))
+      sth=0.
 
       sigma1s=ratqs(ind1,ind2)*zqenv(ind1)
+      sigma2s=0.
 
       xenv=senv/(sqrt2*sigma1s)
@@ -1166 +1169 @@
  
       endif       ! From the separation (thermal/envrionnement) et (environnement) only, l.335 et l.492
+      ! Outputs used to check the PDFs
+      cloudth_senv(ind1,ind2) = senv
+      cloudth_sth(ind1,ind2) = sth
+      cloudth_sigmaenv(ind1,ind2) = sigma1s
+      cloudth_sigmath(ind1,ind2) = sigma2s
+
       enddo       ! from the loop on ngrid l.333