Changeset 1549 for trunk/LMDZ.COMMON/libf/dyn3d_common
- Timestamp:
- May 6, 2016, 12:30:29 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified trunk/LMDZ.COMMON/libf/dyn3d_common/infotrac.F90 ¶
r1508 r1549 728 728 729 729 do iq=nqo+1,nqtot 730 write(lunout,*) 'infotrac 569: iq,tnom_0(iq)=',iq,tnom_0(iq)730 ! write(lunout,*) 'infotrac 569: iq,tnom_0(iq)=',iq,tnom_0(iq) 731 731 do phase=1,nqo 732 732 do ixt= 1,niso_possibles 733 733 tnom_trac=trim(tnom_0(phase))//'_' 734 734 tnom_trac=trim(tnom_trac)//trim(tnom_iso(ixt)) 735 write(*,*) 'phase,ixt,tnom_trac=',phase,ixt,tnom_trac735 ! write(*,*) 'phase,ixt,tnom_trac=',phase,ixt,tnom_trac 736 736 IF (tnom_0(iq) == tnom_trac) then 737 write(lunout,*) 'Ce traceur est un isotope'737 ! write(lunout,*) 'Ce traceur est un isotope' 738 738 nb_iso(ixt,phase)=nb_iso(ixt,phase)+1 739 739 nb_isoind(phase)=nb_isoind(phase)+1 … … 742 742 indnum_fn_num(ixt)=iso_indnum(iq) 743 743 phase_num(iq)=phase 744 write(lunout,*) 'iso_num(iq)=',iso_num(iq)745 write(lunout,*) 'iso_indnum(iq)=',iso_indnum(iq)746 write(lunout,*) 'indnum_fn_num(ixt)=',indnum_fn_num(ixt)747 write(lunout,*) 'phase_num(iq)=',phase_num(iq)744 ! write(lunout,*) 'iso_num(iq)=',iso_num(iq) 745 ! write(lunout,*) 'iso_indnum(iq)=',iso_indnum(iq) 746 ! write(lunout,*) 'indnum_fn_num(ixt)=',indnum_fn_num(ixt) 747 ! write(lunout,*) 'phase_num(iq)=',phase_num(iq) 748 748 goto 20 749 749 else if (iqpere(iq).gt.0) then 750 750 if (tnom_0(iqpere(iq)) == tnom_trac) then 751 write(lunout,*) 'Ce traceur est le fils d''un isotope'751 ! write(lunout,*) 'Ce traceur est le fils d''un isotope' 752 752 ! c'est un traceur d'isotope 753 753 nb_traciso(ixt,phase)=nb_traciso(ixt,phase)+1 … … 756 756 zone_num(iq)=nb_traciso(ixt,phase) 757 757 phase_num(iq)=phase 758 write(lunout,*) 'iso_num(iq)=',iso_num(iq)759 write(lunout,*) 'phase_num(iq)=',phase_num(iq)760 write(lunout,*) 'zone_num(iq)=',zone_num(iq)758 ! write(lunout,*) 'iso_num(iq)=',iso_num(iq) 759 ! write(lunout,*) 'phase_num(iq)=',phase_num(iq) 760 ! write(lunout,*) 'zone_num(iq)=',zone_num(iq) 761 761 goto 20 762 762 endif !if (tnom_0(iqpere(iq)) == trim(tnom_0(phase))//trim(tnom_iso(ixt))) then … … 767 767 enddo !do iq=1,nqtot 768 768 769 write(lunout,*) 'iso_num=',iso_num770 write(lunout,*) 'iso_indnum=',iso_indnum771 write(lunout,*) 'zone_num=',zone_num772 write(lunout,*) 'phase_num=',phase_num773 write(lunout,*) 'indnum_fn_num=',indnum_fn_num769 ! write(lunout,*) 'iso_num=',iso_num 770 ! write(lunout,*) 'iso_indnum=',iso_indnum 771 ! write(lunout,*) 'zone_num=',zone_num 772 ! write(lunout,*) 'phase_num=',phase_num 773 ! write(lunout,*) 'indnum_fn_num=',indnum_fn_num 774 774 775 775 do ixt= 1,niso_possibles … … 824 824 ! dimensions isotopique: 825 825 ntraciso=niso*(ntraceurs_zone+1) 826 WRITE(lunout,*) 'niso=',niso827 WRITE(lunout,*) 'ntraceurs_zone,ntraciso=',ntraceurs_zone,ntraciso826 ! WRITE(lunout,*) 'niso=',niso 827 ! WRITE(lunout,*) 'ntraceurs_zone,ntraciso=',ntraceurs_zone,ntraciso 828 828 829 829 ! flags isotopiques: … … 833 833 ok_isotopes=.false. 834 834 endif 835 WRITE(lunout,*) 'ok_isotopes=',ok_isotopes835 ! WRITE(lunout,*) 'ok_isotopes=',ok_isotopes 836 836 837 837 if (ok_isotopes) then … … 843 843 alpha_ideal=(/1.0,1.01,1.006,1.003,1.0/) 844 844 endif !if (ok_isotopes) then 845 WRITE(lunout,*) 'ok_iso_verif=',ok_iso_verif846 WRITE(lunout,*) 'ok_init_iso=',ok_init_iso845 ! WRITE(lunout,*) 'ok_iso_verif=',ok_iso_verif 846 ! WRITE(lunout,*) 'ok_init_iso=',ok_init_iso 847 847 848 848 if (ntraceurs_zone.gt.0) then … … 851 851 ok_isotrac=.false. 852 852 endif 853 WRITE(lunout,*) 'ok_isotrac=',ok_isotrac853 ! WRITE(lunout,*) 'ok_isotrac=',ok_isotrac 854 854 855 855 ! remplissage du tableau iqiso(ntraciso,phase) … … 875 875 index_trac(:,:)=0.0 876 876 endif !if (ok_isotrac) then 877 write(lunout,*) 'index_trac=',index_trac877 ! write(lunout,*) 'index_trac=',index_trac 878 878 879 879 ! Finalize :
Note: See TracChangeset
for help on using the changeset viewer.