Changeset 2735


Ignore:
Timestamp:
Dec 17, 2016, 3:45:50 PM (7 years ago)
Author:
oboucher
Message:

Swapping the order of ASSO4 and CSSO4.
Preparing the ground for nitrate aerosols.
ALl of this should maintain bit-comparability.

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

Legend:

Unmodified
Added
Removed

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

    r2644 r2735  
    8383  INTEGER, PARAMETER :: la865 = 5
    8484  INTEGER, PARAMETER :: nbre_RH=12
    85   INTEGER, PARAMETER :: naero_soluble=7   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.
     85  INTEGER, PARAMETER :: naero_soluble=9   !  1- BC soluble; 2- POM soluble; 3- SO4 acc.
    8686                                          !  4- SO4 coarse; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
    87   INTEGER, PARAMETER :: naero_insoluble=3 !  1- Dust; 2- BC insoluble; 3- POM insoluble
     87                                          !  8- NO3 coarse;  9- NO3 acc.
     88  INTEGER, PARAMETER :: naero_insoluble=4 !  1- Dust; 2- BC insoluble; 3- POM insoluble; 4- NO3 insoluble
    8889
    8990  REAL :: zrho
    … …  
    101102  REAL :: zdh(klon,klev)
    102103 
    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.
     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.
    104107  REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble)   ! Ext. coeff. ** m2/g
    105    ! Insoluble components 1- Dust: 2- BC; 3- POM
     108   ! Insoluble components 1- Dust: 2- BC; 3- POM; 4- NO3
    106109  REAL :: alpha_aeri_5wv(las,naero_insoluble)         ! Ext. coeff. ** m2/g
    107110
    … …  
    129132! 
    130133 DATA alpha_aers_5wv/ &
    131                                 ! bc soluble
     134                           ! BC soluble
    132135       7.930,7.930,7.930,7.930,7.930,7.930,7.930,7.930,10.893,12.618,14.550,16.613, &
    133136       7.658,7.658,7.658,7.658,7.658,7.658,7.658,7.658,10.351,11.879,13.642,15.510, &
    … …  
    135138       6.736,6.736,6.736,6.736,6.736,6.736,6.736,6.736,8.818,9.938,11.283,12.687,   &
    136139       6.277,6.277,6.277,6.277,6.277,6.277,6.277,6.277,8.123,9.094,10.275,11.501,   &
    137                                 ! pom soluble
     140                           ! POM soluble
    138141       6.676,6.676,6.676,6.676,6.710,6.934,7.141,7.569,8.034,8.529,9.456,10.511,  &
    139142       5.109,5.109,5.109,5.109,5.189,5.535,5.960,6.852,8.008,9.712,12.897,19.676, &
    … …  
    141144       2.849,2.849,2.849,2.849,2.897,3.107,3.365,3.916,4.649,5.760,7.900,12.863,  &
    142145       2.229,2.229,2.229,2.229,2.268,2.437,2.645,3.095,3.692,4.608,6.391,10.633,  &
    143                                 ! Sulfate (Accumulation)
     146                           ! Sulfate (Coarse)
    144147       5.751,6.215,6.690,7.024,7.599,8.195,9.156,10.355,12.660,14.823,18.908,24.508, &
    145148       4.320,4.675,5.052,5.375,5.787,6.274,7.066,8.083,10.088,12.003,15.697,21.133,  &
    … …  
    147150       2.336,2.552,2.781,2.979,3.236,3.540,4.046,4.711,6.056,7.388,10.093,14.313,    &
    148151       1.777,1.949,2.134,2.292,2.503,2.751,3.166,3.712,4.828,5.949,8.264,11.922,     &
    149                                 ! Sulfate (Coarse)
     152                           ! Sulfate (Accumulation)
    150153       5.751,6.215,6.690,7.024,7.599,8.195,9.156,10.355,12.660,14.823,18.908,24.508, &
    151154       4.320,4.675,5.052,5.375,5.787,6.274,7.066,8.083,10.088,12.003,15.697,21.133,  &
    … …  
    153156       2.336,2.552,2.781,2.979,3.236,3.540,4.046,4.711,6.056,7.388,10.093,14.313,    &
    154157       1.777,1.949,2.134,2.292,2.503,2.751,3.166,3.712,4.828,5.949,8.264,11.922,     &
    155                            ! seasalt seasalt Super Coarse Soluble (SS)
     158                           ! Seasalt Super Coarse Soluble (SS)
    156159        0.218, 0.272, 0.293, 0.316, 0.343, 0.380, 0.429, 0.501, 0.636, 0.755, 0.967, 1.495, &
    157160        0.221, 0.275, 0.297, 0.320, 0.348, 0.383, 0.432, 0.509, 0.640, 0.759, 0.972, 1.510, &
    … …  
    159162        0.227, 0.282, 0.303, 0.327, 0.356, 0.392, 0.441, 0.518, 0.652, 0.770, 0.987, 1.529, &
    160163        0.230, 0.285, 0.306, 0.330, 0.359, 0.396, 0.446, 0.522, 0.656, 0.777, 0.993, 1.539, &
    161                            ! seasalt seasalt Coarse Soluble (CS)     
     164                           ! Seasalt Coarse Soluble (CS)     
    162165        0.578, 0.706, 0.756, 0.809, 0.876, 0.964, 1.081, 1.256, 1.577, 1.858, 2.366, 3.613, &
    163166        0.598, 0.725, 0.779, 0.833, 0.898, 0.990, 1.109, 1.290, 1.609, 1.889, 2.398, 3.682, &
    … …  
    165168        0.633, 0.767, 0.820, 0.879, 0.948, 1.044, 1.167, 1.353, 1.683, 1.969, 2.491, 3.785, &
    166169        0.648, 0.785, 0.838, 0.896, 0.967, 1.066, 1.192, 1.381, 1.714, 2.006, 2.531, 3.836, &
    167                            ! seasalt seasalt Accumulation Soluble (AS)
     170                           ! Seasalt Accumulation Soluble (AS)
    168171        4.432, 5.899, 6.505, 7.166, 7.964, 7.962, 9.232,11.257,14.979,18.337,24.223,37.811, &
    169172        3.298, 4.569, 5.110, 5.709, 6.446, 6.268, 7.396, 9.246,12.787,16.113,22.197,37.136, &
    170173        2.340, 3.358, 3.803, 4.303, 4.928, 4.696, 5.629, 7.198,10.308,13.342,19.120,34.296, &
    171174        1.789, 2.626, 2.999, 3.422, 3.955, 3.730, 4.519, 5.864, 8.593,11.319,16.653,31.331, &
    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 /
     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 /
    173188
    174189  DATA alpha_aeri_5wv/ &
    175190                                 ! dust insoluble
    176         0.759, 0.770, 0.775, 0.775, 0.772, &
    177                                  !!jb bc insoluble
     191        0.759, 0.770, 0.775, 0.775, 0.772,  &
     192                                 ! BC insoluble
    178193        11.536,10.033, 8.422, 7.234, 6.270, &
    179                                  ! pom insoluble
    180         5.042, 3.101, 1.890, 1.294, 0.934/
     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 /
    181198
    182199  ! internal mixture
    … …  
    264281     aerosol_name(9) = id_CIDUSTM_phy
    265282     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
    266286  ENDIF
    267287
    … …  
    298318        soluble=.TRUE.
    299319        spsol=2
    300     ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
     320    ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
    301321        soluble=.TRUE.
    302322        spsol=3
    303323        fac=1.375    ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD
    304     ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
     324    ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
    305325        soluble=.TRUE.
    306326        spsol=4
    … …  
    315335        soluble=.TRUE.
    316336        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
    317345    ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN
    318346        soluble=.FALSE.
    … …  
    324352        soluble=.FALSE.
    325353        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
    326358    ELSE
    327359        CYCLE

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

    r2644 r2735  
    33!
    44SUBROUTINE AEROPT_6BANDS_RRTM ( &
    5      pdel, m_allaer, RHcl, &
    6      tau_allaer, piz_allaer, &
    7      cg_allaer, m_allaer_pi, &
    8      flag_aerosol, flag_bc_internal_mixture, zrho )
     5     pdel, m_allaer, RHcl, tau_allaer, piz_allaer, &
     6     cg_allaer, m_allaer_pi, flag_aerosol, flag_bc_internal_mixture, zrho )
    97
    108  USE dimphy
    … …  
    5351  INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA
    5452
    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
     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
    5857
    5958  INTEGER, PARAMETER :: nbre_RH=12
    … …  
    6766  REAL :: fac
    6867  REAL :: zdh(klon,klev)
    69   INTEGER, ALLOCATABLE, DIMENSION(:)   :: aerosol_name
     68  INTEGER, ALLOCATABLE, DIMENSION(:) :: aerosol_name
    7069  INTEGER :: nb_aer
    7170
    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
     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
    7776  REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: piz_ae_pi
    7877  REAL, DIMENSION(klon,klev,id_ASBCM_phy:id_ASBCM_phy,nbands_sw_rrtm) :: cg_ae_pi
    … …  
    8079  ! Proprietes optiques
    8180  !
    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
     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
    8887  !
    8988  ! BC internal mixture
    9089  !
    9190  INTEGER, PARAMETER ::  nbclassbc = 5  ! Added by Rong Wang/OB for the 5 fractions
    92                                        ! of BC in the soluble mode:
    93                                        ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/
     91                                        ! of BC in the soluble mode:
     92                                        ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/
    9493  ! for Maxwell-Garnet internal mixture
    9594  ! Detailed theory can be found in R. Wang Estimation of global black carbon direct
    … …  
    126125  0.377, 0.377, 0.377, 0.377, 0.377, 0.415, 0.456, 0.542, 0.627, 0.636, 0.646, 0.738, &
    127126  0.052, 0.052, 0.052, 0.052, 0.052, 0.057, 0.063, 0.075, 0.087, 0.088, 0.089, 0.102, &
     127       ! sulfate coarse CS
     128  0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
     129  0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
     130  0.715, 0.715, 0.715, 0.788, 0.865, 0.945, 1.029, 1.171, 1.557, 1.936, 1.936, 3.459, &
     131  0.736, 0.736, 0.736, 0.811, 0.891, 0.973, 1.060, 1.206, 1.603, 1.994, 1.994, 3.563, &
     132  0.711, 0.711, 0.711, 0.783, 0.860, 0.939, 1.023, 1.164, 1.548, 1.925, 1.925, 3.439, &
     133  0.602, 0.602, 0.602, 0.664, 0.729, 0.796, 0.867, 0.986, 1.312, 1.631, 1.631, 2.915, &
    128134       ! sulfate AS   
    129135  6.554, 6.554, 6.554, 7.223, 7.931, 8.665, 9.438,10.736,14.275,17.755,17.755,31.722, &
    … …  
    133139  0.312, 0.312, 0.312, 0.344, 0.378, 0.413, 0.450, 0.511, 0.680, 0.846, 0.846, 1.511, &
    134140  0.121, 0.121, 0.121, 0.134, 0.147, 0.161, 0.175, 0.199, 0.264, 0.329, 0.329, 0.588, &
    135        ! sulfate coarse CS
    136   0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
    137   0.693, 0.693, 0.693, 0.764, 0.839, 0.917, 0.999, 1.136, 1.510, 1.879, 1.879, 3.356, &
    138   0.715, 0.715, 0.715, 0.788, 0.865, 0.945, 1.029, 1.171, 1.557, 1.936, 1.936, 3.459, &
    139   0.736, 0.736, 0.736, 0.811, 0.891, 0.973, 1.060, 1.206, 1.603, 1.994, 1.994, 3.563, &
    140   0.711, 0.711, 0.711, 0.783, 0.860, 0.939, 1.023, 1.164, 1.548, 1.925, 1.925, 3.439, &
    141   0.602, 0.602, 0.602, 0.664, 0.729, 0.796, 0.867, 0.986, 1.312, 1.631, 1.631, 2.915, &
    142         ! seasalt seasalt Super Coarse Soluble (SS)
     141        ! Seasalt Super Coarse Soluble (SS)
    143142  0.214, 0.267, 0.287, 0.310, 0.337, 0.373, 0.421, 0.494, 0.625, 0.742, 0.956, 1.480, &
    144143  0.217, 0.270, 0.291, 0.314, 0.341, 0.377, 0.426, 0.499, 0.632, 0.751, 0.963, 1.490, &
    … …  
    147146  0.250, 0.307, 0.330, 0.354, 0.384, 0.424, 0.477, 0.556, 0.696, 0.822, 1.044, 1.592, &
    148147  0.279, 0.347, 0.373, 0.401, 0.434, 0.478, 0.537, 0.625, 0.781, 0.918, 1.158, 1.744, &
    149         ! seasalt seasalt Coarse Soluble (CS)     
     148        ! Seasalt Coarse Soluble (CS)     
    150149  0.550, 0.676, 0.724, 0.779, 0.841, 0.925, 1.040, 1.213, 1.523, 1.803, 2.306, 3.535, &
    151150  0.568, 0.695, 0.744, 0.798, 0.863, 0.950, 1.066, 1.240, 1.556, 1.839, 2.342, 3.588, &
    … …  
    154153  0.663, 0.825, 0.889, 0.957, 1.038, 1.126, 1.268, 1.483, 1.862, 2.196, 2.780, 4.183, &
    155154  0.479, 0.644, 0.711, 0.785, 0.873, 0.904, 1.048, 1.275, 1.692, 2.072, 2.753, 4.430, &
    156         ! seasalt seasalt Accumulation Soluble (AS)
     155        ! Seasalt Accumulation Soluble (AS)
    157156  6.128, 7.564, 8.127, 8.723, 9.421,10.014,11.250,13.113,16.320,19.061,23.675,34.158, &
    158157  5.080, 6.586, 7.197, 7.856, 8.645, 8.820,10.115,12.140,15.775,18.983,24.499,36.993, &
    … …  
    160159  1.389, 2.066, 2.371, 2.719, 3.161, 2.970, 3.623, 4.747, 7.064, 9.419,14.120,27.565, &
    161160  0.309, 0.497, 0.587, 0.693, 0.832, 0.777, 0.986, 1.364, 2.207, 3.136, 5.174,12.065, &
    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  /
     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 /
    163176
    164177  DATA alpha_aeri_6bands/  &
    165178       ! dust insoluble CI
    166179  0.751, 0.751, 0.769, 0.772, 0.672, 0.437, &
    167        ! bc insoluble AI
     180       ! BC insoluble AI
    168181  6.497, 6.497, 5.900, 4.284, 2.163, 0.966, &
    169        ! pom insoluble AI
    170   6.443, 6.443, 4.381, 1.846, 0.377, 0.052 /
     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 /
    171186
    172187  DATA cg_aers_6bands/ &
    … …  
    185200  0.363, 0.363, 0.363, 0.363, 0.363, 0.360, 0.357, 0.352, 0.350, 0.349, 0.349, 0.347, &
    186201  0.142, 0.142, 0.142, 0.142, 0.142, 0.139, 0.137, 0.133, 0.131, 0.131, 0.130, 0.129, &
     202       ! sulfate coarse CS
     203  0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
     204  0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
     205  0.799, 0.799, 0.799, 0.787, 0.777, 0.768, 0.760, 0.747, 0.736, 0.729, 0.725, 0.716, &
     206  0.797, 0.797, 0.797, 0.782, 0.770, 0.760, 0.750, 0.735, 0.722, 0.714, 0.709, 0.698, &
     207  0.810, 0.810, 0.810, 0.794, 0.781, 0.770, 0.759, 0.743, 0.728, 0.719, 0.714, 0.702, &
     208  0.803, 0.803, 0.803, 0.790, 0.779, 0.770, 0.762, 0.748, 0.736, 0.729, 0.725, 0.715, &
    187209       ! sulfate AS
    188210  0.675, 0.675, 0.675, 0.689, 0.701, 0.711, 0.720, 0.735, 0.748, 0.756, 0.760, 0.771, &
    … …  
    192214  0.362, 0.362, 0.362, 0.359, 0.356, 0.354, 0.352, 0.348, 0.345, 0.343, 0.342, 0.340, &
    193215  0.137, 0.137, 0.137, 0.135, 0.133, 0.132, 0.130, 0.128, 0.126, 0.125, 0.124, 0.122, &
    194        ! sulfate coarse CS
    195   0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
    196   0.803, 0.803, 0.803, 0.792, 0.783, 0.776, 0.769, 0.758, 0.747, 0.742, 0.738, 0.730, &
    197   0.799, 0.799, 0.799, 0.787, 0.777, 0.768, 0.760, 0.747, 0.736, 0.729, 0.725, 0.716, &
    198   0.797, 0.797, 0.797, 0.782, 0.770, 0.760, 0.750, 0.735, 0.722, 0.714, 0.709, 0.698, &
    199   0.810, 0.810, 0.810, 0.794, 0.781, 0.770, 0.759, 0.743, 0.728, 0.719, 0.714, 0.702, &
    200   0.803, 0.803, 0.803, 0.790, 0.779, 0.770, 0.762, 0.748, 0.736, 0.729, 0.725, 0.715, &
    201216        ! seasalt seasalt Super Coarse Soluble (SS)
    202217  0.797, 0.800, 0.801, 0.802, 0.804, 0.822, 0.825, 0.828, 0.832, 0.835, 0.838, 0.843, &
    … …  
    219234  0.561, 0.585, 0.593, 0.601, 0.609, 0.637, 0.651, 0.669, 0.691, 0.705, 0.723, 0.745, &
    220235  0.392, 0.427, 0.439, 0.451, 0.464, 0.480, 0.500, 0.526, 0.563, 0.588, 0.621, 0.671, &
    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 /
     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 /
    222251
    223252  DATA cg_aeri_6bands/ &
    224253       ! dust insoluble CI
    225254  0.718, 0.718, 0.699, 0.661, 0.676, 0.670, &
    226        ! bc insoluble AI
     255       ! BC insoluble AI
    227256  0.721, 0.721, 0.643, 0.513, 0.321, 0.153, &
    228        ! pom insoluble AI
    229   0.687, 0.687, 0.658, 0.564, 0.363, 0.142  /
     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 /
    230261
    231262  DATA piz_aers_6bands/&
    232        ! bc soluble AS
     263       ! BC soluble AS
    233264  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
    234265  0.460, 0.460, 0.460, 0.460, 0.460, 0.534, 0.594, 0.688, 0.748, 0.754, 0.760, 0.803, &
    … …  
    244275  0.816, 0.816, 0.816, 0.816, 0.816, 0.839, 0.859, 0.888, 0.908, 0.910, 0.911, 0.925, &
    245276  0.663, 0.663, 0.663, 0.663, 0.663, 0.607, 0.562, 0.492, 0.446, 0.441, 0.437, 0.404, &
     277       ! sulfate coarse CS
     278  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
     279  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
     280  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
     281  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
     282  0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, &
     283  0.877, 0.877, 0.877, 0.873, 0.870, 0.867, 0.864, 0.860, 0.856, 0.854, 0.852, 0.849, &
    246284       ! sulfate AS
    247285  1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    … …  
    251289  0.988, 0.988, 0.988, 0.989, 0.990, 0.990, 0.991, 0.992, 0.993, 0.993, 0.994, 0.994, &
    252290  0.256, 0.256, 0.256, 0.263, 0.268, 0.273, 0.277, 0.284, 0.290, 0.294, 0.296, 0.301, &
    253        ! sulfate coarse CS
    254   1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    255   1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    256   1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    257   1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    258   0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, 0.996, &
    259   0.877, 0.877, 0.877, 0.873, 0.870, 0.867, 0.864, 0.860, 0.856, 0.854, 0.852, 0.849, &
    260291        ! seasalt seasalt Super Coarse Soluble (SS)
    261292  0.983, 0.982, 0.982, 0.982, 0.982, 0.992, 0.993, 0.994, 0.995, 0.996, 0.997, 0.998, &
    … …  
    278309  0.999, 0.999, 0.999, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, &
    279310  0.997, 0.998, 0.998, 0.998, 0.998, 0.999, 0.999, 0.999, 1.000, 1.000, 1.000, 1.000, &
    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 /
     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 /
    281326
    282327  DATA piz_aeri_6bands/ &
    283328       ! dust insoluble CI
    284329  0.954, 0.954, 0.965, 0.981, 0.996, 0.990, &
    285        ! bc insoluble AI
     330       ! BC insoluble AI
    286331  0.460, 0.460, 0.445, 0.394, 0.267, 0.121, &
    287        ! pom insoluble AI
    288   0.973, 0.973, 0.972, 0.940, 0.816, 0.663 /
     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 /
    289336
    290337! Added by R. Wang (July 31 2016)
    … …  
    440487     aerosol_name(1) = id_CIDUSTM_phy
    441488  ELSEIF (flag_aerosol .EQ. 6) THEN
    442      nb_aer = 10
     489     nb_aer = 13
    443490     ALLOCATE (aerosol_name(nb_aer))
    444491     aerosol_name(1) = id_ASSO4M_phy     
    … …  
    452499     aerosol_name(9) = id_CIDUSTM_phy
    453500     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
    454504  ENDIF
    455505
    … …  
    491541        soluble=.TRUE.
    492542        spsol=2
    493      ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
     543     ELSEIF (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
    494544        soluble=.TRUE.
    495545        spsol=3
    496546        fac=1.375    ! (NH4)2-SO4/SO4 132/96 mass conversion factor for AOD
    497      ELSEIF  (aerosol_name(m).EQ.id_CSSO4M_phy) THEN
     547     ELSEIF (aerosol_name(m).EQ.id_ASSO4M_phy) THEN
    498548        soluble=.TRUE.
    499549        spsol=4
    … …  
    508558         soluble=.TRUE.
    509559         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
    510568     ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN
    511569         soluble=.FALSE.
    … …  
    517575         soluble=.FALSE.
    518576         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
    519581     ELSE
    520582         CYCLE
Note: See TracChangeset for help on using the changeset viewer.