Changeset 1254 for LMDZ4/branches


Ignore:
Timestamp:
Oct 30, 2009, 2:41:02 PM (15 years ago)
Author:
yann meurdesoif
Message:

Correction de la regression de la version 1246 sur les datas aerosols.

YM

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r1250 r1254  
    301301!
    302302!
    303 !
    304 ! From here on we look at the optical parameters at 5 wavelengths:
    305 ! 443nm, 550, 670, 765 and 865 nm
    306 !                                   le 12 AVRIL 2006
    307303!
    308   DATA alpha_aers_5wv/ &
    309      ! bc soluble
    310      7.930,7.930,7.930,7.930,7.930,7.930,     &
    311      7.930,7.930,10.893,12.618,14.550,16.613, &
    312      7.658,7.658,7.658,7.658,7.658,7.658,     &
    313      7.658,7.658,10.351,11.879,13.642,15.510, &
    314      7.195,7.195,7.195,7.195,7.195,7.195,     &
    315      7.195,7.195,9.551,10.847,12.381,13.994,  &
    316      6.736,6.736,6.736,6.736,6.736,6.736,     &
    317      6.736,6.736,8.818,9.938,11.283,12.687,   &
    318      6.277,6.277,6.277,6.277,6.277,6.277,     &
    319      6.277,6.277,8.123,9.094,10.275,11.501,   &
    320      ! pom soluble
    321      6.676,6.676,6.676,6.676,6.710,6.934,   &
    322      7.141,7.569,8.034,8.529,9.456,10.511,  &
    323      5.109,5.109,5.109,5.109,5.189,5.535,   &
    324      5.960,6.852,8.008,9.712,12.897,19.676, &
    325      3.718,3.718,3.718,3.718,3.779,4.042,   &
    326      4.364,5.052,5.956,7.314,9.896,15.688,  &
    327      2.849,2.849,2.849,2.849,2.897,3.107,   &
    328      3.365,3.916,4.649,5.760,7.900,12.863,  &
    329      2.229,2.229,2.229,2.229,2.268,2.437,   &
    330      2.645,3.095,3.692,4.608,6.391,10.633,  &
    331      ! Sulfate
    332      5.751,6.215,6.690,7.024,7.599,8.195,      &
    333      9.156,10.355,12.660,14.823,18.908,24.508, &
    334      4.320,4.675,5.052,5.375,5.787,6.274,      &
    335      7.066,8.083,10.088,12.003,15.697,21.133,  &
    336      3.079,3.351,3.639,3.886,4.205,4.584,      &
    337      5.206,6.019,7.648,9.234,12.391,17.220,    &
    338      2.336,2.552,2.781,2.979,3.236,3.540,      &
    339      4.046,4.711,6.056,7.388,10.093,14.313,    &
    340      1.777,1.949,2.134,2.292,2.503,2.751,      &
    341      3.166,3.712,4.828,5.949,8.264,11.922/
     304! 
     305!
     306! From here on we look at the optical parameters at 5 wavelengths: 
     307! 443nm, 550, 670, 765 and 865 nm
     308!                                   le 12 AVRIL 2006
     309! 
     310 DATA alpha_aers_5wv/ &
     311                                ! bc soluble
     312       7.930,7.930,7.930,7.930,7.930,7.930,     &
     313       7.930,7.930,10.893,12.618,14.550,16.613, &
     314       7.658,7.658,7.658,7.658,7.658,7.658,     &
     315       7.658,7.658,10.351,11.879,13.642,15.510, &
     316       7.195,7.195,7.195,7.195,7.195,7.195,     &
     317       7.195,7.195,9.551,10.847,12.381,13.994,  &
     318       6.736,6.736,6.736,6.736,6.736,6.736,     &
     319       6.736,6.736,8.818,9.938,11.283,12.687,   &
     320       6.277,6.277,6.277,6.277,6.277,6.277,     &
     321       6.277,6.277,8.123,9.094,10.275,11.501,   &
     322                                ! pom soluble
     323       6.676,6.676,6.676,6.676,6.710,6.934,   &
     324       7.141,7.569,8.034,8.529,9.456,10.511,  &
     325       5.109,5.109,5.109,5.109,5.189,5.535,   &
     326       5.960,6.852,8.008,9.712,12.897,19.676, &
     327       3.718,3.718,3.718,3.718,3.779,4.042,   &
     328       4.364,5.052,5.956,7.314,9.896,15.688,  &
     329       2.849,2.849,2.849,2.849,2.897,3.107,   &
     330       3.365,3.916,4.649,5.760,7.900,12.863,  &
     331       2.229,2.229,2.229,2.229,2.268,2.437,   &
     332       2.645,3.095,3.692,4.608,6.391,10.633,  &
     333                                ! Sulfate (Accumulation)
     334       5.751,6.215,6.690,7.024,7.599,8.195,      &
     335       9.156,10.355,12.660,14.823,18.908,24.508, &
     336       4.320,4.675,5.052,5.375,5.787,6.274,      &
     337       7.066,8.083,10.088,12.003,15.697,21.133,  &
     338       3.079,3.351,3.639,3.886,4.205,4.584,      &
     339       5.206,6.019,7.648,9.234,12.391,17.220,    &
     340       2.336,2.552,2.781,2.979,3.236,3.540,      &
     341       4.046,4.711,6.056,7.388,10.093,14.313,    &
     342       1.777,1.949,2.134,2.292,2.503,2.751,      &
     343       3.166,3.712,4.828,5.949,8.264,11.922,     &
     344                                ! Sulfate (Coarse)
     345       5.751,6.215,6.690,7.024,7.599,8.195,      &
     346       9.156,10.355,12.660,14.823,18.908,24.508, &
     347       4.320,4.675,5.052,5.375,5.787,6.274,      &
     348       7.066,8.083,10.088,12.003,15.697,21.133,  &
     349       3.079,3.351,3.639,3.886,4.205,4.584,      &
     350       5.206,6.019,7.648,9.234,12.391,17.220,    &
     351       2.336,2.552,2.781,2.979,3.236,3.540,      &
     352       4.046,4.711,6.056,7.388,10.093,14.313,    &
     353       1.777,1.949,2.134,2.292,2.503,2.751,      &
     354       3.166,3.712,4.828,5.949,8.264,11.922,     &
     355                                ! Seasalt soluble super_coarse (computed below for 550nm)
     356       0.50,0.90,1.05,1.21,1.40,2.41, & 
     357       2.66,3.11,3.88,4.52,5.69,8.84, & 
     358       0.000,0.000,0.000,0.000,0.000,0.000, & 
     359       0.000,0.000,0.000,0.000,0.000,0.000, & 
     360     0.52,0.93,1.08,1.24,1.43,2.47, & 
     361     2.73,3.20,3.99,4.64,5.84,9.04, & 
     362     0.52,0.93,1.09,1.25,1.44,2.50, & 
     363     2.76,3.23,4.03,4.68,5.89,9.14, & 
     364     0.52,0.94,1.09,1.26,1.45,2.51, & 
     365     2.78,3.25,4.06,4.72,5.94,9.22, & 
     366                                ! seasalt soluble coarse (computed below for 550nm)
     367       0.50,0.90,1.05,1.21,1.40,2.41, & 
     368       2.66,3.11,3.88,4.52,5.69,8.84, & 
     369       0.000,0.000,0.000,0.000,0.000,0.000, & 
     370       0.000,0.000,0.000,0.000,0.000,0.000, & 
     371     0.52,0.93,1.08,1.24,1.43,2.47, & 
     372     2.73,3.20,3.99,4.64,5.84,9.04, & 
     373     0.52,0.93,1.09,1.25,1.44,2.50, & 
     374     2.76,3.23,4.03,4.68,5.89,9.14, & 
     375     0.52,0.94,1.09,1.26,1.45,2.51, & 
     376     2.78,3.25,4.06,4.72,5.94,9.22, & 
     377                                ! seasalt soluble accumulation (computed below for 550nm)
     378     4.28, 7.17, 8.44, 9.85,11.60,22.44,  & 
     379     25.34,30.54,39.38,46.52,59.33,91.77, & 
     380       0.000,0.000,0.000,0.000,0.000,0.000, & 
     381       0.000,0.000,0.000,0.000,0.000,0.000, & 
     382     2.48, 4.22, 5.02, 5.94, 7.11,15.29,  & 
     383     17.70,22.31,30.73,38.06,52.15,90.59, & 
     384     1.90, 3.29, 3.94, 4.69, 5.65, 12.58, & 
     385     14.68,18.77,26.41,33.25,46.77,85.50, & 
     386     1.47, 2.59, 3.12, 3.74, 4.54, 10.42, & 
     387     12.24,15.82,22.66,28.91,41.54,79.33/
    342388
    343389  DATA alpha_aeri_5wv/ &
    344      ! dust insoluble
    345      0.759, 0.770, 0.775, 0.775, 0.772, &
    346      !!jb bc insoluble
    347      11.536,10.033, 8.422, 7.234, 6.270, &
    348      ! pom insoluble
    349      5.042, 3.101, 1.890, 1.294, 0.934/
    350 
    351   DATA cg_aers_5wv/ &
    352      ! bc soluble
    353      .651, .651, .651, .651, .651, .651, &
    354      .651, .651, .738, .764, .785, .800, &
    355      .597, .597, .597, .597, .597, .597, &
    356      .597, .597, .695, .725, .751, .770, &
    357      .543, .543, .543, .543, .543, .543, &
    358      .543, .543, .650, .684, .714, .736, &
    359      .504, .504, .504, .504, .504, .504, &
    360      .504, .504, .614, .651, .683, .708, &
    361      .469, .469, .469, .469, .469, .469, &
    362      .469, .469, .582, .620, .655, .681, &
    363      ! pom soluble
    364      .679, .679, .679, .679, .683, .691, &
    365      .703, .720, .736, .751, .766, .784, &
    366      .656, .656, .656, .656, .659, .669, &
    367      .681, .699, .717, .735, .750, .779, &
    368      .623, .623, .623, .623, .627, .637, &
    369      .649, .668, .688, .709, .734, .762, &
    370      .592, .592, .592, .592, .595, .605, &
    371      .618, .639, .660, .682, .711, .743, &
    372      .561, .561, .561, .561, .565, .575, &
    373      .588, .609, .632, .656, .688, .724, &
    374      ! sulfate
    375      .671, .684, .697, .704, .714, .723, &
    376      .734, .746, .762, .771, .781, .789, &
    377      .653, .666, .678, .687, .697, .707, &
    378      .719, .732, .751, .762, .775, .789, &
    379      .622, .635, .648, .657, .667, .678, &
    380      .691, .705, .728, .741, .758, .777, &
    381      .591, .604, .617, .627, .638, .650, &
    382      .664, .679, .704, .719, .739, .761, &
    383      .560, .574, .587, .597, .609, .621, &
    384      .637, .653, .680, .697, .719, .745/
    385   !
     390                                 ! dust insoluble
     391        0.759, 0.770, 0.775, 0.775, 0.772, &
     392                                 !!jb bc insoluble
     393        11.536,10.033, 8.422, 7.234, 6.270, &
     394                                 ! pom insoluble
     395        5.042, 3.101, 1.890, 1.294, 0.934/
     396   !
     397  DATA cg_aers_5wv/ & 
     398                                 ! bc soluble
     399      .651, .651, .651, .651, .651, .651, &
     400      .651, .651, .738, .764, .785, .800, &
     401      .597, .597, .597, .597, .597, .597, &
     402      .597, .597, .695, .725, .751, .770, &
     403      .543, .543, .543, .543, .543, .543, &
     404      .543, .543, .650, .684, .714, .736, & 
     405      .504, .504, .504, .504, .504, .504, &
     406      .504, .504, .614, .651, .683, .708, & 
     407      .469, .469, .469, .469, .469, .469, &
     408      .469, .469, .582, .620, .655, .681, &
     409                                 ! pom soluble
     410      .679, .679, .679, .679, .683, .691, &
     411      .703, .720, .736, .751, .766, .784, &
     412      .656, .656, .656, .656, .659, .669, &
     413      .681, .699, .717, .735, .750, .779, & 
     414      .623, .623, .623, .623, .627, .637, &
     415      .649, .668, .688, .709, .734, .762, &
     416      .592, .592, .592, .592, .595, .605, &
     417      .618, .639, .660, .682, .711, .743, &
     418      .561, .561, .561, .561, .565, .575, &
     419      .588, .609, .632, .656, .688, .724, &
     420                                 ! Accumulation sulfate
     421      .671, .684, .697, .704, .714, .723, &
     422      .734, .746, .762, .771, .781, .789, &
     423      .653, .666, .678, .687, .697, .707, &
     424      .719, .732, .751, .762, .775, .789, &
     425      .622, .635, .648, .657, .667, .678, &
     426      .691, .705, .728, .741, .758, .777, &
     427      .591, .604, .617, .627, .638, .650, &
     428      .664, .679, .704, .719, .739, .761, &
     429      .560, .574, .587, .597, .609, .621, & 
     430      .637, .653, .680, .697, .719, .745, &
     431                                 ! Coarse sulfate
     432      .671, .684, .697, .704, .714, .723, &
     433      .734, .746, .762, .771, .781, .789, &
     434      .653, .666, .678, .687, .697, .707, &
     435      .719, .732, .751, .762, .775, .789, &
     436      .622, .635, .648, .657, .667, .678, &
     437      .691, .705, .728, .741, .758, .777, &
     438      .591, .604, .617, .627, .638, .650, &
     439      .664, .679, .704, .719, .739, .761, &
     440      .560, .574, .587, .597, .609, .621, & 
     441      .637, .653, .680, .697, .719, .745, &
     442                                 ! For super coarse seasalt (computed below for 550nm!)
     443      0.730,0.753,0.760,0.766,0.772,0.793, & 
     444      0.797,0.802,0.809,0.813,0.820,0.830, & 
     445      0.000,0.000,0.000,0.000,0.000,0.000, & 
     446      0.000,0.000,0.000,0.000,0.000,0.000, & 
     447      0.721,0.744,0.750,0.756,0.762,0.784, & 
     448      0.787,0.793,0.800,0.804,0.811,0.822, & 
     449      0.717,0.741,0.747,0.753,0.759,0.780, & 
     450      0.784,0.789,0.795,0.800,0.806,0.817, & 
     451      0.715,0.739,0.745,0.751,0.757,0.777, &   
     452      0.781,0.786,0.793,0.797,0.803,0.814, & 
     453                                 ! For coarse-soluble seasalt (computed below for 550nm!)
     454      0.730,0.753,0.760,0.766,0.772,0.793, & 
     455      0.797,0.802,0.809,0.813,0.820,0.830, & 
     456      0.000,0.000,0.000,0.000,0.000,0.000, & 
     457      0.000,0.000,0.000,0.000,0.000,0.000, & 
     458      0.721,0.744,0.750,0.756,0.762,0.784, & 
     459      0.787,0.793,0.800,0.804,0.811,0.822, & 
     460      0.717,0.741,0.747,0.753,0.759,0.780, & 
     461      0.784,0.789,0.795,0.800,0.806,0.817, & 
     462      0.715,0.739,0.745,0.751,0.757,0.777, &   
     463      0.781,0.786,0.793,0.797,0.803,0.814, & 
     464                                 ! accumulation-seasalt soluble (computed below for 550nm!) 
     465      0.698,0.722,0.729,0.736,0.743,0.765, & 
     466      0.768,0.773,0.777,0.779,0.781,0.779, & 
     467      0.000,0.000,0.000,0.000,0.000,0.000, & 
     468      0.000,0.000,0.000,0.000,0.000,0.000, & 
     469      0.658,0.691,0.701,0.710,0.720,0.756, & 
     470      0.763,0.771,0.782,0.788,0.795,0.801, & 
     471      0.632,0.668,0.679,0.690,0.701,0.743, & 
     472      0.750,0.762,0.775,0.783,0.792,0.804, & 
     473      0.605,0.644,0.656,0.669,0.681,0.729, & 
     474      0.737,0.750,0.765,0.775,0.787,0.803/
     475 !
    386476
    387477  DATA cg_aeri_5wv/&
    … …  
    394484  !
    395485  DATA piz_aers_5wv/&
    396      ! bc soluble
    397      .445, .445, .445, .445, .445, .445, &
    398      .445, .445, .470, .487, .508, .531, &
    399      .442, .442, .442, .442, .442, .442, &
    400      .442, .442, .462, .481, .506, .533, &
    401      .427, .427, .427, .427, .427, .427, &
    402      .427, .427, .449, .470, .497, .526, &
    403      .413, .413, .413, .413, .413, .413, &
    404      .413, .413, .437, .458, .486, .516, &
    405      .399, .399, .399, .399, .399, .399, &
    406      .399, .399, .423, .445, .473, .506, &
    407      ! pom soluble
    408      .975, .975, .975, .975, .975, .977, &
    409      .979, .982, .984, .987, .990, .994, &
    410      .972, .972, .972, .972, .973, .974, &
    411      .977, .980, .983, .986, .989, .993, &
    412      .963, .963, .963, .963, .964, .966, &
    413      .969, .974, .977, .982, .986, .991, &
    414      .955, .955, .955, .955, .955, .958, &
    415      .962, .967, .972, .977, .983, .989, &
    416      .944, .944, .944, .944, .944, .948, &
    417      .952, .959, .962, .972, .979, .987, &
    418      ! sulfate
    419      1.000,1.000,1.000,1.000,1.000,1.000, &
    420      1.000,1.000,1.000,1.000,1.000,1.000, &
    421      1.000,1.000,1.000,1.000,1.000,1.000, &
    422      1.000,1.000,1.000,1.000,1.000,1.000, &
    423      1.000,1.000,1.000,1.000,1.000,1.000, &
    424      1.000,1.000,1.000,1.000,1.000,1.000, &
    425      1.000,1.000,1.000,1.000,1.000,1.000, &
    426      1.000,1.000,1.000,1.000,1.000,1.000, &
    427      1.000,1.000,1.000,1.000,1.000,1.000, &
    428      1.000,1.000,1.000,1.000,1.000,1.000/
    429   !
     486                           ! bc soluble
     487  .445, .445, .445, .445, .445, .445, &
     488  .445, .445, .470, .487, .508, .531, &
     489  .442, .442, .442, .442, .442, .442, &
     490  .442, .442, .462, .481, .506, .533, &
     491  .427, .427, .427, .427, .427, .427, &
     492  .427, .427, .449, .470, .497, .526, &
     493  .413, .413, .413, .413, .413, .413, &
     494  .413, .413, .437, .458, .486, .516, &
     495  .399, .399, .399, .399, .399, .399, &
     496  .399, .399, .423, .445, .473, .506, &
     497                           ! pom soluble
     498  .975, .975, .975, .975, .975, .977, &
     499  .979, .982, .984, .987, .990, .994, &
     500  .972, .972, .972, .972, .973, .974, &
     501  .977, .980, .983, .986, .989, .993, &
     502  .963, .963, .963, .963, .964, .966, &
     503  .969, .974, .977, .982, .986, .991, &
     504  .955, .955, .955, .955, .955, .958, &
     505  .962, .967, .972, .977, .983, .989, &
     506  .944, .944, .944, .944, .944, .948, &
     507  .952, .959, .962, .972, .979, .987, &
     508                           ! sulfate soluble accumulation
     509  1.000,1.000,1.000,1.000,1.000,1.000, &
     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                           ! sulfate soluble coarse
     520  1.000,1.000,1.000,1.000,1.000,1.000, & 
     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                           ! seasalt super coarse (computed below for 550nm)
     531  1.000,1.000,1.000,1.000,1.000,1.000, & 
     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                           ! seasalt coarse (computed below for 550nm)
     542  1.000,1.000,1.000,1.000,1.000,1.000, & 
     543  1.000,1.000,1.000,1.000,1.000,1.000, & 
     544  1.000,1.000,1.000,1.000,1.000,1.000, & 
     545  1.000,1.000,1.000,1.000,1.000,1.000, & 
     546  1.000,1.000,1.000,1.000,1.000,1.000, & 
     547  1.000,1.000,1.000,1.000,1.000,1.000, & 
     548  1.000,1.000,1.000,1.000,1.000,1.000, & 
     549  1.000,1.000,1.000,1.000,1.000,1.000, & 
     550  1.000,1.000,1.000,1.000,1.000,1.000, & 
     551  1.000,1.000,1.000,1.000,1.000,1.000, & 
     552                           ! seasalt soluble accumulation (computed below for 550nm)
     553  1.000,1.000,1.000,1.000,1.000,1.000, & 
     554  1.000,1.000,1.000,1.000,1.000,1.000, & 
     555  1.000,1.000,1.000,1.000,1.000,1.000, & 
     556  1.000,1.000,1.000,1.000,1.000,1.000, & 
     557  1.000,1.000,1.000,1.000,1.000,1.000, & 
     558  1.000,1.000,1.000,1.000,1.000,1.000, & 
     559  1.000,1.000,1.000,1.000,1.000,1.000, & 
     560  1.000,1.000,1.000,1.000,1.000,1.000, & 
     561  1.000,1.000,1.000,1.000,1.000,1.000, & 
     562  1.000,1.000,1.000,1.000,1.000,1.000/
     563
     564 !
    430565  DATA piz_aeri_5wv/&
    431566     ! dust insoluble
Note: See TracChangeset for help on using the changeset viewer.