Changeset 3412 for LMDZ6/trunk


Ignore:
Timestamp:
Nov 8, 2018, 10:39:16 PM (6 years ago)
Author:
oboucher
Message:

Externalizing a controlling parameter for activating aerosol radiative feedback
Need to pass the info all the way to recmwf_aero. Default is True as before.

Location:
LMDZ6/trunk/libf
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90

    r3338 r3412  
    113113  INTEGER :: flag_aerosol
    114114  INTEGER :: flag_aerosol_strat
     115  LOGICAL :: flag_aer_feedback
    115116  LOGICAL :: flag_bc_internal_mixture
    116117  LOGICAL :: new_aod
    … …  
    133134                   iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs,            &
    134135                   ok_ade, ok_aie, ok_alw, ok_cdnc, aerosol_couple,     &
    135                    chemistry_couple, flag_aerosol, flag_aerosol_strat,  &
     136                   chemistry_couple, flag_aerosol, flag_aerosol_strat,  &
     137                   flag_aer_feedback,                                   &
    136138                   new_aod, flag_bc_internal_mixture, bl95_b0, bl95_b1, &
    137139                   read_climoz, alp_offset)

  • LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

    r3384 r3412  
    1818       iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
    1919       ok_ade, ok_aie, ok_alw, ok_cdnc, aerosol_couple, chemistry_couple, &
    20        flag_aerosol, flag_aerosol_strat, new_aod, &
     20       flag_aerosol, flag_aerosol_strat, flag_aer_feedback, new_aod, &
    2121       flag_bc_internal_mixture, bl95_b0, bl95_b1,&
    2222       read_climoz, &
    … …  
    7474    INTEGER              :: flag_aerosol
    7575    INTEGER              :: flag_aerosol_strat
     76    LOGICAL              :: flag_aer_feedback
    7677    LOGICAL              :: flag_bc_internal_mixture
    7778    LOGICAL              :: new_aod
    … …  
    9596    INTEGER, SAVE       :: flag_aerosol_omp
    9697    INTEGER, SAVE       :: flag_aerosol_strat_omp
     98    LOGICAL, SAVE       :: flag_aer_feedback_omp
    9799    LOGICAL, SAVE       :: flag_bc_internal_mixture_omp
    98100    LOGICAL, SAVE       :: new_aod_omp
    … …  
    10221024      CALL getin('ok_suntime_rrtm',ok_suntime_rrtm_omp)
    10231025    ENDIF
    1024     !
     1026   
    10251027    !Config Key  = flag_aerosol_strat
    10261028    !Config Desc = use stratospheric aerosols 0, 1, 2
    … …  
    10321034    !Config Help = Used in physiq.F
    10331035    !
    1034     !
    10351036    flag_aerosol_strat_omp = 0
    10361037    CALL getin('flag_aerosol_strat',flag_aerosol_strat_omp)
    10371038
    1038     !
     1039    !Config Key  = flag_aer_feedback
     1040    !Config Desc = (des)activate aerosol radiative feedback
     1041    ! - F = no aerosol radiative feedback
     1042    ! - T = aerosol radiative feedback
     1043    !Config Def  = T
     1044    !Config Help = Used in physiq.F
     1045    !
     1046    flag_aer_feedback_omp = .TRUE.
     1047    IF (iflag_rrtm_omp==1) THEN
     1048       CALL getin('flag_aer_feedback',flag_aer_feedback_omp)
     1049    ENDIF
     1050
    10391051    !Config Key  = iflag_cld_th
    10401052    !Config Desc = 
    … …  
    22972309    flag_aerosol=flag_aerosol_omp
    22982310    flag_aerosol_strat=flag_aerosol_strat_omp
     2311    flag_aer_feedback=flag_aer_feedback_omp
    22992312    flag_bc_internal_mixture=flag_bc_internal_mixture_omp
    23002313    new_aod=new_aod_omp
    … …  
    26352648    WRITE(lunout,*) ' flag_aerosol = ', flag_aerosol
    26362649    WRITE(lunout,*) ' flag_aerosol_strat= ', flag_aerosol_strat
     2650    WRITE(lunout,*) ' flag_aer_feedback= ', flag_aer_feedback
    26372651    WRITE(lunout,*) ' new_aod = ', new_aod
    26382652    WRITE(lunout,*) ' aer_type = ',aer_type

  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r3387 r3412  
    10611061    INTEGER, SAVE :: flag_aerosol_strat
    10621062    !$OMP THREADPRIVATE(flag_aerosol_strat)
     1063    !
     1064    !--INTERACTIVE AEROSOL FEEDBACK ON RADIATION
     1065    LOGICAL, SAVE :: flag_aer_feedback
     1066    !$OMP THREADPRIVATE(flag_aer_feedback)
     1067
    10631068    !c-fin STRAT AEROSOL
    10641069    !
    … …  
    12171222            iflag_cld_th,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
    12181223            ok_ade, ok_aie, ok_alw, ok_cdnc, aerosol_couple,  chemistry_couple, &
    1219             flag_aerosol, flag_aerosol_strat, new_aod, &
     1224            flag_aerosol, flag_aerosol_strat, flag_aer_feedback, new_aod, &
    12201225            flag_bc_internal_mixture, bl95_b0, bl95_b1, &
    12211226                                ! nv flags pour la convection et les
    … …  
    38823887               cldfrarad, cldemirad, cldtaurad, &
    38833888               ok_ade.OR.flag_aerosol_strat.GT.0, ok_aie, flag_aerosol, &
    3884                flag_aerosol_strat, &
     3889               flag_aerosol_strat, flag_aer_feedback, &
    38853890               tau_aero, piz_aero, cg_aero, &
    38863891               tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm, &
    … …  
    39673972                     cldfrarad, cldemirad, cldtaurad, &
    39683973                     ok_ade.OR.flag_aerosol_strat.GT.0, ok_aie, flag_aerosol, &
    3969                      flag_aerosol_strat, &
     3974                     flag_aerosol_strat, flag_aer_feedback, &
    39703975                     tau_aero, piz_aero, cg_aero, &
    39713976                     tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm, &

  • LMDZ6/trunk/libf/phylmd/radlwsw_m.F90

    r3117 r3412  
    1717   cldfra, cldemi, cldtaupd,&
    1818   ok_ade, ok_aie, flag_aerosol,&
    19    flag_aerosol_strat,&
     19   flag_aerosol_strat, flag_aer_feedback, &
    2020   tau_aero, piz_aero, cg_aero,&
    2121   tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm,& ! rajoute par OB pour RRTM
    … …  
    102102  ! flag_aerosol-input-I- aerosol flag from 0 to 6
    103103  ! flag_aerosol_strat-input-I- use stratospheric aerosols flag (0, 1, 2)
     104  ! flag_aer_feedback-input-I- activate aerosol radiative feedback (T, F)
    104105  ! tau_ae, piz_ae, cg_ae-input-R- aerosol optical properties (calculated in aeropt.F)
    105106  ! cldtaupi-input-R- epaisseur optique des nuages dans le visible
    … …  
    195196  INTEGER, INTENT(in)  :: flag_aerosol                                   ! takes value 0 (no aerosol) or 1 to 6 (aerosols)
    196197  INTEGER, INTENT(in)  :: flag_aerosol_strat                             ! use stratospheric aerosols
     198  LOGICAL, INTENT(in)  :: flag_aer_feedback                              ! activate aerosol radiative feedback
    197199  REAL,    INTENT(in)  :: cldfra(KLON,KLEV), cldemi(KLON,KLEV), cldtaupd(KLON,KLEV)
    198200  REAL,    INTENT(in)  :: tau_aero(KLON,KLEV,naero_grp,2)                        ! aerosol optical properties (see aeropt.F)
    … …  
    855857         ZTOPLWAD0AERO,ZSOLLWAD0AERO,&
    856858         ZTOPLWAIAERO,ZSOLLWAIAERO, &
    857          ok_ade, ok_aie, flag_aerosol,flag_aerosol_strat) ! flags aerosols
     859         ok_ade, ok_aie, flag_aerosol,flag_aerosol_strat, flag_aer_feedback) ! flags aerosols
    858860           
    859861!        print *,'RADLWSW: apres RECMWF'

  • LMDZ6/trunk/libf/phylmd/rrtm/recmwf_aero.F90

    r3333 r3412  
    3535 & PTOPLWAIAERO,PSOLLWAIAERO,&
    3636!..end
    37  & ok_ade, ok_aie, flag_aerosol,flag_aerosol_strat)
     37 & ok_ade, ok_aie, flag_aerosol,flag_aerosol_strat,flag_aer_feedback)
    3838!--fin
    3939
    … …  
    8484! flag_aerosol-input-I- aerosol flag from 0 to 7
    8585! flag_aerosol_strat-input-I- use stratospheric aerosols flag (T/F)
     86! flag_aer_feedback-input-I- use aerosols radiative effect flag (T/F)
    8687! PPIZA_NAT  : (KPROMA,KLEV,NSW); Single scattering albedo of natural aerosol
    8788! PCGA_NAT   : (KPROMA,KLEV,NSW); Assymetry factor for natural aerosol
    … …  
    214215INTEGER, INTENT(in)  :: flag_aerosol           ! takes value 0 (no aerosol) or 1 to 6 (aerosols)
    215216LOGICAL, INTENT(in)  :: flag_aerosol_strat     ! use stratospheric aerosols
     217LOGICAL, INTENT(in)  :: flag_aer_feedback      ! use aerosols radiative feedback
    216218REAL(KIND=JPRB)   ,INTENT(out)   :: PTOPSWADAERO(KPROMA), PSOLSWADAERO(KPROMA)       ! Aerosol direct forcing at TOA and surface
    217219REAL(KIND=JPRB)   ,INTENT(OUT)   :: PTOPSWAD0AERO(KPROMA), PSOLSWAD0AERO(KPROMA)     ! Aerosol direct forcing at TOA and surface
    … …  
    298300! FALSE: fluxes use no aerosols (case 1)
    299301! to be used only for maintaining bit reproducibility with aerosol diagnostics activated
    300 LOGICAL :: AEROSOLFEEDBACK_ACTIVE = .TRUE.
     302 LOGICAL :: AEROSOLFEEDBACK_ACTIVE ! now externalized from .def files
    301303
    302304!OB - Fluxes including aerosol effects
    … …  
    335337IBEG=KST
    336338IEND=KEND
     339
     340AEROSOLFEEDBACK_ACTIVE = flag_aer_feedback !NL: externalize aer feedback
     341
    337342
    338343!*       1.    PREPARATORY WORK
Note: See TracChangeset for help on using the changeset viewer.