42 |
|
|
43 |
use dimens_m, only: iim, jjm |
use dimens_m, only: iim, jjm |
44 |
use comconst, only: pi |
use comconst, only: pi |
45 |
use numer_rec, only: assert |
use nr_util, only: assert |
46 |
|
|
47 |
REAL, intent(in):: xdata(:), ydata(:) ! coordinates of input field |
REAL, intent(in):: xdata(:), ydata(:) ! coordinates of input field |
48 |
REAL, intent(in):: zdata(:, :) ! input field |
REAL, intent(in):: zdata(:, :) ! input field |
410 |
integer, PARAMETER:: ISMo=300, JSMo=200 |
integer, PARAMETER:: ISMo=300, JSMo=200 |
411 |
real XF(ISMo, JSMo) |
real XF(ISMo, JSMo) |
412 |
real WEIGHTpb(-1:1, -1:1) |
real WEIGHTpb(-1:1, -1:1) |
413 |
real sum |
real my_sum |
414 |
integer i, is, js, j |
integer i, is, js, j |
415 |
|
|
416 |
if(imar>ismo) stop 'surdimensionner ismo dans mva9 (grid_noro)' |
if(imar>ismo) stop 'surdimensionner ismo dans mva9 (grid_noro)' |
417 |
if(jmar>jsmo) stop 'surdimensionner jsmo dans mva9 (grid_noro)' |
if(jmar>jsmo) stop 'surdimensionner jsmo dans mva9 (grid_noro)' |
418 |
|
|
419 |
SUM=0. |
MY_SUM=0. |
420 |
DO IS=-1, 1 |
DO IS=-1, 1 |
421 |
DO JS=-1, 1 |
DO JS=-1, 1 |
422 |
WEIGHTpb(IS, JS)=1./FLOAT((1+IS**2)*(1+JS**2)) |
WEIGHTpb(IS, JS)=1./FLOAT((1+IS**2)*(1+JS**2)) |
423 |
SUM=SUM+WEIGHTpb(IS, JS) |
MY_SUM=MY_SUM+WEIGHTpb(IS, JS) |
424 |
ENDDO |
ENDDO |
425 |
ENDDO |
ENDDO |
426 |
|
|
427 |
DO IS=-1, 1 |
DO IS=-1, 1 |
428 |
DO JS=-1, 1 |
DO JS=-1, 1 |
429 |
WEIGHTpb(IS, JS)=WEIGHTpb(IS, JS)/SUM |
WEIGHTpb(IS, JS)=WEIGHTpb(IS, JS)/MY_SUM |
430 |
ENDDO |
ENDDO |
431 |
ENDDO |
ENDDO |
432 |
|
|