Changeset 3356

Timestamp:

Jun 4, 2024, 5:56:14 PM (7 months ago)

Author:

afalco

Message:

Pluto PCM:
sedimentation uses pphi. More precise than generic sedimentation?
molrad and visc adapted to Pluto.
AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callsedim.F (modified) (3 diffs)
• newsedim.F (modified) (5 diffs)
• newsedim_pluto.F (added)
• physiq_mod.F90 (modified) (1 diff)

trunk/LMDZ.PLUTO/libf/phypluto/callsedim.F

@@ -1 +1 @@
       SUBROUTINE callsedim(ngrid,nlay, ptimestep,
      &                pplev,zlev, pt, pdt,
-     &                pq, pdqfi, pdqsed,pdqs_sed,nq)
+     &                pq, pdqfi, pdqsed,pdqs_sed,nq,pphi)
 
       use radinc_h, only : naerkind
@@ -39 +39 @@
       real,intent(in):: pdt(ngrid,nlay) ! tendency on temperature
       real,intent(in):: zlev(ngrid,nlay+1)  ! altitude at layer boundaries
+      real,intent(in):: pphi(ngrid,nlay)      ! geopotential
       integer,intent(in) :: nq ! number of tracers
       real,intent(in) :: pq(ngrid,nlay,nq)  ! tracers (kg/kg)
@@ -120 +121 @@
              call newsedim(ngrid,nlay,1,ptimestep,
      &            pplev,masse,epaisseur,zt,radius(iq),rho_q(iq),
-     &            zqi(1,1,iq),wq,iq)
+     &            zqi(1,1,iq),wq,iq,pphi)
       !  endif
 

trunk/LMDZ.PLUTO/libf/phypluto/newsedim.F

@@ -1 +1 @@
       SUBROUTINE newsedim(ngrid,nlay,naersize,ptimestep,
-     &  pplev,masse,epaisseur,pt,rd,rho,pqi,wq,iq)
+     &  pplev,masse,epaisseur,pt,rd,rho,pqi,wq,iq,pphi)
       
       use ioipsl_getin_p_mod, only: getin_p
-      use comcstfi_mod, only: r, g
+      use comcstfi_mod, only: r, g, rad
       use gases_h
       ! use tracer_h, only : igcm_h2o_ice
@@ -41 +41 @@
       real,intent(out) :: wq(ngrid,nlay+1)  ! flux of tracer during timestep (?/m-2)
       integer,intent(in) :: iq ! tracer index
+      real,intent(in):: pphi(ngrid,nlay)      ! geopotential
 
 c   local:
@@ -90 +91 @@
          do igas=1, ngasmx
            if(gfrac(igas).ge.0.0) then
-             if(igas.eq.igas_N2) then
-               molrad = molrad + gfrac(igas)*2.2e-10                              ! N2 
-               visc(:,:) = visc(:,:) + gfrac(igas)*1.0e-5                         ! N2
+             if(igas.eq.igas_CO2) then
+               molrad = molrad + gfrac(igas)*2.2e-10                              ! CO2 
+               visc(:,:) = visc(:,:) + gfrac(igas)*1.0e-5                         ! CO2
              elseif(igas.eq.igas_N2) then
                molrad = molrad + gfrac(igas)*1.8e-10                              ! N2 (Kunze et al. 2022)
@@ -110 +111 @@
                visc(:,:) = visc(:,:) + gfrac(igas)*1.0e-5                         ! CH4 
              else
-               molrad = molrad + gfrac(igas)*2.2e-10                              ! N2 by default
-               visc(:,:) = visc(:,:) + gfrac(igas)*1.e-5                          ! N2 by default
+               molrad = molrad + gfrac(igas)*1.93e-10                              ! N2 by default
+               visc(:,:) = visc(:,:) + 6.67e-6                          ! N2 by default
                write(*,*) trim(gnom(igas))," is not included in"
      &              ," newsedim, N2 is used by default"
@@ -181 +182 @@
               rsurf=rfall
             ! endif
+
+          !  b = 2./9. * g
+            b = 2./9. * ((g*rad-pphi(ig,l))**2/(g*(rad**2))) ! AF24: from Pluto.old 
 
             vstokes(ig,l) = b / visc(ig,l) * rho * rfall**3 / rsurf *

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -1512 +1512 @@
                call callsedim(ngrid,nlayer,ptimestep,       &
                           pplev,zzlev,pt,pdt,pq,pdq,        &
-                          zdqsed,zdqssed,nq)
+                          zdqsed,zdqssed,nq,pphi)
             endif