Ignore:
Timestamp:
Jan 2, 2017, 11:24:30 PM (8 years ago)
Author:
oboucher
Message:

Swapping the order of CSSO4 and ASSO4 aerosols (and fixing an issue on aerindex). Preparing the ground for nitrate aerosols (coarse soluble, accumulation soluble, coarse insoluble). Modifying the LW aeropt routine so that it is compatible with both INCA and climatological aerosols (for dust only). Adding a new flag ok_alw for activating aerosol direct LW effect (for dust only). This change is bit comparable for flag_aerosol=6, flag_rrtm=1, NSW=6.

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r2737 r2738  
    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.
    86                                           !  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
     85  INTEGER, PARAMETER :: naero_soluble=9   !  1- BC soluble; 2- POM soluble; 3- SO4 coarse
     86                                          !  4- SO4 acc; 5 seasalt super-C; 6 seasalt coarse; 7 seasalt acc.
     87                                          !  8- NO3 coarse; 9- NO3 accumulation
     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.
    104104  REAL :: alpha_aers_5wv(nbre_RH,las,naero_soluble)   ! Ext. coeff. ** m2/g
    105    ! Insoluble components 1- Dust: 2- BC; 3- POM
    106105  REAL :: alpha_aeri_5wv(las,naero_insoluble)         ! Ext. coeff. ** m2/g
    107106
    … …  
    110109  !
    111110  INTEGER, PARAMETER ::  nbclassbc = 5  ! Added by Rong Wang/OB for the 5 fractions
    112                                        ! of BC in the soluble mode:
    113                                        ! bc_content/0.001, 0.01, 0.02, 0.05, ! 0.1/
     111                                        ! of BC in the soluble mode:
     112                                        ! bc_content/0.001, 0.01, 0.02, 0.05, 0.1/
    114113  ! for Maxwell-Garnet internal mixture
    115114  ! Detailed theory can be found in R. Wang Estimation of global black carbon ! direct
    … …  
    124123  INTEGER :: n, classbc
    125124
    126 ! From here on we look at the optical parameters at 5 wavelengths: 
    127 ! 443nm, 550, 670, 765 and 865 nm
    128 !                                   le 12 AVRIL 2006
    129 ! 
     125! From here on we look at the optical parameters at 5 wavelengths 443, 550, 670, 765 and 865 nm
     126
    130127 DATA alpha_aers_5wv/ &
    131128                                ! bc soluble
    … …  
    170167        2.340, 3.358, 3.803, 4.303, 4.928, 4.696, 5.629, 7.198,10.308,13.342,19.120,34.296, &
    171168        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 /
     169        1.359, 2.037, 2.343, 2.693, 3.139, 2.940, 3.596, 4.729, 7.076, 9.469,14.266,28.043, &
     170                           ! NO3 coarse (CS)
     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,  &
     176                           ! NO3 accumulation (AS)
     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 /
    173182
    174183  DATA alpha_aeri_5wv/ &
    175184                                 ! dust insoluble
    176         0.759, 0.770, 0.775, 0.775, 0.772, &
     185        0.759, 0.770, 0.775, 0.775, 0.772,  &
    177186                                 !!jb bc insoluble
    178187        11.536,10.033, 8.422, 7.234, 6.270, &
    179188                                 ! pom insoluble
    180         5.042, 3.101, 1.890, 1.294, 0.934/
     189        5.042, 3.101, 1.890, 1.294, 0.934,  &
     190                                 ! NO3 insoluble
     191        0.0, 0.0, 0.0, 0.0, 0.0 /
    181192
    182193  ! internal mixture
    … …  
    252263     aerosol_name(1) = id_CIDUSTM_phy
    253264  ELSEIF (flag_aerosol .EQ. 6) THEN
    254      nb_aer = 10
     265     nb_aer = 13
    255266     ALLOCATE (aerosol_name(nb_aer))
    256267     aerosol_name(1) = id_ASSO4M_phy     
    … …  
    264275     aerosol_name(9) = id_CIDUSTM_phy
    265276     aerosol_name(10)= id_CSSO4M_phy
     277     aerosol_name(11)= id_CSNO3M_phy
     278     aerosol_name(12)= id_ASNO3M_phy
     279     aerosol_name(13)= id_CINO3M_phy
    266280  ENDIF
    267281
    … …  
    315329        soluble=.TRUE.
    316330        spsol=7
     331    ELSEIF (aerosol_name(m).EQ.id_CSNO3M_phy) THEN
     332        soluble=.TRUE.
     333        spsol=8
     334        fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
     335    ELSEIF (aerosol_name(m).EQ.id_ASNO3M_phy) THEN
     336        soluble=.TRUE.
     337        spsol=9
     338        fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
    317339    ELSEIF (aerosol_name(m).EQ.id_CIDUSTM_phy) THEN
    318340        soluble=.FALSE.
    … …  
    324346        soluble=.FALSE.
    325347        spinsol=3
     348    ELSEIF (aerosol_name(m).EQ.id_CINO3M_phy) THEN
     349        soluble=.FALSE.
     350        spinsol=4
     351        fac=1.2903    ! NO3NH4/NO3 / mass conversion factor for AOD
    326352    ELSE
    327353        CYCLE
Note: See TracChangeset for help on using the changeset viewer.