Changeset 3323


Ignore:
Timestamp:
May 7, 2024, 9:39:36 AM (8 months ago)
Author:
flefevre
Message:

1) revised cloud microphysical parameters (this changes the results)
2) the number of modes is passed as an argument to cloud routines (this does not change the results)
3) some cosmetics to chemparam_mod.F90 (this does not change the results)

Location:
trunk/LMDZ.VENUS/libf/phyvenus
Files:
5 edited

Legend:

Unmodified
Added
Removed

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

    r2836 r3323  
    1 MODULE chemparam_mod
    2 
    3 !MODULE qui definit les indices des traceurs et leurs masses molaires.
    4 ! utilise aussi pour variables communes nuages/photochimie
    5 
    6 IMPLICIT NONE
    7 
    8 !---------------------------------------------------------------------------- 
    9 INTEGER, SAVE :: i_co2, i_co, i_h2, i_h2o, i_o1d,        &
     1!============================================================================
     2
     3module chemparam_mod
     4
     5!============================================================================
     6
     7! 1) cloud microphysical parameters for the simplified scheme (cl_scheme = 1)
     8! 2) indexes and molecular mass of chemical species
     9
     10implicit none
     11
     12!----------------------------------------------------------------------------
     13!     chemical tracers
     14!----------------------------------------------------------------------------
     15
     16integer, save :: i_co2, i_co, i_h2, i_h2o, i_o1d,        &
    1017                 i_o, i_o2, i_o2dg, i_o3, i_h,           &
    1118                 i_oh, i_ho2, i_h2o2, i_cl, i_clo,       &
    … …  
    1724                 i_co2plus, i_coplus, i_oplus, i_o2plus, &
    1825                 i_n2plus, i_hplus, i_h2oplus, i_nplus,  &
    19                  i_ohplus, i_cplus, i_noplus, i_h3oplus, & 
    20                  i_hcoplus, i_hco2plus, i_elec     
    21                  
    22 INTEGER, SAVE :: i_h2oliq, i_h2so4liq
    23 
    24 INTEGER, SAVE :: i_m0_aer, i_m3_aer,                       &
     26                 i_ohplus, i_cplus, i_noplus, i_h3oplus, &
     27                 i_hcoplus, i_hco2plus, i_elec
     28
     29integer, save :: i_h2oliq, i_h2so4liq
     30
     31integer, save :: i_m0_aer, i_m3_aer,                       &
    2532                 i_m0_mode1drop, i_m0_mode1ccn,            &
    26                  i_m3_mode1sa, i_m3_mode1w, i_m3_mode1ccn, &
    27                  i_m0_mode2drop, i_m0_mode2ccn,            &
    28                  i_m3_mode2sa, i_m3_mode2w, i_m3_mode2ccn
    29 
    30 INTEGER, SAVE :: nmicro  ! number of microphysical tracers
    31 
    32 REAL, DIMENSION(:), SAVE, ALLOCATABLE :: M_tr     ! Molecular Mass of tracers
    33 REAL, DIMENSION(:), SAVE, ALLOCATABLE :: type_tr  ! type of tracer
    34  
    35 REAL, DIMENSION(:,:),SAVE, ALLOCATABLE :: no_emission
    36 REAL, DIMENSION(:,:),SAVE, ALLOCATABLE :: o2_emission
     33                 i_m3_mode1sa, i_m3_mode1w, i_m3_mode1ccn, &
     34                 i_m0_mode2drop, i_m0_mode2ccn,            &
     35                 i_m3_mode2sa, i_m3_mode2w, i_m3_mode2ccn
     36
     37integer, save :: nmicro  ! number of species in the liquid phase
     38
     39real, dimension(:), save, allocatable :: m_tr           ! molecular mass of tracers
     40real, dimension(:), save, allocatable :: type_tr        ! type of tracer
     41
     42real, dimension(:,:), save, allocatable :: no_emission
     43real, dimension(:,:), save, allocatable :: o2_emission
     44
    3745!----------------------------------------------------------------------------
    38 ! DEF FOR CL_SCHEME = 1 (AURELIEN)
    39 
    40 !     number of clouds mode modelized
    41       INTEGER, PARAMETER :: nbr_mode = 3
    42       INTEGER :: i_cloud
    43       INTEGER, SAVE :: cloudmax
    44       INTEGER, SAVE :: cloudmin
    45       REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: R_MEDIAN
    46       REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: STDDEV
    47      
    48 !     K_MASS coefficient correspondant à la partie condensee de chaque mode
    49       REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: K_MASS
    50 
    51       REAL, SAVE, DIMENSION(:,:,:), ALLOCATABLE :: NBRTOT
    52       REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: WH2SO4
    53       REAL, SAVE, DIMENSION(:,:), ALLOCATABLE :: rho_droplet
     46!     cloud parameters
    5447!----------------------------------------------------------------------------
    55 ! DEF FOR CL_SCHEME = 2 (FULL MICROPHYS)
    56 
    57 !----------------------------------------------------------------------------
    58 
    59 !----------------------------------------------------------------------------
    60   CONTAINS
    61 !----------------------------------------------------------------------------
    62 
    63   SUBROUTINE cloud_ini(nbr_lon,nbr_lev)
    64 
     48
     49integer, save :: cloudmin, cloudmax
     50
     51!     qrad : ratio radius shell model of mode 3 (cimino, icarus, 1982)
     52!     if qrad = 0, fully liquid, if qrad = 1 fully solid
     53
     54real, save :: qrad
     55
     56!     median radius and standard deviation in each mode
     57
     58real, save, dimension(:,:,:), allocatable :: r_median, stddev
     59
     60!     k_mass : defines how the condensed phase is distributed in each mode.
     61!              sum of k_mass = 1
     62
     63real, save, dimension(:,:,:), allocatable :: k_mass
     64
     65real, save, dimension(:,:,:), allocatable :: nbrtot
     66real, save, dimension(:,:), allocatable :: wh2so4
     67real, save, dimension(:,:), allocatable :: rho_droplet
     68
     69contains
     70
     71!============================================================================
     72
     73subroutine cloud_ini(nbr_lon, nbr_lev, nbr_mode)
     74
     75!============================================================================
     76
     77!     sets cloud microphysical parameters for each mode:
     78!     radius, standard deviation, mass distribution
     79
     80integer :: nbr_lon, nbr_lev, nbr_mode
     81integer :: i_lev, ilon
     82
     83allocate(nbrtot(nbr_lon,nbr_lev,nbr_mode))
     84allocate(r_median(nbr_lon,nbr_lev,nbr_mode))
     85allocate(k_mass(nbr_lon,nbr_lev,nbr_mode))
     86allocate(stddev(nbr_lon,nbr_lev,nbr_mode))
     87allocate(wh2so4(nbr_lon,nbr_lev))
     88allocate(rho_droplet(nbr_lon,nbr_lev))
     89
     90!     initialisation
     91
     92r_median(:,:,:)  = 0.     ! median radius
     93stddev(:,:,:)    = 0.     ! geometric std deviation
     94k_mass(:,:,:)    = 0.     ! coeff mass multimodal
     95nbrtot(:,:,:)    = 0.
     96wh2so4(:,:)      = 0.
     97rho_droplet(:,:) = 0.
     98
     99!     minimum and maximum levels for the clouds
     100
     101cloudmin = 20   ! 20: 38 km
     102cloudmax = 50   ! 50: 95 km
     103
     104print*,'================================'
     105print*,'start initialisation cloud layer'
     106print*,'================================'
     107
     108!       ===============================================
     109!       knollenberg & hunten, 1980 and james et al 1997
     110!       ===============================================
     111!       initialisation unimodale
     112!       ===============================================
     113
     114!     lower haze: mode 1
     115!      do i_lev=cloudmin,20
     116!      r_median(:,i_lev,1)=0.2e-6
     117!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     118!      stddev(:,i_lev,1)=1.56
     119!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     120!      k_mass(:,i_lev,1)=1.0
     121!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     122!      end do
     123
     124!     lower cloud: mode 3
     125!      do i_lev=21,23
     126!      r_median(:,i_lev,1)=3.65e-6
     127!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     128!      stddev(:,i_lev,1)=1.28
     129!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     130!      k_mass(:,i_lev,1)=1.0
     131!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     132!      end do
     133
     134!     middle cloud: mode 2 prime
     135!      do i_lev=24,28
     136!      r_median(:,i_lev,1)=1.4e-6
     137!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     138!      stddev(:,i_lev,1)=1.23
     139!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     140!      k_mass(:,i_lev,1)=1.0
     141!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     142!      end do
     143
     144!     upper cloud: mode 2
     145!      do i_lev=29,35
     146!      r_median(:,i_lev,1)=1.0e-6
     147!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     148!      stddev(:,i_lev,1)=1.29
     149!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     150!      k_mass(:,i_lev,1)=1.0
     151!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     152!      end do
     153
     154!     upper haze: mode 1
     155!      do i_lev=36, cloudmax
     156!      r_median(:,i_lev,1)=0.2e-6
     157!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     158!      stddev(:,i_lev,1)=2.16
     159!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     160!      k_mass(:,i_lev,1)=1.0
     161!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     162!      end do
     163
     164!       ===============================================
     165!       initialisation trimodale
     166!       ===============================================
     167
     168!     lower haze: mode 1
     169!      do i_lev=cloudmin,20
     170!      r_median(:,i_lev,1)=0.3e-6
     171!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     172!      stddev(:,i_lev,1)=1.56
     173!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     174!      k_mass(:,i_lev,1)=1.0
     175!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     176!      end do
     177
     178!     lower haze: mode 2
     179   !   do i_lev=cloudmin,20
     180   !   r_median(:,i_lev,2)=1.4e-6
     181   !   print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     182   !   stddev(:,i_lev,2)=1.23
     183   !   print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     184   !   k_mass(:,i_lev,2)=0.0
     185   !   print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     186   !   end do
     187
     188!     lower haze: mode 3
     189!      do i_lev=cloudmin,20
     190!      r_median(:,i_lev,3)=3.65e-6
     191!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     192!      stddev(:,i_lev,3)=1.28
     193!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     194!      k_mass(:,i_lev,3)=0.
     195!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     196!      end do
     197
     198!     lower cloud: mode 1
     199!      do i_lev=21,23
     200!      r_median(:,i_lev,1)=0.3e-6
     201!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     202!      stddev(:,i_lev,1)=1.56
     203!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     204!      k_mass(:,i_lev,1)=0.1
     205!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     206!      end do
     207
     208!     lower cloud: mode 2 prime
     209!      do i_lev=21,23
     210!      r_median(:,i_lev,2)=1.4e-6
     211!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     212!      stddev(:,i_lev,2)=1.23
     213!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     214!      k_mass(:,i_lev,2)=0.4
     215!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     216!      end do
     217
     218!     lower cloud: mode 3
     219!      do i_lev=21,23
     220!      r_median(:,i_lev,3)=3.65e-6
     221!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     222!      stddev(:,i_lev,3)=1.28
     223!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     224!      k_mass(:,i_lev,3)=0.5
     225!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     226!      end do
     227
     228!     middle cloud: mode 1
     229!      do i_lev=24,28
     230!      r_median(:,i_lev,1)=0.3e-6
     231!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     232!      stddev(:,i_lev,1)=1.56
     233!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     234!      k_mass(:,i_lev,1)=0.0
     235!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     236!      end do
     237
     238!     middle cloud: mode 2 prime
     239!      do i_lev=24,28
     240!      r_median(:,i_lev,2)=1.4e-6
     241!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     242!      stddev(:,i_lev,2)=1.23
     243!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     244!      k_mass(:,i_lev,2)=0.8
     245!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     246!      end do
     247
     248!     middle cloud: mode 3
     249!      do i_lev=24,28
     250!      r_median(:,i_lev,3)=3.65e-6
     251!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     252!      stddev(:,i_lev,3)=1.28
     253!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     254!      k_mass(:,i_lev,3)=0.2
     255!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     256!      end do
     257
     258
     259!     upper cloud: mode 1
     260!      do i_lev=29,35
     261!      r_median(:,i_lev,1)=0.3e-6
     262!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     263!      stddev(:,i_lev,1)=1.56
     264!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     265!      k_mass(:,i_lev,1)=0.15
     266!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     267!      end do
     268
     269!     upper cloud: mode 2
     270!      do i_lev=29,35
     271!      r_median(:,i_lev,2)=1.0e-6
     272!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     273!      stddev(:,i_lev,2)=1.29
     274!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     275!      k_mass(:,i_lev,2)=0.85
     276!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     277!      end do
     278
     279!     upper cloud: mode 3
     280!      do i_lev=29,35
     281!      r_median(:,i_lev,3)=3.65e-6
     282!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     283!      stddev(:,i_lev,3)=1.28
     284!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     285!      k_mass(:,i_lev,3)=0.0
     286!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     287!      end do
     288
     289!     upper haze: mode 1
     290!      do i_lev=36, cloudmax
     291!      r_median(:,i_lev,1)=0.3e-6
     292!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     293!      stddev(:,i_lev,1)=1.56
     294!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     295!      k_mass(:,i_lev,1)=1.0
     296!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     297!      end do
     298
     299!     upper haze: mode 2
     300!      do i_lev=36, cloudmax
     301!      r_median(:,i_lev,2)=1.e-6
     302!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     303!      stddev(:,i_lev,2)=1.29
     304!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     305!      k_mass(:,i_lev,2)=0.0
     306!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     307!      end do
     308
     309!     upper haze: mode 3
     310!      do i_lev=36, cloudmax
     311!      r_median(:,i_lev,3)=3.65e-6
     312!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     313!      stddev(:,i_lev,3)=2.16
     314!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     315!      k_mass(:,i_lev,3)=0.0
     316!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     317!      end do
    65318!=============================================================
    66 !       cloud_ini definit le champ 3D des caracteristiques du nuage
    67  
    68       INTEGER :: nbr_lon,nbr_lev,i_lev
    69    
    70       ALLOCATE(NBRTOT(nbr_lon,nbr_lev,nbr_mode))
    71       ALLOCATE(R_MEDIAN(nbr_lon,nbr_lev,nbr_mode))
    72       ALLOCATE(K_MASS(nbr_lon,nbr_lev,nbr_mode))
    73       ALLOCATE(STDDEV(nbr_lon,nbr_lev,nbr_mode))
    74       ALLOCATE(WH2SO4(nbr_lon,nbr_lev))
    75       ALLOCATE(rho_droplet(nbr_lon,nbr_lev))
    76            
    77       PRINT*,'=========================='
    78       PRINT*,'Initialisation cloud layer'
    79       PRINT*,'=========================='
    80       PRINT*,'nbr_lon',nbr_lon
    81       PRINT*,'nbr_lev',nbr_lev
    82       PRINT*,'nbr_mode',nbr_mode
    83        
    84       NBRTOT(:,:,:)    = 0.0E+0
    85       WH2SO4(:,:)      = 0.0E+0
    86       rho_droplet(:,:) = 0.0E+0
    87            
    88 !=============================================================
    89 !                     Initialisation cloud layer 1
    90 !=============================================================
    91 !     cloudmin et cloudmax niveaux du GCM
    92       cloudmin= 18
    93       cloudmax= 50
    94 
    95 !     radius R_MEDIAN en m (donc *e-6 pour microns)
    96        
    97         R_MEDIAN(:,:,:)=0.0E+0           ! Geometric Average Radius
    98         STDDEV(:,:,:)=0.0E+0             ! Geometric Std Deviation
    99         K_MASS(:,:,:)=0.0E+0             ! Coeff multimodal
    100319
    101320!       ===============================================
    102 !       Knollenberg & Hunten, 1980 and James et al 1997
     321!       initialisation trimodale knollenberg
    103322!       ===============================================
    104323
    105 !       ===============================================
    106 !       Initialisation UNIMODALE
    107 !       ===============================================
    108 
    109 !     Lower Haze: mode 1
    110 !      DO i_lev=cloudmin,20
    111 !      R_MEDIAN(:,i_lev,1)=0.2e-6
    112 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    113 !      STDDEV(:,i_lev,1)=1.56
    114 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    115 !      K_MASS(:,i_lev,1)=1.0
    116 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    117 !      END DO
    118 
    119 !     Lower Cloud: mode 3
    120 !      DO i_lev=21,23
    121 !      R_MEDIAN(:,i_lev,1)=3.65e-6
    122 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    123 !      STDDEV(:,i_lev,1)=1.28
    124 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    125 !      K_MASS(:,i_lev,1)=1.0
    126 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    127 !      END DO
    128 
    129 !     Middle Cloud: mode 2 prime
    130 !      DO i_lev=24,28
    131 !      R_MEDIAN(:,i_lev,1)=1.4e-6
    132 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    133 !      STDDEV(:,i_lev,1)=1.23
    134 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    135 !      K_MASS(:,i_lev,1)=1.0
    136 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    137 !      END DO
    138 
    139 !     Upper Cloud: mode 2
    140 !      DO i_lev=29,35
    141 !      R_MEDIAN(:,i_lev,1)=1.0e-6
    142 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    143 !      STDDEV(:,i_lev,1)=1.29
    144 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    145 !      K_MASS(:,i_lev,1)=1.0
    146 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    147 !      END DO
    148 
    149 !     Upper Haze: mode 1
    150 !      DO i_lev=36, cloudmax
    151 !      R_MEDIAN(:,i_lev,1)=0.2e-6
    152 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    153 !      STDDEV(:,i_lev,1)=2.16
    154 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    155 !      K_MASS(:,i_lev,1)=1.0
    156 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    157 !      END DO
    158 
    159 !       ===============================================
    160 !       Initialisation TRIMODALE
    161 !       ===============================================
    162 
    163 !     Lower Haze: mode 1
    164 !      DO i_lev=cloudmin,20
    165 !      R_MEDIAN(:,i_lev,1)=0.3e-6
    166 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    167 !      STDDEV(:,i_lev,1)=1.56
    168 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    169 !      K_MASS(:,i_lev,1)=1.0
    170 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    171 !      END DO
    172                              
    173 !     Lower Haze: mode 2
    174 !      DO i_lev=cloudmin,20
    175 !      R_MEDIAN(:,i_lev,2)=1.4e-6
    176 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    177 !      STDDEV(:,i_lev,2)=1.23
    178 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    179 !      K_MASS(:,i_lev,2)=0.0
    180 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    181 !      END DO
    182 
    183 !     Lower Haze: mode 3
    184 !      DO i_lev=cloudmin,20
    185 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    186 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    187 !      STDDEV(:,i_lev,3)=1.28
    188 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    189 !      K_MASS(:,i_lev,3)=0.
    190 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    191 !      END DO 
    192 
    193 !     Lower Cloud: mode 1
    194 !      DO i_lev=21,23
    195 !      R_MEDIAN(:,i_lev,1)=0.3e-6
    196 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    197 !      STDDEV(:,i_lev,1)=1.56
    198 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    199 !      K_MASS(:,i_lev,1)=0.1
    200 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    201 !      END DO
    202                              
    203 !     Lower Cloud: mode 2 prime
    204 !      DO i_lev=21,23
    205 !      R_MEDIAN(:,i_lev,2)=1.4e-6
    206 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    207 !      STDDEV(:,i_lev,2)=1.23
    208 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    209 !      K_MASS(:,i_lev,2)=0.4
    210 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    211 !      END DO
    212 
    213 !     Lower Cloud: mode 3
    214 !      DO i_lev=21,23
    215 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    216 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    217 !      STDDEV(:,i_lev,3)=1.28
    218 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    219 !      K_MASS(:,i_lev,3)=0.5
    220 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    221 !      END DO
    222 
    223 !     Middle Cloud: mode 1
    224 !      DO i_lev=24,28
    225 !      R_MEDIAN(:,i_lev,1)=0.3e-6
    226 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    227 !      STDDEV(:,i_lev,1)=1.56
    228 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    229 !      K_MASS(:,i_lev,1)=0.0
    230 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    231 !      END DO
    232          
    233 !     Middle Cloud: mode 2 prime
    234 !      DO i_lev=24,28
    235 !      R_MEDIAN(:,i_lev,2)=1.4e-6
    236 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    237 !      STDDEV(:,i_lev,2)=1.23
    238 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    239 !      K_MASS(:,i_lev,2)=0.8
    240 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    241 !      END DO
    242    
    243 !     Middle Cloud: mode 3
    244 !      DO i_lev=24,28
    245 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    246 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    247 !      STDDEV(:,i_lev,3)=1.28
    248 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    249 !      K_MASS(:,i_lev,3)=0.2
    250 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    251 !      END DO
    252 
    253 
    254 !     Upper Cloud: mode 1
    255 !      DO i_lev=29,35
    256 !      R_MEDIAN(:,i_lev,1)=0.3e-6
    257 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    258 !      STDDEV(:,i_lev,1)=1.56
    259 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    260 !      K_MASS(:,i_lev,1)=0.15
    261 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    262 !      END DO
    263 
    264 !     Upper Cloud: mode 2
    265 !      DO i_lev=29,35
    266 !      R_MEDIAN(:,i_lev,2)=1.0e-6
    267 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    268 !      STDDEV(:,i_lev,2)=1.29
    269 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    270 !      K_MASS(:,i_lev,2)=0.85
    271 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    272 !      END DO
    273        
    274 !     Upper Cloud: mode 3
    275 !      DO i_lev=29,35
    276 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    277 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    278 !      STDDEV(:,i_lev,3)=1.28
    279 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    280 !      K_MASS(:,i_lev,3)=0.0
    281 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    282 !      END DO
    283 
    284 !     Upper Haze: mode 1
    285 !      DO i_lev=36, cloudmax
    286 !      R_MEDIAN(:,i_lev,1)=0.3e-6
    287 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    288 !      STDDEV(:,i_lev,1)=1.56
    289 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    290 !      K_MASS(:,i_lev,1)=1.0
    291 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    292 !      END DO
    293 
    294 !     Upper Haze: mode 2
    295 !      DO i_lev=36, cloudmax
    296 !      R_MEDIAN(:,i_lev,2)=1.e-6
    297 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    298 !      STDDEV(:,i_lev,2)=1.29
    299 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    300 !      K_MASS(:,i_lev,2)=0.0
    301 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    302 !      END DO
    303      
    304 !     Upper Haze: mode 3
    305 !      DO i_lev=36, cloudmax
    306 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    307 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    308 !      STDDEV(:,i_lev,3)=2.16
    309 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    310 !      K_MASS(:,i_lev,3)=0.0
    311 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    312 !      END DO       
    313 !=============================================================
    314 
    315 !       ===============================================
    316 !       Initialisation TRIMODALE Knollenberg
    317 !       ===============================================
    318 
    319 !     Lower Haze: mode 1
    320       DO i_lev=cloudmin,22
    321       R_MEDIAN(:,i_lev,1)=0.1e-6
    322       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    323       STDDEV(:,i_lev,1)=1.57
    324       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    325       K_MASS(:,i_lev,1)=1.0
    326       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    327       END DO
    328                              
    329 !     Lower Haze: mode 2
    330       DO i_lev=cloudmin,22
    331       R_MEDIAN(:,i_lev,2)=1.4e-6
    332       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    333       STDDEV(:,i_lev,2)=1.23
    334       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    335       K_MASS(:,i_lev,2)=0.0
    336       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    337       END DO
    338 
    339 !     Lower Haze: mode 3
    340       DO i_lev=cloudmin,22
    341       R_MEDIAN(:,i_lev,3)=3.65e-6
    342       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    343       STDDEV(:,i_lev,3)=1.28
    344       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    345       K_MASS(:,i_lev,3)=0.0
    346       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    347       END DO 
    348 
    349 !     Pre Cloud: mode 1
    350       DO i_lev=23,23
    351       R_MEDIAN(:,i_lev,1)=0.15e-6
    352       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    353       STDDEV(:,i_lev,1)=1.8
    354       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    355       K_MASS(:,i_lev,1)=0.04
    356       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    357       END DO
    358                              
    359 !     Pre Cloud: mode 2
    360       DO i_lev=23,23
    361       R_MEDIAN(:,i_lev,2)=1.0e-6
    362       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    363       STDDEV(:,i_lev,2)=1.29
    364       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    365       K_MASS(:,i_lev,2)=0.96
    366       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    367       END DO
    368 
    369 !     Pre Cloud: mode 3
    370       DO i_lev=23,23
    371       R_MEDIAN(:,i_lev,3)=3.65e-6
    372       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    373       STDDEV(:,i_lev,3)=1.28
    374       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    375       K_MASS(:,i_lev,3)=0.0
    376       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    377       END DO
    378                  
    379 !      Lower Cloud: mode 1
    380       DO i_lev=24,24
    381       R_MEDIAN(:,i_lev,1)=0.2e-6
    382       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    383       STDDEV(:,i_lev,1)=1.8
    384       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    385       K_MASS(:,i_lev,1)=0.014
    386       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    387       END DO
    388                
    389 !     Lower Cloud: mode 2
    390       DO i_lev=24,24
    391       R_MEDIAN(:,i_lev,2)=1.0e-6
    392       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    393       STDDEV(:,i_lev,2)=1.29
    394       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    395       K_MASS(:,i_lev,2)=0.02
    396       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    397       END DO
    398 
    399 !     Lower Cloud: mode 3
    400       DO i_lev=24,24
    401       R_MEDIAN(:,i_lev,3)=3.65e-6
    402       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    403       STDDEV(:,i_lev,3)=1.28
    404       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    405       K_MASS(:,i_lev,3)=0.966
    406       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    407       END DO
    408 
    409 !     Middle Cloud: mode 1
    410       DO i_lev=25,28
    411       R_MEDIAN(:,i_lev,1)=0.15e-6
    412       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    413       STDDEV(:,i_lev,1)=1.9
    414       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    415       K_MASS(:,i_lev,1)=0.0084
    416       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    417       END DO
    418          
    419 !     Middle Cloud: mode 2 prime
    420       DO i_lev=25,28
    421       R_MEDIAN(:,i_lev,2)=1.4e-6
    422       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    423       STDDEV(:,i_lev,2)=1.23
    424       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    425       K_MASS(:,i_lev,2)=0.21
    426       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    427       END DO
    428    
    429 !     Middle Cloud: mode 3
    430       DO i_lev=25,28
    431       R_MEDIAN(:,i_lev,3)=3.65e-6
    432       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    433       STDDEV(:,i_lev,3)=1.28
    434       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    435       K_MASS(:,i_lev,3)=0.7816
    436       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    437       END DO
     324!     lower haze: mode 1
     325      ! do i_lev=cloudmin,22
     326      ! r_median(:,i_lev,1)=0.1e-6
     327      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     328      ! stddev(:,i_lev,1)=1.57
     329      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     330      ! k_mass(:,i_lev,1)=1.0
     331      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     332      ! end do
     333
     334!     lower haze: mode 2
     335      ! do i_lev=cloudmin,22
     336      ! r_median(:,i_lev,2)=1.4e-6
     337      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     338      ! stddev(:,i_lev,2)=1.23
     339      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     340      ! k_mass(:,i_lev,2)=0.0
     341      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     342      ! end do
     343
     344!     lower haze: mode 3
     345      ! do i_lev=cloudmin,22
     346      ! r_median(:,i_lev,3)=3.65e-6
     347      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     348      ! stddev(:,i_lev,3)=1.28
     349      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     350      ! k_mass(:,i_lev,3)=0.0
     351      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     352      ! end do
     353
     354!     pre cloud: mode 1
     355      ! do i_lev=23,23
     356      ! r_median(:,i_lev,1)=0.15e-6
     357      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     358      ! stddev(:,i_lev,1)=1.8
     359      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     360      ! k_mass(:,i_lev,1)=0.04
     361      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     362      ! end do
     363
     364!     pre cloud: mode 2
     365      ! do i_lev=23,23
     366      ! r_median(:,i_lev,2)=1.0e-6
     367      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     368      ! stddev(:,i_lev,2)=1.29
     369      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     370      ! k_mass(:,i_lev,2)=0.96
     371      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     372      ! end do
     373
     374!     pre cloud: mode 3
     375      ! do i_lev=23,23
     376      ! r_median(:,i_lev,3)=3.65e-6
     377      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     378      ! stddev(:,i_lev,3)=1.28
     379      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     380      ! k_mass(:,i_lev,3)=0.0
     381      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     382      ! end do
     383
     384!      lower cloud: mode 1
     385      ! do i_lev=24,24
     386      ! r_median(:,i_lev,1)=0.2e-6
     387      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     388      ! stddev(:,i_lev,1)=1.8
     389      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     390      ! k_mass(:,i_lev,1)=0.014
     391      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     392      ! end do
     393
     394!     lower cloud: mode 2
     395      ! do i_lev=24,24
     396      ! r_median(:,i_lev,2)=1.0e-6
     397      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     398      ! stddev(:,i_lev,2)=1.29
     399      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     400      ! k_mass(:,i_lev,2)=0.02
     401      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     402      ! end do
     403
     404!     lower cloud: mode 3
     405      ! do i_lev=24,24
     406      ! r_median(:,i_lev,3)=3.65e-6
     407      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     408      ! stddev(:,i_lev,3)=1.28
     409      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     410      ! k_mass(:,i_lev,3)=0.966
     411      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     412      ! end do
     413
     414!     middle cloud: mode 1
     415      ! do i_lev=25,28
     416      ! r_median(:,i_lev,1)=0.15e-6
     417      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     418      ! stddev(:,i_lev,1)=1.9
     419      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     420      ! k_mass(:,i_lev,1)=0.0084
     421      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     422      ! end do
     423
     424!     middle cloud: mode 2 prime
     425      ! do i_lev=25,28
     426      ! r_median(:,i_lev,2)=1.4e-6
     427      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     428      ! stddev(:,i_lev,2)=1.23
     429      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     430      ! k_mass(:,i_lev,2)=0.21
     431      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     432      ! end do
     433
     434!     middle cloud: mode 3
     435      ! do i_lev=25,28
     436      ! r_median(:,i_lev,3)=3.65e-6
     437      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     438      ! stddev(:,i_lev,3)=1.28
     439      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     440      ! k_mass(:,i_lev,3)=0.7816
     441      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     442      ! end do
    438443
    439444!      option: upper haze remplacee par extension upper cloud
    … …  
    441446!       ===============================================
    442447
    443 !     Upper Cloud: mode 1
    444       DO i_lev=29,35 !cloudmax
    445       R_MEDIAN(:,i_lev,1)=0.2e-6
    446       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    447       STDDEV(:,i_lev,1)=2.16
    448       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    449       K_MASS(:,i_lev,1)=0.72
    450       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    451       END DO
    452 
    453 !     Upper Cloud: mode 2
    454       DO i_lev=29,35 !cloudmax
    455       R_MEDIAN(:,i_lev,2)=1.0e-6
    456       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    457       STDDEV(:,i_lev,2)=1.29
    458       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    459       K_MASS(:,i_lev,2)=0.28
    460       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    461       END DO
    462        
    463 !     Upper Cloud: mode 3
    464       DO i_lev=29,35 !cloudmax
    465       R_MEDIAN(:,i_lev,3)=3.65e-6
    466       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    467       STDDEV(:,i_lev,3)=1.28
    468       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    469       K_MASS(:,i_lev,3)=0.0
    470       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    471       END DO
    472 
    473 !     Upper Haze: mode 1
    474       DO i_lev=36, cloudmax
    475       R_MEDIAN(:,i_lev,1)=0.2e-6
    476       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    477       STDDEV(:,i_lev,1)=2.16
    478       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    479       K_MASS(:,i_lev,1)=1.0
    480       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    481       END DO
    482 
    483 !     Upper Haze: mode 2
    484       DO i_lev=36, cloudmax
    485       R_MEDIAN(:,i_lev,2)=1.e-6
    486       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    487       STDDEV(:,i_lev,2)=1.29
    488       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    489       K_MASS(:,i_lev,2)=0.0
    490       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    491       END DO
    492      
    493 !     Upper Haze: mode 3
    494       DO i_lev=36, cloudmax
    495       R_MEDIAN(:,i_lev,3)=3.65e-6
    496       PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    497       STDDEV(:,i_lev,3)=2.16
    498       PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    499       K_MASS(:,i_lev,3)=0.0
    500       PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    501       END DO       
     448!     upper cloud: mode 1
     449      ! do i_lev=29,35 !cloudmax
     450      ! r_median(:,i_lev,1)=0.2e-6
     451      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     452      ! stddev(:,i_lev,1)=2.16
     453      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     454      ! k_mass(:,i_lev,1)=0.72
     455      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     456      ! end do
     457
     458!     upper cloud: mode 2
     459      ! do i_lev=29,35 !cloudmax
     460      ! r_median(:,i_lev,2)=1.0e-6
     461      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     462      ! stddev(:,i_lev,2)=1.29
     463      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     464      ! k_mass(:,i_lev,2)=0.28
     465      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     466      ! end do
     467
     468!     upper cloud: mode 3
     469      ! do i_lev=29,35 !cloudmax
     470      ! r_median(:,i_lev,3)=3.65e-6
     471      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     472      ! stddev(:,i_lev,3)=1.28
     473      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     474      ! k_mass(:,i_lev,3)=0.0
     475      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     476      ! end do
     477
     478!     upper haze: mode 1
     479      ! do i_lev=36, cloudmax
     480      ! r_median(:,i_lev,1)=0.2e-6
     481      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     482      ! stddev(:,i_lev,1)=2.16
     483      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     484      ! k_mass(:,i_lev,1)=1.0
     485      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     486      ! end do
     487
     488!     upper haze: mode 2
     489      ! do i_lev=36, cloudmax
     490      ! r_median(:,i_lev,2)=1.e-6
     491      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     492      ! stddev(:,i_lev,2)=1.29
     493      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     494      ! k_mass(:,i_lev,2)=0.0
     495      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     496      ! end do
     497
     498!     upper haze: mode 3
     499      ! do i_lev=36, cloudmax
     500      ! r_median(:,i_lev,3)=3.65e-6
     501      ! print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     502      ! stddev(:,i_lev,3)=2.16
     503      ! print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     504      ! k_mass(:,i_lev,3)=0.0
     505      ! print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     506      ! end do
    502507
    503508!=============================================================
    504509
    505510!       ===============================================================
    506 !       Initialisation TRIMODALE "Knollenberg" sans Mode3, Mode2 etendu
     511!       initialisation trimodale "knollenberg" sans mode3, mode2 etendu
    507512!       ===============================================================
    508513
    509 !     Lower Haze: mode 1
    510 !      DO i_lev=cloudmin,22
    511 !      R_MEDIAN(:,i_lev,1)=0.1e-6
    512 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    513 !      STDDEV(:,i_lev,1)=1.57
    514 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    515 !      K_MASS(:,i_lev,1)=1.0
    516 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    517 !      END DO
    518                              
    519 !     Lower Haze: mode 2
    520 !      DO i_lev=cloudmin,22
    521 !      R_MEDIAN(:,i_lev,2)=1.4e-6
    522 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    523 !      STDDEV(:,i_lev,2)=1.23
    524 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    525 !      K_MASS(:,i_lev,2)=0.0
    526 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    527 !      END DO
    528 
    529 !     Lower Haze: mode 3
    530 !      DO i_lev=cloudmin,22
    531 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    532 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    533 !      STDDEV(:,i_lev,3)=1.28
    534 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    535 !      K_MASS(:,i_lev,3)=0.0
    536 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    537 !      END DO 
    538 
    539 !     Pre Cloud: mode 1
    540 !      DO i_lev=23,23
    541 !      R_MEDIAN(:,i_lev,1)=0.15e-6
    542 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    543 !      STDDEV(:,i_lev,1)=1.8
    544 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    545 !      K_MASS(:,i_lev,1)=0.04
    546 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    547 !      END DO
    548                              
    549 !     Pre Cloud: mode 2
    550 !      DO i_lev=23,23
    551 !      R_MEDIAN(:,i_lev,2)=1.0e-6
    552 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    553 !      STDDEV(:,i_lev,2)=1.29
    554 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    555 !      K_MASS(:,i_lev,2)=0.96
    556 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    557 !      END DO
    558 
    559 !     Pre Cloud: mode 3
    560 !      DO i_lev=23,23
    561 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    562 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    563 !      STDDEV(:,i_lev,3)=1.28
    564 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    565 !      K_MASS(:,i_lev,3)=0.0
    566 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    567 !      END DO
    568                  
    569 !      Lower Cloud: mode 1
    570 !      DO i_lev=24,24
    571 !      R_MEDIAN(:,i_lev,1)=0.2e-6
    572 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    573 !      STDDEV(:,i_lev,1)=1.8
    574 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    575 !      K_MASS(:,i_lev,1)=0.014
    576 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    577 !      END DO
    578                
    579 !     Lower Cloud: mode 2
    580 !      DO i_lev=24,24
    581 !      R_MEDIAN(:,i_lev,2)=1.0e-6
    582 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    583 !      STDDEV(:,i_lev,2)=1.6
    584 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    585 !      K_MASS(:,i_lev,2)=0.986
    586 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    587 !      END DO
    588 
    589 !     Lower Cloud: mode 3
    590 !      DO i_lev=24,24
    591 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    592 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    593 !      STDDEV(:,i_lev,3)=1.28
    594 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    595 !      K_MASS(:,i_lev,3)=0.
    596 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    597 !      END DO
    598 
    599 !     Middle Cloud: mode 1
    600 !      DO i_lev=25,28
    601 !      R_MEDIAN(:,i_lev,1)=0.15e-6
    602 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    603 !      STDDEV(:,i_lev,1)=1.9
    604 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    605 !      K_MASS(:,i_lev,1)=0.0084
    606 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    607 !      END DO
    608          
    609 !     Middle Cloud: mode 2 prime
    610 !      DO i_lev=25,28
    611 !      R_MEDIAN(:,i_lev,2)=1.4e-6
    612 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    613 !      STDDEV(:,i_lev,2)=1.6
    614 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    615 !      K_MASS(:,i_lev,2)=0.9916
    616 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    617 !      END DO
    618    
    619 !     Middle Cloud: mode 3
    620 !      DO i_lev=25,28
    621 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    622 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    623 !      STDDEV(:,i_lev,3)=1.28
    624 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    625 !      K_MASS(:,i_lev,3)=0.0
    626 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    627 !      END DO
    628 
    629 
    630 !     Upper Cloud: mode 1
    631 !      DO i_lev=29,35
    632 !      R_MEDIAN(:,i_lev,1)=0.2e-6
    633 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    634 !      STDDEV(:,i_lev,1)=2.16
    635 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    636 !      K_MASS(:,i_lev,1)=0.72
    637 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    638 !      END DO
    639 
    640 !     Upper Cloud: mode 2
    641 !      DO i_lev=29,35
    642 !      R_MEDIAN(:,i_lev,2)=1.0e-6
    643 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    644 !      STDDEV(:,i_lev,2)=1.29
    645 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    646 !      K_MASS(:,i_lev,2)=0.28
    647 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    648 !      END DO
    649        
    650 !     Upper Cloud: mode 3
    651 !      DO i_lev=29,35
    652 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    653 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    654 !      STDDEV(:,i_lev,3)=1.28
    655 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    656 !      K_MASS(:,i_lev,3)=0.0
    657 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    658 !      END DO
    659 
    660 !     Upper Haze: mode 1
    661 !      DO i_lev=36, cloudmax
    662 !      R_MEDIAN(:,i_lev,1)=0.2e-6
    663 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,1)
    664 !      STDDEV(:,i_lev,1)=2.16
    665 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,1)
    666 !      K_MASS(:,i_lev,1)=1.0
    667 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,1)
    668 !      END DO
    669 
    670 !     Upper Haze: mode 2
    671 !      DO i_lev=36, cloudmax
    672 !      R_MEDIAN(:,i_lev,2)=1.e-6
    673 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,2)
    674 !      STDDEV(:,i_lev,2)=1.29
    675 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,2)
    676 !      K_MASS(:,i_lev,2)=0.0
    677 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,2)
    678 !      END DO
    679      
    680 !     Upper Haze: mode 3
    681 !      DO i_lev=36, cloudmax
    682 !      R_MEDIAN(:,i_lev,3)=3.65e-6
    683 !      PRINT*,'level',i_lev,'R R_MEDIAN',R_MEDIAN(1,i_lev,3)
    684 !      STDDEV(:,i_lev,3)=2.16
    685 !      PRINT*,'level',i_lev,'Dev Std',STDDEV(1,i_lev,3)
    686 !      K_MASS(:,i_lev,3)=0.0
    687 !      PRINT*,'level',i_lev,'Coeff Mass: k_mass',K_MASS(1,i_lev,3)
    688 !      END DO       
    689 
    690 !=============================================================
    691       PRINT*,'==============================='
    692       PRINT*,'FIN Initialisation cloud layer'
    693       PRINT*,'==============================='
    694      
    695   END SUBROUTINE cloud_ini
    696  
    697 ! ===========================================================
    698 
    699   SUBROUTINE chemparam_ini
    700   USE infotrac_phy, ONLY: nqtot, tname
    701   IMPLICIT NONE
    702   INTEGER :: i
    703  
    704   ALLOCATE(M_tr(nqtot))    ! Molecular Mass of Chemistry Tracers
    705   ALLOCATE(type_tr(nqtot)) ! Type of chemistry tracers 1: neutral, 2: ion, 3: liquide
    706 
    707 ! Initialisation index chemistry tracers
    708 
    709   ! Neutral gas Tracer
    710   i_co2      = 0
    711   i_co       = 0
    712   i_h2       = 0
    713   i_h2o      = 0
    714   i_o1d      = 0
    715   i_o        = 0
    716   i_o2       = 0
    717   i_o2dg     = 0
    718   i_o3       = 0
    719   i_h        = 0
    720   i_oh       = 0
    721   i_ho2      = 0
    722   i_h2o2     = 0
    723   i_cl       = 0
    724   i_clo      = 0
    725   i_cl2      = 0
    726   i_hcl      = 0
    727   i_hocl     = 0
    728   i_clco     = 0
    729   i_clco3    = 0
    730   i_cocl2    = 0
    731   i_s        = 0
    732   i_so       = 0
    733   i_so2      = 0
    734   i_so3      = 0
    735   i_s2o2     = 0
    736   i_ocs      = 0
    737   i_hso3     = 0
    738   i_h2so4    = 0
    739   i_s2       = 0
    740   i_clso2    = 0
    741   i_oscl     = 0
    742   i_n2       = 0
    743   i_he       = 0
    744   i_n2d      = 0
    745   i_n        = 0
    746   i_no       = 0
    747   i_no2      = 0
    748   ! ION TRACERS FOR IONCHEM = TRUE
    749   i_co2plus  = 0
    750   i_coplus   = 0
    751   i_oplus    = 0
    752   i_o2plus   = 0
    753   i_n2plus   = 0
    754   i_hplus    = 0
    755   i_h2oplus  = 0
    756   i_nplus    = 0
    757   i_ohplus   = 0
    758   i_cplus    = 0
    759   i_noplus   = 0
    760   i_h3oplus  = 0
    761   i_hcoplus  = 0
    762   i_hco2plus = 0
    763   i_elec     = 0 
    764   ! MICROPHYSICAL TRACERS FOR CL_SCHEME=1
    765   i_h2oliq   = 0
    766   i_h2so4liq = 0 
    767 
    768   DO i=1, nqtot                                 
    769        
    770         PRINT*,'i',i
    771         PRINT*,'tname(i)',tname(i)
    772        
    773         SELECT CASE(tname(i))
    774 ! NEUTRAL TRACERS
    775                 CASE('co2')
    776                 i_co2=i
    777                 PRINT*,'co2',i_co2
    778                 M_tr(i_co2) = 44.0095
    779       type_tr(i_co2) = 1 
    780                 CASE('co')
    781                 i_co=i
    782                 PRINT*,'co',i_co
    783                 M_tr(i_co)=28.0101
    784       type_tr(i_co) = 1
    785                 CASE('h2')
    786                 i_h2=i
    787                 PRINT*,'h2',i_h2
    788                 M_tr(i_h2)= 2.01588
    789       type_tr(i_h2) = 1
    790                 CASE('h2o')
    791                 i_h2o=i
    792                 PRINT*,'h2o',i_h2o
    793                 M_tr(i_h2o)=18.0153
    794       type_tr(i_h2o) = 1
    795                 CASE('o1d')
    796                 i_o1d=i
    797                 PRINT*,'o1d',i_o1d
    798                 M_tr(i_o1d)=15.994
    799       type_tr(i_o1d) = 1
    800                 CASE('o')
    801                 i_o=i
    802                 PRINT*,'o',i_o
    803                 M_tr(i_o)=15.994
    804       type_tr(i_o) = 1
    805                 CASE('o2')
    806                 i_o2=i
    807                 PRINT*,'o2',i_o2
    808                 M_tr(i_o2)=31.9988
    809       type_tr(i_o2) = 1
    810                 CASE('o2dg')
    811                 i_o2dg=i
    812                 PRINT*,'o2dg',i_o2dg
    813                 M_tr(i_o2dg)=31.9988
    814       type_tr(i_o2dg) = 1
    815                 CASE('o3')
    816                 i_o3=i
    817                 PRINT*,'o3',i_o3
    818                 M_tr(i_o3)= 47.9982
    819       type_tr(i_o3) = 1
    820                 CASE('h')
    821                 i_h=i
    822                 PRINT*,'h',i_h
    823                 M_tr(i_h)= 1.00794
    824       type_tr(i_h) = 1 
    825                 CASE('oh')
    826                 i_oh=i
    827                 PRINT*,'oh',i_oh
    828                 M_tr(i_oh)=17.0073
    829       type_tr(i_oh) = 1
    830                 CASE('ho2')
    831                 i_ho2=i
    832                 PRINT*,'ho2',i_ho2
    833                 M_tr(i_ho2)=33.0067
    834       type_tr(i_ho2) = 1
    835                 CASE('h2o2')
    836                 i_h2o2=i
    837                 PRINT*,'h2o2',i_h2o2
    838                 M_tr(i_h2o2)=34.0147
    839       type_tr(i_h2o2) = 1
    840                 CASE('cl')
    841                 i_cl=i
    842                 PRINT*,'cl',i_cl
    843                 M_tr(i_cl)=35.453
    844       type_tr(i_cl) = 1
    845                 CASE('clo')
    846                 i_clo=i
    847                 PRINT*,'clo',i_clo
    848                 M_tr(i_clo)=51.452
    849       type_tr(i_clo) = 1
    850                 CASE('cl2')
    851                 i_cl2=i
    852                 PRINT*,'cl2',i_cl2
    853                 M_tr(i_cl2)=70.906
    854       type_tr(i_cl2) = 1
    855                 CASE('hcl')
    856                 i_hcl=i
    857                 PRINT*,'hcl',i_hcl
    858                 M_tr(i_hcl)=36.461
    859       type_tr(i_hcl) = 1
    860                 CASE('hocl')
    861                 i_hocl=i
    862                 PRINT*,'hocl',i_hocl
    863                 M_tr(i_hocl)=52.46
    864                 type_tr(i_hocl) = 1
    865       CASE('clco')
    866                 i_clco=i
    867                 PRINT*,'clco',i_clco
    868                 M_tr(i_clco)=63.463
    869                 type_tr(i_clco) = 1
    870       CASE('clco3')
    871                 i_clco3=i
    872                 PRINT*,'clco3',i_clco3
    873                 M_tr(i_clco3)=95.462
    874                 type_tr(i_clco3) = 1
    875       CASE('cocl2')
    876                 i_cocl2=i
    877                 PRINT*,'cocl2',i_cocl2
    878                 M_tr(i_cocl2)=98.916
    879                 type_tr(i_cocl2) = 1
    880       CASE('s')
    881                 i_s=i
    882                 PRINT*,'s',i_s
    883                 M_tr(i_s)=32.065
    884                 type_tr(i_s) = 1
    885       CASE('so')
    886                 i_so=i
    887                 PRINT*,'so',i_so
    888                 M_tr(i_so)=48.0644
    889                 type_tr(i_so) = 1
    890       CASE('so2')
    891                 i_so2=i
    892                 PRINT*,'so2',i_so2
    893                 M_tr(i_so2)=64.064
    894                 type_tr(i_so2) = 1
    895       CASE('so3')
    896                 i_so3=i
    897                 PRINT*,'so3',i_so3
    898                 M_tr(i_so3)=80.063
    899                 type_tr(i_so3) = 1
    900       CASE('s2o2')
    901                 i_s2o2=i
    902                 PRINT*,'s2o2',i_s2o2
    903                 M_tr(i_s2o2)= 96.1288
    904                 type_tr(i_s2o2) = 1
    905       CASE('ocs')
    906                 i_ocs=i
    907                 PRINT*,'ocs',i_ocs
    908                 M_tr(i_ocs)=60.0751
    909                 type_tr(i_ocs) = 1
    910       CASE('hso3')
    911                 i_hso3=i
    912                 PRINT*,'hso3',i_hso3
    913                 M_tr(i_hso3)=81.071
    914                 type_tr(i_hso3) = 1
    915       CASE('h2so4')
    916                 i_h2so4=i
    917                 PRINT*,'h2so4',i_h2so4
    918                 M_tr(i_h2so4)=98.078
    919                 type_tr(i_h2so4) = 1
    920       CASE('s2')
    921                 i_s2=i
    922                 PRINT*,'s2',i_s2
    923                 M_tr(i_s2)=64.13
    924                 type_tr(i_s2) = 1
    925       CASE('clso2')
    926                 i_clso2=i
    927                 PRINT*,'clso2',i_clso2
    928                 M_tr(i_clso2)=99.517
    929                 type_tr(i_clso2) = 1
    930       CASE('oscl')
    931                 i_oscl=i
    932                 PRINT*,'oscl',i_oscl
    933                 M_tr(i_oscl)=83.517
    934                   type_tr(i_oscl) = 1
    935       CASE('n2')
    936                   i_n2=i
    937       PRINT*,'n2',i_n2
    938                   M_tr(i_n2)=28.013
    939                   type_tr(i_n2) = 1
    940       CASE('he')
    941                   i_he=i
    942       PRINT*,'he',i_he
    943                   M_tr(i_he)=4.0026
    944                   type_tr(i_he) = 1
    945       CASE('n2d')
    946                   i_n2d=i
    947       PRINT*,'n2d',i_n2d
    948                   M_tr(i_n2d)=14.0067
    949                   type_tr(i_n2d) = 1
    950       CASE('n')
    951                   i_n=i
    952       PRINT*,'n',i_n
    953                   M_tr(i_n)=14.0067
    954                   type_tr(i_n) = 1
    955       CASE('no')
    956                   i_no=i
    957       PRINT*,'no',i_no
    958                   M_tr(i_no)=30.0061   
    959                   type_tr(i_no) = 1
    960       CASE('no2')
    961                   i_no2=i
    962       PRINT*,'no2',i_no2
    963                   M_tr(i_no2)=46.0055
    964       type_tr(i_no2) = 1         
    965 ! ION TRACERS FOR OK_IONCHEM = TRUE
    966       CASE('co2plus')
    967       i_co2plus=i
    968       PRINT*,'co2plus',i_co2plus
    969       M_tr(i_co2plus)=44.0095
    970       type_tr(i_co2plus) = 2
    971       CASE('coplus')
    972       i_coplus=i
    973       PRINT*,'coplus',i_coplus
    974       M_tr(i_coplus)=28.0101
    975       type_tr(i_coplus) = 2
    976       CASE('oplus')
    977       i_oplus=i
    978       PRINT*,'oplus',i_oplus
    979       M_tr(i_oplus)=15.994
    980       type_tr(i_oplus) = 2
    981       CASE('o2plus')
    982       i_o2plus=i
    983       PRINT*,'o2plus',i_o2plus
    984       M_tr(i_o2plus)=31.9988
    985       type_tr(i_o2plus) = 2
    986       CASE('n2plus')
    987       i_n2plus=i
    988       PRINT*,'n2plus',i_n2plus
    989       M_tr(i_n2plus)=28.013
    990       type_tr(i_n2plus) = 2
    991       CASE('hplus')
    992       i_hplus=i
    993       PRINT*,'hplus',i_hplus
    994       M_tr(i_hplus)=1.00794
    995                   type_tr(i_hplus) = 2
    996       CASE('h2oplus')
    997       i_h2oplus=i
    998       PRINT*,'h2oplus',i_h2oplus
    999       M_tr(i_h2oplus)=18.0153
    1000                   type_tr(i_h2oplus) = 2     
    1001       CASE('nplus')
    1002                   i_nplus=i
    1003       PRINT*,'nplus',i_nplus
    1004                   M_tr(i_nplus)=14.0067
    1005                   type_tr(i_nplus) = 2
    1006       CASE('ohplus')
    1007                   i_ohplus=i
    1008       PRINT*,'ohplus',i_ohplus
    1009                   M_tr(i_ohplus)=17.0073
    1010                   type_tr(i_ohplus) = 2
    1011       CASE('cplus')
    1012                   i_cplus=i
    1013       PRINT*,'cplus',i_cplus
    1014                   M_tr(i_cplus)=12.011
    1015                   type_tr(i_cplus) = 2
    1016       CASE('noplus')
    1017                   i_noplus=i
    1018       PRINT*,'noplus',i_noplus
    1019                   M_tr(i_noplus)=30.0061
    1020                   type_tr(i_noplus) = 2
    1021       CASE('h3oplus')
    1022                   i_h3oplus=i
    1023       PRINT*,'h3oplus',i_h3oplus
    1024                   M_tr(i_h3oplus)=19.0232
    1025                   type_tr(i_h3oplus) = 2
    1026       CASE('hcoplus')
    1027                   i_hcoplus=i
    1028       PRINT*,'hcoplus',i_hcoplus
    1029                   M_tr(i_hcoplus)=29.0180
    1030                   type_tr(i_hcoplus) = 2
    1031       CASE('hco2plus')
    1032                   i_hco2plus=i
    1033       PRINT*,'hco2plus',i_hco2plus
    1034                   M_tr(i_hco2plus)=45.
    1035                   type_tr(i_hco2plus) = 2
    1036       CASE('elec')
    1037                   i_elec=i
    1038       PRINT*,'elec',i_elec
    1039                   M_tr(i_elec)=1./1822.89
    1040       type_tr(i_elec) = 2       
    1041 ! MICROPHYSICAL TRACERS FOR CL_SCHEME=1
    1042                 CASE('h2oliq')
    1043                 i_h2oliq=i
    1044                 PRINT*,'h2oliq',i_h2oliq
    1045                 M_tr(i_h2oliq)=18.0153
    1046       type_tr(i_h2oliq) = 3
    1047                 CASE('h2so4liq')
    1048                 i_h2so4liq=i
    1049                 PRINT*,'h2so4liq',i_h2so4liq
    1050                 M_tr(i_h2so4liq)=98.078
    1051       type_tr(i_h2so4liq) = 3
    1052 ! MICROPHYSICAL TRACERS FOR CL_SCHEME=2
    1053     CASE('M0_aer')
    1054                 i_m0_aer=i
    1055     type_tr(i_m0_aer) = 10
    1056                 PRINT*,'M0_aer',i_m0_aer
    1057                 CASE('M3_aer')
    1058                 i_m3_aer=i
    1059     type_tr(i_m3_aer) = 10
    1060                 PRINT*,'M3_aer',i_m3_aer
    1061     CASE('M0_m1drop')
    1062                 i_m0_mode1drop=i
    1063     type_tr(i_m0_mode1drop) = 10
    1064                 PRINT*,'M0_m1drop',i_m0_mode1drop
    1065     CASE('M0_m1ccn')
    1066                 i_m0_mode1ccn=i
    1067     type_tr(i_m0_mode1ccn) = 10
    1068                 PRINT*,'M0_m1ccn',i_m0_mode1ccn
    1069                 CASE('M3_m1sa')
    1070                 i_m3_mode1sa=i
    1071     type_tr(i_m3_mode1sa) = 10
    1072                 PRINT*,'M3_m1sa',i_m3_mode1sa
    1073                 CASE('M3_m1w')
    1074                 i_m3_mode1w=i
    1075     type_tr(i_m3_mode1w) = 10
    1076                 PRINT*,'M3_m1w',i_m3_mode1w
    1077                 CASE('M3_m1ccn')
    1078                 i_m3_mode1ccn=i
    1079     type_tr(i_m3_mode1ccn) = 10
    1080                 PRINT*,'M3_m1ccn',i_m3_mode1ccn
    1081     CASE('M0_m2drop')
    1082                 i_m0_mode2drop=i
    1083     type_tr(i_m0_mode2drop) = 10
    1084                 PRINT*,'M0_m2drop',i_m0_mode2drop
    1085     CASE('M0_m2ccn')
    1086                 i_m0_mode2ccn=i
    1087     type_tr(i_m0_mode2ccn) = 10
    1088                 PRINT*,'M0_m2ccn',i_m0_mode2ccn
    1089                 CASE('M3_m2sa')
    1090                 i_m3_mode2sa=i
    1091     type_tr(i_m3_mode2sa) = 10
    1092                 PRINT*,'M3_m2sa',i_m3_mode2sa
    1093                 CASE('M3_m2w')
    1094                 i_m3_mode2w=i
    1095     type_tr(i_m3_mode2w) = 10
    1096                 PRINT*,'M3_m2w',i_m3_mode2w
    1097                 CASE('M3_m2ccn')
    1098                 i_m3_mode2ccn=i
    1099     type_tr(i_m3_mode2ccn) = 10
    1100                 PRINT*,'M3_m2ccn',i_m3_mode2ccn
    1101         END SELECT
    1102        
    1103 !       PRINT*,'M_tr(i)',M_tr(i)
    1104   END DO
    1105  
    1106   END SUBROUTINE chemparam_ini
    1107 
    1108 ! ===========================================================
    1109 
    1110   SUBROUTINE vapors4muphy_ini(nlon,nlev,trac)
    1111   USE infotrac_phy, ONLY: nqtot, tname
    1112   IMPLICIT NONE
     514!     lower haze: mode 1
     515!      do i_lev=cloudmin,22
     516!      r_median(:,i_lev,1)=0.1e-6
     517!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     518!      stddev(:,i_lev,1)=1.57
     519!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     520!      k_mass(:,i_lev,1)=1.0
     521!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     522!      end do
     523
     524!     lower haze: mode 2
     525!      do i_lev=cloudmin,22
     526!      r_median(:,i_lev,2)=1.4e-6
     527!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     528!      stddev(:,i_lev,2)=1.23
     529!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     530!      k_mass(:,i_lev,2)=0.0
     531!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     532!      end do
     533
     534!     lower haze: mode 3
     535!      do i_lev=cloudmin,22
     536!      r_median(:,i_lev,3)=3.65e-6
     537!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     538!      stddev(:,i_lev,3)=1.28
     539!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     540!      k_mass(:,i_lev,3)=0.0
     541!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     542!      end do
     543
     544!     pre cloud: mode 1
     545!      do i_lev=23,23
     546!      r_median(:,i_lev,1)=0.15e-6
     547!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     548!      stddev(:,i_lev,1)=1.8
     549!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     550!      k_mass(:,i_lev,1)=0.04
     551!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     552!      end do
     553
     554!     pre cloud: mode 2
     555!      do i_lev=23,23
     556!      r_median(:,i_lev,2)=1.0e-6
     557!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     558!      stddev(:,i_lev,2)=1.29
     559!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     560!      k_mass(:,i_lev,2)=0.96
     561!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     562!      end do
     563
     564!     pre cloud: mode 3
     565!      do i_lev=23,23
     566!      r_median(:,i_lev,3)=3.65e-6
     567!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     568!      stddev(:,i_lev,3)=1.28
     569!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     570!      k_mass(:,i_lev,3)=0.0
     571!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     572!      end do
     573
     574!      lower cloud: mode 1
     575!      do i_lev=24,24
     576!      r_median(:,i_lev,1)=0.2e-6
     577!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     578!      stddev(:,i_lev,1)=1.8
     579!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     580!      k_mass(:,i_lev,1)=0.014
     581!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     582!      end do
     583
     584!     lower cloud: mode 2
     585!      do i_lev=24,24
     586!      r_median(:,i_lev,2)=1.0e-6
     587!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     588!      stddev(:,i_lev,2)=1.6
     589!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     590!      k_mass(:,i_lev,2)=0.986
     591!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     592!      end do
     593
     594!     lower cloud: mode 3
     595!      do i_lev=24,24
     596!      r_median(:,i_lev,3)=3.65e-6
     597!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     598!      stddev(:,i_lev,3)=1.28
     599!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     600!      k_mass(:,i_lev,3)=0.
     601!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     602!      end do
     603
     604!     middle cloud: mode 1
     605!      do i_lev=25,28
     606!      r_median(:,i_lev,1)=0.15e-6
     607!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     608!      stddev(:,i_lev,1)=1.9
     609!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     610!      k_mass(:,i_lev,1)=0.0084
     611!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     612!      end do
     613
     614!     middle cloud: mode 2 prime
     615!      do i_lev=25,28
     616!      r_median(:,i_lev,2)=1.4e-6
     617!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     618!      stddev(:,i_lev,2)=1.6
     619!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     620!      k_mass(:,i_lev,2)=0.9916
     621!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     622!      end do
     623
     624!     middle cloud: mode 3
     625!      do i_lev=25,28
     626!      r_median(:,i_lev,3)=3.65e-6
     627!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     628!      stddev(:,i_lev,3)=1.28
     629!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     630!      k_mass(:,i_lev,3)=0.0
     631!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     632!      end do
     633
     634
     635!     upper cloud: mode 1
     636!      do i_lev=29,35
     637!      r_median(:,i_lev,1)=0.2e-6
     638!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     639!      stddev(:,i_lev,1)=2.16
     640!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     641!      k_mass(:,i_lev,1)=0.72
     642!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     643!      end do
     644
     645!     upper cloud: mode 2
     646!      do i_lev=29,35
     647!      r_median(:,i_lev,2)=1.0e-6
     648!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     649!      stddev(:,i_lev,2)=1.29
     650!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     651!      k_mass(:,i_lev,2)=0.28
     652!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     653!      end do
     654
     655!     upper cloud: mode 3
     656!      do i_lev=29,35
     657!      r_median(:,i_lev,3)=3.65e-6
     658!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     659!      stddev(:,i_lev,3)=1.28
     660!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     661!      k_mass(:,i_lev,3)=0.0
     662!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     663!      end do
     664
     665!     upper haze: mode 1
     666!      do i_lev=36, cloudmax
     667!      r_median(:,i_lev,1)=0.2e-6
     668!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     669!      stddev(:,i_lev,1)=2.16
     670!      print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     671!      k_mass(:,i_lev,1)=1.0
     672!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     673!      end do
     674
     675!     upper haze: mode 2
     676!      do i_lev=36, cloudmax
     677!      r_median(:,i_lev,2)=1.e-6
     678!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     679!      stddev(:,i_lev,2)=1.29
     680!      print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     681!      k_mass(:,i_lev,2)=0.0
     682!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     683!      end do
     684
     685!     upper haze: mode 3
     686!      do i_lev=36, cloudmax
     687!      r_median(:,i_lev,3)=3.65e-6
     688!      print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     689!      stddev(:,i_lev,3)=2.16
     690!      print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     691!      k_mass(:,i_lev,3)=0.0
     692!      print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     693!      end do
     694
     695   ! ========================================================
     696   ! initialisation bimodale k&h 1980 with mode 3 fully solid
     697   ! ========================================================
     698   !    ! mode 3 fully solid
     699   !    qrad=1
     700   !    ! normally nb_mode=2 in physiq.def !!!
     701   !    do ilon=1,nbr_lon
     702   ! !     mode 1
     703   !       do i_lev=cloudmin,20
     704   !          r_median(ilon,i_lev,1)=0.125e-6
     705   !          print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     706   !          stddev(ilon,i_lev,1)=1.57
     707   !          print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     708   !          k_mass(ilon,i_lev,1)=1.0
     709   !          print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     710   !       end do
     711
     712   !       r_median(ilon,21,1)=0.125e-6
     713   !       print*,'level',21,'r r_median',r_median(1,21,1)
     714   !       stddev(ilon,21,1)=1.57
     715   !       print*,'level',21,'dev std',stddev(1,21,1)
     716   !       k_mass(ilon,21,1)=0.02
     717   !       print*,'level',21,'coeff mass: k_mass',k_mass(1,21,1)
     718
     719   !       r_median(ilon,22,1)=0.2e-6
     720   !       print*,'level',22,'r r_median',r_median(1,22,1)
     721   !       stddev(ilon,22,1)=1.8
     722   !       print*,'level',22,'dev std',stddev(1,22,1)
     723   !       k_mass(ilon,22,1)=0.02
     724   !       print*,'level',22,'coeff mass: k_mass',k_mass(1,22,1)
     725
     726   !       r_median(ilon,23,1)=0.15e-6
     727   !       print*,'level',23,'r r_median',r_median(1,23,1)
     728   !       stddev(ilon,23,1)=1.8
     729   !       print*,'level',23,'dev std',stddev(1,23,1)
     730   !       k_mass(ilon,23,1)=0.02
     731   !       print*,'level',23,'coeff mass: k_mass',k_mass(1,23,1)
     732
     733   !       do i_lev=24,25
     734   !          r_median(ilon,i_lev,1)=0.15e-6
     735   !          print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     736   !          stddev(ilon,i_lev,1)=1.9
     737   !          print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     738   !          k_mass(ilon,i_lev,1)=0.02
     739   !          print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     740   !       end do
     741
     742   !       r_median(ilon,26,1)=0.175e-6
     743   !       print*,'level',26,'r r_median',r_median(1,26,1)
     744   !       stddev(ilon,26,1)=2.16
     745   !       print*,'level',26,'dev std',stddev(1,26,1)
     746   !       k_mass(ilon,26,1)=0.175
     747   !       print*,'level',26,'coeff mass: k_mass',k_mass(1,26,1)
     748
     749   !       do i_lev=27,33
     750   !          r_median(ilon,i_lev,1)=0.175e-6
     751   !          print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     752   !          stddev(ilon,i_lev,1)=2.16
     753   !          print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     754   !          k_mass(ilon,i_lev,1)=0.25
     755   !          print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     756   !       end do
     757
     758   !       do i_lev=34,cloudmax
     759   !          r_median(ilon,i_lev,1)=0.175e-6
     760   !          print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     761   !          stddev(ilon,i_lev,1)=2.16
     762   !          print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     763   !          k_mass(ilon,i_lev,1)=1.0
     764   !          print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     765   !       end do
     766
     767   !    !     mode 2
     768   !       do i_lev=cloudmin,20
     769   !          r_median(ilon,i_lev,2)=1.4e-6
     770   !          print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     771   !          stddev(ilon,i_lev,2)=1.35
     772   !          print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     773   !          k_mass(ilon,i_lev,2)=0.0
     774   !          print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     775   !       end do
     776
     777   !       do i_lev=21,22
     778   !          r_median(ilon,i_lev,2)=1.4e-6
     779   !          print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     780   !          stddev(ilon,i_lev,2)=1.35
     781   !          print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     782   !          k_mass(ilon,i_lev,2)=0.98
     783   !          print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     784   !       end do
     785
     786   !       r_median(ilon,23,2)=1.35e-6
     787   !       print*,'level',23,'r r_median',r_median(1,23,2)
     788   !       stddev(ilon,23,2)=1.25
     789   !       print*,'level',23,'dev std',stddev(1,23,2)
     790   !       k_mass(ilon,23,2)=0.98
     791   !       print*,'level',23,'coeff mass: k_mass',k_mass(1,23,2)
     792
     793
     794   !       r_median(ilon,24,2)=1.375e-6
     795   !       print*,'level',24,'r r_median',r_median(1,24,2)
     796   !       stddev(ilon,24,2)=1.2
     797   !       print*,'level',24,'dev std',stddev(1,24,2)
     798   !       k_mass(ilon,24,2)=0.98
     799   !       print*,'level',24,'coeff mass: k_mass',k_mass(1,24,2)
     800
     801
     802   !       r_median(ilon,25,2)=1.4e-6
     803   !       print*,'level',25,'r r_median',r_median(1,25,2)
     804   !       stddev(ilon,25,2)=1.16
     805   !       print*,'level',25,'dev std',stddev(1,25,2)
     806   !       k_mass(ilon,25,2)=0.98
     807   !       print*,'level',25,'coeff mass: k_mass',k_mass(1,25,2)
     808
     809
     810   !       r_median(ilon,26,2)=1.15e-6
     811   !       print*,'level',26,'r r_median',r_median(1,26,2)
     812   !       stddev(ilon,26,2)=1.34
     813   !       print*,'level',26,'dev std',stddev(1,26,2)
     814   !       k_mass(ilon,26,2)=0.825
     815   !       print*,'level',26,'coeff mass: k_mass',k_mass(1,26,2)
     816
     817
     818   !       r_median(ilon,27,2)=1.14e-6
     819   !       print*,'level',27,'r r_median',r_median(1,27,2)
     820   !       stddev(ilon,27,2)=1.33
     821   !       print*,'level',27,'dev std',stddev(1,27,2)
     822   !       k_mass(ilon,27,2)=0.75
     823   !       print*,'level',27,'coeff mass: k_mass',k_mass(1,27,2)
     824
     825
     826   !       r_median(ilon,28,2)=1.35e-6
     827   !       print*,'level',28,'r r_median',r_median(1,28,2)
     828   !       stddev(ilon,28,2)=1.32
     829   !       print*,'level',28,'dev std',stddev(1,28,2)
     830   !       k_mass(ilon,28,2)=0.75
     831   !       print*,'level',28,'coeff mass: k_mass',k_mass(1,28,2)
     832
     833
     834   !       r_median(ilon,29,2)=1.125e-6
     835   !       print*,'level',29,'r r_median',r_median(1,29,2)
     836   !       stddev(ilon,29,2)=1.31
     837   !       print*,'level',29,'dev std',stddev(1,29,2)
     838   !       k_mass(ilon,29,2)=0.75
     839   !       print*,'level',29,'coeff mass: k_mass',k_mass(1,29,2)
     840
     841
     842   !       r_median(ilon,30,2)=1.118e-6
     843   !       print*,'level',30,'r r_median',r_median(1,30,2)
     844   !       stddev(ilon,30,2)=1.30
     845   !       print*,'level',30,'dev std',stddev(1,30,2)
     846   !       k_mass(ilon,30,2)=0.75
     847   !       print*,'level',30,'coeff mass: k_mass',k_mass(1,30,2)
     848
     849
     850   !       r_median(ilon,31,2)=1.11e-6
     851   !       print*,'level',31,'r r_median',r_median(1,31,2)
     852   !       stddev(ilon,31,2)=1.29
     853   !       print*,'level',31,'dev std',stddev(1,31,2)
     854   !       k_mass(ilon,31,2)=0.75
     855   !       print*,'level',31,'coeff mass: k_mass',k_mass(1,31,2)
     856
     857   !       do i_lev=32,33
     858   !          r_median(ilon,i_lev,2)=1.1e-6
     859   !          print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     860   !          stddev(ilon,i_lev,2)=1.28
     861   !          print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     862   !          k_mass(ilon,i_lev,2)=0.75
     863   !          print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     864   !       end do
     865
     866   !       do i_lev=34,cloudmax
     867   !          r_median(ilon,i_lev,2)=1.1e-6
     868   !          print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     869   !          stddev(ilon,i_lev,2)=1.28
     870   !          print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     871   !          k_mass(ilon,i_lev,2)=0.0
     872   !          print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     873   !       end do
     874   !    end do
     875
     876! ======================================================
     877! knollenberg & hunten (1980) with mode 3 97% solid
     878! ======================================================
     879
     880! mode 3 97% solid => qrad=0.97
     881
     882qrad = 0.97
     883
     884do ilon = 1,nbr_lon
     885
     886   ! mode 1
     887
     888   do i_lev = cloudmin,21
     889      r_median(ilon,i_lev,1) = 0.125e-6
     890      stddev(ilon,i_lev,1)   = 1.57
     891      k_mass(ilon,i_lev,1)   = 1.0
     892   end do
     893
     894   r_median(ilon,22,1)       = 0.125e-6
     895   stddev(ilon,22,1)         = 1.57
     896   k_mass(ilon,22,1)         = 1.0
     897
     898   r_median(ilon,23,1)       = 0.2e-6
     899   stddev(ilon,23,1)         = 1.8
     900   k_mass(ilon,23,1)         = 0.01
     901
     902   r_median(ilon,24,1)       = 0.15e-6
     903   stddev(ilon,24,1)         = 1.8
     904   k_mass(ilon,24,1)         = 0.01
     905
     906   do i_lev = 25,26
     907      r_median(ilon,i_lev,1) = 0.15e-6
     908      stddev(ilon,i_lev,1)   = 1.9
     909      k_mass(ilon,i_lev,1)   = 0.01
     910   end do
     911
     912   r_median(ilon,27,1)       = 0.175e-6
     913   stddev(ilon,27,1)         = 2.16
     914   k_mass(ilon,27,1)         = 0.175
     915
     916   do i_lev = 28,34
     917      r_median(ilon,i_lev,1) = 0.175e-6
     918      stddev(ilon,i_lev,1)   = 2.16
     919      k_mass(ilon,i_lev,1)   = 0.25
     920   end do
     921
     922   do i_lev = 35,cloudmax
     923      r_median(ilon,i_lev,1) = 0.175e-6
     924      stddev(ilon,i_lev,1)   = 2.16
     925      k_mass(ilon,i_lev,1)   = 0.25
     926   end do
     927
     928   !  mode 2
     929
     930   do i_lev = cloudmin,21
     931      r_median(ilon,i_lev,2) = 1.4e-6
     932      stddev(ilon,i_lev,2)   = 1.35
     933      k_mass(ilon,i_lev,2)   = 0.0
     934   end do
     935
     936   r_median(ilon,22,2)       = 1.4e-6
     937   stddev(ilon,22,2)         = 1.25
     938   k_mass(ilon,22,2)         = 0.0
     939
     940   r_median(ilon,23,2)       = 1.4e-6
     941   stddev(ilon,23,2)         = 1.25
     942   k_mass(ilon,23,2)         = 0.35
     943
     944   r_median(ilon,24,2)       = 1.35e-6
     945   stddev(ilon,24,2)         = 1.25
     946   k_mass(ilon,24,2)         = 0.35
     947
     948   r_median(ilon,25,2)       = 1.375e-6
     949   stddev(ilon,25,2)         = 1.2
     950   k_mass(ilon,25,2)         = 0.35
     951
     952   r_median(ilon,26,2)       = 1.4e-6
     953   stddev(ilon,26,2)         = 1.16
     954   k_mass(ilon,26,2)         = 0.35
     955
     956   r_median(ilon,27,2)       = 1.15e-6
     957   stddev(ilon,27,2)         = 1.34
     958   k_mass(ilon,27,2)         = 0.825
     959
     960   r_median(ilon,28,2)       = 1.14e-6
     961   stddev(ilon,28,2)         = 1.33
     962   k_mass(ilon,28,2)         = 0.75
     963
     964   r_median(ilon,29,2)       = 1.35e-6
     965   stddev(ilon,29,2)         = 1.32
     966   k_mass(ilon,29,2)         = 0.75
     967
     968   r_median(ilon,30,2)       = 1.125e-6
     969   stddev(ilon,30,2)         = 1.31
     970   k_mass(ilon,30,2)         = 0.75
     971
     972   r_median(ilon,31,2)       = 1.118e-6
     973   stddev(ilon,31,2)         = 1.31
     974   k_mass(ilon,31,2)         = 0.75
     975
     976   r_median(ilon,32,2)       = 1.11e-6
     977   stddev(ilon,32,2)         = 1.29
     978   k_mass(ilon,32,2)         = 0.75
     979
     980   do i_lev = 33,34
     981      r_median(ilon,i_lev,2) = 1.1e-6
     982      stddev(ilon,i_lev,2)   = 1.28
     983      k_mass(ilon,i_lev,2)   = 0.75
     984   end do
     985
     986   do i_lev = 35,cloudmax
     987      r_median(ilon,i_lev,2) = 1.1e-6
     988      stddev(ilon,i_lev,2)   = 1.28
     989      k_mass(ilon,i_lev,2)   = 0.75
     990   end do
     991
     992   ! mode 3
     993
     994   do i_lev = cloudmin,22
     995      r_median(ilon,i_lev,3) = 3.65e-6
     996      stddev(ilon,i_lev,3)   = 1.28
     997      k_mass(ilon,i_lev,3)   = 0.0
     998   end do
     999
     1000   do i_lev = 23,26
     1001      r_median(ilon,i_lev,3) = 3.65e-6
     1002      stddev(ilon,i_lev,3)   = 1.28
     1003      k_mass(ilon,i_lev,3)   = 0.64
     1004   end do
     1005
     1006   do i_lev = 27,cloudmax
     1007      r_median(ilon,i_lev,3) = 3.65e-6
     1008      stddev(ilon,i_lev,3)   = 1.28
     1009      k_mass(ilon,i_lev,3)   = 0.0
     1010   end do
     1011end do
     1012
     1013      ! ! ==================================================================
     1014      ! ! initialisation bimodale k&h 1980 with mode 3 0% solid fully liquid
     1015      ! ! ==================================================================
     1016      ! ! mode 3 0% solid, fully liquid
     1017      ! qrad=0.0
     1018      ! ! normally nb_mode=3 in physiq.def !!!
     1019      ! do ilon=1,nbr_lon
     1020      !    ! mode 1
     1021      !    do i_lev=cloudmin,20
     1022      !       r_median(ilon,i_lev,1)=0.125e-6
     1023      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     1024      !       stddev(ilon,i_lev,1)=1.57
     1025      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     1026      !       k_mass(ilon,i_lev,1)=1.0
     1027      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     1028      !    end do
     1029
     1030      !    r_median(ilon,21,1)=0.125e-6
     1031      !    print*,'level',21,'r r_median',r_median(1,21,1)
     1032      !    stddev(ilon,21,1)=1.57
     1033      !    print*,'level',21,'dev std',stddev(1,21,1)
     1034      !    k_mass(ilon,21,1)=1.0
     1035      !    print*,'level',21,'coeff mass: k_mass',k_mass(1,21,1)
     1036
     1037      !    r_median(ilon,22,1)=0.2e-6
     1038      !    print*,'level',22,'r r_median',r_median(1,22,1)
     1039      !    stddev(ilon,22,1)=1.8
     1040      !    print*,'level',22,'dev std',stddev(1,22,1)
     1041      !    k_mass(ilon,22,1)=0.01
     1042      !    print*,'level',22,'coeff mass: k_mass',k_mass(1,22,1)
     1043
     1044      !    r_median(ilon,23,1)=0.15e-6
     1045      !    print*,'level',23,'r r_median',r_median(1,23,1)
     1046      !    stddev(ilon,23,1)=1.8
     1047      !    print*,'level',23,'dev std',stddev(1,23,1)
     1048      !    k_mass(ilon,23,1)=0.01
     1049      !    print*,'level',23,'coeff mass: k_mass',k_mass(1,23,1)
     1050
     1051      !    do i_lev=24,25
     1052      !       r_median(ilon,i_lev,1)=0.15e-6
     1053      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     1054      !       stddev(ilon,i_lev,1)=1.9
     1055      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     1056      !       k_mass(ilon,i_lev,1)=0.01
     1057      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     1058      !    end do
     1059
     1060      !    r_median(ilon,26,1)=0.175e-6
     1061      !    print*,'level',26,'r r_median',r_median(1,26,1)
     1062      !    stddev(ilon,26,1)=2.16
     1063      !    print*,'level',26,'dev std',stddev(1,26,1)
     1064      !    k_mass(ilon,26,1)=0.175
     1065      !    print*,'level',26,'coeff mass: k_mass',k_mass(1,26,1)
     1066
     1067      !    do i_lev=27,33
     1068      !       r_median(ilon,i_lev,1)=0.175e-6
     1069      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     1070      !       stddev(ilon,i_lev,1)=2.16
     1071      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     1072      !       k_mass(ilon,i_lev,1)=0.25
     1073      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     1074      !    end do
     1075
     1076      !    do i_lev=34,cloudmax
     1077      !       r_median(ilon,i_lev,1)=0.175e-6
     1078      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     1079      !       stddev(ilon,i_lev,1)=2.16
     1080      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     1081      !       k_mass(ilon,i_lev,1)=0.25
     1082      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     1083      !    end do
     1084
     1085      !    !  mode 2
     1086      !    do i_lev=cloudmin,20
     1087      !       r_median(ilon,i_lev,2)=1.4e-6
     1088      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     1089      !       stddev(ilon,i_lev,2)=1.35
     1090      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     1091      !       k_mass(ilon,i_lev,2)=0.0
     1092      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     1093      !    end do
     1094
     1095      !    r_median(ilon,21,2)=1.4e-6
     1096      !    print*,'level',21,'r r_median',r_median(1,21,2)
     1097      !    stddev(ilon,21,2)=1.25
     1098      !    print*,'level',21,'dev std',stddev(1,21,2)
     1099      !    k_mass(ilon,21,2)=0.0
     1100      !    print*,'level',21,'coeff mass: k_mass',k_mass(1,21,2)
     1101
     1102      !    r_median(ilon,22,2)=1.4e-6
     1103      !    print*,'level',22,'r r_median',r_median(1,22,2)
     1104      !    stddev(ilon,22,2)=1.25
     1105      !    print*,'level',22,'dev std',stddev(1,22,2)
     1106      !    k_mass(ilon,22,2)=0.04
     1107      !    print*,'level',22,'coeff mass: k_mass',k_mass(1,22,2)
     1108
     1109      !    r_median(ilon,23,2)=1.35e-6
     1110      !    print*,'level',23,'r r_median',r_median(1,23,2)
     1111      !    stddev(ilon,23,2)=1.25
     1112      !    print*,'level',23,'dev std',stddev(1,23,2)
     1113      !    k_mass(ilon,23,2)=0.04
     1114      !    print*,'level',23,'coeff mass: k_mass',k_mass(1,23,2)
     1115
     1116
     1117      !    r_median(ilon,24,2)=1.375e-6
     1118      !    print*,'level',24,'r r_median',r_median(1,24,2)
     1119      !    stddev(ilon,24,2)=1.2
     1120      !    print*,'level',24,'dev std',stddev(1,24,2)
     1121      !    k_mass(ilon,24,2)=0.04
     1122      !    print*,'level',24,'coeff mass: k_mass',k_mass(1,24,2)
     1123
     1124
     1125      !    r_median(ilon,25,2)=1.4e-6
     1126      !    print*,'level',25,'r r_median',r_median(1,25,2)
     1127      !    stddev(ilon,25,2)=1.16
     1128      !    print*,'level',25,'dev std',stddev(1,25,2)
     1129      !    k_mass(ilon,25,2)=0.04
     1130      !    print*,'level',25,'coeff mass: k_mass',k_mass(1,25,2)
     1131
     1132
     1133      !    r_median(ilon,26,2)=1.15e-6
     1134      !    print*,'level',26,'r r_median',r_median(1,26,2)
     1135      !    stddev(ilon,26,2)=1.34
     1136      !    print*,'level',26,'dev std',stddev(1,26,2)
     1137      !    k_mass(ilon,26,2)=0.825
     1138      !    print*,'level',26,'coeff mass: k_mass',k_mass(1,26,2)
     1139
     1140
     1141      !    r_median(ilon,27,2)=1.14e-6
     1142      !    print*,'level',27,'r r_median',r_median(1,27,2)
     1143      !    stddev(ilon,27,2)=1.33
     1144      !    print*,'level',27,'dev std',stddev(1,27,2)
     1145      !    k_mass(ilon,27,2)=0.75
     1146      !    print*,'level',27,'coeff mass: k_mass',k_mass(1,27,2)
     1147
     1148
     1149      !    r_median(ilon,28,2)=1.35e-6
     1150      !    print*,'level',28,'r r_median',r_median(1,28,2)
     1151      !    stddev(ilon,28,2)=1.32
     1152      !    print*,'level',28,'dev std',stddev(1,28,2)
     1153      !    k_mass(ilon,28,2)=0.75
     1154      !    print*,'level',28,'coeff mass: k_mass',k_mass(1,28,2)
     1155
     1156
     1157      !    r_median(ilon,29,2)=1.125e-6
     1158      !    print*,'level',29,'r r_median',r_median(1,29,2)
     1159      !    stddev(ilon,29,2)=1.31
     1160      !    print*,'level',29,'dev std',stddev(1,29,2)
     1161      !    k_mass(ilon,29,2)=0.75
     1162      !    print*,'level',29,'coeff mass: k_mass',k_mass(1,29,2)
     1163
     1164
     1165      !    r_median(ilon,30,2)=1.118e-6
     1166      !    print*,'level',30,'r r_median',r_median(1,30,2)
     1167      !    stddev(ilon,30,2)=1.30
     1168      !    print*,'level',30,'dev std',stddev(1,30,2)
     1169      !    k_mass(ilon,30,2)=0.75
     1170      !    print*,'level',30,'coeff mass: k_mass',k_mass(1,30,2)
     1171
     1172
     1173      !    r_median(ilon,31,2)=1.11e-6
     1174      !    print*,'level',31,'r r_median',r_median(1,31,2)
     1175      !    stddev(ilon,31,2)=1.29
     1176      !    print*,'level',31,'dev std',stddev(1,31,2)
     1177      !    k_mass(ilon,31,2)=0.75
     1178      !    print*,'level',31,'coeff mass: k_mass',k_mass(1,31,2)
     1179
     1180      !    do i_lev=32,33
     1181      !       r_median(ilon,i_lev,2)=1.1e-6
     1182      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     1183      !       stddev(ilon,i_lev,2)=1.28
     1184      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     1185      !       k_mass(ilon,i_lev,2)=0.75
     1186      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     1187      !    end do
     1188
     1189      !    ! if k_mass > 0 it means we have a bimodal upper haze.
     1190      !    do i_lev=34,cloudmax
     1191      !       r_median(ilon,i_lev,2)=1.1e-6
     1192      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     1193      !       stddev(ilon,i_lev,2)=1.28
     1194      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     1195      !       k_mass(ilon,i_lev,2)=0.75
     1196      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     1197      !    end do
     1198
     1199      !    ! mode 3
     1200      !    do i_lev=cloudmin,21
     1201      !       r_median(ilon,i_lev,3)=3.65e-6
     1202      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     1203      !       stddev(ilon,i_lev,3)=1.28
     1204      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     1205      !       k_mass(ilon,i_lev,3)=0.0
     1206      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     1207      !    end do
     1208
     1209      !    do i_lev=22,25
     1210      !       r_median(ilon,i_lev,3)=3.65e-6
     1211      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     1212      !       stddev(ilon,i_lev,3)=1.28
     1213      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     1214      !       k_mass(ilon,i_lev,3)=0.95
     1215      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     1216      !    end do
     1217
     1218      !    do i_lev=26,cloudmax
     1219      !       r_median(ilon,i_lev,3)=3.65e-6
     1220      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     1221      !       stddev(ilon,i_lev,3)=1.28
     1222      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     1223      !       k_mass(ilon,i_lev,3)=0.0
     1224      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     1225      !    end do
     1226      ! end do
     1227
     1228      ! print*,'r_median size',size(r_median)
     1229
     1230      ! print*,'r_median(1,28,3)',r_median(1,28,3)
     1231      ! print*,'r_median(1,28,3) ne devrait pas exister puisque les dim sont (285,78,2)'
     1232      ! print*,'r_median(1,28,3) a la valeur de k_mass(1,28,1)', k_mass(1,28,1)
     1233      ! print*,'r_median(1,28,1) just to see',r_median(1,28,1)
     1234      ! print*,'r_median size after adding some shenanigans',size(r_median)
     1235
     1236      ! ! ==================================================================
     1237      ! ! initialisation stupid with mode 3 100 um
     1238      ! ! ==================================================================
     1239      ! ! mode 3 0% solid, fully liquid
     1240      ! qrad=0.0
     1241      ! ! normally nb_mode=3 in physiq.def !!!
     1242      ! do ilon=1,nbr_lon
     1243      !    ! mode 1
     1244      !    do i_lev=cloudmin,20
     1245      !       r_median(ilon,i_lev,1)=0.125e-6
     1246      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     1247      !       stddev(ilon,i_lev,1)=1.57
     1248      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     1249      !       k_mass(ilon,i_lev,1)=1.0
     1250      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     1251      !    end do
     1252
     1253      !    r_median(ilon,21,1)=0.125e-6
     1254      !    print*,'level',21,'r r_median',r_median(1,21,1)
     1255      !    stddev(ilon,21,1)=1.57
     1256      !    print*,'level',21,'dev std',stddev(1,21,1)
     1257      !    k_mass(ilon,21,1)=1.0
     1258      !    print*,'level',21,'coeff mass: k_mass',k_mass(1,21,1)
     1259
     1260      !    r_median(ilon,22,1)=0.2e-6
     1261      !    print*,'level',22,'r r_median',r_median(1,22,1)
     1262      !    stddev(ilon,22,1)=1.8
     1263      !    print*,'level',22,'dev std',stddev(1,22,1)
     1264      !    k_mass(ilon,22,1)=0.0
     1265      !    print*,'level',22,'coeff mass: k_mass',k_mass(1,22,1)
     1266
     1267      !    r_median(ilon,23,1)=0.15e-6
     1268      !    print*,'level',23,'r r_median',r_median(1,23,1)
     1269      !    stddev(ilon,23,1)=1.8
     1270      !    print*,'level',23,'dev std',stddev(1,23,1)
     1271      !    k_mass(ilon,23,1)=0.0
     1272      !    print*,'level',23,'coeff mass: k_mass',k_mass(1,23,1)
     1273
     1274      !    do i_lev=24,25
     1275      !       r_median(ilon,i_lev,1)=0.15e-6
     1276      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     1277      !       stddev(ilon,i_lev,1)=1.9
     1278      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     1279      !       k_mass(ilon,i_lev,1)=0.0
     1280      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     1281      !    end do
     1282
     1283      !    r_median(ilon,26,1)=0.175e-6
     1284      !    print*,'level',26,'r r_median',r_median(1,26,1)
     1285      !    stddev(ilon,26,1)=2.16
     1286      !    print*,'level',26,'dev std',stddev(1,26,1)
     1287      !    k_mass(ilon,26,1)=0.0
     1288      !    print*,'level',26,'coeff mass: k_mass',k_mass(1,26,1)
     1289
     1290      !    do i_lev=27,33
     1291      !       r_median(ilon,i_lev,1)=0.175e-6
     1292      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     1293      !       stddev(ilon,i_lev,1)=2.16
     1294      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     1295      !       k_mass(ilon,i_lev,1)=0.0
     1296      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     1297      !    end do
     1298
     1299      !    do i_lev=34,cloudmax
     1300      !       r_median(ilon,i_lev,1)=0.175e-6
     1301      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,1)
     1302      !       stddev(ilon,i_lev,1)=2.16
     1303      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,1)
     1304      !       k_mass(ilon,i_lev,1)=0.0
     1305      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,1)
     1306      !    end do
     1307
     1308      !    !  mode 2
     1309      !    do i_lev=cloudmin,20
     1310      !       r_median(ilon,i_lev,2)=1.4e-6
     1311      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     1312      !       stddev(ilon,i_lev,2)=1.35
     1313      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     1314      !       k_mass(ilon,i_lev,2)=0.0
     1315      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     1316      !    end do
     1317
     1318      !    r_median(ilon,21,2)=1.4e-6
     1319      !    print*,'level',21,'r r_median',r_median(1,21,2)
     1320      !    stddev(ilon,21,2)=1.25
     1321      !    print*,'level',21,'dev std',stddev(1,21,2)
     1322      !    k_mass(ilon,21,2)=0.0
     1323      !    print*,'level',21,'coeff mass: k_mass',k_mass(1,21,2)
     1324
     1325      !    r_median(ilon,22,2)=1.4e-6
     1326      !    print*,'level',22,'r r_median',r_median(1,22,2)
     1327      !    stddev(ilon,22,2)=1.25
     1328      !    print*,'level',22,'dev std',stddev(1,22,2)
     1329      !    k_mass(ilon,22,2)=0.0
     1330      !    print*,'level',22,'coeff mass: k_mass',k_mass(1,22,2)
     1331
     1332      !    r_median(ilon,23,2)=1.35e-6
     1333      !    print*,'level',23,'r r_median',r_median(1,23,2)
     1334      !    stddev(ilon,23,2)=1.25
     1335      !    print*,'level',23,'dev std',stddev(1,23,2)
     1336      !    k_mass(ilon,23,2)=0.0
     1337      !    print*,'level',23,'coeff mass: k_mass',k_mass(1,23,2)
     1338
     1339
     1340      !    r_median(ilon,24,2)=1.375e-6
     1341      !    print*,'level',24,'r r_median',r_median(1,24,2)
     1342      !    stddev(ilon,24,2)=1.2
     1343      !    print*,'level',24,'dev std',stddev(1,24,2)
     1344      !    k_mass(ilon,24,2)=0.0
     1345      !    print*,'level',24,'coeff mass: k_mass',k_mass(1,24,2)
     1346
     1347
     1348      !    r_median(ilon,25,2)=1.4e-6
     1349      !    print*,'level',25,'r r_median',r_median(1,25,2)
     1350      !    stddev(ilon,25,2)=1.16
     1351      !    print*,'level',25,'dev std',stddev(1,25,2)
     1352      !    k_mass(ilon,25,2)=0.0
     1353      !    print*,'level',25,'coeff mass: k_mass',k_mass(1,25,2)
     1354
     1355
     1356      !    r_median(ilon,26,2)=1.15e-6
     1357      !    print*,'level',26,'r r_median',r_median(1,26,2)
     1358      !    stddev(ilon,26,2)=1.34
     1359      !    print*,'level',26,'dev std',stddev(1,26,2)
     1360      !    k_mass(ilon,26,2)=0.0
     1361      !    print*,'level',26,'coeff mass: k_mass',k_mass(1,26,2)
     1362
     1363
     1364      !    r_median(ilon,27,2)=1.14e-6
     1365      !    print*,'level',27,'r r_median',r_median(1,27,2)
     1366      !    stddev(ilon,27,2)=1.33
     1367      !    print*,'level',27,'dev std',stddev(1,27,2)
     1368      !    k_mass(ilon,27,2)=0.0
     1369      !    print*,'level',27,'coeff mass: k_mass',k_mass(1,27,2)
     1370
     1371
     1372      !    r_median(ilon,28,2)=1.35e-6
     1373      !    print*,'level',28,'r r_median',r_median(1,28,2)
     1374      !    stddev(ilon,28,2)=1.32
     1375      !    print*,'level',28,'dev std',stddev(1,28,2)
     1376      !    k_mass(ilon,28,2)=0.0
     1377      !    print*,'level',28,'coeff mass: k_mass',k_mass(1,28,2)
     1378
     1379
     1380      !    r_median(ilon,29,2)=1.125e-6
     1381      !    print*,'level',29,'r r_median',r_median(1,29,2)
     1382      !    stddev(ilon,29,2)=1.31
     1383      !    print*,'level',29,'dev std',stddev(1,29,2)
     1384      !    k_mass(ilon,29,2)=0.0
     1385      !    print*,'level',29,'coeff mass: k_mass',k_mass(1,29,2)
     1386
     1387
     1388      !    r_median(ilon,30,2)=1.118e-6
     1389      !    print*,'level',30,'r r_median',r_median(1,30,2)
     1390      !    stddev(ilon,30,2)=1.30
     1391      !    print*,'level',30,'dev std',stddev(1,30,2)
     1392      !    k_mass(ilon,30,2)=0.0
     1393      !    print*,'level',30,'coeff mass: k_mass',k_mass(1,30,2)
     1394
     1395
     1396      !    r_median(ilon,31,2)=1.11e-6
     1397      !    print*,'level',31,'r r_median',r_median(1,31,2)
     1398      !    stddev(ilon,31,2)=1.29
     1399      !    print*,'level',31,'dev std',stddev(1,31,2)
     1400      !    k_mass(ilon,31,2)=0.0
     1401      !    print*,'level',31,'coeff mass: k_mass',k_mass(1,31,2)
     1402
     1403      !    do i_lev=32,33
     1404      !       r_median(ilon,i_lev,2)=1.1e-6
     1405      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     1406      !       stddev(ilon,i_lev,2)=1.28
     1407      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     1408      !       k_mass(ilon,i_lev,2)=0.0
     1409      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     1410      !    end do
     1411
     1412      !    ! if k_mass > 0 it means we have a bimodal upper haze.
     1413      !    do i_lev=34,cloudmax
     1414      !       r_median(ilon,i_lev,2)=1.1e-6
     1415      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,2)
     1416      !       stddev(ilon,i_lev,2)=1.28
     1417      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,2)
     1418      !       k_mass(ilon,i_lev,2)=0.0
     1419      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,2)
     1420      !    end do
     1421
     1422      !    ! mode 3
     1423      !    do i_lev=cloudmin,21
     1424      !       r_median(ilon,i_lev,3)=100.0e-6
     1425      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     1426      !       stddev(ilon,i_lev,3)=1.28
     1427      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     1428      !       k_mass(ilon,i_lev,3)=1.0
     1429      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     1430      !    end do
     1431
     1432      !    do i_lev=22,25
     1433      !       r_median(ilon,i_lev,3)=100.0e-6
     1434      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     1435      !       stddev(ilon,i_lev,3)=1.28
     1436      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     1437      !       k_mass(ilon,i_lev,3)=1.0
     1438      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     1439      !    end do
     1440
     1441      !    do i_lev=26,cloudmax
     1442      !       r_median(ilon,i_lev,3)=100.0e-6
     1443      !       print*,'level',i_lev,'r r_median',r_median(1,i_lev,3)
     1444      !       stddev(ilon,i_lev,3)=1.28
     1445      !       print*,'level',i_lev,'dev std',stddev(1,i_lev,3)
     1446      !       k_mass(ilon,i_lev,3)=1.0
     1447      !       print*,'level',i_lev,'coeff mass: k_mass',k_mass(1,i_lev,3)
     1448      !    end do
     1449      ! end do
     1450
     1451      ! print*,'r_median size',size(r_median)
     1452
     1453      ! print*,'r_median(1,28,3)',r_median(1,28,3)
     1454      ! print*,'r_median(1,28,3) ne devrait pas exister puisque les dim sont (285,78,2)'
     1455      ! print*,'r_median(1,28,3) a la valeur de k_mass(1,28,1)', k_mass(1,28,1)
     1456      ! print*,'r_median(1,28,1) just to see',r_median(1,28,1)
     1457      ! print*,'r_median size after adding some shenanigans',size(r_median)
     1458
     1459! check if sum of kmass = 1 at each level
     1460
     1461print*, 'chemparam_mod: start checking k_mass=1'
     1462
     1463do i_lev = cloudmin,cloudmax
     1464   if (nbr_mode == 3) then
     1465      if ((k_mass(1,i_lev,1) + k_mass(1,i_lev,2) + k_mass(1,i_lev,3)) /= 1.0) then
     1466         print*, 'kmass total is not 1.0'
     1467         print*, k_mass(1,i_lev,1) + k_mass(1,i_lev,2) + k_mass(1,i_lev,3)
     1468         print*, 'at level: ',i_lev
     1469         print*, 'check chemparam_mod cloud structure definition'
     1470         stop
     1471      end if
     1472   else
     1473      if ((k_mass(1,i_lev,1) + k_mass(1,i_lev,2)) /= 1.0) then
     1474         print*, 'kmass total is not 1.0'
     1475         print*, k_mass(1,i_lev,1) + k_mass(1,i_lev,2)
     1476         print*, 'at level: ',i_lev
     1477         print*, 'check chemparam_mod cloud structure definition'
     1478      end if
     1479   end if
     1480end do
     1481
     1482print*, 'chemparam_mod: end checking k_mass=1'
     1483print*, 'k_mass is fine'
     1484print*
     1485print*,'==============================='
     1486print*,'end initialisation cloud layer'
     1487print*,'==============================='
     1488
     1489end subroutine cloud_ini
     1490
     1491!============================================================================
     1492
     1493  subroutine chemparam_ini
     1494
     1495!============================================================================
     1496
     1497      use infotrac_phy, only: nqtot, tname
     1498
     1499      implicit none
     1500
     1501      integer :: i
     1502
     1503      allocate(m_tr(nqtot))    ! molecular mass of tracers
     1504      allocate(type_tr(nqtot)) ! type of chemical tracers 1: neutral, 2: ion, 3: liquid
     1505
     1506      ! index for chemical tracers
     1507
     1508      ! neutrals
     1509
     1510      i_co2      = 0
     1511      i_co       = 0
     1512      i_h2       = 0
     1513      i_h2o      = 0
     1514      i_o1d      = 0
     1515      i_o        = 0
     1516      i_o2       = 0
     1517      i_o2dg     = 0
     1518      i_o3       = 0
     1519      i_h        = 0
     1520      i_oh       = 0
     1521      i_ho2      = 0
     1522      i_h2o2     = 0
     1523      i_cl       = 0
     1524      i_clo      = 0
     1525      i_cl2      = 0
     1526      i_hcl      = 0
     1527      i_hocl     = 0
     1528      i_clco     = 0
     1529      i_clco3    = 0
     1530      i_cocl2    = 0
     1531      i_s        = 0
     1532      i_so       = 0
     1533      i_so2      = 0
     1534      i_so3      = 0
     1535      i_s2o2     = 0
     1536      i_ocs      = 0
     1537      i_hso3     = 0
     1538      i_h2so4    = 0
     1539      i_s2       = 0
     1540      i_clso2    = 0
     1541      i_oscl     = 0
     1542      i_n2       = 0
     1543      i_he       = 0
     1544      i_n2d      = 0
     1545      i_n        = 0
     1546      i_no       = 0
     1547      i_no2      = 0
     1548
     1549      ! ions
     1550
     1551      i_co2plus  = 0
     1552      i_coplus   = 0
     1553      i_oplus    = 0
     1554      i_o2plus   = 0
     1555      i_n2plus   = 0
     1556      i_hplus    = 0
     1557      i_h2oplus  = 0
     1558      i_nplus    = 0
     1559      i_ohplus   = 0
     1560      i_cplus    = 0
     1561      i_noplus   = 0
     1562      i_h3oplus  = 0
     1563      i_hcoplus  = 0
     1564      i_hco2plus = 0
     1565      i_elec     = 0
     1566
     1567      ! liquid
     1568
     1569      i_h2oliq   = 0
     1570      i_h2so4liq = 0
     1571
     1572      do i = 1,nqtot
     1573         print*,'tname(i)',tname(i)
     1574         select case(tname(i))
     1575
     1576            ! neutrals
     1577
     1578            case('co2')
     1579               i_co2 = i
     1580               print*,'co2',i_co2
     1581               m_tr(i_co2) = 44.0095
     1582               type_tr(i_co2) = 1
     1583            case('co')
     1584               i_co = i
     1585               print*,'co',i_co
     1586               m_tr(i_co) = 28.0101
     1587               type_tr(i_co) = 1
     1588            case('h2')
     1589               i_h2 = i
     1590               print*,'h2',i_h2
     1591               m_tr(i_h2) = 2.01588
     1592               type_tr(i_h2) = 1
     1593            case('h2o')
     1594               i_h2o = i
     1595               print*,'h2o',i_h2o
     1596               m_tr(i_h2o) = 18.0153
     1597               type_tr(i_h2o) = 1
     1598            case('o1d')
     1599               i_o1d = i
     1600               print*,'o1d',i_o1d
     1601               m_tr(i_o1d) = 15.994
     1602               type_tr(i_o1d) = 1
     1603            case('o')
     1604               i_o = i
     1605               print*,'o',i_o
     1606               m_tr(i_o) = 15.994
     1607               type_tr(i_o) = 1
     1608            case('o2')
     1609               i_o2 = i
     1610               print*,'o2',i_o2
     1611               m_tr(i_o2) = 31.9988
     1612               type_tr(i_o2) = 1
     1613            case('o2dg')
     1614               i_o2dg = i
     1615               print*,'o2dg',i_o2dg
     1616               m_tr(i_o2dg) = 31.9988
     1617               type_tr(i_o2dg) = 1
     1618            case('o3')
     1619               i_o3 = i
     1620               print*,'o3',i_o3
     1621               m_tr(i_o3) = 47.9982
     1622               type_tr(i_o3) = 1
     1623            case('h')
     1624               i_h = i
     1625               print*,'h',i_h
     1626               m_tr(i_h) = 1.00794
     1627               type_tr(i_h) = 1
     1628            case('oh')
     1629               i_oh = i
     1630               print*,'oh',i_oh
     1631               m_tr(i_oh) = 17.0073
     1632               type_tr(i_oh) = 1
     1633            case('ho2')
     1634               i_ho2 = i
     1635               print*,'ho2',i_ho2
     1636               m_tr(i_ho2) = 33.0067
     1637               type_tr(i_ho2) = 1
     1638            case('h2o2')
     1639               i_h2o2 = i
     1640               print*,'h2o2',i_h2o2
     1641               m_tr(i_h2o2) = 34.0147
     1642               type_tr(i_h2o2) = 1
     1643            case('cl')
     1644               i_cl = i
     1645               print*,'cl',i_cl
     1646               m_tr(i_cl) = 35.453
     1647               type_tr(i_cl) = 1
     1648            case('clo')
     1649               i_clo = i
     1650               print*,'clo',i_clo
     1651               m_tr(i_clo) = 51.452
     1652               type_tr(i_clo) = 1
     1653            case('cl2')
     1654               i_cl2 = i
     1655               print*,'cl2',i_cl2
     1656               m_tr(i_cl2) = 70.906
     1657               type_tr(i_cl2) = 1
     1658            case('hcl')
     1659               i_hcl = i
     1660               print*,'hcl',i_hcl
     1661               m_tr(i_hcl) = 36.461
     1662               type_tr(i_hcl) = 1
     1663            case('hocl')
     1664               i_hocl = i
     1665               print*,'hocl',i_hocl
     1666               m_tr(i_hocl) = 52.46
     1667               type_tr(i_hocl) = 1
     1668            case('clco')
     1669               i_clco = i
     1670               print*,'clco',i_clco
     1671               m_tr(i_clco) = 63.463
     1672               type_tr(i_clco) = 1
     1673            case('clco3')
     1674               i_clco3 = i
     1675               print*,'clco3',i_clco3
     1676               m_tr(i_clco3) = 95.462
     1677               type_tr(i_clco3) = 1
     1678            case('cocl2')
     1679               i_cocl2 = i
     1680               print*,'cocl2',i_cocl2
     1681               m_tr(i_cocl2) = 98.916
     1682               type_tr(i_cocl2) = 1
     1683            case('s')
     1684               i_s = i
     1685               print*,'s',i_s
     1686               m_tr(i_s) = 32.065
     1687               type_tr(i_s) = 1
     1688            case('so')
     1689               i_so = i
     1690               print*,'so',i_so
     1691               m_tr(i_so) = 48.0644
     1692               type_tr(i_so) = 1
     1693            case('so2')
     1694               i_so2 = i
     1695               print*,'so2',i_so2
     1696               m_tr(i_so2) = 64.064
     1697               type_tr(i_so2) = 1
     1698            case('so3')
     1699               i_so3 = i
     1700               print*,'so3',i_so3
     1701               m_tr(i_so3) = 80.063
     1702               type_tr(i_so3) = 1
     1703            case('s2o2')
     1704               i_s2o2 = i
     1705               print*,'s2o2',i_s2o2
     1706               m_tr(i_s2o2)= 96.1288
     1707               type_tr(i_s2o2) = 1
     1708            case('ocs')
     1709               i_ocs = i
     1710               print*,'ocs',i_ocs
     1711               m_tr(i_ocs) = 60.0751
     1712               type_tr(i_ocs) = 1
     1713            case('hso3')
     1714               i_hso3 = i
     1715               print*,'hso3',i_hso3
     1716               m_tr(i_hso3) = 81.071
     1717               type_tr(i_hso3) = 1
     1718            case('h2so4')
     1719               i_h2so4 = i
     1720               print*,'h2so4',i_h2so4
     1721               m_tr(i_h2so4) = 98.078
     1722               type_tr(i_h2so4) = 1
     1723            case('s2')
     1724               i_s2 = i
     1725               print*,'s2',i_s2
     1726               m_tr(i_s2) = 64.13
     1727               type_tr(i_s2) = 1
     1728            case('clso2')
     1729               i_clso2 = i
     1730               print*,'clso2',i_clso2
     1731               m_tr(i_clso2) = 99.517
     1732               type_tr(i_clso2) = 1
     1733            case('oscl')
     1734               i_oscl = i
     1735               print*,'oscl',i_oscl
     1736               m_tr(i_oscl) = 83.517
     1737               type_tr(i_oscl) = 1
     1738            case('n2')
     1739               i_n2 = i
     1740               print*,'n2',i_n2
     1741               m_tr(i_n2) = 28.013
     1742               type_tr(i_n2) = 1
     1743            case('he')
     1744               i_he = i
     1745               print*,'he',i_he
     1746               m_tr(i_he) = 4.0026
     1747               type_tr(i_he) = 1
     1748            case('n2d')
     1749               i_n2d = i
     1750               print*,'n2d',i_n2d
     1751               m_tr(i_n2d) = 14.0067
     1752               type_tr(i_n2d) = 1
     1753            case('n')
     1754               i_n = i
     1755               print*,'n',i_n
     1756               m_tr(i_n) = 14.0067
     1757               type_tr(i_n) = 1
     1758            case('no')
     1759               i_no = i
     1760               print*,'no',i_no
     1761               m_tr(i_no) = 30.0061
     1762               type_tr(i_no) = 1
     1763            case('no2')
     1764               i_no2 = i
     1765               print*,'no2',i_no2
     1766               m_tr(i_no2) = 46.0055
     1767               type_tr(i_no2) = 1
     1768
     1769            ! ions
     1770
     1771            case('co2plus')
     1772               i_co2plus = i
     1773               print*,'co2plus',i_co2plus
     1774               m_tr(i_co2plus) = 44.0095
     1775               type_tr(i_co2plus) = 2
     1776            case('coplus')
     1777               i_coplus = i
     1778               print*,'coplus',i_coplus
     1779               m_tr(i_coplus) = 28.0101
     1780               type_tr(i_coplus) = 2
     1781            case('oplus')
     1782               i_oplus = i
     1783               print*,'oplus',i_oplus
     1784               m_tr(i_oplus) = 15.994
     1785               type_tr(i_oplus) = 2
     1786            case('o2plus')
     1787               i_o2plus = i
     1788               print*,'o2plus',i_o2plus
     1789               m_tr(i_o2plus) = 31.9988
     1790               type_tr(i_o2plus) = 2
     1791            case('n2plus')
     1792               i_n2plus = i
     1793               print*,'n2plus',i_n2plus
     1794               m_tr(i_n2plus) = 28.013
     1795               type_tr(i_n2plus) = 2
     1796            case('hplus')
     1797               i_hplus = i
     1798               print*,'hplus',i_hplus
     1799               m_tr(i_hplus) = 1.00794
     1800               type_tr(i_hplus) = 2
     1801            case('h2oplus')
     1802               i_h2oplus = i
     1803               print*,'h2oplus',i_h2oplus
     1804               m_tr(i_h2oplus) = 18.0153
     1805               type_tr(i_h2oplus) = 2
     1806            case('nplus')
     1807               i_nplus = i
     1808               print*,'nplus',i_nplus
     1809               m_tr(i_nplus) = 14.0067
     1810               type_tr(i_nplus) = 2
     1811            case('ohplus')
     1812               i_ohplus = i
     1813               print*,'ohplus',i_ohplus
     1814               m_tr(i_ohplus) = 17.0073
     1815               type_tr(i_ohplus) = 2
     1816            case('cplus')
     1817               i_cplus = i
     1818               print*,'cplus',i_cplus
     1819               m_tr(i_cplus) = 12.011
     1820               type_tr(i_cplus) = 2
     1821            case('noplus')
     1822               i_noplus = i
     1823               print*,'noplus',i_noplus
     1824               m_tr(i_noplus) = 30.0061
     1825               type_tr(i_noplus) = 2
     1826            case('h3oplus')
     1827               i_h3oplus = i
     1828               print*,'h3oplus',i_h3oplus
     1829               m_tr(i_h3oplus) = 19.0232
     1830               type_tr(i_h3oplus) = 2
     1831            case('hcoplus')
     1832               i_hcoplus = i
     1833               print*,'hcoplus',i_hcoplus
     1834               m_tr(i_hcoplus) = 29.0180
     1835               type_tr(i_hcoplus) = 2
     1836            case('hco2plus')
     1837               i_hco2plus = i
     1838               print*,'hco2plus',i_hco2plus
     1839               m_tr(i_hco2plus) = 45.
     1840               type_tr(i_hco2plus) = 2
     1841            case('elec')
     1842               i_elec = i
     1843               print*,'elec',i_elec
     1844               m_tr(i_elec) = 1./1822.89
     1845               type_tr(i_elec) = 2
     1846
     1847            ! liquid tracers (cl_scheme = 1)
     1848
     1849            case('h2oliq')
     1850               i_h2oliq = i
     1851               print*,'h2oliq',i_h2oliq
     1852               m_tr(i_h2oliq) = 18.0153
     1853               type_tr(i_h2oliq) = 3
     1854            case('h2so4liq')
     1855               i_h2so4liq = i
     1856               print*,'h2so4liq',i_h2so4liq
     1857               m_tr(i_h2so4liq) = 98.078
     1858               type_tr(i_h2so4liq) = 3
     1859
     1860            ! liquid tracers (cl_scheme = 2)
     1861
     1862            case('m0_aer')
     1863               i_m0_aer = i
     1864               print*,'m0_aer',i_m0_aer
     1865               type_tr(i_m0_aer) = 10
     1866            case('m3_aer')
     1867               i_m3_aer = i
     1868               print*,'m3_aer',i_m3_aer
     1869               type_tr(i_m3_aer) = 10
     1870            case('m0_m1drop')
     1871               i_m0_mode1drop = i
     1872               print*,'m0_m1drop',i_m0_mode1drop
     1873               type_tr(i_m0_mode1drop) = 10
     1874            case('m0_m1ccn')
     1875               i_m0_mode1ccn = i
     1876               print*,'m0_m1ccn',i_m0_mode1ccn
     1877               type_tr(i_m0_mode1ccn) = 10
     1878            case('m3_m1sa')
     1879               i_m3_mode1sa = i
     1880               print*,'m3_m1sa',i_m3_mode1sa
     1881               type_tr(i_m3_mode1sa) = 10
     1882            case('m3_m1w')
     1883               i_m3_mode1w = i
     1884               print*,'m3_m1w',i_m3_mode1w
     1885               type_tr(i_m3_mode1w) = 10
     1886            case('m3_m1ccn')
     1887               i_m3_mode1ccn = i
     1888               print*,'m3_m1ccn',i_m3_mode1ccn
     1889               type_tr(i_m3_mode1ccn) = 10
     1890            case('m0_m2drop')
     1891               i_m0_mode2drop = i
     1892               print*,'m0_m2drop',i_m0_mode2drop
     1893               type_tr(i_m0_mode2drop) = 10
     1894            case('m0_m2ccn')
     1895               i_m0_mode2ccn = i
     1896               print*,'m0_m2ccn',i_m0_mode2ccn
     1897               type_tr(i_m0_mode2ccn) = 10
     1898            case('m3_m2sa')
     1899               i_m3_mode2sa = i
     1900               print*,'m3_m2sa',i_m3_mode2sa
     1901               type_tr(i_m3_mode2sa) = 10
     1902            case('m3_m2w')
     1903               i_m3_mode2w = i
     1904               print*,'m3_m2w',i_m3_mode2w
     1905               type_tr(i_m3_mode2w) = 10
     1906            case('m3_m2ccn')
     1907               i_m3_mode2ccn = i
     1908               print*,'m3_m2ccn',i_m3_mode2ccn
     1909               type_tr(i_m3_mode2ccn) = 10
     1910         end select
     1911      end do
     1912
     1913  end subroutine chemparam_ini
     1914
     1915!============================================================================
     1916
     1917  subroutine vapors4muphy_ini(nlon,nlev,trac)
     1918
     1919!============================================================================
     1920
     1921  use infotrac_phy, only: nqtot, tname
    11131922
    11141923  integer :: nlon, nlev
    … …  
    11171926!  integer :: i
    11181927!  real    :: trac1d(nlev,2) ! traceur lu ( en vmr)
    1119  
    1120 ! lecture d'un fichier texte contenant les profils de trac1d(:1) = H2O et trac1d(:,2) = H2SO4
    1121 !  DO i=1,nlon
     1928
     1929! lecture d'un fichier texte contenant les profils de trac1d(:1) = h2o et trac1d(:,2) = h2so4
     1930!  do i=1,nlon
    11221931!     trac(i,:,i_h2o) = trac1d(:,1)
    11231932!     trac(i,:,i_h2so4) = trac1d(:,2)
    1124 !  ENDDO
    1125 
    1126 !  intitialisation profils altitude H2O et H2SO4
    1127 !  profil H2O initial vap+liq == que vap
    1128    trac(:,1:24,i_h2o) = 30.E-6 !
    1129    trac(:,25:50,i_h2o) = 1.E-6 !
    1130 
    1131    trac(:,:,i_h2so4) = 3.E-9 ! Limite sup Sandor 2012
    1132    trac(:,23:50,i_h2so4) = 2.E-6 ! Profil H2SO4 initial => vap+liq
    1133 
    1134   END SUBROUTINE vapors4muphy_ini
    1135 
    1136 END MODULE chemparam_mod
    1137 
     1933!  enddo
     1934
     1935!  intitialisation profils altitude h2o et h2so4
     1936!  profil h2o initial vap+liq == que vap
     1937   trac(:,1:24,i_h2o) = 30.e-6 !
     1938   trac(:,25:50,i_h2o) = 1.e-6 !
     1939
     1940   trac(:,:,i_h2so4) = 3.e-9 ! limite sup sandor 2012
     1941   trac(:,23:50,i_h2so4) = 2.e-6 ! profil h2so4 initial => vap+liq
     1942
     1943  end subroutine vapors4muphy_ini
     1944
     1945end module chemparam_mod

  • trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/new_cloud_sedim.F

    r2464 r3323  
    1       subroutine new_cloud_sedim(n_lon, n_lev, ptimestep,
     1      subroutine new_cloud_sedim(nbr_mode, n_lon, n_lev, ptimestep,
    22     $                           pmidlay, pbndlay, pt, pq,
    33     $                           d_tr_chem, pdqsed,
    … …  
    2323      INTEGER n_lon                 ! number of horizontal grid points
    2424      INTEGER n_lev                 ! number of atmospheric layers
     25      integer nbr_mode
    2526      REAL ptimestep                ! physics time step (s)
    2627      REAL pmidlay(n_lon,n_lev)     ! pressure at middle layers (Pa)
    … …  
    5657      real,parameter :: molrad=2.2e-10   ! CO2
    5758     
    58 c     Ratio radius shell model du mode 3
    59 c     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    60 c       Ce ratio correspond aux mesures effectuées par J. Cimino (1982), Icarus
    61 c     Fixer ce parametre a 0 revient a une gouttelette pure en liquide acide sulfurique
    62 c     ATTENTION ! DOIT ETRE COHERENT AVEC new_cloud_venus !
    63       REAL, PARAMETER :: qrad = 0.97
    6459      REAL :: qmass
    6560c       masse volumique du coeur (kg.m-3)

  • trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/new_cloud_venus.F

    r1687 r3323  
    1111!*
    1212!----------------------------------------------------------------------------
    13       SUBROUTINE new_cloud_venus(
     13      SUBROUTINE new_cloud_venus(nbr_mode,
    1414     + nblev, nblon,
    1515     + TT,PP,
    … …  
    2222#include "YOMCST.h"
    2323
     24      integer :: nbr_mode
    2425      INTEGER, INTENT(IN) :: nblon  ! nombre de points horizontaux
    2526      INTEGER, INTENT(IN) :: nblev  ! nombre de couches verticales
    … …  
    5758
    5859      INTEGER :: FLAG
    59 !----------------------------------------------------------------------------
    60 c     Ratio radius shell model du mode 3
    61 c     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    62 c       Ce ratio correspond aux mesures effectuées par J. Cimino (1982), Icarus
    63 c     Fixer ce parametre a 0 revient a une gouttelette pure en liquide acide sulfurique
    64 c     ATTENTION ! DOIT ETRE COHERENT AVEC new_cloud_sedim !
    65       REAL, PARAMETER :: qrad = 0.97
     60     
    6661      REAL :: qmass
    6762c       masse volumique du coeur (kg.m-3)

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

    r3035 r3323  
    746746
    747747      if (ok_cloud .and. (cl_scheme == 1)) then
    748          call cloud_ini(nlon,nlev)
     748         call cloud_ini(nlon,nlev,nb_mode)
    749749      end if
    750750
    … …  
    10961096
    10971097#ifndef MESOSCALE
    1098                call new_cloud_sedim(klon,
     1098               call new_cloud_sedim(nb_mode,
     1099     $                              klon,
    10991100     $                              nlev,
    11001101     $                              zctime,

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

    r3035 r3323  
    177177                 trac(:,:,iq) = trac(:,:,iq)*mmean(:,:)/m_tr(iq)
    178178              end do
     179
    179180              print*, "SO2 is re-initialised"
     181
    180182              if (i_so2 /= 0) then
     183                 trac(:,:,i_so2) = 0.
    181184                 trac(:,1:22,i_so2) = 10.e-6
    182185 
    … …  
    293296!     call microphysics
    294297
    295          call new_cloud_venus(nlev, nlon, temp, pplay,
     298         call new_cloud_venus(nb_mode, nlev, nlon, temp, pplay,
    296299     $                        mrtwv, mrtsa, mrwv, mrsa)
    297300
Note: See TracChangeset for help on using the changeset viewer.