Changeset 1267
- Timestamp:
- Nov 23, 2009, 5:37:51 PM (14 years ago)
- Location:
- LMDZ4/branches/LMDZ4-dev/libf/phylmd
- Files:
-
- 2 edited
-
aeropt_2bands.F90 (modified) (9 diffs)
-
aeropt_5wv.F90 (modified) (10 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_2bands.F90
r1250 r1267 944 944 tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+ & 945 945 tau_ae(i,k,id_SSSSM,inu)+ tau_ae(i,k,id_CIDUSTM,inu) 946 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e- 20)946 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 947 947 948 948 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ & … … 957 957 tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)) & 958 958 /tau_allaer(i,k,mrfspecies,inu) 959 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu), 1e-20)959 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 960 960 961 961 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu)+ & … … 988 988 tau_ae_pi(i,k,id_SSSSM,inu)+ & 989 989 tau_ae_pi(i,k,id_CIDUSTM,inu) 990 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e- 20)990 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 991 991 992 992 piz_allaer(i,k,mrfspecies,inu)=(tau_ae_pi(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ & … … 1001 1001 tau_ae_pi(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)) & 1002 1002 /tau_allaer(i,k,mrfspecies,inu) 1003 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu), 1e-20)1003 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 1004 1004 1005 1005 cg_allaer(i,k,mrfspecies,inu)=(& … … 1023 1023 DO i=1, KLON 1024 1024 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu) 1025 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e- 20)1025 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 1026 1026 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) & 1027 1027 +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu))/ & 1028 1028 tau_allaer(i,k,mrfspecies,inu) 1029 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu), 1e-20)1029 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 1030 1030 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) *cg_ae(i,k,id_ASBCM,inu)& 1031 1031 +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)*cg_ae(i,k,id_AIBCM,inu))/ & … … 1039 1039 DO i=1, KLON 1040 1040 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu) 1041 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e- 20)1041 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 1042 1042 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) & 1043 1043 +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu))/ & 1044 1044 tau_allaer(i,k,mrfspecies,inu) 1045 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu), 1e-20)1045 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 1046 1046 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) *cg_ae(i,k,id_CSSO4M,inu)& 1047 1047 +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu))/ & … … 1055 1055 DO i=1, KLON 1056 1056 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu) 1057 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e- 20)1057 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 1058 1058 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) & 1059 1059 +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu))/ & 1060 1060 tau_allaer(i,k,mrfspecies,inu) 1061 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu), 1e-20)1061 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 1062 1062 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) *cg_ae(i,k,id_ASPOMM,inu)& 1063 1063 +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)*cg_ae(i,k,id_AIPOMM,inu))/ & … … 1071 1071 DO i=1, KLON 1072 1072 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM,inu) 1073 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e- 20)1073 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 1074 1074 piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM,inu) 1075 1075 cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM,inu) … … 1082 1082 DO i=1, KLON 1083 1083 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu) 1084 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e- 20)1084 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-5) 1085 1085 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) & 1086 1086 +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) & 1087 1087 +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ & 1088 1088 tau_allaer(i,k,mrfspecies,inu) 1089 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu), 1e-20)1089 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),0.1) 1090 1090 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) *cg_ae(i,k,id_ASSSM,inu)& 1091 1091 +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu) & -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90
r1266 r1267 76 76 INTEGER :: i, k, ierr, m 77 77 INTEGER :: spsol, spinsol, spss, la 78 INTEGER :: RH_num 78 INTEGER :: RH_num(klon,klev) 79 79 INTEGER, PARAMETER :: la443 = 1 80 80 INTEGER, PARAMETER :: la550 = 2 … … 144 144 145 145 REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble) ! ext. coeff. Soluble comp. units *** m2/g 146 REAL :: alpha_aeri_5wv(las,naero_insoluble) ! ext. coeff. Insoluble comp. 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 147 148 REAL :: cg_aers_5wv(nbre_RH,las,naero_soluble) ! Asym. param. soluble comp. 148 REAL :: cg_aeri_5wv(las,naero_insoluble) ! Asym. param. insoluble 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 149 151 REAL :: piz_aers_5wv(nbre_RH,las,naero_soluble) 150 REAL :: piz_aeri_5wv(las,naero_insoluble) ! Insoluble comp. 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 151 154 152 155 REAL, DIMENSION(klon,klev,naero_spc) :: mass_temp … … 159 162 REAL :: fact_RH(nbre_RH) 160 163 LOGICAL :: used_tau(naero_spc) 164 INTEGER :: n 161 165 162 166 DATA presnivs_19/& … … 305 309 ! le 12 AVRIL 2006 306 310 ! 307 311 DATA alpha_aers_5wv/ & 308 312 ! bc soluble 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 soluble320 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, &313 7.930,7.930,7.930,7.930,7.930,7.930, & 314 7.930,7.930,10.893,12.618,14.550,16.613, & 315 7.658,7.658,7.658,7.658,7.658,7.658, & 316 7.658,7.658,10.351,11.879,13.642,15.510, & 317 7.195,7.195,7.195,7.195,7.195,7.195, & 318 7.195,7.195,9.551,10.847,12.381,13.994, & 319 6.736,6.736,6.736,6.736,6.736,6.736, & 320 6.736,6.736,8.818,9.938,11.283,12.687, & 321 6.277,6.277,6.277,6.277,6.277,6.277, & 322 6.277,6.277,8.123,9.094,10.275,11.501, & 323 ! pom soluble 324 6.676,6.676,6.676,6.676,6.710,6.934, & 325 7.141,7.569,8.034,8.529,9.456,10.511, & 326 5.109,5.109,5.109,5.109,5.189,5.535, & 327 5.960,6.852,8.008,9.712,12.897,19.676, & 328 3.718,3.718,3.718,3.718,3.779,4.042, & 329 4.364,5.052,5.956,7.314,9.896,15.688, & 330 2.849,2.849,2.849,2.849,2.897,3.107, & 331 3.365,3.916,4.649,5.760,7.900,12.863, & 332 2.229,2.229,2.229,2.229,2.268,2.437, & 333 2.645,3.095,3.692,4.608,6.391,10.633, & 334 ! Sulfate (Accumulation) 335 5.751,6.215,6.690,7.024,7.599,8.195, & 336 9.156,10.355,12.660,14.823,18.908,24.508, & 337 4.320,4.675,5.052,5.375,5.787,6.274, & 338 7.066,8.083,10.088,12.003,15.697,21.133, & 339 3.079,3.351,3.639,3.886,4.205,4.584, & 340 5.206,6.019,7.648,9.234,12.391,17.220, & 341 2.336,2.552,2.781,2.979,3.236,3.540, & 342 4.046,4.711,6.056,7.388,10.093,14.313, & 343 1.777,1.949,2.134,2.292,2.503,2.751, & 344 3.166,3.712,4.828,5.949,8.264,11.922, & 345 ! Sulfate (Coarse) 346 5.751,6.215,6.690,7.024,7.599,8.195, & 347 9.156,10.355,12.660,14.823,18.908,24.508, & 348 4.320,4.675,5.052,5.375,5.787,6.274, & 349 7.066,8.083,10.088,12.003,15.697,21.133, & 350 3.079,3.351,3.639,3.886,4.205,4.584, & 351 5.206,6.019,7.648,9.234,12.391,17.220, & 352 2.336,2.552,2.781,2.979,3.236,3.540, & 353 4.046,4.711,6.056,7.388,10.093,14.313, & 354 1.777,1.949,2.134,2.292,2.503,2.751, & 355 3.166,3.712,4.828,5.949,8.264,11.922, & 356 ! Seasalt soluble super_coarse (computed below for 550nm) 357 0.50,0.90,1.05,1.21,1.40,2.41, & 358 2.66,3.11,3.88,4.52,5.69,8.84, & 359 0.000,0.000,0.000,0.000,0.000,0.000, & 360 0.000,0.000,0.000,0.000,0.000,0.000, & 357 361 0.52,0.93,1.08,1.24,1.43,2.47, & 358 362 2.73,3.20,3.99,4.64,5.84,9.04, & … … 361 365 0.52,0.94,1.09,1.26,1.45,2.51, & 362 366 2.78,3.25,4.06,4.72,5.94,9.22, & 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, &367 ! seasalt soluble coarse (computed below for 550nm) 368 0.50,0.90,1.05,1.21,1.40,2.41, & 369 2.66,3.11,3.88,4.52,5.69,8.84, & 370 0.000,0.000,0.000,0.000,0.000,0.000, & 371 0.000,0.000,0.000,0.000,0.000,0.000, & 368 372 0.52,0.93,1.08,1.24,1.43,2.47, & 369 373 2.73,3.20,3.99,4.64,5.84,9.04, & … … 372 376 0.52,0.94,1.09,1.26,1.45,2.51, & 373 377 2.78,3.25,4.06,4.72,5.94,9.22, & 374 ! seasalt soluble accumulation (computed below for 550nm)378 ! seasalt soluble accumulation (computed below for 550nm) 375 379 4.28, 7.17, 8.44, 9.85,11.60,22.44, & 376 380 25.34,30.54,39.38,46.52,59.33,91.77, & 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, &381 0.000,0.000,0.000,0.000,0.000,0.000, & 382 0.000,0.000,0.000,0.000,0.000,0.000, & 379 383 2.48, 4.22, 5.02, 5.94, 7.11,15.29, & 380 384 17.70,22.31,30.73,38.06,52.15,90.59, & … … 383 387 1.47, 2.59, 3.12, 3.74, 4.54, 10.42, & 384 388 12.24,15.82,22.66,28.91,41.54,79.33/ 385 389 386 390 DATA alpha_aeri_5wv/ & 387 ! dust insoluble388 0.759, 0.770, 0.775, 0.775, 0.772, &389 !!jb bc insoluble390 11.536,10.033, 8.422, 7.234, 6.270, &391 ! pom insoluble392 5.042, 3.101, 1.890, 1.294, 0.934/393 !391 ! dust insoluble 392 0.759, 0.770, 0.775, 0.775, 0.772, & 393 !!jb bc insoluble 394 11.536,10.033, 8.422, 7.234, 6.270, & 395 ! pom insoluble 396 5.042, 3.101, 1.890, 1.294, 0.934/ 397 ! 394 398 DATA cg_aers_5wv/ & 395 399 ! bc soluble 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 soluble407 .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 sulfate418 .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 sulfate429 .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 400 .651, .651, .651, .651, .651, .651, & 401 .651, .651, .738, .764, .785, .800, & 402 .597, .597, .597, .597, .597, .597, & 403 .597, .597, .695, .725, .751, .770, & 404 .543, .543, .543, .543, .543, .543, & 405 .543, .543, .650, .684, .714, .736, & 406 .504, .504, .504, .504, .504, .504, & 407 .504, .504, .614, .651, .683, .708, & 408 .469, .469, .469, .469, .469, .469, & 409 .469, .469, .582, .620, .655, .681, & 410 ! pom soluble 411 .679, .679, .679, .679, .683, .691, & 412 .703, .720, .736, .751, .766, .784, & 413 .656, .656, .656, .656, .659, .669, & 414 .681, .699, .717, .735, .750, .779, & 415 .623, .623, .623, .623, .627, .637, & 416 .649, .668, .688, .709, .734, .762, & 417 .592, .592, .592, .592, .595, .605, & 418 .618, .639, .660, .682, .711, .743, & 419 .561, .561, .561, .561, .565, .575, & 420 .588, .609, .632, .656, .688, .724, & 421 ! Accumulation sulfate 422 .671, .684, .697, .704, .714, .723, & 423 .734, .746, .762, .771, .781, .789, & 424 .653, .666, .678, .687, .697, .707, & 425 .719, .732, .751, .762, .775, .789, & 426 .622, .635, .648, .657, .667, .678, & 427 .691, .705, .728, .741, .758, .777, & 428 .591, .604, .617, .627, .638, .650, & 429 .664, .679, .704, .719, .739, .761, & 430 .560, .574, .587, .597, .609, .621, & 431 .637, .653, .680, .697, .719, .745, & 432 ! Coarse sulfate 433 .671, .684, .697, .704, .714, .723, & 434 .734, .746, .762, .771, .781, .789, & 435 .653, .666, .678, .687, .697, .707, & 436 .719, .732, .751, .762, .775, .789, & 437 .622, .635, .648, .657, .667, .678, & 438 .691, .705, .728, .741, .758, .777, & 439 .591, .604, .617, .627, .638, .650, & 440 .664, .679, .704, .719, .739, .761, & 441 .560, .574, .587, .597, .609, .621, & 442 .637, .653, .680, .697, .719, .745, & 443 ! For super coarse seasalt (computed below for 550nm!) 444 0.730,0.753,0.760,0.766,0.772,0.793, & 445 0.797,0.802,0.809,0.813,0.820,0.830, & 446 0.000,0.000,0.000,0.000,0.000,0.000, & 447 0.000,0.000,0.000,0.000,0.000,0.000, & 448 0.721,0.744,0.750,0.756,0.762,0.784, & 449 0.787,0.793,0.800,0.804,0.811,0.822, & 450 0.717,0.741,0.747,0.753,0.759,0.780, & 451 0.784,0.789,0.795,0.800,0.806,0.817, & 452 0.715,0.739,0.745,0.751,0.757,0.777, & 453 0.781,0.786,0.793,0.797,0.803,0.814, & 454 ! For coarse-soluble seasalt (computed below for 550nm!) 455 0.730,0.753,0.760,0.766,0.772,0.793, & 456 0.797,0.802,0.809,0.813,0.820,0.830, & 457 0.000,0.000,0.000,0.000,0.000,0.000, & 458 0.000,0.000,0.000,0.000,0.000,0.000, & 459 0.721,0.744,0.750,0.756,0.762,0.784, & 460 0.787,0.793,0.800,0.804,0.811,0.822, & 461 0.717,0.741,0.747,0.753,0.759,0.780, & 462 0.784,0.789,0.795,0.800,0.806,0.817, & 463 0.715,0.739,0.745,0.751,0.757,0.777, & 464 0.781,0.786,0.793,0.797,0.803,0.814, & 465 ! accumulation-seasalt soluble (computed below for 550nm!) 466 0.698,0.722,0.729,0.736,0.743,0.765, & 467 0.768,0.773,0.777,0.779,0.781,0.779, & 468 0.000,0.000,0.000,0.000,0.000,0.000, & 469 0.000,0.000,0.000,0.000,0.000,0.000, & 470 0.658,0.691,0.701,0.710,0.720,0.756, & 471 0.763,0.771,0.782,0.788,0.795,0.801, & 472 0.632,0.668,0.679,0.690,0.701,0.743, & 473 0.750,0.762,0.775,0.783,0.792,0.804, & 474 0.605,0.644,0.656,0.669,0.681,0.729, & 475 0.737,0.750,0.765,0.775,0.787,0.803/ 476 ! 477 474 478 DATA cg_aeri_5wv/& 475 479 ! dust insoluble … … 482 486 DATA piz_aers_5wv/& 483 487 ! bc soluble 484 485 486 487 488 489 490 491 492 493 494 ! pom soluble495 496 497 498 499 500 501 502 503 504 505 ! sulfate soluble accumulation506 507 508 509 510 511 512 513 514 515 516 ! sulfate soluble coarse517 518 519 520 521 522 523 524 525 526 527 ! seasalt super coarse (computed below for 550nm)528 529 530 531 532 533 534 535 536 537 538 ! seasalt coarse (computed below for 550nm)539 540 541 542 543 544 545 546 547 548 549 ! seasalt soluble accumulation (computed below for 550nm)550 551 552 553 554 555 556 557 558 559 560 561 488 .445, .445, .445, .445, .445, .445, & 489 .445, .445, .470, .487, .508, .531, & 490 .442, .442, .442, .442, .442, .442, & 491 .442, .442, .462, .481, .506, .533, & 492 .427, .427, .427, .427, .427, .427, & 493 .427, .427, .449, .470, .497, .526, & 494 .413, .413, .413, .413, .413, .413, & 495 .413, .413, .437, .458, .486, .516, & 496 .399, .399, .399, .399, .399, .399, & 497 .399, .399, .423, .445, .473, .506, & 498 ! pom soluble 499 .975, .975, .975, .975, .975, .977, & 500 .979, .982, .984, .987, .990, .994, & 501 .972, .972, .972, .972, .973, .974, & 502 .977, .980, .983, .986, .989, .993, & 503 .963, .963, .963, .963, .964, .966, & 504 .969, .974, .977, .982, .986, .991, & 505 .955, .955, .955, .955, .955, .958, & 506 .962, .967, .972, .977, .983, .989, & 507 .944, .944, .944, .944, .944, .948, & 508 .952, .959, .962, .972, .979, .987, & 509 ! sulfate soluble accumulation 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 ! sulfate soluble coarse 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 super coarse (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, & 542 ! seasalt coarse (computed below for 550nm) 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 1.000,1.000,1.000,1.000,1.000,1.000, & 553 ! seasalt soluble accumulation (computed below for 550nm) 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 1.000,1.000,1.000,1.000,1.000,1.000/ 564 565 ! 562 566 DATA piz_aeri_5wv/& 563 567 ! dust insoluble … … 680 684 ! 681 685 682 ! !CDIR ON_ADB(RH_tab)686 !CDIR ON_ADB(RH_tab) 683 687 !CDIR ON_ADB(fact_RH) 684 688 !CDIR NOVECTOR 685 DO RH_Num=1,nbre_RH-1686 fact_RH( RH_num)=1./(RH_tab(RH_num+1)-RH_tab(RH_num))689 DO n=1,nbre_RH-1 690 fact_RH(n)=1./(RH_tab(n+1)-RH_tab(n)) 687 691 ENDDO 688 692 689 693 DO k=1, KLEV 690 ! !CDIR ON_ADB(RH_tab)694 !CDIR ON_ADB(RH_tab) 691 695 !CDIR ON_ADB(fact_RH) 692 696 DO i=1, KLON 693 697 rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX) 694 RH_num = INT( rh(i,k)/10. + 1.)695 IF (rh(i,k).GT.85.) RH_num =10696 IF (rh(i,k).GT.90.) RH_num =11697 DELTA(i,k)=(rh(i,k)-RH_tab(RH_num ))*fact_RH(RH_num)698 RH_num(i,k) = INT( rh(i,k)/10. + 1.) 699 IF (rh(i,k).GT.85.) RH_num(i,k)=10 700 IF (rh(i,k).GT.90.) RH_num(i,k)=11 701 DELTA(i,k)=(rh(i,k)-RH_tab(RH_num(i,k)))*fact_RH(RH_num(i,k)) 698 702 ENDDO 699 703 ENDDO … … 805 809 !CDIR ON_ADB(alpha_aers_5wv) 806 810 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 812 tau_ae5wv_int = alpha_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* & 813 (alpha_aers_5wv(RH_num+1,la,spsol) - & 814 alpha_aers_5wv(RH_num,la,spsol)) 811 tau_ae5wv_int = alpha_aers_5wv(RH_num(i,k),la,spsol)+DELTA(i,k)* & 812 (alpha_aers_5wv(RH_num(i,k)+1,la,spsol) - & 813 alpha_aers_5wv(RH_num(i,k),la,spsol)) 815 814 816 815 tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k) &
Note: See TracChangeset
for help on using the changeset viewer.
