Changeset 3421 for trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto

Timestamp:

Aug 26, 2024, 6:11:23 PM (17 months ago)

Author:

afalco

Message:

Pluto PCM: Added vertical mixing in no gcm, ch4/co surf, no_n2frost.
AF

File:

• trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/leapfrog_nogcm.F (modified) (8 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/leapfrog_nogcm.F

@@ -15 +15 @@
       USE guide_mod, ONLY : guide_main
       USE write_field, ONLY: writefield
+      USE callkeys_mod, only: tau_n2, tau_ch4, tau_co
       USE control_mod, ONLY: planet_type,nday,day_step,iperiod,iphysiq,
      &                       less1day,fractday,ndynstep,iconser,
@@ -109 +110 @@
 !       ED18 nogcm
       REAL tau_ps
-      REAL tau_n2
+      REAL tau_x
+      REAL tau_def
       REAL tau_teta
       REAL tetadpmean
@@ -199 +201 @@
 
       REAL psmean                    ! pression moyenne 
-      REAL pqn2mean                 ! moyenne globale ps*qco
+      REAL pqxmean                 ! moyenne globale ps*qco
       REAL p0                        ! pression de reference
       REAL p00d                        ! globalaverage(kpd)
-      REAL qmean_n2,qmean_n2_vert ! mass mean mixing ratio vap n2
-      REAL pqn2(ip1jmp1)           ! average n2 mass index : ps*q_n2
+      REAL qmean_x,qmean_x_vert ! mass mean mixing ratio vap n2
+      REAL pqx(ip1jmp1)           ! average n2 mass index : ps*q_n2
       REAL oldps(ip1jmp1)           ! saving old pressure ps to calculate qch4
 
@@ -243 +245 @@
       real :: n2mass(iip1,jjp1)
       real :: n2ice_ij(iip1,jjp1)
-      integer,save :: igcm_n2=0 ! index of CO2 tracer (if any)
+      integer,save :: igcm_n2=0 ! index of N2 tracer (if any)
+      integer,save :: igcm_ch4=0 ! index of CH4 tracer (if any)
+      integer,save :: igcm_co=0 ! index of CO tracer (if any)
       integer :: i,j,ig
       integer, parameter :: ngrid = 2+(jjm-1)*iim
@@ -339 +343 @@
       p00d=globaverage2d(kpd) ! mean pres at ref level
       tau_ps = 1. ! constante de rappel for pressure  (s) 
-      tau_n2 = 1.E5 !E5 ! constante de rappel for mix ratio qn2 (s) 
+      tau_n2 = 1 ! constante de rappel for mix ratio qn2 (s) 
+      tau_def = 1.E7 ! default constante de rappel for mix ratio qX (s) 
       tau_teta = 1.E7 !constante de rappel for potentiel temperature
 
-      PRINT*,'igcm_n2 = ',igcm_n2
       do iq=1,nqtot
         if (tname(iq)=="n2") then
           igcm_n2=iq
-          exit
+      !     exit
+        else if (tname(iq)=="ch4_gas") then
+            igcm_ch4=iq
+        else if (tname(iq)=="co_gas") then
+            igcm_co=iq
         endif
       enddo
+      PRINT*,'igcm_n2 = ',igcm_n2
+      PRINT*,'igcm_ch4 = ',igcm_ch4
+      PRINT*,'igcm_co = ',igcm_co
 
 c-----------------------------------------------------------------------
@@ -658 +669 @@
 !          ENDIF
 
-c          Mixing CO2 vertically
+c          Mixing N2 vertically ! not used for pluto ?
 c          --------------------------
-           if (igcm_n2.ne.0) then
-            DO ij=1,ip1jmp1
-               qmean_n2_vert=0.
-               DO l=1, llm 
-                 qmean_n2_vert= qmean_n2_vert
-     &          + q(ij,l,igcm_n2)*( p(ij,l) - p(ij,l+1))
-               END DO
-               qmean_n2_vert= qmean_n2_vert/ps(ij)
-               DO l=1, llm 
-                 q(ij,l,igcm_n2)= qmean_n2_vert
-               END DO
-            END DO
-           end if
+!            if (igcm_n2.ne.0) then
+!             DO ij=1,ip1jmp1
+!                qmean_x_vert=0.
+!                DO l=1, llm 
+!                  qmean_x_vert= qmean_x_vert
+!      &          + q(ij,l,igcm_n2)*( p(ij,l) - p(ij,l+1))
+!                END DO
+!                qmean_x_vert= qmean_x_vert/ps(ij)
+!                DO l=1, llm 
+!                  q(ij,l,igcm_n2)= qmean_x_vert
+!                END DO
+!             END DO
+!            end if
 
 
@@ -684 +695 @@
 c        ---------------------------------------------------------------------------------  
 
-         ! Simulate redistribution by dynamics for qn2 
-           if (igcm_n2.ne.0) then
+         ! Simulate redistribution by dynamics for qX
+         DO iq=1,nqtot
+           if ((iq.eq.igcm_n2).or.(iq.eq.igcm_ch4).or.
+     &         (iq.eq.igcm_co)) then
 
               DO ij=1,ip1jmp1
-                 pqn2(ij)= ps(ij) * q(ij,1,igcm_n2)
+                 pqx(ij)= ps(ij) * q(ij,1,iq)
               ENDDO
-              pqn2mean=globaverage2d(pqn2) 
-
-         !    Rappel newtonien vers qn2_mean
-              qmean_n2= pqn2mean / psmean
+              pqxmean=globaverage2d(pqx) 
+
+         !    Rappel newtonien vers qx_mean
+              qmean_x= pqxmean / psmean
+              
+              tau_x = tau_def
+              if (iq.eq.igcm_n2) then
+                  tau_x = tau_n2
+              else if (iq.eq.igcm_ch4) then
+                  tau_x = tau_ch4
+              else if (iq.eq.igcm_co) then
+                  tau_x = tau_co
+              end if
+              PRINT*,' tau_x ',iq,tau_x
 
               DO ij=1,ip1jmp1
-                  q(ij,1,igcm_n2)=q(ij,1,igcm_n2)+
-     &                  (qmean_n2-q(ij,1,igcm_n2))*
-     &                  (1.-exp(-dtphys/tau_n2))
+                  q(ij,1,iq)=q(ij,1,iq)+
+     &                  (qmean_x-q(ij,1,iq))*
+     &                  (1.-exp(-dtphys/tau_x))
               ENDDO
 
               DO l=2, llm 
                  DO ij=1,ip1jmp1
-                     q(ij,l,igcm_n2)=q(ij,1,igcm_n2)
+                     q(ij,l,iq)=q(ij,1,iq)
                  END DO
               END DO
@@ -709 +732 @@
 !             TEMPORAIRE (ED)
 !             PRINT*,'psmean = ',psmean
-!             PRINT*,'qmean_n2 = ',qmean_n2
-!             PRINT*,'pqn2mean = ',pqn2mean
-!             PRINT*,'q(50,1,igcm_n2) = ',q(50,1,igcm_n2)
-!             PRINT*,'q(50,2,igcm_n2) = ',q(50,2,igcm_n2)
-!             PRINT*,'q(50,3,igcm_n2) = ',q(50,3,igcm_n2)
+!             PRINT*,'qmean_x = ',qmean_x
+!             PRINT*,'pqxmean = ',pqxmean
+            ! PRINT*,' q(50,1) = ',iq,q(50,1,iq)
+            ! PRINT*,' q(50,2) = ',iq,q(50,2,iq)
+            ! PRINT*,' q(50,3) = ',iq,q(50,3,iq)
 
            endif ! igcm_n2.ne.0
-
-
-!       ********************************************************
+      enddo
+
+
+!       *****************************************s***************
 
 c        Horizontal mixing of pressure