Changeset 4625


Ignore:
Timestamp:
Jul 12, 2023, 11:31:50 PM (10 months ago)
Author:
dcugnet
Message:

StratAer? commit, N. Lebas

Location:
LMDZ6/trunk/libf
Files:
2 added
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.F90

    r4601 r4625  
    4343  REAL paprs_glo(klon_glo,klev+1)
    4444
    45   REAL, POINTER:: latitude(:)
     45  REAL, ALLOCATABLE :: latitude(:)
    4646! (of input data sorted in strictly ascending order)
    4747
    48   REAL, POINTER:: longitude(:)
     48  REAL, ALLOCATABLE :: longitude(:)
    4949! (of input data sorted in strictly ascending order)
    5050
    51   REAL, POINTER:: time(:)
     51  REAL, ALLOCATABLE :: time(:)
    5252! (of input data sorted in strictly ascending order)
    5353
    54   REAL, POINTER:: lev(:)
     54  REAL, ALLOCATABLE :: lev(:)
    5555! levels of input data
    5656

  • LMDZ6/trunk/libf/phylmd/StratAer/strataer_emiss_mod.F90

    r4601 r4625  
    9999       CALL getin_p('xlat_sai',xlat_sai)
    100100       CALL getin_p('xlon_sai',xlon_sai)
     101       CALL getin_p('year_emit_sai_start',year_emit_sai_start)
     102       CALL getin_p('year_emit_sai_end',year_emit_sai_end)
     103       CALL getin_p('mth_emit_sai_start',mth_emit_sai_start)
     104       CALL getin_p('mth_emit_sai_end',mth_emit_sai_end)
     105       CALL getin_p('day_emit_sai_start',day_emit_sai_start)
     106       CALL getin_p('day_emit_sai_end',day_emit_sai_end)
    101107       
    102108    ELSEIF (flag_emit == 3) THEN ! SAI between latitudes
    … …  
    169175          WRITE(lunout,*) 'IN STRATAER : altemiss_sai',altemiss_sai
    170176          WRITE(lunout,*) 'IN STRATAER : sigma_alt_sai',sigma_alt_sai
     177          WRITE(lunout,*) 'IN STRATAER : year_emit_sai start/end',year_emit_sai_start,year_emit_sai_end
     178          WRITE(lunout,*) 'IN STRATAER : mth_emit_sai start/end',mth_emit_sai_start,mth_emit_sai_end
     179          WRITE(lunout,*) 'IN STRATAER : day_emit_sai start/end',day_emit_sai_start,day_emit_sai_end
    171180          WRITE(lunout,*) 'IN STRATAER : xlat_min_sai',xlat_min_sai
    172181          WRITE(lunout,*) 'IN STRATAER : xlat_max_sai',xlat_max_sai

  • LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.F90

    r4601 r4625  
    9999  REAL,SAVE    :: sigma_alt_sai          ! standard deviation of emission altitude in m
    100100  !$OMP THREADPRIVATE(m_aer_emiss_sai, altemiss_sai, sigma_alt_sai)
     101  INTEGER,SAVE    :: year_emit_sai_start
     102  INTEGER,SAVE    :: year_emit_sai_end
     103  INTEGER,SAVE    :: mth_emit_sai_start
     104  INTEGER,SAVE    :: mth_emit_sai_end
     105  INTEGER,SAVE    :: day_emit_sai_start
     106  INTEGER,SAVE    :: day_emit_sai_end
     107  !$OMP THREADPRIVATE(year_emit_sai_start, year_emit_sai_end)
     108  !$OMP THREADPRIVATE(mth_emit_sai_start, mth_emit_sai_end)
     109  !$OMP THREADPRIVATE(day_emit_sai_start, day_emit_sai_end)
    101110  REAL,SAVE    :: xlat_sai               ! latitude of SAI in degree
    102111  REAL,SAVE    :: xlon_sai               ! longitude of SAI in degree

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

    r4514 r4625  
    9999!     Droplet size shrinks by evap
    100100      frac_aer=0.5
     101      debut=.FALSE.
    101102!
    102103      OPEN(99,file='lsc_scav_param.data',status='old', &

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

    r4619 r4625  
    12321232
    12331233    REAL pi
    1234 
     1234    INTEGER ieru
    12351235
    12361236    !======================================================================!
    … …  
    15451545       piz_aero(:,:,:,:) = 1.
    15461546       cg_aero(:,:,:,:)  = 0.
     1547       d_q_ch4(:,:) = 0.
    15471548
    15481549       IF (aerosol_couple .AND. (config_inca /= "aero" &
    … …  
    48464847    !DC Calcul de la tendance due au methane
    48474848    IF (ok_qch4) THEN
    4848        CALL METHOX(1,klon,klon,klev,q_seri,d_q_ch4,pplay)
     4849       CALL stratH2O_methox(debut,paprs,d_q_ch4)
     4850!      d_q_ch4: H2O source in ppmv/sec - ecmwf routine METHOX
     4851!       CALL METHOX(1,klon,klon,klev,q_seri,d_q_ch4,pplay)
    48494852       ! ajout de la tendance d'humidite due au methane
    48504853       d_q_ch4_dtime(:,:) = d_q_ch4(:,:)*phys_tstep

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

    r4619 r4625  
    13251325   
    13261326    REAL pi
     1327    INTEGER ieru
    13271328
    13281329    pi = 4. * ATAN(1.)
    … …  
    61436144    !DC Calcul de la tendance due au methane
    61446145    IF (ok_qch4) THEN
    6145        CALL METHOX(1,klon,klon,klev,q_seri,d_q_ch4,pplay)
     6146       CALL stratH2O_methox(debut,paprs,d_q_ch4)
     6147!      d_q_ch4: H2O source in ppmv/sec - ecmwf routine METHOX
     6148!       CALL METHOX(1,klon,klon,klev,q_seri,d_q_ch4,pplay)
    61466149       ! ajout de la tendance d'humidite due au methane
    61476150       d_q_ch4_dtime(:,:) = d_q_ch4(:,:)*phys_tstep
Note: See TracChangeset for help on using the changeset viewer.