Changeset 2320 for trunk/LMDZ.VENUS

Timestamp:

May 12, 2020, 7:40:51 PM (5 years ago)

Author:

flefevre

Message:

1) correction de bug sur la photolyse de S2
2) ajustement de la valeur de SO2 en cas de reinitialisation de la chimie
3) passage des latitude et longitude vraies dans la chimie (utile pour certains tests)

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• photochemistry_venus.F90 (modified) (6 diffs)
• phytrac_chimie.F (modified) (2 diffs)

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

@@ -1 @@
-subroutine photochemistry_venus(nz, n_lon, ptimestep, p, t, tr, mumean, sza_input, nesp, iter)
+subroutine photochemistry_venus(nz, n_lon, ptimestep, p, t, tr, mumean, sza_input, lon, lat, nesp, iter)
 
 use chemparam_mod
@@ -17 +17 @@
 real :: ptimestep                 ! physics timestep (s)
 real :: sza_input                 ! solar zenith angle (degrees)
+real :: lon, lat                  ! longitude and latitude (degrees)
 
 integer :: n_lon                  ! for 1D test
@@ -176 +177 @@
 
    call define_dt(nesp, dt_corrected, dt_guess, ctimestep, cold(:), c(iz,:), &
-                  mat1, prod, loss, conc(iz))
+                  mat1, prod, loss, conc(iz), lon, lat)
 
    if (time + dt_corrected > ptimestep) then
@@ -2649 +2650 @@
 !      v_phot(:,nb_phot) = g023(:)*conc(:)
 
-!     Changement pour g023 ==> s2 + hv -> s + s
-!     Pas encore inclu dans la table jphot
-
-!     Pas de photodissociation sous le nuage photochimique 
-      g023(1:28)  = 0.0E+0
-
-!     Dependance en sza pour la photodissociation
-!     moins de W.m-2      
-      IF (sza_input.LT.90.0) THEN
-        g023(29:50) = 6.5E-3*COS(sza_input*pi/180.0)
-      ELSE
-        g023(29:50) = 0.0E+0
-      END IF
+!---  temporary treatment of s2 photolysis
+
+!     s2 + hv -> s + s
+
+!     under the clouds: no photolysis
+
+      g023(1:28)  = 0.
+
+!     above the clouds: value derived from Mills, 1998
+
+      if (sza_input < 90.) then
+        g023(29:nz) = 6.5e-3*cos(sza_input*pi/180.)
+      else
+        g023(29:nz) = 0.
+      end if
            
       nb_phot = nb_phot + 1
@@ -2980 +2983 @@
 
  subroutine define_dt(nesp, dtnew, dtold, ctimestep, cold, ccur, mat1, &
-                      prod, loss, dens)
+                      prod, loss, dens, lon, lat)
 
 !================================================================
@@ -3000 +3003 @@
 real, dimension(nesp)      :: prod, loss
 real                       :: dens
+real                       :: lon, lat
 
 ! output

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

@@ -83 +83 @@
                trac(:,1:22,i_co)  = 25.e-6
                trac(:,:,i_hcl)    = 0.4e-6
-               trac(:,1:22,i_so2) = 10.e-6
+               trac(:,1:22,i_so2) = 7.e-6
                trac(:,1:22,i_h2o) = 30.e-6
                trac(:,:,i_n2)     = 0.35e-1
@@ -257 +257 @@
      $                                ztrac(ilon,:,:),
      $                                mmean(ilon,:),
-     $                                sza_local, nqmax, iter(ilon,:))
+     $                                sza_local, 
+     $                                lon(ilon), lat(ilon),
+     $                                nqmax, iter(ilon,:))
 
          end do