Commit 6297df40 authored by agebhard's avatar agebhard

write(*,*) removed

parent a97072ac
......@@ -219,7 +219,7 @@ c write(*,*) "wrong y dimension of tiles"
c IERR=1
c RETURN
c END IF
c write(*,*)"dx,dy:",dx,dy
DO 1 I=1,NY
YG(I)=YSW+DY*(I-1)
1 CONTINUE
......@@ -264,8 +264,8 @@ c loop over all points and pass them to KRIGE:
DO 10 J=1,NY
if (retlm.eq.1) then
PCNT=PCNT+1
BITS(I+NY*(J-1)) = PCNT
end if
BITS(I+NY*(J-1)) = PCNT
X0=XG(I)
Y0=YG(J)
DO0(1)=DOG(I,J)
......
......@@ -12,12 +12,14 @@ c x (LON,LAT)
c
INTEGER I,J,IERR
c write(*,*)"ndo:",ndo
c determine search neighbourhood
DO 1 I=1,N
INDSNW(I)=0
1 CONTINUE
DO 5 I=1,NDO
c write(*,*)"snb for (",LON0(INDDO(I)),":",LAT0(INDDO(I)),")"
c determine search neighbourhood for x0_i
c calculate and sort distances
CALL DVEC(LON0(INDDO(I)),LAT0(INDDO(I)),LON,LAT,DIST,N)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment