Changeset 4625 for LMDZ6/trunk
- Timestamp:
- Jul 12, 2023, 11:31:50 PM (17 months ago)
- Location:
- LMDZ6/trunk/libf
- Files:
-
- 2 added
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.F90
r4601 r4625 43 43 REAL paprs_glo(klon_glo,klev+1) 44 44 45 REAL, POINTER:: latitude(:)45 REAL, ALLOCATABLE :: latitude(:) 46 46 ! (of input data sorted in strictly ascending order) 47 47 48 REAL, POINTER:: longitude(:)48 REAL, ALLOCATABLE :: longitude(:) 49 49 ! (of input data sorted in strictly ascending order) 50 50 51 REAL, POINTER:: time(:)51 REAL, ALLOCATABLE :: time(:) 52 52 ! (of input data sorted in strictly ascending order) 53 53 54 REAL, POINTER:: lev(:)54 REAL, ALLOCATABLE :: lev(:) 55 55 ! levels of input data 56 56 -
LMDZ6/trunk/libf/phylmd/StratAer/strataer_emiss_mod.F90
r4601 r4625 99 99 CALL getin_p('xlat_sai',xlat_sai) 100 100 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) 101 107 102 108 ELSEIF (flag_emit == 3) THEN ! SAI between latitudes … … 169 175 WRITE(lunout,*) 'IN STRATAER : altemiss_sai',altemiss_sai 170 176 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 171 180 WRITE(lunout,*) 'IN STRATAER : xlat_min_sai',xlat_min_sai 172 181 WRITE(lunout,*) 'IN STRATAER : xlat_max_sai',xlat_max_sai -
LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.F90
r4601 r4625 99 99 REAL,SAVE :: sigma_alt_sai ! standard deviation of emission altitude in m 100 100 !$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) 101 110 REAL,SAVE :: xlat_sai ! latitude of SAI in degree 102 111 REAL,SAVE :: xlon_sai ! longitude of SAI in degree -
LMDZ6/trunk/libf/phylmd/lsc_scav.F90
r4514 r4625 99 99 ! Droplet size shrinks by evap 100 100 frac_aer=0.5 101 debut=.FALSE. 101 102 ! 102 103 OPEN(99,file='lsc_scav_param.data',status='old', & -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r4619 r4625 1232 1232 1233 1233 REAL pi 1234 1234 INTEGER ieru 1235 1235 1236 1236 !======================================================================! … … 1545 1545 piz_aero(:,:,:,:) = 1. 1546 1546 cg_aero(:,:,:,:) = 0. 1547 d_q_ch4(:,:) = 0. 1547 1548 1548 1549 IF (aerosol_couple .AND. (config_inca /= "aero" & … … 4846 4847 !DC Calcul de la tendance due au methane 4847 4848 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) 4849 4852 ! ajout de la tendance d'humidite due au methane 4850 4853 d_q_ch4_dtime(:,:) = d_q_ch4(:,:)*phys_tstep -
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r4619 r4625 1325 1325 1326 1326 REAL pi 1327 INTEGER ieru 1327 1328 1328 1329 pi = 4. * ATAN(1.) … … 6143 6144 !DC Calcul de la tendance due au methane 6144 6145 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) 6146 6149 ! ajout de la tendance d'humidite due au methane 6147 6150 d_q_ch4_dtime(:,:) = d_q_ch4(:,:)*phys_tstep
Note: See TracChangeset
for help on using the changeset viewer.