Changeset 3954

Timestamp:

Nov 6, 2025, 4:06:36 PM (3 months ago)

Author:

slebonnois

Message:

SL: VCD 2.4 version. Updates in photochemistry (just tuning adjustment for O production), new molecular diffusion scheme, and adjustments of the non-orographic gravity waves (in particular change in direction of propagation at latitudes above 55 degrees).

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• moldiff_MPF.F90 (added)
• moldiff_red.old (moved) (moved from trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90)
• nonoro_gwd_ran_mod.F90 (modified) (5 diffs)
• photochemistry_venus.F90 (modified) (1 diff)
• physiq_mod.F (modified) (2 diffs)
• phytrac_chimie.F (modified) (2 diffs)

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

@@ -84 +84 @@
      ! 0.3.1 Gravity waves specifications
      INTEGER, parameter:: nk = 2 ! number of horizontal wavenumbers
-     INTEGER, parameter:: np = 2 ! directions (eastward and westward) phase speed
+     !----------------------------------------
+     ! VCD 2.0 tuning
+     ! INTEGER, parameter:: np = 2 ! directions (eastward and westward) phase speed
+     !----------------------------------------  
+     ! VCD 2.4 tuning
+     INTEGER, parameter:: np = 4 ! directions (east, north south, westward) phase speed
      INTEGER, parameter:: no = 2 ! absolute values of phase speed
      INTEGER, parameter:: na = 5 ! Number of realizations to get the phase speed
@@ -140 +145 @@
 !----------------------------------------
 ! VCD 2.0 tuning
-     REAL, parameter:: sat   = 0.6    ! saturation parameter
+!    REAL, parameter:: sat   = 0.6    ! saturation parameter
+!    REAL, parameter:: rdiss = 8.e-4  ! coefficient of dissipation
+!----------------------------------------
+! VCD 2.4 tuning
+     REAL, parameter:: sat   = 0.60   ! saturation parameter
      REAL, parameter:: rdiss = 8.e-4  ! coefficient of dissipation
 !----------------------------------------
@@ -168 +177 @@
      REAL href(nlayer + 1)           ! Reference altitude for launching altitude
 
-     REAL ran_num_1, ran_num_2, ran_num_3
+     REAL ran_num_1, ran_num_2, ran_num_3, ran_num_4
    
      LOGICAL, save :: firstcall = .true.
@@ -260 +269 @@
            ran_num_1 = mod(tt(ii, jw) * 10., 1.)
            ran_num_2 = mod(tt(ii, jw) * 100., 1.)
-
+           
 !! OPTIONS GENERIC DIFF VENUS !! 
+           !----------------------------------------
+           ! VCD 2.0
            ! angle (random) - reference 
+           !   zp(jw, ii) = (sign(1., 0.5 - ran_num_1) &
+           !            + 1.) * RPI / 2.
+           !
+           !
+           !!   zp(jw, ii) = atan2(4.*vh(ii, launch),uh(ii, launch)) &
+           !!            + (sign(1., 0.5 - ran_num_1) + 1.) * RPI / 2.
+           !! --------- TRY 00 -----------------------
+           !!IF((abs(latitude_deg(ii)) .le. 55.)) THEN
+           !!  zp(jw, ii) = 0. + 2. * ( ran_num_1 - 0.5 ) * RPI /180. &
+           !!                       * (10.) ! +/- 10 deg par rapport equateur
+           !!ELSE IF ((abs(latitude_deg(ii)) .le. 75.)) THEN
+           !!  zp(jw, ii) = 0. + 2. * ( ran_num_1 - 0.5 ) * RPI /180. &
+           !!                       * (10. +(90.-10.)                 &
+           !!                       * (latitude_deg(ii)-55.)/20.)
+           !!ELSE 
+           !!  zp(jw, ii) = ran_num_1 * RPI
+           !!ENDIF
+           !!zp(jw, ii) = mod(zp(jw, ii),2.*RPI)
+           !----------------------------------------
+           ! VCD 2.4
+           IF((abs(latitude_deg(ii)) .le. 43.125)) THEN ! EASTWARD WESTWARD
              zp(jw, ii) = (sign(1., 0.5 - ran_num_1) &
                         + 1.) * RPI / 2.
-           !   zp(jw, ii) = atan2(4.*vh(ii, launch),uh(ii, launch)) &
-           !            + (sign(1., 0.5 - ran_num_1) + 1.) * RPI / 2.
-           ! --------- TRY 00 -----------------------
-           !IF((abs(latitude_deg(ii)) .le. 55.)) THEN
-           !  zp(jw, ii) = 0. + 2. * ( ran_num_1 - 0.5 ) * RPI /180. &
-           !                       * (10.) ! +/- 10 deg par rapport equateur
-           !ELSE IF ((abs(latitude_deg(ii)) .le. 75.)) THEN
-           !  zp(jw, ii) = 0. + 2. * ( ran_num_1 - 0.5 ) * RPI /180. &
-           !                       * (10. +(90.-10.)                 &
-           !                       * (latitude_deg(ii)-55.)/20.)
-           !ELSE 
-           !  zp(jw, ii) = ran_num_1 * RPI
-           !ENDIF
-           !zp(jw, ii) = mod(zp(jw, ii),2.*RPI)
-           ! ------ TRY 01-------------------------------
-           !IF((abs(latitude_deg(ii)) .le. 55)) THEN
-           !  zp(jw, ii) = (sign(1., 0.5 - ran_num_1) &
-           !             + 1.) * RPI / 2.
-           !ELSE
-           !  zp(jw, ii) = ran_num_1 * RPI
-           !ENDIF
+           ELSE IF((abs(latitude_deg(ii)) .le. 56.25)) THEN ! TRANSITION BETWEEN ALL DIRECTION AND WEST/EASTWARD
+             zp(jw, ii) =              (0.5 - ran_num_1)            &
+                        * RPI      / (56.25 - 43.125)               &
+                        * (abs(latitude_deg(ii)) - 43.125)
+             IF (zp(jw, ii) < 0.) THEN
+               zp(jw, ii) = zp(jw, ii) + 2.*RPI
+             ENDIF
+           ELSE ! all direction if lat > 56.25 
+
+             zp(jw, ii) = ran_num_1 * RPI ! Pas besoin de plus comme cpha change l'orientation zp de 180 deg
+           
+           ENDIF
            ! ---- angle (0 or RPI) -----
            !!zp(jw, ii) = RPI*mod(jw,2)
@@ -308 +331 @@
               cpha = - 1. * cpha
               zp (jw, ii) = zp(jw, ii) + RPI
+           ENDIF
+           IF (zp(jw, ii) > 2.*RPI) THEN
+             zp(jw, ii) = zp(jw, ii) - 2.*RPI
            ENDIF
            !IF (cpha < 1.) THEN

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

@@ -4736 +4736 @@
       !--
       ! VCD 2.1 tuning
-          v_phot(65:nz,4) = v_phot(65:nz,4)*8. ! CO2 + hv ==> O(1D) + CO
+      !   v_phot(65:nz,4) = v_phot(65:nz,4)*8. ! CO2 + hv ==> O(1D) + CO
+      !--
+      ! VCD 2.4 tuning
+          v_phot(65:nz,4) = v_phot(65:nz,4)*7. ! CO2 + hv ==> O(1D) + CO
       !--
       ! TEST TUNEH

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

@@ -234 +234 @@
       REAL Fsedim(klon,klev+1)  ! Flux de sedimentation (kg.m-2)
 
+      EXTERNAL moldiff_MPF
 
 c
@@ -1670 +1671 @@
          IF (callthermos .and. ok_chem) THEN
 
-             call moldiff_red(klon, klev, nqmax,
+         !----------------------------------------
+         ! VCD 2.0
+
+!             call moldiff_red(klon, klev, nqmax,
+!     &                   pplay,paprs,t_seri, tr_seri, pdtphys,
+!     &                   d_t_euv,d_t_conduc,d_q_moldif)
+
+         !----------------------------------------
+         ! VCD 2.4
+
+             call moldiff_MPF(klon, klev, nqmax,
      &                   pplay,paprs,t_seri, tr_seri, pdtphys,
      &                   d_t_euv,d_t_conduc,d_q_moldif)
-
 
 ! --- update tendencies tracers ---

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

@@ -182 +182 @@
               if (i_so2 /= 0) then 
                  trac(:,25:,i_so2)  = 100.e-9
-                 trac(:,1:24,i_so2) = 10.e-6
+                 trac(:,1:24,i_so2) = 15.e-6
                  trac(:,25:,i_h2o)  = 1.e-6
                  trac(:,1:24,i_h2o) = 30.e-6
@@ -189 +189 @@
                  trac(:,:,i_osso_trans) = 0.
                  trac(:,:,i_s2o2_cyc)   = 0.
+                 trac(:,:,i_cl2so2)     = 0.
   
 ! AL TRACERS!!!!!!!!!!!!!!!!!!!!!!!!!!!!