Changeset 1267 for LMDZ4/branches


Ignore:
Timestamp:
Nov 23, 2009, 5:37:51 PM (15 years ago)
Author:
yann meurdesoif
Message:

Correction bug aerosol (division / 0 + réécriture de la correction d'A. Cozic)

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

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_2bands.F90

    r1250 r1267  
    944944                                           tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+   &
    945945                                           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)
    947947                 
    948948             piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &
     
    957957                                             tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)) &
    958958                                            /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)
    960960
    961961             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)+ &
     
    988988                      tau_ae_pi(i,k,id_SSSSM,inu)+ &
    989989                      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)
    991991
    992992                 piz_allaer(i,k,mrfspecies,inu)=(tau_ae_pi(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &
     
    10011001                      tau_ae_pi(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)) &
    10021002                      /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)
    10041004
    10051005                 cg_allaer(i,k,mrfspecies,inu)=(&
     
    10231023          DO i=1, KLON
    10241024            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)
    10261026            piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) &
    10271027                      +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu))/ &
    10281028                      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)
    10301030            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)&
    10311031                      +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)*cg_ae(i,k,id_AIBCM,inu))/ &
     
    10391039          DO i=1, KLON
    10401040            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)
    10421042            piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) &
    10431043                      +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu))/ &
    10441044                      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)
    10461046            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)&
    10471047                      +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu))/ &
     
    10551055          DO i=1, KLON
    10561056            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)
    10581058            piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) &
    10591059                      +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu))/ &
    10601060                      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)
    10621062            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)&
    10631063                      +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)*cg_ae(i,k,id_AIPOMM,inu))/ &
     
    10711071          DO i=1, KLON
    10721072            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)
    10741074            piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM,inu)
    10751075            cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM,inu)
     
    10821082          DO i=1, KLON
    10831083            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)
    10851085            piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) &
    10861086                    +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) &
    10871087                    +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ &
    10881088                    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)
    10901090            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)&
    10911091                    +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  
    7676  INTEGER :: i, k, ierr, m
    7777  INTEGER :: spsol, spinsol, spss, la
    78   INTEGER :: RH_num
     78  INTEGER :: RH_num(klon,klev)
    7979  INTEGER, PARAMETER :: la443 = 1
    8080  INTEGER, PARAMETER :: la550 = 2
     
    144144 
    145145  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
    147148  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
    149151  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
    151154
    152155  REAL, DIMENSION(klon,klev,naero_spc) :: mass_temp
     
    159162  REAL :: fact_RH(nbre_RH)
    160163  LOGICAL :: used_tau(naero_spc)
     164  INTEGER :: n
    161165 
    162166  DATA presnivs_19/&
     
    305309!                                   le 12 AVRIL 2006
    306310
    307   DATA alpha_aers_5wv/ &
     311 DATA alpha_aers_5wv/ &
    308312                                ! 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 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, & 
     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, & 
    357361     0.52,0.93,1.08,1.24,1.43,2.47, & 
    358362     2.73,3.20,3.99,4.64,5.84,9.04, & 
     
    361365     0.52,0.94,1.09,1.26,1.45,2.51, & 
    362366     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, & 
    368372     0.52,0.93,1.08,1.24,1.43,2.47, & 
    369373     2.73,3.20,3.99,4.64,5.84,9.04, & 
     
    372376     0.52,0.94,1.09,1.26,1.45,2.51, & 
    373377     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)
    375379     4.28, 7.17, 8.44, 9.85,11.60,22.44,  & 
    376380     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, & 
    379383     2.48, 4.22, 5.02, 5.94, 7.11,15.29,  & 
    380384     17.70,22.31,30.73,38.06,52.15,90.59, & 
     
    383387     1.47, 2.59, 3.12, 3.74, 4.54, 10.42, & 
    384388     12.24,15.82,22.66,28.91,41.54,79.33/
    385  
     389
    386390  DATA alpha_aeri_5wv/ &
    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   !
     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   !
    394398  DATA cg_aers_5wv/ & 
    395399                                 ! 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 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  
     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
    474478  DATA cg_aeri_5wv/&
    475479     ! dust insoluble
     
    482486  DATA piz_aers_5wv/&
    483487                           ! bc soluble
    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   !
     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 !
    562566  DATA piz_aeri_5wv/&
    563567     ! dust insoluble
     
    680684  !     
    681685
    682 !!CDIR ON_ADB(RH_tab)
     686!CDIR ON_ADB(RH_tab)
    683687!CDIR ON_ADB(fact_RH)
    684688!CDIR NOVECTOR
    685   DO RH_Num=1,nbre_RH-1
    686     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))
    687691  ENDDO
    688692   
    689693  DO k=1, KLEV
    690 !!CDIR ON_ADB(RH_tab)
     694!CDIR ON_ADB(RH_tab)
    691695!CDIR ON_ADB(fact_RH)
    692696    DO i=1, KLON
    693697      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=10
    696       IF (rh(i,k).GT.90.) RH_num=11
    697       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))
    698702    ENDDO
    699703  ENDDO
     
    805809!CDIR ON_ADB(alpha_aers_5wv)
    806810            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))
    815814
    816815              tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k)   &
Note: See TracChangeset for help on using the changeset viewer.