41 |
! (d) |
! (d) |
42 |
|
|
43 |
use dimens_m, only: iim, jjm |
use dimens_m, only: iim, jjm |
44 |
use comconst, only: pi |
use nr_util, only: assert, pi |
|
use numer_rec, only: assert |
|
45 |
|
|
46 |
REAL, intent(in):: xdata(:), ydata(:) ! coordinates of input field |
REAL, intent(in):: xdata(:), ydata(:) ! coordinates of input field |
47 |
REAL, intent(in):: zdata(:, :) ! input field |
REAL, intent(in):: zdata(:, :) ! input field |
254 |
zxtzy(ii, jj)=zxtzy(ii, jj)/weight(ii, jj) |
zxtzy(ii, jj)=zxtzy(ii, jj)/weight(ii, jj) |
255 |
ztz(ii, jj) =ztz(ii, jj)/weight(ii, jj) |
ztz(ii, jj) =ztz(ii, jj)/weight(ii, jj) |
256 |
! Standard deviation: |
! Standard deviation: |
257 |
zstd(ii, jj)=sqrt(AMAX1(0., ztz(ii, jj)-zmea(ii, jj)**2)) |
zstd(ii, jj)=sqrt(MAX(0., ztz(ii, jj) - zmea(ii, jj)**2)) |
258 |
ENDDO |
ENDDO |
259 |
ENDDO |
ENDDO |
260 |
|
|