Changeset 2523

Timestamp:

May 17, 2021, 12:16:36 PM (4 years ago)

Author:

slebonnois

Message:

SL: Venus GCM outputs for VCD + bug correction in moldiff_red

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• moldiff_red.F90 (modified) (1 diff)
• moldiffcoeff_red.F (modified) (3 diffs)
• physiq_mod.F (modified) (6 diffs)
• phytrac_chimie.F (modified) (1 diff)

trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90

@@ -176 +176 @@
         allocate(dij(ncompdiff,ncompdiff))
 
-        call moldiffcoeff_red(dij,indic_diff,gcmind,ncompdiff,nq)
+        call moldiffcoeff_red(dij,indic_diff,gcmind,ncompdiff)
         print*,'MOLDIFF  EXO'
 

trunk/LMDZ.VENUS/libf/phyvenus/moldiffcoeff_red.F

@@ -84 +84 @@
       real dnh
       logical,save :: firstcall=.true.
+!     logical,parameter :: outputcoeffs=.true. ! to output 'coeffs.dat' file,
       logical,parameter :: outputcoeffs=.false. ! to output 'coeffs.dat' file,
                                                 ! set outputcoeffs=.true.
@@ -282 +283 @@
       enddo 
       
+!!! TEST !!! Augmentation des coeffs CO2/CO et CO2/O
+!     do n=2,ncompdiff2
+!     if ((tname(gcmind(n)).eq.'co').or.(tname(gcmind(n)).eq.'o')) then
+!      dij(1,n)=dij(1,n)*1000.
+!      dij(n,1)=dij(1,n)
+!     endif
+!     enddo
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
       if (outputcoeffs) then
        ! output coefficients in 'coeffs.dat' file
@@ -287 +297 @@
        do n=1,ncompdiff2
         do nn=n,ncompdiff2
-          write(56,*) n,nn,dij(n,nn)    !*1.e5/1.381e-23/(273**1.75)
+          write(56,*) tname(gcmind(n)),tname(gcmind(nn)),dij(n,nn)      !*1.e5/1.381e-23/(273**1.75)
         enddo
        enddo

trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F

@@ -178 +178 @@
       integer,save :: itap        ! physics counter
       REAL delp(klon,klev)        ! epaisseur d'une couche
-      REAL omega(klon,klev)       ! vitesse verticale en Pa/s
+      REAL omega(klon,klev)       ! vitesse verticale en Pa/s (+ downward)
+      REAL vertwind(klon,klev)    ! vitesse verticale en m/s (+ upward)
       
       INTEGER igwd,idx(klon),itest(klon)
@@ -359 +360 @@
       REAL :: tr_hedin(klon,klev,nqmax)
       REAL :: d_tr(klon,klev,nqmax)
+c pour sorties
+      REAL :: col_dens_tr(klon,nqmax)
 
 c pour ioipsl
@@ -861 +864 @@
 c====================================================================
 
-c Calcule de vitesse verticale a partir de flux de masse verticale
-      DO k = 1, klev
-       DO i = 1, klon
-        omega(i,k) = RG*flxmw(i,k) / cell_area(i)
-       END DO
-      END DO
+! Compute vertical velocity (Pa/s) from vertical mass flux 
+! Need to linearly interpolate mass flux to mid-layers 
+      do k=1,klev-1
+         omega(1:klon,k) = 0.5*RG*(flxmw(1:klon,k)+flxmw(1:klon,k+1)) 
+     .                      / cell_area(1:klon)
+      enddo
+      omega(1:klon,klev) = 0.5*RG*flxmw(1:klon,klev) / cell_area(1:klon)
 
 c======
@@ -1930 +1934 @@
       endif
 
+c calcul vitesse verticale en m/s
+      DO k = 1, klev
+       DO i = 1, klon
+        vertwind(i,k) = -omega(i,k)/(rho(i,k)*RG)
+       END DO
+      END DO
+
 c------------------------
 c Calcul moment cinetique
@@ -2013 +2024 @@
       CALL send_xios_field("vitv",-1.*v_seri)
       CALL send_xios_field("vitw",omega)
+      CALL send_xios_field("vitwz",vertwind)
       CALL send_xios_field("Kz",ycoefh)
       CALL send_xios_field("mmean",mmean)
@@ -2059 +2071 @@
      $                           tr_seri(:,:,iq)*mmean(:,:)/m_tr(iq))
          end do
+
+! tracers in gas phase, column densities
+         do iq = 1,nqmax - nmicro
+            col_dens_tr(:,iq)=0.
+            do k = 1, klev
+              col_dens_tr(:,iq) = col_dens_tr(:,iq) + 
+     $              tr_seri(:,k,iq)* (paprs(:,k)-paprs(:,k+1)) / RG
+            end do
+            call send_xios_field("col_"//tname(iq),col_dens_tr(:,iq))
+         end do
+
 
 ! tracers in liquid phase, volume mixing ratio

trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F

@@ -68 +68 @@
       if (debutphy) then
       
+!!--- Adjustment of Helium amount
+!        if (i_he/=0) then
+!           trac(:,:,i_he)=trac(:,:,i_he)*20.
+!        endif
+!!---
+
 !-------------------------------------------------------------------
 !        case of tracers re-initialisation with chemistry