Changeset 2714

Timestamp:

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

Author:

aslmd

Message:

Print tendencies in diagfi.nc and small correction

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• condensation_generic_mod.F90 (modified) (1 diff)
• dyn1d/rcm1d.F (modified) (1 diff)
• generic_cloud_common_h.F90 (modified) (1 diff)
• phys_state_var_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (3 diffs)

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

@@ -21 +21 @@
 !     Authors
 !     -------
-!     Adapted from largescale.F90 by Lucas Teinturier (2022)
+!     Adapted from largescale.F90 by Lucas Teinturier & Noé Clément (2022)
 !     largescale.F90 adapted from the LMDTERRE code by R. Wordsworth (2009)
 !     Original author Z. X. Li (1993)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -687 +687 @@
                   q(:,iq)=0. !kg/kg of atmosphere
                   qsurf(iq) = 0. !kg/kg of atmosphere
-                  ! Look for a "profile_Fe_vap" input file*
+                  ! Look for a "profile_...._vap" input file
                   tracer_profile_file_name=""
                   tracer_profile_file_name='profile_'//txt

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

@@ -1 +1 @@
 module generic_cloud_common_h
+
+    !============================================================================
+    ! Purpose
+    ! -------
+    ! Set up relevant constants and parameters for generic condensable tracers cycles, 
+    ! and generic condensable tracers cloud properties
+    !
+    ! Authors
+    ! -------
+    ! Adapted from watercommon_h.F90
+    ! by Lucas Teinturier & Noé Clément (2022)
+    !============================================================================
+
     use comcstfi_mod, only: r, cpp, mugaz
     implicit none 

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

@@ -115 +115 @@
 !$OMP THREADPRIVATE(madjdE,madjdEz,lscaledE,lscaledEz)
       
+      real,dimension(:),allocatable,save :: genericconddE ! Heat from generic condensation (W.m-2)
+!$OMP THREADPRIVATE(genericconddE)
+
       real,dimension(:),allocatable,save :: H2Omaxcol ! Maximum possible H2O column amount (at 100% saturation) (kg.m-2).
 !$OMP THREADPRIVATE(H2Omaxcol)
@@ -202 +205 @@
         ALLOCATE(lscaledE(klon))
         ALLOCATE(lscaledEz(klon,klev))
+        ALLOCATE(genericconddE(klon))
         ALLOCATE(H2Omaxcol(klon))
         ! This is defined in comsaison_h
@@ -293 +297 @@
         DEALLOCATE(lscaledE)
         DEALLOCATE(lscaledEz)
+        DEALLOCATE(genericconddE)
         DEALLOCATE(H2Omaxcol)        
         DEALLOCATE(mu0)

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

@@ -1556 +1556 @@
             pdq(1:ngrid,1:nlayer,1:nq) = pdq(1:ngrid,1:nlayer,1:nq)+dqvaplscale_generic(1:ngrid,1:nlayer,1:nq)
             pdq(1:ngrid,1:nlayer,1:nq) = pdq(1:ngrid,1:nlayer,1:nq)+dqcldlscale_generic(1:ngrid,1:nlayer,1:nq)
+
+            if(enertest)then
+               do ig=1,ngrid
+                  genericconddE(ig) = cpp*SUM(mass(:,:)*dtlscale_generic(:,:))
+               enddo
+
+               call planetwide_sumval(cpp*massarea(:,:)*dtlscale_generic(:,:)/totarea_planet,dEtot)
+
+               if (is_master) print*,'In generic condensation atmospheric energy change =',dEtot,' W m-2'
+            end if
          endif !generic_condensation
 
@@ -1605 +1615 @@
 
          end if ! end of 'sedimentation'
+
+         write(*,*) 'end of sedimentation pdt', pdt
 
 
@@ -2266 +2278 @@
 !             call writediagfi(ngrid,"h2o_max_col","maximum H2O column amount","kg.m^-2",2,H2Omaxcol)
 
+         endif
+
+         if (generic_condensation) then
+
+            call writediagfi(ngrid,"genericconddE","heat from generic condensation","W m-2",2,genericconddE)
+         
          endif