Ignore:
Timestamp:
Sep 3, 2009, 2:03:33 PM (15 years ago)
Author:
Laurent Fairhead
Message:

Des modifications sur la lecture des aerosols par Michael
Correction du test sur le jour de lecture des aerosols qui ne marchait
pas avec le nouveau calendrier (a revoir?)
Menage sur quelques prints
SD/MAF

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90

    r1224 r1237  
    1 
     1!
     2! $Id$
     3!
    24
    35SUBROUTINE AEROPT_5WV(&
    … …  
    911  USE DIMPHY
    1012  USE aero_mod
    11 
     13  USE mod_phys_lmdz_para, ONLY : mpi_rank
    1214  !
    1315  !    Yves Balkanski le 12 avril 2006
    … …  
    7274  LOGICAL :: soluble
    7375 
    74   INTEGER :: i, k, m
     76  INTEGER :: i, k, ierr, m
    7577  INTEGER :: spsol, spinsol, spss, la
    7678  INTEGER :: RH_num
    … …  
    8183  INTEGER, PARAMETER :: la865 = 5
    8284  INTEGER, PARAMETER :: nbre_RH=12
    83   INTEGER, PARAMETER :: nbsol_compaer=3   !  1- BC soluble; 2- POM soluble; 3- SO4.
    84   INTEGER, PARAMETER :: nbinsol_compaer=3 !  1- Dust; 2- BC insoluble; 3- POM insoluble
     85  INTEGER, PARAMETER :: naero_soluble=7   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.
     86                                          ! 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
     87  INTEGER, PARAMETER :: naero_insoluble=3 !  1- Dust; 2- BC insoluble; 3- POM insoluble
    8588  INTEGER, PARAMETER :: nb_level = 19 ! number of vertical levels
    8689  LOGICAL, SAVE :: firstcall=.TRUE.
    … …  
    98101
    99102  ! Coefficient optiques interpole sur le nombre de niveau du modele
    100   REAL, DIMENSION(klev) :: A1_ASSSM, A2_ASSSM, A3_ASSSM,&
     103  REAL, ALLOCATABLE,  DIMENSION(:), SAVE :: &
     104          A1_ASSSM, A2_ASSSM, A3_ASSSM,&
    101105          B1_ASSSM, B2_ASSSM, C1_ASSSM, C2_ASSSM,&
    102106          A1_CSSSM, A2_CSSSM, A3_CSSSM,&
    … …  
    124128
    125129 
    126   REAL :: alpha_aers_5wv(nbre_RH,las,nbsol_compaer)   ! ext. coeff. Soluble comp. units *** m2/g
    127                                                       ! 1- Seasalt AS: 2- Sesalt CS; 3- BC; 4- POM; 5- SO4.
    128   REAL :: alpha_aeri_5wv(las,nbinsol_compaer)         ! ext. coeff. Insoluble comp. 1- Dust: 2- BC; 3- POM
    129   REAL :: cg_aers_5wv(nbre_RH,las,nbsol_compaer)      ! Asym. param. soluble comp.
    130                                                       ! 1- Seasalt AS: 2- Sesalt CS; 3- BC; 4- POM; 5- SO4.
    131   REAL :: cg_aeri_5wv(las,nbinsol_compaer)            ! Asym. param. insoluble comp. 1- Dust: 2- BC; 3- POM
    132   REAL :: piz_aers_5wv(nbre_RH,las,nbsol_compaer)   
    133                                                       ! 1- Seasalt AS: 2- Sesalt CS; 3- BC; 4- POM; 5- SO4.
    134   REAL :: piz_aeri_5wv(las,nbinsol_compaer)           ! Insoluble comp. 1- Dust: 2- BC; 3- POM
     130  REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble)   ! ext. coeff. Soluble comp. units *** m2/g
     131   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
     132  REAL :: alpha_aeri_5wv(las,naero_insoluble)         ! ext. coeff. Insoluble comp. 1- Dust: 2- BC; 3- POM
     133  REAL :: cg_aers_5wv(nbre_RH,las,naero_soluble)      ! Asym. param. soluble comp.
     134   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
     135  REAL :: cg_aeri_5wv(las,naero_insoluble)            ! Asym. param. insoluble comp. 1- Dust: 2- BC; 3- POM
     136  REAL :: piz_aers_5wv(nbre_RH,las,naero_soluble)   
     137   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
     138  REAL :: piz_aeri_5wv(las,naero_insoluble)           ! Insoluble comp. 1- Dust: 2- BC; 3- POM
    135139
    136140  REAL, DIMENSION(klon,klev,naero_spc) :: mass_temp
    … …  
    284288  !                                   le 12 AVRIL 2006
    285289  !
    286   DATA alpha_aers_5wv/ &
    287      ! bc soluble
    288      7.930,7.930,7.930,7.930,7.930,7.930,     &
    289      7.930,7.930,10.893,12.618,14.550,16.613, &
    290      7.658,7.658,7.658,7.658,7.658,7.658,     &
    291      7.658,7.658,10.351,11.879,13.642,15.510, &
    292      7.195,7.195,7.195,7.195,7.195,7.195,     &
    293      7.195,7.195,9.551,10.847,12.381,13.994,  &
    294      6.736,6.736,6.736,6.736,6.736,6.736,     &
    295      6.736,6.736,8.818,9.938,11.283,12.687,   &
    296      6.277,6.277,6.277,6.277,6.277,6.277,     &
    297      6.277,6.277,8.123,9.094,10.275,11.501,   &
    298      ! pom soluble
    299      6.676,6.676,6.676,6.676,6.710,6.934,   &
    300      7.141,7.569,8.034,8.529,9.456,10.511,  &
    301      5.109,5.109,5.109,5.109,5.189,5.535,   &
    302      5.960,6.852,8.008,9.712,12.897,19.676, &
    303      3.718,3.718,3.718,3.718,3.779,4.042,   &
    304      4.364,5.052,5.956,7.314,9.896,15.688,  &
    305      2.849,2.849,2.849,2.849,2.897,3.107,   &
    306      3.365,3.916,4.649,5.760,7.900,12.863,  &
    307      2.229,2.229,2.229,2.229,2.268,2.437,   &
    308      2.645,3.095,3.692,4.608,6.391,10.633,  &
    309      ! Sulfate
    310      5.751,6.215,6.690,7.024,7.599,8.195,      &
    311      9.156,10.355,12.660,14.823,18.908,24.508, &
    312      4.320,4.675,5.052,5.375,5.787,6.274,      &
    313      7.066,8.083,10.088,12.003,15.697,21.133,  &
    314      3.079,3.351,3.639,3.886,4.205,4.584,      &
    315      5.206,6.019,7.648,9.234,12.391,17.220,    &
    316      2.336,2.552,2.781,2.979,3.236,3.540,      &
    317      4.046,4.711,6.056,7.388,10.093,14.313,    &
    318      1.777,1.949,2.134,2.292,2.503,2.751,      &
    319      3.166,3.712,4.828,5.949,8.264,11.922/
     290 DATA alpha_aers_5wv/ &
     291                                ! bc soluble
     292       7.930,7.930,7.930,7.930,7.930,7.930,     &
     293       7.930,7.930,10.893,12.618,14.550,16.613, &
     294       7.658,7.658,7.658,7.658,7.658,7.658,     &
     295       7.658,7.658,10.351,11.879,13.642,15.510, &
     296       7.195,7.195,7.195,7.195,7.195,7.195,     &
     297       7.195,7.195,9.551,10.847,12.381,13.994,  &
     298       6.736,6.736,6.736,6.736,6.736,6.736,     &
     299       6.736,6.736,8.818,9.938,11.283,12.687,   &
     300       6.277,6.277,6.277,6.277,6.277,6.277,     &
     301       6.277,6.277,8.123,9.094,10.275,11.501,   &
     302                                ! pom soluble
     303       6.676,6.676,6.676,6.676,6.710,6.934,   &
     304       7.141,7.569,8.034,8.529,9.456,10.511,  &
     305       5.109,5.109,5.109,5.109,5.189,5.535,   &
     306       5.960,6.852,8.008,9.712,12.897,19.676, &
     307       3.718,3.718,3.718,3.718,3.779,4.042,   &
     308       4.364,5.052,5.956,7.314,9.896,15.688,  &
     309       2.849,2.849,2.849,2.849,2.897,3.107,   &
     310       3.365,3.916,4.649,5.760,7.900,12.863,  &
     311       2.229,2.229,2.229,2.229,2.268,2.437,   &
     312       2.645,3.095,3.692,4.608,6.391,10.633,  &
     313                                ! Sulfate (Accumulation)
     314       5.751,6.215,6.690,7.024,7.599,8.195,      &
     315       9.156,10.355,12.660,14.823,18.908,24.508, &
     316       4.320,4.675,5.052,5.375,5.787,6.274,      &
     317       7.066,8.083,10.088,12.003,15.697,21.133,  &
     318       3.079,3.351,3.639,3.886,4.205,4.584,      &
     319       5.206,6.019,7.648,9.234,12.391,17.220,    &
     320       2.336,2.552,2.781,2.979,3.236,3.540,      &
     321       4.046,4.711,6.056,7.388,10.093,14.313,    &
     322       1.777,1.949,2.134,2.292,2.503,2.751,      &
     323       3.166,3.712,4.828,5.949,8.264,11.922,     &
     324                                ! Sulfate (Coarse)
     325       5.751,6.215,6.690,7.024,7.599,8.195,      &
     326       9.156,10.355,12.660,14.823,18.908,24.508, &
     327       4.320,4.675,5.052,5.375,5.787,6.274,      &
     328       7.066,8.083,10.088,12.003,15.697,21.133,  &
     329       3.079,3.351,3.639,3.886,4.205,4.584,      &
     330       5.206,6.019,7.648,9.234,12.391,17.220,    &
     331       2.336,2.552,2.781,2.979,3.236,3.540,      &
     332       4.046,4.711,6.056,7.388,10.093,14.313,    &
     333       1.777,1.949,2.134,2.292,2.503,2.751,      &
     334       3.166,3.712,4.828,5.949,8.264,11.922,     &
     335                                ! Seasalt soluble super_coarse (computed below for 550nm)
     336       0.50,0.90,1.05,1.21,1.40,2.41, &
     337       2.66,3.11,3.88,4.52,5.69,8.84, &
     338       0.000,0.000,0.000,0.000,0.000,0.000, &
     339       0.000,0.000,0.000,0.000,0.000,0.000, &
     340     0.52,0.93,1.08,1.24,1.43,2.47, &
     341     2.73,3.20,3.99,4.64,5.84,9.04, &
     342     0.52,0.93,1.09,1.25,1.44,2.50, &
     343     2.76,3.23,4.03,4.68,5.89,9.14, &
     344     0.52,0.94,1.09,1.26,1.45,2.51, &
     345     2.78,3.25,4.06,4.72,5.94,9.22, &
     346                                ! seasalt soluble coarse (computed below for 550nm)
     347       0.50,0.90,1.05,1.21,1.40,2.41, &
     348       2.66,3.11,3.88,4.52,5.69,8.84, &
     349       0.000,0.000,0.000,0.000,0.000,0.000, &
     350       0.000,0.000,0.000,0.000,0.000,0.000, &
     351     0.52,0.93,1.08,1.24,1.43,2.47, &
     352     2.73,3.20,3.99,4.64,5.84,9.04, &
     353     0.52,0.93,1.09,1.25,1.44,2.50, &
     354     2.76,3.23,4.03,4.68,5.89,9.14, &
     355     0.52,0.94,1.09,1.26,1.45,2.51, &
     356     2.78,3.25,4.06,4.72,5.94,9.22, &
     357                                ! seasalt soluble accumulation (computed below for 550nm)
     358     4.28, 7.17, 8.44, 9.85,11.60,22.44,  &
     359     25.34,30.54,39.38,46.52,59.33,91.77, &
     360       0.000,0.000,0.000,0.000,0.000,0.000, &
     361       0.000,0.000,0.000,0.000,0.000,0.000, &
     362     2.48, 4.22, 5.02, 5.94, 7.11,15.29,  &
     363     17.70,22.31,30.73,38.06,52.15,90.59, &
     364     1.90, 3.29, 3.94, 4.69, 5.65, 12.58, &
     365     14.68,18.77,26.41,33.25,46.77,85.50, &
     366     1.47, 2.59, 3.12, 3.74, 4.54, 10.42, &
     367     12.24,15.82,22.66,28.91,41.54,79.33/
    320368
    321369  DATA alpha_aeri_5wv/ &
    322      ! dust insoluble
    323      0.759, 0.770, 0.775, 0.775, 0.772, &
    324      !!jb bc insoluble
    325      11.536,10.033, 8.422, 7.234, 6.270, &
    326      ! pom insoluble
    327      5.042, 3.101, 1.890, 1.294, 0.934/
    328 
    329   DATA cg_aers_5wv/ &
    330      ! bc soluble
    331      .651, .651, .651, .651, .651, .651, &
    332      .651, .651, .738, .764, .785, .800, &
    333      .597, .597, .597, .597, .597, .597, &
    334      .597, .597, .695, .725, .751, .770, &
    335      .543, .543, .543, .543, .543, .543, &
    336      .543, .543, .650, .684, .714, .736, &
    337      .504, .504, .504, .504, .504, .504, &
    338      .504, .504, .614, .651, .683, .708, &
    339      .469, .469, .469, .469, .469, .469, &
    340      .469, .469, .582, .620, .655, .681, &
    341      ! pom soluble
    342      .679, .679, .679, .679, .683, .691, &
    343      .703, .720, .736, .751, .766, .784, &
    344      .656, .656, .656, .656, .659, .669, &
    345      .681, .699, .717, .735, .750, .779, &
    346      .623, .623, .623, .623, .627, .637, &
    347      .649, .668, .688, .709, .734, .762, &
    348      .592, .592, .592, .592, .595, .605, &
    349      .618, .639, .660, .682, .711, .743, &
    350      .561, .561, .561, .561, .565, .575, &
    351      .588, .609, .632, .656, .688, .724, &
    352      ! sulfate
    353      .671, .684, .697, .704, .714, .723, &
    354      .734, .746, .762, .771, .781, .789, &
    355      .653, .666, .678, .687, .697, .707, &
    356      .719, .732, .751, .762, .775, .789, &
    357      .622, .635, .648, .657, .667, .678, &
    358      .691, .705, .728, .741, .758, .777, &
    359      .591, .604, .617, .627, .638, .650, &
    360      .664, .679, .704, .719, .739, .761, &
    361      .560, .574, .587, .597, .609, .621, &
    362      .637, .653, .680, .697, .719, .745/
    363   !
     370                                ! dust insoluble
     371       0.759, 0.770, 0.775, 0.775, 0.772, &
     372                                !!jb bc insoluble
     373       11.536,10.033, 8.422, 7.234, 6.270, &
     374                                ! pom insoluble
     375       5.042, 3.101, 1.890, 1.294, 0.934/
     376  !
     377 DATA cg_aers_5wv/ &
     378                                ! bc soluble
     379       .651, .651, .651, .651, .651, .651, &
     380       .651, .651, .738, .764, .785, .800, &
     381       .597, .597, .597, .597, .597, .597, &
     382       .597, .597, .695, .725, .751, .770, &
     383       .543, .543, .543, .543, .543, .543, &
     384       .543, .543, .650, .684, .714, .736, &
     385       .504, .504, .504, .504, .504, .504, &
     386       .504, .504, .614, .651, .683, .708, &
     387       .469, .469, .469, .469, .469, .469, &
     388       .469, .469, .582, .620, .655, .681, &
     389                                ! pom soluble
     390       .679, .679, .679, .679, .683, .691, &
     391       .703, .720, .736, .751, .766, .784, &
     392       .656, .656, .656, .656, .659, .669, &
     393       .681, .699, .717, .735, .750, .779, &
     394       .623, .623, .623, .623, .627, .637, &
     395       .649, .668, .688, .709, .734, .762, &
     396       .592, .592, .592, .592, .595, .605, &
     397       .618, .639, .660, .682, .711, .743, &
     398       .561, .561, .561, .561, .565, .575, &
     399       .588, .609, .632, .656, .688, .724, &
     400                                ! Accumulation sulfate
     401       .671, .684, .697, .704, .714, .723, &
     402       .734, .746, .762, .771, .781, .789, &
     403       .653, .666, .678, .687, .697, .707, &
     404       .719, .732, .751, .762, .775, .789, &
     405       .622, .635, .648, .657, .667, .678, &
     406       .691, .705, .728, .741, .758, .777, &
     407       .591, .604, .617, .627, .638, .650, &
     408       .664, .679, .704, .719, .739, .761, &
     409       .560, .574, .587, .597, .609, .621, &
     410       .637, .653, .680, .697, .719, .745, &
     411                                ! Coarse sulfate
     412       .671, .684, .697, .704, .714, .723, &
     413       .734, .746, .762, .771, .781, .789, &
     414       .653, .666, .678, .687, .697, .707, &
     415       .719, .732, .751, .762, .775, .789, &
     416       .622, .635, .648, .657, .667, .678, &
     417       .691, .705, .728, .741, .758, .777, &
     418       .591, .604, .617, .627, .638, .650, &
     419       .664, .679, .704, .719, .739, .761, &
     420       .560, .574, .587, .597, .609, .621, &
     421       .637, .653, .680, .697, .719, .745, &
     422                                ! For super coarse seasalt (computed below for 550nm!)
     423      0.730,0.753,0.760,0.766,0.772,0.793, &
     424      0.797,0.802,0.809,0.813,0.820,0.830, &
     425       0.000,0.000,0.000,0.000,0.000,0.000, &
     426       0.000,0.000,0.000,0.000,0.000,0.000, &
     427     0.721,0.744,0.750,0.756,0.762,0.784, &
     428     0.787,0.793,0.800,0.804,0.811,0.822, &
     429     0.717,0.741,0.747,0.753,0.759,0.780, &
     430     0.784,0.789,0.795,0.800,0.806,0.817, &
     431     0.715,0.739,0.745,0.751,0.757,0.777, & 
     432     0.781,0.786,0.793,0.797,0.803,0.814, &
     433                                ! For coarse-soluble seasalt (computed below for 550nm!)
     434     0.730,0.753,0.760,0.766,0.772,0.793, &
     435     0.797,0.802,0.809,0.813,0.820,0.830, &
     436       0.000,0.000,0.000,0.000,0.000,0.000, &
     437       0.000,0.000,0.000,0.000,0.000,0.000, &
     438     0.721,0.744,0.750,0.756,0.762,0.784, &
     439     0.787,0.793,0.800,0.804,0.811,0.822, &
     440     0.717,0.741,0.747,0.753,0.759,0.780, &
     441     0.784,0.789,0.795,0.800,0.806,0.817, &
     442     0.715,0.739,0.745,0.751,0.757,0.777, & 
     443     0.781,0.786,0.793,0.797,0.803,0.814, &
     444                                ! accumulation-seasalt soluble (computed below for 550nm!)
     445     0.698,0.722,0.729,0.736,0.743,0.765, &
     446     0.768,0.773,0.777,0.779,0.781,0.779, &
     447       0.000,0.000,0.000,0.000,0.000,0.000, &
     448       0.000,0.000,0.000,0.000,0.000,0.000, &
     449     0.658,0.691,0.701,0.710,0.720,0.756, &
     450     0.763,0.771,0.782,0.788,0.795,0.801, &
     451     0.632,0.668,0.679,0.690,0.701,0.743, &
     452     0.750,0.762,0.775,0.783,0.792,0.804, &
     453     0.605,0.644,0.656,0.669,0.681,0.729, &
     454     0.737,0.750,0.765,0.775,0.787,0.803/
    364455
    365456  DATA cg_aeri_5wv/&
    … …  
    372463  !
    373464  DATA piz_aers_5wv/&
    374      ! bc soluble
    375      .445, .445, .445, .445, .445, .445, &
    376      .445, .445, .470, .487, .508, .531, &
    377      .442, .442, .442, .442, .442, .442, &
    378      .442, .442, .462, .481, .506, .533, &
    379      .427, .427, .427, .427, .427, .427, &
    380      .427, .427, .449, .470, .497, .526, &
    381      .413, .413, .413, .413, .413, .413, &
    382      .413, .413, .437, .458, .486, .516, &
    383      .399, .399, .399, .399, .399, .399, &
    384      .399, .399, .423, .445, .473, .506, &
    385      ! pom soluble
    386      .975, .975, .975, .975, .975, .977, &
    387      .979, .982, .984, .987, .990, .994, &
    388      .972, .972, .972, .972, .973, .974, &
    389      .977, .980, .983, .986, .989, .993, &
    390      .963, .963, .963, .963, .964, .966, &
    391      .969, .974, .977, .982, .986, .991, &
    392      .955, .955, .955, .955, .955, .958, &
    393      .962, .967, .972, .977, .983, .989, &
    394      .944, .944, .944, .944, .944, .948, &
    395      .952, .959, .962, .972, .979, .987, &
    396      ! sulfate
    397      1.000,1.000,1.000,1.000,1.000,1.000, &
    398      1.000,1.000,1.000,1.000,1.000,1.000, &
    399      1.000,1.000,1.000,1.000,1.000,1.000, &
    400      1.000,1.000,1.000,1.000,1.000,1.000, &
    401      1.000,1.000,1.000,1.000,1.000,1.000, &
    402      1.000,1.000,1.000,1.000,1.000,1.000, &
    403      1.000,1.000,1.000,1.000,1.000,1.000, &
    404      1.000,1.000,1.000,1.000,1.000,1.000, &
    405      1.000,1.000,1.000,1.000,1.000,1.000, &
    406      1.000,1.000,1.000,1.000,1.000,1.000/
     465                                ! bc soluble
     466       .445, .445, .445, .445, .445, .445, &
     467       .445, .445, .470, .487, .508, .531, &
     468       .442, .442, .442, .442, .442, .442, &
     469       .442, .442, .462, .481, .506, .533, &
     470       .427, .427, .427, .427, .427, .427, &
     471       .427, .427, .449, .470, .497, .526, &
     472       .413, .413, .413, .413, .413, .413, &
     473       .413, .413, .437, .458, .486, .516, &
     474       .399, .399, .399, .399, .399, .399, &
     475       .399, .399, .423, .445, .473, .506, &
     476                                ! pom soluble
     477       .975, .975, .975, .975, .975, .977, &
     478       .979, .982, .984, .987, .990, .994, &
     479       .972, .972, .972, .972, .973, .974, &
     480       .977, .980, .983, .986, .989, .993, &
     481       .963, .963, .963, .963, .964, .966, &
     482       .969, .974, .977, .982, .986, .991, &
     483       .955, .955, .955, .955, .955, .958, &
     484       .962, .967, .972, .977, .983, .989, &
     485       .944, .944, .944, .944, .944, .948, &
     486       .952, .959, .962, .972, .979, .987, &
     487                                ! sulfate soluble accumulation
     488       1.000,1.000,1.000,1.000,1.000,1.000, &
     489       1.000,1.000,1.000,1.000,1.000,1.000, &
     490       1.000,1.000,1.000,1.000,1.000,1.000, &
     491       1.000,1.000,1.000,1.000,1.000,1.000, &
     492       1.000,1.000,1.000,1.000,1.000,1.000, &
     493       1.000,1.000,1.000,1.000,1.000,1.000, &
     494       1.000,1.000,1.000,1.000,1.000,1.000, &
     495       1.000,1.000,1.000,1.000,1.000,1.000, &
     496       1.000,1.000,1.000,1.000,1.000,1.000, &
     497       1.000,1.000,1.000,1.000,1.000,1.000, &
     498                                ! sulfate soluble coarse
     499       1.000,1.000,1.000,1.000,1.000,1.000, &
     500       1.000,1.000,1.000,1.000,1.000,1.000, &
     501       1.000,1.000,1.000,1.000,1.000,1.000, &
     502       1.000,1.000,1.000,1.000,1.000,1.000, &
     503       1.000,1.000,1.000,1.000,1.000,1.000, &
     504       1.000,1.000,1.000,1.000,1.000,1.000, &
     505       1.000,1.000,1.000,1.000,1.000,1.000, &
     506       1.000,1.000,1.000,1.000,1.000,1.000, &
     507       1.000,1.000,1.000,1.000,1.000,1.000, &
     508       1.000,1.000,1.000,1.000,1.000,1.000, &
     509                                ! seasalt super coarse (computed below for 550nm)
     510       1.000,1.000,1.000,1.000,1.000,1.000, &
     511       1.000,1.000,1.000,1.000,1.000,1.000, &
     512       1.000,1.000,1.000,1.000,1.000,1.000, &
     513       1.000,1.000,1.000,1.000,1.000,1.000, &
     514       1.000,1.000,1.000,1.000,1.000,1.000, &
     515       1.000,1.000,1.000,1.000,1.000,1.000, &
     516       1.000,1.000,1.000,1.000,1.000,1.000, &
     517       1.000,1.000,1.000,1.000,1.000,1.000, &
     518       1.000,1.000,1.000,1.000,1.000,1.000, &
     519       1.000,1.000,1.000,1.000,1.000,1.000, &
     520                                ! seasalt coarse (computed below for 550nm)
     521       1.000,1.000,1.000,1.000,1.000,1.000, &
     522       1.000,1.000,1.000,1.000,1.000,1.000, &
     523       1.000,1.000,1.000,1.000,1.000,1.000, &
     524       1.000,1.000,1.000,1.000,1.000,1.000, &
     525       1.000,1.000,1.000,1.000,1.000,1.000, &
     526       1.000,1.000,1.000,1.000,1.000,1.000, &
     527       1.000,1.000,1.000,1.000,1.000,1.000, &
     528       1.000,1.000,1.000,1.000,1.000,1.000, &
     529       1.000,1.000,1.000,1.000,1.000,1.000, &
     530       1.000,1.000,1.000,1.000,1.000,1.000, &
     531                                ! seasalt soluble accumulation (computed below for 550nm)
     532       1.000,1.000,1.000,1.000,1.000,1.000, &
     533       1.000,1.000,1.000,1.000,1.000,1.000, &
     534       1.000,1.000,1.000,1.000,1.000,1.000, &
     535       1.000,1.000,1.000,1.000,1.000,1.000, &
     536       1.000,1.000,1.000,1.000,1.000,1.000, &
     537       1.000,1.000,1.000,1.000,1.000,1.000, &
     538       1.000,1.000,1.000,1.000,1.000,1.000, &
     539       1.000,1.000,1.000,1.000,1.000,1.000, &
     540       1.000,1.000,1.000,1.000,1.000,1.000, &
     541       1.000,1.000,1.000,1.000,1.000,1.000/
    407542  !
    408543  DATA piz_aeri_5wv/&
    … …  
    417552  IF (firstcall) THEN
    418553     firstcall=.FALSE.
     554! Allocation
     555    IF (.NOT. ALLOCATED(A1_ASSSM)) THEN
     556        ALLOCATE(A1_ASSSM(klev),A2_ASSSM(klev), A3_ASSSM(klev),&
     557          B1_ASSSM(klev), B2_ASSSM(klev), C1_ASSSM(klev), C2_ASSSM(klev),&
     558          A1_CSSSM(klev), A2_CSSSM(klev), A3_CSSSM(klev),&
     559          B1_CSSSM(klev), B2_CSSSM(klev), C1_CSSSM(klev), C2_CSSSM(klev),&
     560          A1_SSSSM(klev), A2_SSSSM(klev), A3_SSSSM(klev),&
     561          B1_SSSSM(klev), B2_SSSSM(klev), C1_SSSSM(klev), C2_SSSSM(klev), stat=ierr)
     562        IF (ierr /= 0) CALL abort_gcm('aeropt_5mw', 'pb in allocation 1',1)
     563     END IF
     564
    419565!Accumulation mode
    420566     CALL pres2lev(A1_ASSSM_19, A1_ASSSM, nb_level, klev, presnivs_19, presnivs, 1, 1, .FALSE.)
    … …  
    456602  DO k=1, klev
    457603    DO i=1, klon
    458       IF (t_seri(i,k).EQ.0) stop 'stop aeropt_5wv T '
    459       IF (pplay(i,k).EQ.0) stop  'stop aeropt_5wv p '
    460 
     604!      IF (t_seri(i,k).EQ.0) stop 'stop aeropt_5wv T '
     605!      IF (pplay(i,k).EQ.0) stop  'stop aeropt_5wv p '
    461606      zrho=pplay(i,k)/t_seri(i,k)/RD                  ! kg/m3
    462607      mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9
    463 
    464608    ENDDO
    465609  ENDDO
    … …  
    519663  zdp1=pdel(:,:)/(gravit*delt)
    520664 
    521   DO m=1,nb_aer   ! tau is only computed for each mass
    522    
     665  DO m=1,nb_aer   ! tau is only computed for each mass   
    523666    fac=1.0
    524667    IF (aerosol_name(m).EQ.id_ASBCM) THEN
    … …  
    530673        spsol=2
    531674        spss=0
    532     ELSEIF ((aerosol_name(m).EQ.id_ASSO4M) .OR.  (aerosol_name(m).EQ.id_CSSO4M)) THEN
     675    ELSEIF (aerosol_name(m).EQ.id_ASSO4M) THEN
    533676        soluble=.TRUE.
    534677        spsol=3
     678        spss=0
     679        fac=1.375    ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD
     680    ELSEIF (aerosol_name(m).EQ.id_CSSO4M) THEN
     681        soluble=.TRUE.
     682        spsol=4
    535683        spss=0
    536684        fac=1.375    ! (NH4)2-SO4/SO4 132/96 mass conversion factor for OD
    … …  
    581729          IF (soluble) THEN
    582730
    583 
    584731              IF((la.EQ.2).AND.(spss.NE.0)) THEN !la=2 corresponds to 550 nm
    585732                  H=rh/100
    … …  
    620767                 mass_temp(i,k,spsol)*1000.*zdp1(i,k)*tau_ae5wv_int(i,k,la)*delt*fac
    621768
    622           ELSE
     769               tausum(i,la,spsol)=tausum(i,la,spsol)+tau3d(i,k)
     770           tau(i,k,la,spsol)=tau3d(i,k)
     771
     772          ELSE                                           ! Case insoluble aerosol
    623773              tau_ae5wv_int(i,k,la) = alpha_aeri_5wv(la,spinsol)
    624774              piz_ae5wv_int(i,k,la) = piz_aeri_5wv(la,spinsol)
    … …  
    626776
    627777              tau3d(i,k) = &
    628                  mass_temp(i,k,7+spinsol)*1000.*zdp1(i,k)*tau_ae5wv_int(i,k,la)*delt*fac
     778                 mass_temp(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)* &
     779                      tau_ae5wv_int(i,k,la)*delt*fac
     780
     781               tausum(i,la,naero_soluble+spinsol)=tausum(i,la,spsol)+tau3d(i,k)
     782           tau(i,k,la,naero_soluble+spinsol)=tau3d(i,k)
     783
    629784          ENDIF
    630          
    631          
     785
    632786        ENDDO     ! Boucle sur les points géographiques (grille horizontale)
    633787      ENDDO     ! Boucle sur les niveaux verticaux
    634 
    635       IF (soluble) THEN
    636 
    637           tau(:,:,la,spsol)=tau3d(:,:)
    638      
    639           DO k=1, KLEV
    640             DO i=1,KLON
    641               tausum(i,la,spsol)=tausum(i,la,spsol)+tau3d(i,k)
    642             ENDDO
    643           ENDDO
    644       ELSE
    645           tau(:,:,la,spsol)=tau3d(:,:)
    646      
    647           DO k=1, KLEV
    648             DO i=1,KLON
    649               tausum(i,la,5+spinsol)=tausum(i,la,5+spinsol)+tau3d(i,k)
    650             ENDDO
    651           ENDDO
    652       ENDIF
    653 
    654 
    655 
    656788    ENDDO   ! boucle sur les longueurs d'onde
    657789  ENDDO     ! Boucle  sur les masses de traceurs
Note: See TracChangeset for help on using the changeset viewer.