Changeset 2737 for LMDZ5/trunk


Ignore:
Timestamp:
Dec 18, 2016, 4:51:56 PM (8 years ago)
Author:
fhourdin
Message:

Retour en arriere sur les dernieres commissions d'Olivier.

Location:
LMDZ5/trunk/libf/phylmd/rrtm
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90

    r2735 r2737  
    8383  INTEGER, PARAMETER :: la865 = 5
    8484  INTEGER, PARAMETER :: nbre_RH=12
    85   INTEGER, PARAMETER :: naero_soluble=9   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.
     85  INTEGER, PARAMETER :: naero_soluble=7   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.
    8686                                          !  4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
    87                                           !  8- NO3 coarse;  9- NO3 acc.
    88   INTEGER, PARAMETER :: naero_insoluble=4 !  1- Dust; 2- BC insoluble; 3- POM insoluble; 4- NO3 insoluble
     87  INTEGER, PARAMETER :: naero_insoluble=3 !  1- Dust; 2- BC insoluble; 3- POM insoluble
    8988
    9089  REAL :: zrho
    … …  
    102101  REAL :: zdh(klon,klev)
    103102 
    104    ! Soluble components 1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse;
    105    !                    5- seasalt super-coarse; 6- seasalt coarse; 7- seasalt acc.
    106    !                    8- NO3 coarse; 9- NO3 acc.
     103   ! Soluble components 1- BC soluble; 2- POM soluble; 3- SO4 acc.; 4- SO4 coarse; 5 seasalt super-coarse; 6 seasalt coarse; 7 seasalt acc.
    107104  REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble)   ! Ext. coeff. ** m2/g
    108    ! Insoluble components 1- Dust: 2- BC; 3- POM; 4- NO3
     105   ! Insoluble components 1- Dust: 2- BC; 3- POM
    109106  REAL :: alpha_aeri_5wv(las,naero_insoluble)         ! Ext. coeff. ** m2/g
    110107
    … …  
    132129! 
    133130 DATA alpha_aers_5wv/ &
    134                            ! BC soluble
     131                                ! bc soluble
    135132       7.930,7.930,7.930,7.930,7.930,7.930,7.930,7.930,10.893,12.618,14.550,16.613, &
    136133       7.658,7.658,7.658,7.658,7.658,7.658,7.658,7.658,10.351,11.879,13.642,15.510, &
    … …  
    138135       6.736,6.736,6.736,6.736,6.736,6.736,6.736,6.736,8.818,9.938,11.283,12.687,   &
    139136       6.277,6.277,6.277,6.277,6.277,6.277,6.277,6.277,8.123,9.094,10.275,11.501,   &
    140                            ! POM soluble
     137                                ! pom soluble
    141138       6.676,6.676,6.676,6.676,6.710,6.934,7.141,7.569,8.034,8.529,9.456,10.511,  &
    142139       5.109,5.109,5.109,5.109,5.189,5.535,5.960,6.852,8.008,9.712,12.897,19.676, &
    … …  
    144141       2.849,2.849,2.849,2.849,2.897,3.107,3.365,3.916,4.649,5.760,7.900,12.863,  &
    145142       2.229,2.229,2.229,2.229,2.268,2.437,2.645,3.095,3.692,4.608,6.391,10.633,  &
    146                            ! Sulfate (Coarse)
     143                                ! Sulfate (Accumulation)
    147144       5.751,6.215,6.690,7.024,7.599,8.195,9.156,10.355,12.660,14.823,18.908,24.508, &
    148145       4.320,4.675,5.052,5.375,5.787,6.274,7.066,8.083,10.088,12.003,15.697,21.133,  &
    … …  
    150147       2.336,2.552,2.781,2.979,3.236,3.540,4.046,4.711,6.056,7.388,10.093,14.313,    &
    151148       1.777,1.949,2.134,2.292,2.503,2.751,3.166,3.712,4.828,5.949,8.264,11.922,     &
    152                            ! Sulfate (Accumulation)
     149                                ! Sulfate (Coarse)
    153150       5.751,6.215,6.690,7.024,7.599,8.195,9.156,10.355,12.660,14.823,18.908,24.508, &
    154151       4.320,4.675,5.052,5.375,5.787,6.274,7.066,8.083,10.088,12.003,15.697,21.133,  &
    … …  
    156153       2.336,2.552,2.781,2.979,3.236,3.540,4.046,4.711,6.056,7.388,10.093,14.313,    &
    157154       1.777,1.949,2.134,2.292,2.503,2.751,3.166,3.712,4.828,5.949,8.264,11.922,     &
    158                            ! Seasalt Super Coarse Soluble (SS)
     155                           ! seasalt seasalt Super Coarse Soluble (SS)
    159156        0.218, 0.272, 0.293, 0.316, 0.343, 0.380, 0.429, 0.501, 0.636, 0.755, 0.967, 1.495, &
    160157        0.221, 0.275, 0.297, 0.320, 0.348, 0.383, 0.432, 0.509, 0.640, 0.759, 0.972, 1.510, &
    … …  
    162159        0.227, 0.282, 0.303, 0.327, 0.356, 0.392, 0.441, 0.518, 0.652, 0.770, 0.987, 1.529, &
    163160        0.230, 0.285, 0.306, 0.330, 0.359, 0.396, 0.446, 0.522, 0.656, 0.777, 0.993, 1.539, &
    164                            ! Seasalt Coarse Soluble (CS)     
     161                           ! seasalt seasalt Coarse Soluble (CS)     
    165162        0.578, 0.706, 0.756, 0.809, 0.876, 0.964, 1.081, 1.256, 1.577, 1.858, 2.366, 3.613, &
    166163        0.598, 0.725, 0.779, 0.833, 0.898, 0.990, 1.109, 1.290, 1.609, 1.889, 2.398, 3.682, &
    … …  
    168165        0.633, 0.767, 0.820, 0.879, 0.948, 1.044, 1.167, 1.353, 1.683, 1.969, 2.491, 3.785, &
    169166        0.648, 0.785, 0.838, 0.896, 0.967, 1.066, 1.192, 1.381, 1.714, 2.006, 2.531, 3.836, &
    170                            ! Seasalt Accumulation Soluble (AS)
     167                           ! seasalt seasalt Accumulation Soluble (AS)
    171168        4.432, 5.899, 6.505, 7.166, 7.964, 7.962, 9.232,11.257,14.979,18.337,24.223,37.811, &
    172169        3.298, 4.569, 5.110, 5.709, 6.446, 6.268, 7.396, 9.246,12.787,16.113,22.197,37.136, &
    173170        2.340, 3.358, 3.803, 4.303, 4.928, 4.696, 5.629, 7.198,10.308,13.342,19.120,34.296, &
    174171        1.789, 2.626, 2.999, 3.422, 3.955, 3.730, 4.519, 5.864, 8.593,11.319,16.653,31.331, &
    175         1.359, 2.037, 2.343, 2.693, 3.139, 2.940, 3.596, 4.729, 7.076, 9.469,14.266,28.043, &
    176                            ! NO3 coarse (CS)
    177         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    178         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    179         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    180         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    181         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    182                            ! NO3 accumulation (AS)
    183         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    184         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    185         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    186         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    187         0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
     172        1.359, 2.037, 2.343, 2.693, 3.139, 2.940, 3.596, 4.729, 7.076, 9.469,14.266,28.043 /
    188173
    189174  DATA alpha_aeri_5wv/ &
    190175                                 ! dust insoluble
    191         0.759, 0.770, 0.775, 0.775, 0.772,  &
    192                                  ! BC insoluble
     176        0.759, 0.770, 0.775, 0.775, 0.772, &
     177                                 !!jb bc insoluble
    193178        11.536,10.033, 8.422, 7.234, 6.270, &
    194                                  ! POM insoluble
    195         5.042, 3.101, 1.890, 1.294, 0.934,  &
    196                                  ! NO3 insoluble
    197         0.0, 0.0, 0.0, 0.0, 0.0 /
     179                                 ! pom insoluble
     180        5.042, 3.101, 1.890, 1.294, 0.934/
    198181
    199182  ! internal mixture
    … …  
    281264     aerosol_name(9) = id_CIDUSTM_phy
    282265     aerosol_name(10)= id_CSSO4M_phy
    283      aerosol_name(11)= id_CSNO3M_phy
    284      aerosol_name(12)= id_ASNO3M_phy
    285      aerosol_name(13)= id_CINO3M_phy
    286266  ENDIF
    287267
    … …  
    318298        soluble=.TRUE.
    319299        spsol=2
    320     ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
     300    ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
    321301        soluble=.TRUE.
    322302        spsol=3
    323303        fac=1.375    ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD
    324     ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
     304    ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
    325305        soluble=.TRUE.
    326306        spsol=4
    … …  
    335315        soluble=.TRUE.
    336316        spsol=7
    337     ELSEIF (aerosol_name(m).EQ.id_CSNO3M_phy) THEN
    338         soluble=.TRUE.
    339         spsol=8
    340         fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
    341      ELSEIF (aerosol_name(m).EQ.id_ASNO3M_phy) THEN
    342         soluble=.TRUE.
    343         spsol=9
    344         fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
    345317    ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN
    346318        soluble=.FALSE.
    … …  
    352324        soluble=.FALSE.
    353325        spinsol=3
    354     ELSEIF (aerosol_name(m).EQ.id_CINO3M_phy) THEN
    355         soluble=.FALSE.
    356         spinsol=4
    357         fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
    358326    ELSE
    359327        CYCLE

  • LMDZ5/trunk/libf/phylmd/rrtm/aeropt_6bands_rrtm.F90

    r2736 r2737  
    33!
    44SUBROUTINE AEROPT_6BANDS_RRTM ( &
    5      pdel, m_allaer, RHcl, tau_allaer, piz_allaer, &
    6      cg_allaer, m_allaer_pi, flag_aerosol, flag_bc_internal_mixture, zrho )
     5     pdel, m_allaer, RHcl, &
     6     tau_allaer, piz_allaer, &
     7     cg_allaer, m_allaer_pi, &
     8     flag_aerosol, flag_bc_internal_mixture, zrho )
    79
    810  USE dimphy
    … …  
    5153  INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA
    5254
    53   INTEGER, PARAMETER :: naero_soluble=9    ! 1- BC soluble; 2- POM soluble; 3- SO4. acc. 4- SO4 coarse
    54                                            ! 5- seasalt super coarse  6- seasalt coarse  7- seasalt acc.
    55                                            ! 8- NO3 coarse 9- NO3 acc.
    56   INTEGER, PARAMETER :: naero_insoluble=4  ! 1- Dust; 2- BC insoluble; 3- POM insoluble; 4- NO3 insoluble
     55  INTEGER, PARAMETER :: naero_soluble=7    ! 1- BC soluble; 2- POM soluble; 3- SO4. acc. 4- SO4 coarse
     56                                           ! 5- seasalt super coarse  6- seasalt coarse   7- seasalt acc.
     57  INTEGER, PARAMETER :: naero_insoluble=3  ! 1- Dust; 2- BC insoluble; 3- POM insoluble
    5758
    5859  INTEGER, PARAMETER :: nbre_RH=12
    … …  
    6667  REAL :: fac
    6768  REAL :: zdh(klon,klev)
    68   INTEGER, ALLOCATABLE, DIMENSION(:) :: aerosol_name
     69  INTEGER, ALLOCATABLE, DIMENSION(:)   :: aerosol_name
    6970  INTEGER :: nb_aer
    7071
    71   REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: tau_ae
    72   REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: piz_ae
    73   REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: cg_ae
    74 
    75   REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) :: tau_ae_pi
     72  REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) ::  tau_ae
     73  REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) ::  piz_ae
     74  REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) ::  cg_ae
     75
     76  REAL, DIMENSION(klon,klev,naero_tot,nbands_sw_rrtm) ::  tau_ae_pi
    7677  REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: piz_ae_pi
    7778  REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: cg_ae_pi
    … …  
    7980  ! Proprietes optiques
    8081  !
    81   REAL:: alpha_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble) !--unit m2/g
    82   REAL:: alpha_aeri_6bands(nbands_sw_rrtm,naero_insoluble)       !--unit m2/g
    83   REAL:: cg_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)    !--unitless
    84   REAL:: cg_aeri_6bands(nbands_sw_rrtm,naero_insoluble)          !--unitless
    85   REAL:: piz_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)   !-- unitless
    86   REAL:: piz_aeri_6bands(nbands_sw_rrtm,naero_insoluble)         !-- unitless
     82  REAL:: alpha_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)   !--unit m2/g SO4
     83  REAL:: alpha_aeri_6bands(nbands_sw_rrtm,naero_insoluble)
     84  REAL:: cg_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)      !--unit
     85  REAL:: cg_aeri_6bands(nbands_sw_rrtm,naero_insoluble)
     86  REAL:: piz_aers_6bands(nbre_RH,nbands_sw_rrtm,naero_soluble)     !-- unit
     87  REAL:: piz_aeri_6bands(nbands_sw_rrtm,naero_insoluble)        !-- unit
    8788  !
    8889  ! BC internal mixture
    8990  !
    9091  INTEGER, PARAMETER ::  nbclassbc = 5  ! Added by Rong Wang/OB for the 5 fractions
    91                                         ! of BC in the soluble mode:
    92                                         ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/
     92                                       ! of BC in the soluble mode:
     93                                       ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/
    9394  ! for Maxwell-Garnet internal mixture
    9495  ! Detailed theory can be found in R. Wang Estimation of global black carbon direct
    … …  
    99100  REAL :: piz_MG_6bands(nbre_RH,nbands_sw_rrtm,nbclassbc)
    100101  !
    101   INTEGER :: aerindex, classbc, classbc_pi
     102  INTEGER :: id, classbc, classbc_pi
    102103  REAL :: tmp_var, tmp_var_pi, BC_massfra, BC_massfra_pi
    103104
    … …  
    125126  0.377, 0.377, 0.377, 0.377, 0.377, 0.415, 0.456, 0.542, 0.627, 0.636, 0.646, 0.738, &
    126127  0.052, 0.052, 0.052, 0.052, 0.052, 0.057, 0.063, 0.075, 0.087, 0.088, 0.089, 0.102, &
     128       ! sulfate AS   
     129  6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, &
     130  6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, &
     131  4.381, 4.381, 4.381, 4.828, 5.301, 5.792, 6.309, 7.176, 9.542,11.868,11.868,21.204, &
     132  1.727, 1.727, 1.727, 1.903, 2.090, 2.283, 2.487, 2.829, 3.762, 4.679, 4.679, 8.359, &
     133  0.312, 0.312, 0.312, 0.344, 0.378, 0.413, 0.450, 0.511, 0.680, 0.846, 0.846, 1.511, &
     134  0.121, 0.121, 0.121, 0.134, 0.147, 0.161, 0.175, 0.199, 0.264, 0.329, 0.329, 0.588, &
    127135       ! sulfate coarse CS
    128136  0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
    … …  
    132140  0.711, 0.711, 0.711, 0.783, 0.860, 0.939, 1.023, 1.164, 1.548, 1.925, 1.925, 3.439, &
    133141  0.602, 0.602, 0.602, 0.664, 0.729, 0.796, 0.867, 0.986, 1.312, 1.631, 1.631, 2.915, &
    134        ! sulfate AS   
    135   6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, &
    136   6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, &
    137   4.381, 4.381, 4.381, 4.828, 5.301, 5.792, 6.309, 7.176, 9.542,11.868,11.868,21.204, &
    138   1.727, 1.727, 1.727, 1.903, 2.090, 2.283, 2.487, 2.829, 3.762, 4.679, 4.679, 8.359, &
    139   0.312, 0.312, 0.312, 0.344, 0.378, 0.413, 0.450, 0.511, 0.680, 0.846, 0.846, 1.511, &
    140   0.121, 0.121, 0.121, 0.134, 0.147, 0.161, 0.175, 0.199, 0.264, 0.329, 0.329, 0.588, &
    141         ! Seasalt Super Coarse Soluble (SS)
     142        ! seasalt seasalt Super Coarse Soluble (SS)
    142143  0.214, 0.267, 0.287, 0.310, 0.337, 0.373, 0.421, 0.494, 0.625, 0.742, 0.956, 1.480, &
    143144  0.217, 0.270, 0.291, 0.314, 0.341, 0.377, 0.426, 0.499, 0.632, 0.751, 0.963, 1.490, &
    … …  
    146147  0.250, 0.307, 0.330, 0.354, 0.384, 0.424, 0.477, 0.556, 0.696, 0.822, 1.044, 1.592, &
    147148  0.279, 0.347, 0.373, 0.401, 0.434, 0.478, 0.537, 0.625, 0.781, 0.918, 1.158, 1.744, &
    148         ! Seasalt Coarse Soluble (CS)     
     149        ! seasalt seasalt Coarse Soluble (CS)     
    149150  0.550, 0.676, 0.724, 0.779, 0.841, 0.925, 1.040, 1.213, 1.523, 1.803, 2.306, 3.535, &
    150151  0.568, 0.695, 0.744, 0.798, 0.863, 0.950, 1.066, 1.240, 1.556, 1.839, 2.342, 3.588, &
    … …  
    153154  0.663, 0.825, 0.889, 0.957, 1.038, 1.126, 1.268, 1.483, 1.862, 2.196, 2.780, 4.183, &
    154155  0.479, 0.644, 0.711, 0.785, 0.873, 0.904, 1.048, 1.275, 1.692, 2.072, 2.753, 4.430, &
    155         ! Seasalt Accumulation Soluble (AS)
     156        ! seasalt seasalt Accumulation Soluble (AS)
    156157  6.128, 7.564, 8.127, 8.723, 9.421,10.014,11.250,13.113,16.320,19.061,23.675,34.158, &
    157158  5.080, 6.586, 7.197, 7.856, 8.645, 8.820,10.115,12.140,15.775,18.983,24.499,36.993, &
    … …  
    159160  1.389, 2.066, 2.371, 2.719, 3.161, 2.970, 3.623, 4.747, 7.064, 9.419,14.120,27.565, &
    160161  0.309, 0.497, 0.587, 0.693, 0.832, 0.777, 0.986, 1.364, 2.207, 3.136, 5.174,12.065, &
    161   0.021, 0.037, 0.045, 0.054, 0.067, 0.065, 0.087, 0.129, 0.230, 0.353, 0.654, 1.885, &
    162         ! NO3 coarse (CS)
    163   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    164   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    165   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    166   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    167   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    168   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    169         ! NO3 accumulation (AS)
    170   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    171   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    172   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    173   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    174   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    175   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
     162  0.021, 0.037, 0.045, 0.054, 0.067, 0.065, 0.087, 0.129, 0.230, 0.353, 0.654, 1.885  /
    176163
    177164  DATA alpha_aeri_6bands/  &
    178165       ! dust insoluble CI
    179166  0.751, 0.751, 0.769, 0.772, 0.672, 0.437, &
    180        ! BC insoluble AI
     167       ! bc insoluble AI
    181168  6.497, 6.497, 5.900, 4.284, 2.163, 0.966, &
    182        ! POM insoluble AI
    183   6.443, 6.443, 4.381, 1.846, 0.377, 0.052, &
    184        ! NO3 insoluble CI
    185   0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
     169       ! pom insoluble AI
     170  6.443, 6.443, 4.381, 1.846, 0.377, 0.052 /
    186171
    187172  DATA cg_aers_6bands/ &
    … …  
    200185  0.363, 0.363, 0.363, 0.363, 0.363, 0.360, 0.357, 0.352, 0.350, 0.349, 0.349, 0.347, &
    201186  0.142, 0.142, 0.142, 0.142, 0.142, 0.139, 0.137, 0.133, 0.131, 0.131, 0.130, 0.129, &
     187       ! sulfate AS
     188  0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, &
     189  0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, &
     190  0.653, 0.653, 0.653, 0.662, 0.670, 0.676, 0.683, 0.692, 0.701, 0.706, 0.709, 0.716, &
     191  0.563, 0.563, 0.563, 0.565, 0.567, 0.569, 0.570, 0.573, 0.575, 0.576, 0.577, 0.579, &
     192  0.362, 0.362, 0.362, 0.359, 0.356, 0.354, 0.352, 0.348, 0.345, 0.343, 0.342, 0.340, &
     193  0.137, 0.137, 0.137, 0.135, 0.133, 0.132, 0.130, 0.128, 0.126, 0.125, 0.124, 0.122, &
    202194       ! sulfate coarse CS
    203195  0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
    … …  
    207199  0.810, 0.810, 0.810, 0.794, 0.781, 0.770, 0.759, 0.743, 0.728, 0.719, 0.714, 0.702, &
    208200  0.803, 0.803, 0.803, 0.790, 0.779, 0.770, 0.762, 0.748, 0.736, 0.729, 0.725, 0.715, &
    209        ! sulfate AS
    210   0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, &
    211   0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, &
    212   0.653, 0.653, 0.653, 0.662, 0.670, 0.676, 0.683, 0.692, 0.701, 0.706, 0.709, 0.716, &
    213   0.563, 0.563, 0.563, 0.565, 0.567, 0.569, 0.570, 0.573, 0.575, 0.576, 0.577, 0.579, &
    214   0.362, 0.362, 0.362, 0.359, 0.356, 0.354, 0.352, 0.348, 0.345, 0.343, 0.342, 0.340, &
    215   0.137, 0.137, 0.137, 0.135, 0.133, 0.132, 0.130, 0.128, 0.126, 0.125, 0.124, 0.122, &
    216201        ! seasalt seasalt Super Coarse Soluble (SS)
    217202  0.797, 0.800, 0.801, 0.802, 0.804, 0.822, 0.825, 0.828, 0.832, 0.835, 0.838, 0.843, &
    … …  
    234219  0.561, 0.585, 0.593, 0.601, 0.609, 0.637, 0.651, 0.669, 0.691, 0.705, 0.723, 0.745, &
    235220  0.392, 0.427, 0.439, 0.451, 0.464, 0.480, 0.500, 0.526, 0.563, 0.588, 0.621, 0.671, &
    236   0.144, 0.170, 0.179, 0.189, 0.201, 0.207, 0.224, 0.248, 0.285, 0.315, 0.359, 0.439, &
    237         ! NO3 coarse (CS)
    238   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    239   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    240   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    241   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    242   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    243   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    244         ! NO3 accumulation (AS)
    245   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    246   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    247   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    248   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    249   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,  &
    250   0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
     221  0.144, 0.170, 0.179, 0.189, 0.201, 0.207, 0.224, 0.248, 0.285, 0.315, 0.359, 0.439 /
    251222
    252223  DATA cg_aeri_6bands/ &
    253224       ! dust insoluble CI
    254225  0.718, 0.718, 0.699, 0.661, 0.676, 0.670, &
    255        ! BC insoluble AI
     226       ! bc insoluble AI
    256227  0.721, 0.721, 0.643, 0.513, 0.321, 0.153, &
    257        ! POM insoluble AI
    258   0.687, 0.687, 0.658, 0.564, 0.363, 0.142, &
    259        ! NO3 insoluble CI
    260   0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
     228       ! pom insoluble AI
     229  0.687, 0.687, 0.658, 0.564, 0.363, 0.142  /
    261230
    262231  DATA piz_aers_6bands/&
    263        ! BC soluble AS
     232       ! bc soluble AS
    264233  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
    265234  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
    … …  
    275244  0.816, 0.816, 0.816, 0.816, 0.816, 0.839, 0.859, 0.888, 0.908, 0.910, 0.911, 0.925, &
    276245  0.663, 0.663, 0.663, 0.663, 0.663, 0.607, 0.562, 0.492, 0.446, 0.441, 0.437, 0.404, &
     246       ! sulfate AS
     247  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
     248  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
     249  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
     250  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
     251  0.988, 0.988, 0.988, 0.989, 0.990, 0.990, 0.991, 0.992, 0.993, 0.993, 0.994, 0.994, &
     252  0.256, 0.256, 0.256, 0.263, 0.268, 0.273, 0.277, 0.284, 0.290, 0.294, 0.296, 0.301, &
    277253       ! sulfate coarse CS
    278254  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    … …  
    282258  0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, &
    283259  0.877, 0.877, 0.877, 0.873, 0.870, 0.867, 0.864, 0.860, 0.856, 0.854, 0.852, 0.849, &
    284        ! sulfate AS
    285   1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    286   1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    287   1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    288   1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    289   0.988, 0.988, 0.988, 0.989, 0.990, 0.990, 0.991, 0.992, 0.993, 0.993, 0.994, 0.994, &
    290   0.256, 0.256, 0.256, 0.263, 0.268, 0.273, 0.277, 0.284, 0.290, 0.294, 0.296, 0.301, &
    291260        ! seasalt seasalt Super Coarse Soluble (SS)
    292261  0.983, 0.982, 0.982, 0.982, 0.982, 0.992, 0.993, 0.994, 0.995, 0.996, 0.997, 0.998, &
    … …  
    309278  0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    310279  0.997, 0.998, 0.998, 0.998, 0.998, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, &
    311   0.985, 0.989, 0.990, 0.990, 0.991, 0.996, 0.997, 0.998, 0.999, 0.999, 0.999, 1.000,  &
    312         ! NO3 coarse (CS)
    313   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    314   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    315   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    316   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    317   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    318   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    319         ! NO3 accumulation (AS)
    320   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    321   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    322   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    323   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    324   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0,  &
    325   1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /
     280  0.985, 0.989, 0.990, 0.990, 0.991, 0.996, 0.997, 0.998, 0.999, 0.999, 0.999, 1.000 /
    326281
    327282  DATA piz_aeri_6bands/ &
    328283       ! dust insoluble CI
    329284  0.954, 0.954, 0.965, 0.981, 0.996, 0.990, &
    330        ! BC insoluble AI
     285       ! bc insoluble AI
    331286  0.460, 0.460, 0.445, 0.394, 0.267, 0.121, &
    332        ! POM insoluble AI
    333   0.973, 0.973, 0.972, 0.940, 0.816, 0.663, &
    334        ! NO3 insoluble CI
    335   1.0, 1.0, 1.0, 1.0, 1.0, 1.0 /
     287       ! pom insoluble AI
     288  0.973, 0.973, 0.972, 0.940, 0.816, 0.663 /
    336289
    337290! Added by R. Wang (July 31 2016)
    … …  
    487440     aerosol_name(1) = id_CIDUSTM_phy
    488441  ELSEIF (flag_aerosol .EQ. 6) THEN
    489      nb_aer = 13
     442     nb_aer = 10
    490443     ALLOCATE (aerosol_name(nb_aer))
    491444     aerosol_name(1) = id_ASSO4M_phy     
    … …  
    499452     aerosol_name(9) = id_CIDUSTM_phy
    500453     aerosol_name(10)= id_CSSO4M_phy
    501      aerosol_name(11)= id_CSNO3M_phy
    502      aerosol_name(12)= id_ASNO3M_phy
    503      aerosol_name(13)= id_CINO3M_phy
    504454  ENDIF
    505455
    … …  
    541491        soluble=.TRUE.
    542492        spsol=2
    543      ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
     493     ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
    544494        soluble=.TRUE.
    545495        spsol=3
    546496        fac=1.375    ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD
    547      ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
     497     ELSEIF  (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
    548498        soluble=.TRUE.
    549499        spsol=4
    … …  
    558508         soluble=.TRUE.
    559509         spsol=7
    560      ELSEIF (aerosol_name(m).EQ.id_CSNO3M_phy) THEN
    561          soluble=.TRUE.
    562          spsol=8
    563          fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
    564      ELSEIF (aerosol_name(m).EQ.id_ASNO3M_phy) THEN
    565          soluble=.TRUE.
    566          spsol=9
    567          fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
    568510     ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN
    569511         soluble=.FALSE.
    … …  
    575517         soluble=.FALSE.
    576518         spinsol=3
    577      ELSEIF (aerosol_name(m).EQ.id_CINO3M_phy) THEN
    578          soluble=.FALSE.
    579          spinsol=4
    580          fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
    581519     ELSE
    582520         CYCLE
    … …  
    584522
    585523    !--shortname for aerosol index
    586     aerindex=aerosol_name(m)
     524    id=aerosol_name(m)
    587525
    588526    IF (soluble) THEN
    … …  
    594532           DO i=1, klon
    595533
    596              tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
    597              tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
     534             tmp_var=m_allaer(i,k,spsol)/1.e6*zdh(i,k)*fac
     535             tmp_var_pi=m_allaer_pi(i,k,spsol)/1.e6*zdh(i,k)*fac
    598536
    599537             ! Calculate the dry BC/(BC+SUL) mass ratio for all (natural+anthropogenic) aerosols
    … …  
    642580                             cg_MG_6bands(RH_num(i,k),inu,classbc))
    643581
    644                tau_ae(i,k,aerindex,inu)    = tmp_var*tau_ae2b_int
    645                piz_ae(i,k,aerindex,inu)    = piz_ae2b_int
    646                cg_ae(i,k,aerindex,inu)     = cg_ae2b_int
     582               tau_ae(i,k,id,inu)    = tmp_var*tau_ae2b_int
     583               piz_ae(i,k,id,inu)    = piz_ae2b_int
     584               cg_ae(i,k,id,inu)     = cg_ae2b_int
    647585
    648586               !--natural aerosols
    … …  
    659597                             cg_MG_6bands(RH_num(i,k),inu,classbc_pi))
    660598
    661                tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
    662                piz_ae_pi(i,k,aerindex,inu) = piz_ae2b_int
    663                cg_ae_pi(i,k,aerindex,inu)  = cg_ae2b_int
     599               tau_ae_pi(i,k,id,inu) = tmp_var_pi* tau_ae2b_int
     600               piz_ae_pi(i,k,id,inu) = piz_ae2b_int
     601               cg_ae_pi(i,k,id,inu)  = cg_ae2b_int
    664602                       
    665603             ENDDO
    … …  
    672610         DO k=1, klev
    673611           DO i=1, klon
    674              tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
    675              tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
     612             tmp_var=m_allaer(i,k,spsol)/1.e6*zdh(i,k)*fac
     613             tmp_var_pi=m_allaer_pi(i,k,spsol)/1.e6*zdh(i,k)*fac
    676614
    677615             DO inu=1,NSW
    … …  
    689627                             cg_aers_6bands(RH_num(i,k),inu,spsol))
    690628
    691                tau_ae(i,k,aerindex,inu)    = tmp_var*tau_ae2b_int
    692                tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi* tau_ae2b_int
    693                piz_ae(i,k,aerindex,inu)    = piz_ae2b_int
    694                cg_ae(i,k,aerindex,inu)     = cg_ae2b_int
     629               tau_ae(i,k,id,inu)    = tmp_var*tau_ae2b_int
     630               tau_ae_pi(i,k,id,inu) = tmp_var_pi* tau_ae2b_int
     631               piz_ae(i,k,id,inu)    = piz_ae2b_int
     632               cg_ae(i,k,id,inu)     = cg_ae2b_int
    695633                       
    696634             ENDDO
    … …  
    700638         !--external mixture case for soluble BC
    701639         IF (spsol.EQ.1) THEN
    702            piz_ae_pi(:,:,aerindex,:) = piz_ae(:,:,aerindex,:)
    703            cg_ae_pi(:,:,aerindex,:)  = cg_ae(:,:,aerindex,:)
     640           piz_ae_pi(:,:,id,:) = piz_ae(:,:,id,:)
     641           cg_ae_pi(:,:,id,:)  = cg_ae(:,:,id,:)
    704642         ENDIF
    705643
    … …  
    710648       DO k=1, klev
    711649         DO i=1, klon
    712            tmp_var=m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*fac
    713            tmp_var_pi=m_allaer_pi(i,k,aerindex)/1.e6*zdh(i,k)*fac
     650           tmp_var=m_allaer(i,k,naero_soluble+spinsol)/1.e6*zdh(i,k)*fac
     651           tmp_var_pi=m_allaer_pi(i,k,naero_soluble+spinsol)/1.e6*zdh(i,k)*fac
    714652
    715653           DO inu=1,NSW
    … …  
    718656             cg_ae2b_int = cg_aeri_6bands(inu,spinsol)
    719657
    720              tau_ae(i,k,aerindex,inu) = tmp_var*tau_ae2b_int
    721              tau_ae_pi(i,k,aerindex,inu) = tmp_var_pi*tau_ae2b_int
    722              piz_ae(i,k,aerindex,inu) = piz_ae2b_int
    723              cg_ae(i,k,aerindex,inu)= cg_ae2b_int
     658             tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int
     659             tau_ae_pi(i,k,id,inu) = tmp_var_pi*tau_ae2b_int
     660             piz_ae(i,k,id,inu) = piz_ae2b_int
     661             cg_ae(i,k,id,inu)= cg_ae2b_int
    724662           ENDDO
    725663         ENDDO
Note: See TracChangeset for help on using the changeset viewer.