Changeset 1691 for trunk/LMDZ.VENUS

Timestamp:

Apr 10, 2017, 11:14:59 AM (8 years ago)

Author:

slebonnois

Message:

SL: corrections after testing of cloud microphysics in 1D

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• chemparam_mod.F90 (modified) (2 diffs)
• dyn1d/rcm1d.F (modified) (4 diffs)
• physiq_mod.F (modified) (1 diff)

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

@@ -88 +88 @@
 !     radius R_MEDIAN en m (donc *e-6 pour microns)
         
-        R_MEDIAN(:,:,:)=0.0E+0             ! Geometric Average Radius
+        R_MEDIAN(:,:,:)=0.0E+0           ! Geometric Average Radius
         STDDEV(:,:,:)=0.0E+0             ! Geometric Std Deviation
         K_MASS(:,:,:)=0.0E+0             ! Coeff multimodal
@@ -892 +892 @@
   real    :: trac(nlon,nlev,nqtot) ! traceur ( en vmr)
 
-  integer :: i
-  real    :: trac1d(nlev,2) ! traceur lu ( en vmr)
+!  integer :: i
+!  real    :: trac1d(nlev,2) ! traceur lu ( en vmr)
  
 ! lecture d'un fichier texte contenant les profils de trac1d(:1) = H2O et trac1d(:,2) = H2SO4
-
-  DO i=1,nlon
-     trac(i,:,i_h2o) = trac1d(:,1)
-     trac(i,:,i_h2so4) = trac1d(:,2)
-  ENDDO
-
+!  DO i=1,nlon
+!     trac(i,:,i_h2o) = trac1d(:,1)
+!     trac(i,:,i_h2so4) = trac1d(:,2)
+!  ENDDO
+
+!  intitialisation profils altitude H2O et H2SO4
+!  profil H2O initial vap+liq == que vap
+   trac(:,1:24,i_h2o) = 30.E-6 !
+   trac(:,25:50,i_h2o) = 1.E-6 !
+
+   trac(:,:,i_h2so4) = 3.E-9 ! Limite sup Sandor 2012
+   trac(:,23:50,i_h2so4) = 2.E-6 ! Profil H2SO4 initial => vap+liq
 
   END SUBROUTINE vapors4muphy_ini

trunk/LMDZ.VENUS/libf/phyvenus/dyn1d/rcm1d.F

@@ -194 +194 @@
       READ(unit,*) lati(1)
       PRINT *,lati(1)
+      lati(1)=lati(1)*pi/180.  ! must be in radians.
       long(1)=0.E+0
 
@@ -206 +207 @@
       cufi(1)=1.E+0
       cvfi(1)=1.E+0
+
+      call ini_cpdet
 
 c Ehouarn: iniphysiq requires arrays related to (3D) dynamics grid,
@@ -220 +223 @@
      &            rad,g,r,cpp,1)
 
-      call ini_cpdet
-
 c   le geopotentiel au sol est inutile en 1D car tout est controle
 c   par la pression de surface --->
@@ -348 +349 @@
       solsw(1)    = 0.
       sollw(1)    = 0.
-      fder(1)      = 0.
+      fder(1)     = 0.
       radsol(1)   = 0.
      

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

@@ -1870 +1870 @@
       CALL send_xios_field("fluxec",flux_ec)
 
-c plusieurs traceurs  !!!outputs in [vmr]
+c When using tracers
       IF (iflag_trac.eq.1) THEN
-         DO iq=1,nqmax
+c photochemical compounds  !!!outputs in [vmr]
+         DO iq=1,nqmax-nmicro
        CALL send_xios_field(tname(iq),qx(:,:,iq)*mmean(:,:)/M_tr(iq))
          ENDDO