Changeset 2394 for LMDZ5


Ignore:
Timestamp:
Nov 18, 2015, 12:41:49 PM (9 years ago)
Author:
acozic
Message:

to fit with inca we need to increase dimension of naero_grp in aero_mod.
When lmdz is running without inca dim 10 to 12 will be empty and there is no change in results

these modifications only impact the old physics coupled with inca

Location:
LMDZ5/trunk/libf/phylmd
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/aero_mod.F90

    r2146 r2394  
    6262
    6363! 3/ Number of aerosol groups
    64   INTEGER, PARAMETER :: naero_grp = 9
     64  INTEGER, PARAMETER :: naero_grp = 12
     65  ! if info_trac = inca
     66  ! 1 = ZERO   
     67  ! 2 = AER total   
     68  ! 3 = NAT   
     69  ! 4 = BC   
     70  ! 5 = SO4   
     71  ! 6 = POM   
     72  ! 7 = DUST   
     73  ! 8 = SS   
     74  ! 9 = FNO3   
     75  ! 10 = DNO3
     76  ! 11 = SNO3
     77  ! 12 = empty actually
     78  ! else
    6579  ! 1 = ZERO   
    6680  ! 2 = AER total   
  • LMDZ5/trunk/libf/phylmd/phytrac_mod.F90

    r2351 r2394  
    100100    USE mod_phys_lmdz_mpi_data, ONLY :  is_mpi_root
    101101    USE print_control_mod, ONLY: lunout
     102    USE aero_mod, ONLY : naero_grp
    102103
    103104    IMPLICIT NONE
     
    168169    LOGICAL,INTENT(IN)                       :: aerosol_couple
    169170    REAL,DIMENSION(klon,klev),INTENT(IN)     :: flxmass_w
    170     REAL,DIMENSION(klon,klev,9,2),INTENT(IN) :: tau_aero
    171     REAL,DIMENSION(klon,klev,9,2),INTENT(IN) :: piz_aero
    172     REAL,DIMENSION(klon,klev,9,2),INTENT(IN) :: cg_aero
    173     CHARACTER(len=4),DIMENSION(9),INTENT(IN) :: rfname
     171    REAL,DIMENSION(klon,klev,naero_grp,2),INTENT(IN) :: tau_aero
     172    REAL,DIMENSION(klon,klev,naero_grp,2),INTENT(IN) :: piz_aero
     173    REAL,DIMENSION(klon,klev,naero_grp,2),INTENT(IN) :: cg_aero
     174    CHARACTER(len=4),DIMENSION(naero_grp),INTENT(IN) :: rfname
    174175    REAL,DIMENSION(klon,klev,2),INTENT(IN)   :: ccm
    175176    !... K.Emanuel
  • LMDZ5/trunk/libf/phylmd/radlwsw_m.F90

    r2366 r2394  
    7777      USE YOERRTWN , ONLY : DELWAVE   ,TOTPLNK     
    7878      USE YOMPHY3  , ONLY : RII0
    79 #else
    80       USE aero_mod, ONLY : nbands_lw_rrtm
    8179#endif
     80      USE aero_mod
    8281
    8382  !======================================================================
     
    193192  LOGICAL, INTENT(in)  :: flag_aerosol_strat                             ! use stratospheric aerosols
    194193  REAL,    INTENT(in)  :: cldfra(KLON,KLEV), cldemi(KLON,KLEV), cldtaupd(KLON,KLEV)
    195   REAL,    INTENT(in)  :: tau_aero(KLON,KLEV,9,2)                        ! aerosol optical properties (see aeropt.F)
    196   REAL,    INTENT(in)  :: piz_aero(KLON,KLEV,9,2)                        ! aerosol optical properties (see aeropt.F)
    197   REAL,    INTENT(in)  :: cg_aero(KLON,KLEV,9,2)                         ! aerosol optical properties (see aeropt.F)
     194  REAL,    INTENT(in)  :: tau_aero(KLON,KLEV,naero_grp,2)                        ! aerosol optical properties (see aeropt.F)
     195  REAL,    INTENT(in)  :: piz_aero(KLON,KLEV,naero_grp,2)                        ! aerosol optical properties (see aeropt.F)
     196  REAL,    INTENT(in)  :: cg_aero(KLON,KLEV,naero_grp,2)                         ! aerosol optical properties (see aeropt.F)
    198197!--OB
    199198  REAL,    INTENT(in)  :: tau_aero_sw_rrtm(KLON,KLEV,2,NSW)                 ! aerosol optical properties RRTM
     
    291290  REAL(KIND=8) zsolsw0(kdlon), zsollw0(kdlon)
    292291  REAL(KIND=8) zznormcp
    293   REAL(KIND=8) tauaero(kdlon,kflev,9,2)                     ! aer opt properties
    294   REAL(KIND=8) pizaero(kdlon,kflev,9,2)
    295   REAL(KIND=8) cgaero(kdlon,kflev,9,2)
     292  REAL(KIND=8) tauaero(kdlon,kflev,naero_grp,2)                     ! aer opt properties
     293  REAL(KIND=8) pizaero(kdlon,kflev,naero_grp,2)
     294  REAL(KIND=8) cgaero(kdlon,kflev,naero_grp,2)
    296295  REAL(KIND=8) PTAUA(kdlon,2,kflev)                         ! present-day value of cloud opt thickness (PTAU is pre-industrial value), local use
    297296  REAL(KIND=8) POMEGAA(kdlon,2,kflev)                       ! dito for single scatt albedo
  • LMDZ5/trunk/libf/phylmd/sw_aeroAR4.F90

    r2311 r2394  
    2323  USE phys_output_mod, ONLY : swaero_diag
    2424  USE print_control_mod, ONLY: lunout
     25  USE aero_mod, ONLY : naero_grp
    2526  IMPLICIT NONE
    2627
     
    140141  LOGICAL flag_aerosol_strat ! use stratospehric aerosols
    141142  INTEGER flag_aerosol      ! global flag for aerosol 0 (no aerosol) or 1-5 (aerosols)
    142   REAL(KIND=8) tauaero(kdlon,kflev,9,2)  ! aerosol optical properties
    143   REAL(KIND=8) pizaero(kdlon,kflev,9,2)  ! (see aeropt.F)
    144   REAL(KIND=8) cgaero(kdlon,kflev,9,2)   ! -"-
     143  REAL(KIND=8) tauaero(kdlon,kflev,naero_grp,2)  ! aerosol optical properties
     144  REAL(KIND=8) pizaero(kdlon,kflev,naero_grp,2)  ! (see aeropt.F)
     145  REAL(KIND=8) cgaero(kdlon,kflev,naero_grp,2)   ! -"-
    145146  REAL(KIND=8) PTAUA(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS (present-day value)
    146147  REAL(KIND=8) POMEGAA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
  • LMDZ5/trunk/libf/phylmd/tracinca_mod.F90

    r2372 r2394  
    5454    USE geometry_mod, ONLY: cell_area
    5555    USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat
     56    USE aero_mod, ONLY : naero_grp
    5657    IMPLICIT NONE
    5758   
     
    99100    LOGICAL,INTENT(IN)                       :: aerosol_couple
    100101    REAL,DIMENSION(klon,klev),INTENT(IN)     :: flxmass_w
    101     REAL,DIMENSION(klon,klev,9,2),INTENT(IN) :: tau_aero
    102     REAL,DIMENSION(klon,klev,9,2),INTENT(IN) :: piz_aero
    103     REAL,DIMENSION(klon,klev,9,2),INTENT(IN) :: cg_aero
    104     CHARACTER(len=4),DIMENSION(9),INTENT(IN) :: rfname
     102    REAL,DIMENSION(klon,klev,naero_grp,2),INTENT(IN) :: tau_aero
     103    REAL,DIMENSION(klon,klev,naero_grp,2),INTENT(IN) :: piz_aero
     104    REAL,DIMENSION(klon,klev,naero_grp,2),INTENT(IN) :: cg_aero
     105    CHARACTER(len=4),DIMENSION(naero_grp),INTENT(IN) :: rfname
    105106    REAL,DIMENSION(klon,klev,2),INTENT(IN)   :: ccm
    106107
Note: See TracChangeset for help on using the changeset viewer.