1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
|
do j=1,nby
utmy(j) = yorigr + (j-1)*dgrid
do i=1,nbx
utmx(i) = xorigr + (i-1)*dgrid
call utm_ll(iutmzn,utmx(i),utmy(j),lat,lon)
call ll_xy(lat,lon,hpolelat,hpolelon,cal_xmn(i,j),cal_ymn(i,j)) ! 2 fonctions validées
! ici cal_xmn et cal_ymn sont corrects :
print*,cal_xmn(i,j),cal_ymn(i,j)
enddo
! ici cal_xmn et cal_ymn sont toujours corrects :
do i=1,nbx
print*,cal_xmn(i,j),cal_ymn(i,j)
enddo
enddo
!ici, ils ne sont plus bons! les valeurs semblent 'mélangées', par
!exemple, cal_xmn(1,1) ici est différent de cal_xmn(1,1) dans la boucle!
do j=1,nby
do i=1,nbx
print*,cal_xmn(i,j),cal_ymn(i,j)
enddo
enddo |
Partager