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

Timestamp:

Jan 19, 2026, 2:31:03 PM (6 days ago)

Author:

debatzbr

Message:

Titan PCM: Minor merge and clean.
BBT

File:

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

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

@@ -220 +220 @@
       real zh(ngrid,nlayer)               ! Potential temperature (K).
       real pw(ngrid,nlayer)               ! Vertical velocity (m/s). (NOTE : >0 WHEN DOWNWARDS !!)
+      real omega(ngrid,nlayer)            ! omega velocity (Pa/s, >0 when downward)
 
       integer l,ig,ierr,iq,nw,isoil,ilat,lat_idx,i,j
@@ -814 +815 @@
          pw(:,l)=(pw(:,l)*r*pt(:,l)) / (pplay(:,l)*cell_area(:))
       enddo
+      ! omega in Pa/s
+      do l=1,nlayer-1
+         omega(1:ngrid,l)=0.5*(flxw(1:ngrid,l)+flxw(1:ngrid,l+1))
+       enddo
+       omega(1:ngrid,nlayer)=0.5*flxw(1:ngrid,nlayer) ! since flxw(nlayer+1)=0
+       do l=1,nlayer
+         omega(1:ngrid,l)=g*omega(1:ngrid,l)/cell_area(1:ngrid)
+       enddo
 
 !---------------------------------
@@ -1437 +1446 @@
       enddo
 
-      ! [Forcage de la photochimie pour les nuages]
+      ! [Temporary: adapted photochemistry for clouds]
       if (callclouds) then
          do ig = 1, ngrid
@@ -1444 +1453 @@
                !-------
                if(trim(nameOfTracer(gazs_indx(iq))) .eq. "C2H2") then
-                  pdq(ig,nlayer-3:,gazs_indx(iq)) = (4.0e-5 * rat_mmol(gazs_indx(iq)) - pq(ig,nlayer-3:,gazs_indx(iq))) / ptimestep
+                  pdq(ig,nlayer-3:,gazs_indx(iq)) = (3.0e-5 * rat_mmol(gazs_indx(iq)) - pq(ig,nlayer-3:,gazs_indx(iq))) / ptimestep
                endif
                ! C2H6 :
                !-------
                if(trim(nameOfTracer(gazs_indx(iq))) .eq. "C2H6") then
-                  pdq(ig,nlayer-3:,gazs_indx(iq)) = (8.0e-5 * rat_mmol(gazs_indx(iq)) - pq(ig,nlayer-3:,gazs_indx(iq))) / ptimestep
+                  pdq(ig,nlayer-3:,gazs_indx(iq)) = (6.0e-5 * rat_mmol(gazs_indx(iq)) - pq(ig,nlayer-3:,gazs_indx(iq))) / ptimestep
+               endif
+               ! AC6H6 :
+               !--------
+               if(trim(nameOfTracer(gazs_indx(iq))) .eq. "AC6H6") then
+                  pdq(ig,nlayer-5:,gazs_indx(iq)) = (5.0e-11 * rat_mmol(gazs_indx(iq)) - pq(ig,nlayer-5:,gazs_indx(iq))) / ptimestep
                endif
                ! HCN :
                !------
                if(trim(nameOfTracer(gazs_indx(iq))) .eq. "HCN") then
-                  pdq(ig,nlayer-3:,gazs_indx(iq)) = (2.0e-5 * rat_mmol(gazs_indx(iq)) - pq(ig,nlayer-3:,gazs_indx(iq))) / ptimestep
+                  pdq(ig,nlayer-3:,gazs_indx(iq)) = (5.0e-6 * rat_mmol(gazs_indx(iq)) - pq(ig,nlayer-3:,gazs_indx(iq))) / ptimestep
                endif
-               ! AC6H6 :
-               !--------
-               if(trim(nameOfTracer(gazs_indx(iq))) .eq. "AC6H6") then
-                  pdq(ig,nlayer-3:,gazs_indx(iq)) = (2.0e-5 * rat_mmol(gazs_indx(iq)) - pq(ig,nlayer-3:,gazs_indx(iq))) / ptimestep
+               ! HC3N :
+               !-------
+               if(trim(nameOfTracer(gazs_indx(iq))) .eq. "HC3N") then
+                  pdq(ig,nlayer-3:,gazs_indx(iq)) = (5.0e-8 * rat_mmol(gazs_indx(iq)) - pq(ig,nlayer-3:,gazs_indx(iq))) / ptimestep
                endif
             enddo
@@ -1737 +1751 @@
       CALL send_xios_field("v",zv)
       CALL send_xios_field("w",pw)
+      CALL send_xios_field("omega",omega)
 
       CALL send_xios_field("area",cell_area)
@@ -1964 +1979 @@
          CALL send_xios_field("evapCH4",dycevapCH4(:)) ! Pseudo-evaporation flux (mol/mol/s)
          CALL send_xios_field("tankCH4",tankCH4(:))    ! CH4 tank at the surface (m)
+         CALL send_xios_field("sphumCH4",zq(:,:,7+nmicro)/(1+zq(:,:,7+nmicro))) ! CH4 Specific Humidity (kg/kg)
 
          ! Atmosphere (3D) :