From 1bfaf942e2b546ce2a78f177d00617ad0a40424e Mon Sep 17 00:00:00 2001 From: Lauri Tuppi Date: Tue, 28 May 2024 08:36:24 +0300 Subject: [PATCH] Test: Junk removed from SATELLITE/STATST/fortran-programs/plots_for_one_location.f95 --- .../plots_for_one_location.f95 | 56 ------------------- 1 file changed, 56 deletions(-) diff --git a/SATELLITE/STATST/fortran-programs/plots_for_one_location.f95 b/SATELLITE/STATST/fortran-programs/plots_for_one_location.f95 index 4d5e293ce..e8e1eeff0 100644 --- a/SATELLITE/STATST/fortran-programs/plots_for_one_location.f95 +++ b/SATELLITE/STATST/fortran-programs/plots_for_one_location.f95 @@ -256,11 +256,7 @@ program plots_for_one_station print*,' Reading sim_file: ',sim_file open(unit=1,form='formatted',file=sim_file,status='old') do while(.true.) - !do i=1,nlines 11 read(1,*,err=11,end=12) yyear,mmonth,dday,hhour,f1 -!11 read(1,*,err=11,end=12) station_char,yyear,mmonth,dday,f1,hhour -! read(1,*) station_char,yyear,mmonth,dday,f1,hhour - !print*,yyear,mmonth,dday,f1,hhour if(l_round_6h)then call round_6h(yyear,mmonth,dday,hhour,year,month,day,hour) ! Rounding to nearest 6-h? else @@ -269,18 +265,12 @@ program plots_for_one_station day=dday hour=12 !hhour endif - !print*,i,yyear,mmonth,dday,f1,hhour -! print*,'hep1' - !print*,1+hour/6 if((year.gt.year1.or.(year.eq.year1.and.month.ge.month1)).and. & ! Is the time within the (year.lt.year2.or.(year.eq.year2.and.month.le.month2))) then ! wanted range? - !print*,'hep' if(.not.(month.eq.2.and.day.eq.29))then ! leap days are skipped from plotting f(year-year1+1,day_of_year(month,day),1)=f1 ! hard-coding to 6-hourly resolution endif endif - !print*,'hep2',f1, f(year-year1+1,day_of_year(month,day),1) - !print*,i,yyear,mmonth,dday,f1,hhour enddo 12 continue close(1) @@ -293,7 +283,6 @@ program plots_for_one_station do while(.true.) if(l_code_in_char)then 21 read(1,*,err=21,end=23)station_char,lon,lat,month,day,hour,(quant1(i),i=1,nquant) - !print*,station_char,lon,lat,month,day,hour else 22 read(1,*,err=22,end=23)station_int,lon,lat,month,day,hour,(quant1(i),i=1,nquant) endif @@ -309,40 +298,25 @@ program plots_for_one_station ! ! 3) quantile bin frequencies and related MSD and p-value statistics ! -! open(unit=2,form='formatted',status='old',file=rank_histogram_file) -! nlines = 0 -! DO -! READ(2,*,iostat=io) -! IF (io/=0) EXIT -! nlines = nlines + 1 -! ENDDO -! close(2) freq=miss if(include_rank_histograms)then print*,' Reading rank histograms: ',rank_histogram_file open(unit=1,form='formatted',file=rank_histogram_file,status='old') do while(.true.) - !do j=2,nlines if(l_code_in_char)then - !print*,'hep0' 31 read(1,*,err=31,end=33)station_char,lon,lat,hour,ntot,msd1,p_value1,(freq1(i),i=1,nquant+1) else 32 read(1,*,err=32,end=33)station_int,lon,lat,hour,ntot,msd1,p_value1,(freq1(i),i=1,nquant+1) read(1,*)station_int,lon,lat,hour,ntot,msd1,p_value1,(freq1(i),i=1,nquant+1) endif - !if((hour.eq.0.or.hour.eq.6.or.hour.eq.12.or.hour.eq.18.or.hour.eq.24))then - !print*,'hep1' if((hour.eq.12))then - !print*,'hep1',hour,msd1,p_value1 msd(1)=msd1 p_value(1)=p_value1 do i=1,nquant+1 freq(i,1)=freq1(i) enddo - !print*,'hep2',msd(1),p_value(1) endif - !print*,'hep2' enddo 33 continue close(1) @@ -399,11 +373,9 @@ program plots_for_one_station if(nyear.eq.1)then num_int_arr=nint(1.+n_colors)/2. write(1,'(A11,I2)')'set ccolor ',color(num_int_arr) - !write(1,'(A11,I2)')'set ccolor ',color(nint(1.+n_colors)/2.) else num_int_arr=nint(1+(i-1.)/(nyear-1.)*(n_colors-1.)) write(1,'(A11,I2)')'set ccolor ',color(num_int_arr) - !write(1,'(A11,I2)')'set ccolor ',color(nint(1+(i-1.)/(nyear-1.)*(n_colors-1.))) endif write(1,'(A7,I2,A1)')'d &0(z=',i,')' enddo @@ -415,30 +387,6 @@ program plots_for_one_station close(1) if(include_rank_histograms)then -! open(unit=1,form='formatted',file='msd_and_p-value_00') -! write(1,'(A15)')'set strsiz 0.17' -! write(1,'(A18)')'set string 1 l 4 0' -! write(1,'(A25,F5.3)')'draw string 1.3 &0 MSD=',msd(1) -! write(1,'(A29,F5.3)')'draw string 3.0 &0 p-value=',p_value(1) -! close(1) -! open(unit=1,form='formatted',file='msd_and_p-value_06') -! write(1,'(A15)')'set strsiz 0.17' -! write(1,'(A18)')'set string 1 l 4 0' -! write(1,'(A25,F5.3)')'draw string 1.3 &0 MSD=',msd(2) -! write(1,'(A29,F5.3)')'draw string 3.0 &0 p-value=',p_value(2) -! close(1) -! open(unit=1,form='formatted',file='msd_and_p-value_12') -! write(1,'(A15)')'set strsiz 0.17' -! write(1,'(A18)')'set string 1 l 4 0' -! write(1,'(A25,F5.3)')'draw string 1.3 &0 MSD=',msd(3) -! write(1,'(A29,F5.3)')'draw string 3.0 &0 p-value=',p_value(3) -! close(1) -! open(unit=1,form='formatted',file='msd_and_p-value_18') -! write(1,'(A15)')'set strsiz 0.17' -! write(1,'(A18)')'set string 1 l 4 0' -! write(1,'(A25,F5.3)')'draw string 1.3 &0 MSD=',msd(4) -! write(1,'(A29,F5.3)')'draw string 3.0 &0 p-value=',p_value(4) -! close(1) open(unit=1,form='formatted',file='msd_and_p-value') write(1,'(A15)')'set strsiz 0.17' write(1,'(A18)')'set string 1 l 4 0' @@ -508,10 +456,6 @@ program plots_for_one_station do j=1,365 !do hour=0,18,6 do hour=12,12 - !print*,'###########################################' - !print*,year-year1+1,j,1+hour/6 !f(1,1,1) - !print*,f(year-year1+1,j,1+hour/6),miss - !print*,f(:,j,1) if(abs(f(year-year1+1,j,1)).lt.abs(miss))then if(l_code_in_char)then write(1,'(A3,2F16.6,1X,I4,1X,I3,1X,I2,F16.6)')station_char,lon,lat,year,j,hour,f(year-year1+1,j,1) -- GitLab