- Timestamp:
- Oct 17, 2024, 5:48:37 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/aeronomars/photochemistry.F90
r3290 r3462 161 161 !integer,parameter :: i_hdo2 = 38 162 162 163 integer , save:: i_last164 !$OMP THREADPRIVATE(i_last) 163 integer :: i_last 164 165 165 !Tracer indexes for photionization coeffs 166 166 … … 180 180 181 181 integer :: ilay 182 integer , save:: ind_norec183 integer , save:: ind_orec184 !$OMP THREADPRIVATE(ind_norec, ind_orec) 182 integer :: ind_norec 183 integer :: ind_orec 184 185 185 real :: ctimestep ! standard timestep for the chemistry (s) 186 186 real :: dt_guess ! first-guess timestep (s) … … 503 503 em_no(:)=c(:,i_o)*c(:,i_n)*v_4(:,ind_norec) !2.8e-17*(300./temp(:)))**0.5 504 504 em_o2(:)=0.75*c(:,i_o)*c(:,i_o)*c(:,i_co2)*v_3(:,ind_orec) !2.5*9.46e-34*exp(485./temp(:))*dens(:) 505 !print*, "ind_norec= ", ind_norec 506 !print*, "ind_orec= ", ind_orec 505 507 506 !=================================================================== 508 507 ! save chemical species for the gcm … … 660 659 661 660 use comcstfi_h 662 !use photolysis_mod, only : nphot661 use photolysis_mod, only : nphot 663 662 664 663 implicit none … … 676 675 logical, intent(in) :: deutchem 677 676 integer :: lswitch ! interface level between lower 678 integer,parameter :: nphot0 =13679 677 ! atmosphere and thermosphere chemistries 680 678 real, dimension(nlayer) :: dens ! total number density (molecule.cm-3) … … 697 695 real (kind = 8), dimension(nlayer,nb_reaction_3_max) :: v_3 698 696 real (kind = 8), dimension(nlayer,nb_reaction_4_max) :: v_4 699 integer , intent(out):: ind_norec700 integer , intent(out):: ind_orec697 integer :: ind_norec 698 integer :: ind_orec 701 699 702 700 !---------------------------------------------------------------------- … … 705 703 706 704 integer :: ilev 707 integer, save :: nb_phot, nb_reaction_3, nb_reaction_4 708 !$OMP THREADPRIVATE(nb_phot, nb_reaction_3, nb_reaction_4) 705 integer :: nb_phot, nb_reaction_3, nb_reaction_4 709 706 real :: ak0, ak1, xpo, rate, rate1, rate2 710 707 real :: k1a0, k1b0, k1ainf, k1a, k1b, fc, fx, x, y, gam … … 740 737 ! initialisation 741 738 !---------------------------------------------------------------------- 742 nb_phot = 0 743 nb_phot = nphot 0+ nphotion ! initialised to the number of photolysis + number of photoionization rates739 740 nb_phot = nphot + nphotion ! initialised to the number of photolysis + number of photoionization rates 744 741 nb_reaction_3 = 0 745 742 nb_reaction_4 = 0 … … 779 776 v_3(:,nb_reaction_3) = a002(:) 780 777 ind_orec=nb_reaction_3 781 ! print*,"ind_orec",ind_orec 778 782 779 !--- a003: o + o3 -> o2 + o2 783 780 … … 1203 1200 v_4(:,nb_reaction_4) = d007(:) 1204 1201 ind_norec = nb_reaction_4 1205 ! print*,"ind_norec=",ind_norec 1202 1206 1203 !--- d008: n + ho2 -> no + oh 1207 1204 … … 2461 2458 ! local 2462 2459 2463 integer, save :: nb_phot, nb_reaction_3, nb_reaction_4 2464 !$OMP THREADPRIVATE(nb_phot, nb_reaction_3, nb_reaction_4) 2460 integer :: nb_phot, nb_reaction_3, nb_reaction_4 2465 2461 integer :: i_dummy 2466 2462
Note: See TracChangeset
for help on using the changeset viewer.