--- trunk/libf/phylmd/Orography/grid_noro_m.f90 2009/12/14 15:25:16 23 +++ trunk/libf/phylmd/Orography/grid_noro_m.f90 2010/12/02 17:11:04 36 @@ -42,7 +42,7 @@ use dimens_m, only: iim, jjm use comconst, only: pi - use numer_rec, only: assert + use nr_util, only: assert REAL, intent(in):: xdata(:), ydata(:) ! coordinates of input field REAL, intent(in):: zdata(:, :) ! input field @@ -410,23 +410,23 @@ integer, PARAMETER:: ISMo=300, JSMo=200 real XF(ISMo, JSMo) real WEIGHTpb(-1:1, -1:1) - real sum + real my_sum integer i, is, js, j if(imar>ismo) stop 'surdimensionner ismo dans mva9 (grid_noro)' if(jmar>jsmo) stop 'surdimensionner jsmo dans mva9 (grid_noro)' - SUM=0. + MY_SUM=0. DO IS=-1, 1 DO JS=-1, 1 WEIGHTpb(IS, JS)=1./FLOAT((1+IS**2)*(1+JS**2)) - SUM=SUM+WEIGHTpb(IS, JS) + MY_SUM=MY_SUM+WEIGHTpb(IS, JS) ENDDO ENDDO DO IS=-1, 1 DO JS=-1, 1 - WEIGHTpb(IS, JS)=WEIGHTpb(IS, JS)/SUM + WEIGHTpb(IS, JS)=WEIGHTpb(IS, JS)/MY_SUM ENDDO ENDDO