/[lmdze]/trunk/libf/phylmd/Orography/grid_noro_m.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/Orography/grid_noro_m.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 23 by guez, Mon Dec 14 15:25:16 2009 UTC revision 40 by guez, Tue Feb 22 13:49:36 2011 UTC
# Line 41  contains Line 41  contains
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
# Line 255  contains Line 254  contains
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    
# Line 410  contains Line 409  contains
409      integer, PARAMETER:: ISMo=300, JSMo=200      integer, PARAMETER:: ISMo=300, JSMo=200
410      real XF(ISMo, JSMo)      real XF(ISMo, JSMo)
411      real WEIGHTpb(-1:1, -1:1)      real WEIGHTpb(-1:1, -1:1)
412      real sum      real my_sum
413      integer i, is, js, j      integer i, is, js, j
414    
415      if(imar>ismo) stop 'surdimensionner ismo dans mva9 (grid_noro)'      if(imar>ismo) stop 'surdimensionner ismo dans mva9 (grid_noro)'
416      if(jmar>jsmo) stop 'surdimensionner jsmo dans mva9 (grid_noro)'      if(jmar>jsmo) stop 'surdimensionner jsmo dans mva9 (grid_noro)'
417    
418      SUM=0.      MY_SUM=0.
419      DO IS=-1, 1      DO IS=-1, 1
420         DO JS=-1, 1         DO JS=-1, 1
421            WEIGHTpb(IS, JS)=1./FLOAT((1+IS**2)*(1+JS**2))            WEIGHTpb(IS, JS)=1./FLOAT((1+IS**2)*(1+JS**2))
422            SUM=SUM+WEIGHTpb(IS, JS)            MY_SUM=MY_SUM+WEIGHTpb(IS, JS)
423         ENDDO         ENDDO
424      ENDDO      ENDDO
425    
426      DO IS=-1, 1      DO IS=-1, 1
427         DO JS=-1, 1         DO JS=-1, 1
428            WEIGHTpb(IS, JS)=WEIGHTpb(IS, JS)/SUM            WEIGHTpb(IS, JS)=WEIGHTpb(IS, JS)/MY_SUM
429         ENDDO         ENDDO
430      ENDDO      ENDDO
431    
Legend:
Removed from v.23  
changed lines
  Added in v.40
  ViewVC Help
Powered by ViewVC 1.1.21