Changeset 1956 for trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

Timestamp:

Jul 1, 2018, 5:07:48 PM (7 years ago)

Author:

jvatant

Message:

+ Add surface methane CH4 flux pseudo-evap diagnostic and include it in bottom layer zdqcond as >0.
+ minor cosmetic changes in gptitan.c
--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 (modified) (4 diffs)

trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

@@ -387 +387 @@
       real temp_eq(nlayer), press_eq(nlayer) ! Planetary averages for the init. of saturation profiles (K,mbar)
 
-      ! Surface methane tank
-      real,dimension(:),allocatable,save :: tankCH4 ! Depth of surface methane tank (m)
-!$OMP THREADPRIVATE(tankCH4)
+      ! Surface methane
+      real, dimension(:), allocatable, save :: tankCH4    ! Depth of surface methane tank (m)
+      real, dimension(:), allocatable, save :: zdyevapCH4 ! Surface pseudo-evaporation flux (chemistry keeping constant surface humidity) (mol/mol/s).
+!$OMP THREADPRIVATE(tankCH4,zdyevapCH4)
 
       ! -----******----- FOR MUPHYS OPTICS -----******----- 
@@ -522 +523 @@
             allocate(dycchi(ngrid,nlayer,nkim)) ! only for chemical tracers
             allocate(qysat(nlayer,nkim))
+            allocate(zdyevapCH4(ngrid))
             
             ! Chemistry timestep
@@ -1199 +1201 @@
                  ! Here we send zonal average fields ( corrected with cond ) from dynamics to chem. module
                  call calchim(ngrid,ychimbar,declin,ctimestep,ztfibar,zphibar,  &
-                              zplaybar,zplevbar,zzlaybar,zzlevbar,dycchi)
+                              zplaybar,zplevbar,zzlaybar,zzlevbar,dycchi,zdyevapCH4)
                else ! 3D chemistry (or 1D run)
                  call calchim(ngrid,ychim,declin,ctimestep,pt,pphi,  &
-                              pplay,pplev,zzlay,zzlev,dycchi)
+                              pplay,pplev,zzlay,zzlev,dycchi,zdyevapCH4)
                endif ! if moyzon
 
             endif
+            
+            ! Add diagnostic-only surface pseudo-evapoaration in condensation tendency for bottom layer
+            zdqcond(:,1,chimi_indx(7)) = zdyevapCH4(:)*rat_mmol(chimi_indx(7))
             
             do iq=1,nkim
@@ -1559 +1564 @@
              call writediagfi(ngrid,cnames(iq),cnames(iq),'mol/mol',3,zq(:,:,iq+nmicro)/rat_mmol(iq+nmicro))
            enddo
+           call writediagfi(ngrid,"fluxCH4","Surface CH4 pseudo-evaporation",'mol/mol/s',2,zdyevapCH4)
          endif