Changeset 1266 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90
- Timestamp:
- Nov 20, 2009, 5:27:40 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90
r1254 r1266 144 144 145 145 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. 148 147 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. 151 149 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. 154 151 155 152 REAL, DIMENSION(klon,klev,naero_spc) :: mass_temp … … 308 305 ! le 12 AVRIL 2006 309 306 ! 310 DATA alpha_aers_5wv/ &307 DATA alpha_aers_5wv/ & 311 308 ! bc soluble 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 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, & 360 357 0.52,0.93,1.08,1.24,1.43,2.47, & 361 358 2.73,3.20,3.99,4.64,5.84,9.04, & … … 364 361 0.52,0.94,1.09,1.26,1.45,2.51, & 365 362 2.78,3.25,4.06,4.72,5.94,9.22, & 366 367 368 369 370 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, & 371 368 0.52,0.93,1.08,1.24,1.43,2.47, & 372 369 2.73,3.20,3.99,4.64,5.84,9.04, & … … 375 372 0.52,0.94,1.09,1.26,1.45,2.51, & 376 373 2.78,3.25,4.06,4.72,5.94,9.22, & 377 374 ! seasalt soluble accumulation (computed below for 550nm) 378 375 4.28, 7.17, 8.44, 9.85,11.60,22.44, & 379 376 25.34,30.54,39.38,46.52,59.33,91.77, & 380 381 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, & 382 379 2.48, 4.22, 5.02, 5.94, 7.11,15.29, & 383 380 17.70,22.31,30.73,38.06,52.15,90.59, & … … 386 383 1.47, 2.59, 3.12, 3.74, 4.54, 10.42, & 387 384 12.24,15.82,22.66,28.91,41.54,79.33/ 388 385 389 386 DATA alpha_aeri_5wv/ & 390 391 392 393 394 395 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 ! 397 394 DATA cg_aers_5wv/ & 398 395 ! bc soluble 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 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 477 474 DATA cg_aeri_5wv/& 478 475 ! dust insoluble … … 485 482 DATA piz_aers_5wv/& 486 483 ! 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 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 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 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 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 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 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 ! 565 562 DATA piz_aeri_5wv/& 566 563 ! dust insoluble … … 808 805 !CDIR ON_ADB(alpha_aers_5wv) 809 806 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 810 812 tau_ae5wv_int = alpha_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* & 811 813 (alpha_aers_5wv(RH_num+1,la,spsol) - &
Note: See TracChangeset
for help on using the changeset viewer.