Changeset 1266


Ignore:
Timestamp:
Nov 20, 2009, 5:27:40 PM (15 years ago)
Author:
Laurent Fairhead
Message:

Correction d'un bug dans aeropt_5wv.F90 : une variable était utilisée sans
être calculée
Ajout de tests sur l'écriture de variables dans les fichiers de sortie
Changement des arguments d'appels à INCA
ACo

Location:
LMDZ4/branches/LMDZ4-dev/libf/phylmd
Files:
4 edited

Legend:

Unmodified
Added
Removed

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

    r1254 r1266  
    144144 
    145145  REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble)   ! ext. coeff. Soluble comp. units *** m2/g
    146    !  1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
    147   REAL :: alpha_aeri_5wv(las,naero_insoluble)         ! ext. coeff. Insoluble comp. 1- Dust: 2- BC; 3- POM
     146  REAL :: alpha_aeri_5wv(las,naero_insoluble)         ! ext. coeff. Insoluble comp.
    148147  REAL :: cg_aers_5wv(nbre_RH,las,naero_soluble)      ! Asym. param. soluble comp.
    149    !  1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
    150   REAL :: cg_aeri_5wv(las,naero_insoluble)            ! Asym. param. insoluble comp. 1- Dust: 2- BC; 3- POM
     148  REAL :: cg_aeri_5wv(las,naero_insoluble)            ! Asym. param. insoluble comp.
    151149  REAL :: piz_aers_5wv(nbre_RH,las,naero_soluble)   
    152    !  1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
    153   REAL :: piz_aeri_5wv(las,naero_insoluble)           ! Insoluble comp. 1- Dust: 2- BC; 3- POM
     150  REAL :: piz_aeri_5wv(las,naero_insoluble)           ! Insoluble comp.
    154151
    155152  REAL, DIMENSION(klon,klev,naero_spc) :: mass_temp
    … …  
    308305!                                   le 12 AVRIL 2006
    309306! 
    310  DATA alpha_aers_5wv/ &
     307  DATA alpha_aers_5wv/ &
    311308                                ! 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, & 
     309     7.930,7.930,7.930,7.930,7.930,7.930,     &
     310     7.930,7.930,10.893,12.618,14.550,16.613, &
     311     7.658,7.658,7.658,7.658,7.658,7.658,     &
     312     7.658,7.658,10.351,11.879,13.642,15.510, &
     313     7.195,7.195,7.195,7.195,7.195,7.195,     &
     314     7.195,7.195,9.551,10.847,12.381,13.994,  &
     315     6.736,6.736,6.736,6.736,6.736,6.736,     &
     316     6.736,6.736,8.818,9.938,11.283,12.687,   &
     317     6.277,6.277,6.277,6.277,6.277,6.277,     &
     318     6.277,6.277,8.123,9.094,10.275,11.501,   &
     319     ! pom soluble
     320     6.676,6.676,6.676,6.676,6.710,6.934,   &
     321     7.141,7.569,8.034,8.529,9.456,10.511,  &
     322     5.109,5.109,5.109,5.109,5.189,5.535,   &
     323     5.960,6.852,8.008,9.712,12.897,19.676, &
     324     3.718,3.718,3.718,3.718,3.779,4.042,   &
     325     4.364,5.052,5.956,7.314,9.896,15.688,  &
     326     2.849,2.849,2.849,2.849,2.897,3.107,   &
     327     3.365,3.916,4.649,5.760,7.900,12.863,  &
     328     2.229,2.229,2.229,2.229,2.268,2.437,   &
     329     2.645,3.095,3.692,4.608,6.391,10.633,  &
     330     ! Sulfate (Accumulation)
     331     5.751,6.215,6.690,7.024,7.599,8.195,      &
     332     9.156,10.355,12.660,14.823,18.908,24.508, &
     333     4.320,4.675,5.052,5.375,5.787,6.274,      &
     334     7.066,8.083,10.088,12.003,15.697,21.133,  &
     335     3.079,3.351,3.639,3.886,4.205,4.584,      &
     336     5.206,6.019,7.648,9.234,12.391,17.220,    &
     337     2.336,2.552,2.781,2.979,3.236,3.540,      &
     338     4.046,4.711,6.056,7.388,10.093,14.313,    &
     339     1.777,1.949,2.134,2.292,2.503,2.751,      &
     340     3.166,3.712,4.828,5.949,8.264,11.922,     &
     341     ! Sulfate (Coarse)
     342     5.751,6.215,6.690,7.024,7.599,8.195,      &
     343     9.156,10.355,12.660,14.823,18.908,24.508, &
     344     4.320,4.675,5.052,5.375,5.787,6.274,      &
     345     7.066,8.083,10.088,12.003,15.697,21.133,  &
     346     3.079,3.351,3.639,3.886,4.205,4.584,      &
     347     5.206,6.019,7.648,9.234,12.391,17.220,    &
     348     2.336,2.552,2.781,2.979,3.236,3.540,      &
     349     4.046,4.711,6.056,7.388,10.093,14.313,    &
     350     1.777,1.949,2.134,2.292,2.503,2.751,      &
     351     3.166,3.712,4.828,5.949,8.264,11.922,     &
     352     ! Seasalt soluble super_coarse (computed below for 550nm)
     353     0.50,0.90,1.05,1.21,1.40,2.41, & 
     354     2.66,3.11,3.88,4.52,5.69,8.84, & 
     355     0.000,0.000,0.000,0.000,0.000,0.000, & 
     356     0.000,0.000,0.000,0.000,0.000,0.000, & 
    360357     0.52,0.93,1.08,1.24,1.43,2.47, & 
    361358     2.73,3.20,3.99,4.64,5.84,9.04, & 
    … …  
    364361     0.52,0.94,1.09,1.26,1.45,2.51, & 
    365362     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, & 
     363     ! seasalt soluble coarse (computed below for 550nm)
     364     0.50,0.90,1.05,1.21,1.40,2.41, & 
     365     2.66,3.11,3.88,4.52,5.69,8.84, & 
     366     0.000,0.000,0.000,0.000,0.000,0.000, & 
     367     0.000,0.000,0.000,0.000,0.000,0.000, & 
    371368     0.52,0.93,1.08,1.24,1.43,2.47, & 
    372369     2.73,3.20,3.99,4.64,5.84,9.04, & 
    … …  
    375372     0.52,0.94,1.09,1.26,1.45,2.51, & 
    376373     2.78,3.25,4.06,4.72,5.94,9.22, & 
    377                                 ! seasalt soluble accumulation (computed below for 550nm)
     374     ! seasalt soluble accumulation (computed below for 550nm)
    378375     4.28, 7.17, 8.44, 9.85,11.60,22.44,  & 
    379376     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, & 
     377     0.000,0.000,0.000,0.000,0.000,0.000, & 
     378     0.000,0.000,0.000,0.000,0.000,0.000, & 
    382379     2.48, 4.22, 5.02, 5.94, 7.11,15.29,  & 
    383380     17.70,22.31,30.73,38.06,52.15,90.59, & 
    … …  
    386383     1.47, 2.59, 3.12, 3.74, 4.54, 10.42, & 
    387384     12.24,15.82,22.66,28.91,41.54,79.33/
    388 
     385 
    389386  DATA alpha_aeri_5wv/ &
    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    !
     387     ! dust insoluble
     388     0.759, 0.770, 0.775, 0.775, 0.772, &
     389     !!jb bc insoluble
     390     11.536,10.033, 8.422, 7.234, 6.270, &
     391     ! pom insoluble
     392     5.042, 3.101, 1.890, 1.294, 0.934/
     393  !
    397394  DATA cg_aers_5wv/ & 
    398395                                 ! 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  !
    476 
     396     .651, .651, .651, .651, .651, .651, &
     397     .651, .651, .738, .764, .785, .800, &
     398     .597, .597, .597, .597, .597, .597, &
     399     .597, .597, .695, .725, .751, .770, &
     400     .543, .543, .543, .543, .543, .543, &
     401     .543, .543, .650, .684, .714, .736, & 
     402     .504, .504, .504, .504, .504, .504, &
     403     .504, .504, .614, .651, .683, .708, & 
     404     .469, .469, .469, .469, .469, .469, &
     405     .469, .469, .582, .620, .655, .681, &
     406     ! pom soluble
     407     .679, .679, .679, .679, .683, .691, &
     408     .703, .720, .736, .751, .766, .784, &
     409     .656, .656, .656, .656, .659, .669, &
     410     .681, .699, .717, .735, .750, .779, & 
     411     .623, .623, .623, .623, .627, .637, &
     412     .649, .668, .688, .709, .734, .762, &
     413     .592, .592, .592, .592, .595, .605, &
     414     .618, .639, .660, .682, .711, .743, &
     415     .561, .561, .561, .561, .565, .575, &
     416     .588, .609, .632, .656, .688, .724, &
     417     ! Accumulation sulfate
     418     .671, .684, .697, .704, .714, .723, &
     419     .734, .746, .762, .771, .781, .789, &
     420     .653, .666, .678, .687, .697, .707, &
     421     .719, .732, .751, .762, .775, .789, &
     422     .622, .635, .648, .657, .667, .678, &
     423     .691, .705, .728, .741, .758, .777, &
     424     .591, .604, .617, .627, .638, .650, &
     425     .664, .679, .704, .719, .739, .761, &
     426     .560, .574, .587, .597, .609, .621, & 
     427     .637, .653, .680, .697, .719, .745, &
     428     ! Coarse sulfate
     429     .671, .684, .697, .704, .714, .723, &
     430     .734, .746, .762, .771, .781, .789, &
     431     .653, .666, .678, .687, .697, .707, &
     432     .719, .732, .751, .762, .775, .789, &
     433     .622, .635, .648, .657, .667, .678, &
     434     .691, .705, .728, .741, .758, .777, &
     435     .591, .604, .617, .627, .638, .650, &
     436     .664, .679, .704, .719, .739, .761, &
     437     .560, .574, .587, .597, .609, .621, & 
     438     .637, .653, .680, .697, .719, .745, &
     439     ! For super coarse seasalt (computed below for 550nm!)
     440     0.730,0.753,0.760,0.766,0.772,0.793, & 
     441     0.797,0.802,0.809,0.813,0.820,0.830, & 
     442     0.000,0.000,0.000,0.000,0.000,0.000, & 
     443     0.000,0.000,0.000,0.000,0.000,0.000, & 
     444     0.721,0.744,0.750,0.756,0.762,0.784, & 
     445     0.787,0.793,0.800,0.804,0.811,0.822, & 
     446     0.717,0.741,0.747,0.753,0.759,0.780, & 
     447     0.784,0.789,0.795,0.800,0.806,0.817, & 
     448     0.715,0.739,0.745,0.751,0.757,0.777, &   
     449     0.781,0.786,0.793,0.797,0.803,0.814, & 
     450     ! For coarse-soluble seasalt (computed below for 550nm!)
     451     0.730,0.753,0.760,0.766,0.772,0.793, & 
     452     0.797,0.802,0.809,0.813,0.820,0.830, & 
     453     0.000,0.000,0.000,0.000,0.000,0.000, & 
     454     0.000,0.000,0.000,0.000,0.000,0.000, & 
     455     0.721,0.744,0.750,0.756,0.762,0.784, & 
     456     0.787,0.793,0.800,0.804,0.811,0.822, & 
     457     0.717,0.741,0.747,0.753,0.759,0.780, & 
     458     0.784,0.789,0.795,0.800,0.806,0.817, & 
     459     0.715,0.739,0.745,0.751,0.757,0.777, &   
     460     0.781,0.786,0.793,0.797,0.803,0.814, & 
     461     ! accumulation-seasalt soluble (computed below for 550nm!) 
     462     0.698,0.722,0.729,0.736,0.743,0.765, & 
     463     0.768,0.773,0.777,0.779,0.781,0.779, & 
     464     0.000,0.000,0.000,0.000,0.000,0.000, & 
     465     0.000,0.000,0.000,0.000,0.000,0.000, & 
     466     0.658,0.691,0.701,0.710,0.720,0.756, & 
     467     0.763,0.771,0.782,0.788,0.795,0.801, & 
     468     0.632,0.668,0.679,0.690,0.701,0.743, & 
     469     0.750,0.762,0.775,0.783,0.792,0.804, & 
     470     0.605,0.644,0.656,0.669,0.681,0.729, & 
     471     0.737,0.750,0.765,0.775,0.787,0.803/
     472  !
     473 
    477474  DATA cg_aeri_5wv/&
    478475     ! dust insoluble
    … …  
    485482  DATA piz_aers_5wv/&
    486483                           ! 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  !
     484     .445, .445, .445, .445, .445, .445, &
     485     .445, .445, .470, .487, .508, .531, &
     486     .442, .442, .442, .442, .442, .442, &
     487     .442, .442, .462, .481, .506, .533, &
     488     .427, .427, .427, .427, .427, .427, &
     489     .427, .427, .449, .470, .497, .526, &
     490     .413, .413, .413, .413, .413, .413, &
     491     .413, .413, .437, .458, .486, .516, &
     492     .399, .399, .399, .399, .399, .399, &
     493     .399, .399, .423, .445, .473, .506, &
     494     ! pom soluble
     495     .975, .975, .975, .975, .975, .977, &
     496     .979, .982, .984, .987, .990, .994, &
     497     .972, .972, .972, .972, .973, .974, &
     498     .977, .980, .983, .986, .989, .993, &
     499     .963, .963, .963, .963, .964, .966, &
     500     .969, .974, .977, .982, .986, .991, &
     501     .955, .955, .955, .955, .955, .958, &
     502     .962, .967, .972, .977, .983, .989, &
     503     .944, .944, .944, .944, .944, .948, &
     504     .952, .959, .962, .972, .979, .987, &
     505     ! sulfate soluble accumulation
     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     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     ! sulfate soluble coarse
     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     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     ! seasalt super coarse (computed below for 550nm)
     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     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     ! seasalt coarse (computed below for 550nm)
     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, & 
     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     ! seasalt soluble accumulation (computed below for 550nm)
     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     1.000,1.000,1.000,1.000,1.000,1.000, & 
     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 
     561  !
    565562  DATA piz_aeri_5wv/&
    566563     ! dust insoluble
    … …  
    808805!CDIR ON_ADB(alpha_aers_5wv)
    809806            DO i=1, KLON
     807              rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
     808              RH_num = INT( rh(i,k)/10. + 1.)
     809              IF (rh(i,k).GT.85.) RH_num=10
     810              IF (rh(i,k).GT.90.) RH_num=11
     811
    810812              tau_ae5wv_int = alpha_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* &
    811813                             (alpha_aers_5wv(RH_num+1,la,spsol) - &

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

    r1263 r1266  
    486486       ctetaSTD,dtime, ok_veget, &
    487487       type_ocean, iflag_pbl,ok_mensuel,ok_journe, &
    488        ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz)
     488       ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz, &
     489       new_aod, aerosol_couple)   
     490
    489491
    490492  USE iophy
    … …  
    507509  logical                               :: ok_mensuel, ok_journe, ok_hf, ok_instan
    508510  logical                               :: ok_LES,ok_ade,ok_aie
    509 
     511  logical                               :: new_aod, aerosol_couple
    510512  integer, intent(in)::  read_climoz ! read ozone climatology
    511513  !     Allowed values are 0, 1 and 2
    … …  
    833835 CALL histdef2d(iff,o_rugs_srf(nsrf)%flag,o_rugs_srf(nsrf)%name,"Latent heat flux "//clnsurf(nsrf),"W/m2")
    834836 CALL histdef2d(iff,o_ages_srf(nsrf)%flag,o_ages_srf(nsrf)%name,"Snow age", "day")
    835      END DO
    836 
    837      DO naero = 1, naero_spc
    838 CALL histdef2d(iff,o_tausumaero(naero)%flag,o_tausumaero(naero)%name,"Aerosol Optical depth at 550 nm "//name_aero(naero),"1")
    839      END DO
     837END DO
     838
     839IF (new_aod .AND. (.NOT. aerosol_couple)) THEN
     840  DO naero = 1, naero_spc
     841  CALL histdef2d(iff,o_tausumaero(naero)%flag,o_tausumaero(naero)%name,"Aerosol Optical depth at 550 nm "//name_aero(naero),"1")
     842  END DO
     843ENDIF
    840844
    841845 IF (ok_ade) THEN
    … …  
    853857 CALL histdef2d(iff,o_swsrfcs_ant%flag,o_swsrfcs_ant%name, "Anthropogenic aerosol radiative forcing clear-sky at SRF", "W/m2")
    854858
     859 IF (.NOT. aerosol_couple) THEN
    855860 CALL histdef2d(iff,o_swtoacf_nat%flag,o_swtoacf_nat%name, "Natural aerosol impact on cloud radiative forcing at TOA", "W/m2")
    856861 CALL histdef2d(iff,o_swsrfcf_nat%flag,o_swsrfcf_nat%name, "Natural aerosol impact on cloud radiative forcing  at SRF", "W/m2")
    … …  
    859864 CALL histdef2d(iff,o_swtoacf_zero%flag,o_swtoacf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at TOA", "W/m2")
    860865 CALL histdef2d(iff,o_swsrfcf_zero%flag,o_swsrfcf_zero%name, "Cloud radiative forcing (allsky-clearsky fluxes) at SRF", "W/m2")
     866 ENDIF
    861867
    862868 ENDIF

  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h

    r1263 r1266  
    834834
    835835! OD550 per species
    836       IF (new_aod) THEN
     836      IF (new_aod .and. (.not. aerosol_couple)) THEN
    837837      DO naero = 1, naero_spc
    838838          IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN
    … …  
    898898
    899899!cf
     900
     901        if (.not. aerosol_couple) then
    900902        IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN
    901903        CALL histwrite_phy(nid_files(iff),o_swtoacf_nat%name,itau_w,
    … …  
    927929     $      solswcf_aero(:,3))
    928930        ENDIF
     931        endif
    929932
    930933        endif ! new_aod

  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/physiq.F

    r1263 r1266  
    14571457     &                        ctetaSTD,dtime,ok_veget,
    14581458     &                        type_ocean,iflag_pbl,ok_mensuel,ok_journe,
    1459      &              ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz)
     1459     &                        ok_hf,ok_instan,ok_LES,ok_ade,ok_aie,
     1460     &                        read_climoz, new_aod, aerosol_couple)
    14601461c$OMP END MASTER
    14611462c$OMP BARRIER
    … …  
    15331534     $                   pdtphys,
    15341535     $                   annee_ref,
    1535      $                   day_ini)
     1536     $                   day_ref,
     1537     $                   itau_phy)
    15361538
    15371539         CALL VTe(VTinca)
    … …  
    28262828            CALL AEROSOL_METEO_CALC(
    28272829     $           calday,pdtphys,pplay,paprs,t,pmflxr,pmflxs,
    2828      $           prfl,psfl,pctsrf,airephy,xjour,rlat,rlon,u10m,v10m)
     2830     $           prfl,psfl,pctsrf,airephy,rlat,rlon,u10m,v10m)
    28292831         END IF
    28302832
    … …  
    33853387     $                        day_ini,
    33863388     $                        airephy,
    3387      $                        xjour,
    33883389     $                        pphi,
    33893390     $                        pphis,
Note: See TracChangeset for help on using the changeset viewer.