/[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 27 by guez, Thu Mar 25 14:29:07 2010 UTC
# Line 410  contains Line 410  contains
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    
Legend:
Removed from v.23  
changed lines
  Added in v.27
  ViewVC Help
Powered by ViewVC 1.1.21