开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 并行计算 > 查看源码
dxm2d_12.f
资源名称：wave_forward.rar　[点击查看]
上传用户：ls4004p
上传日期：2007-08-05
资源大小：2314k
文件大小：2k
源码类别：

并行计算

开发平台：

Matlab

dxm2d_12.f：源码内容
  1. C**********************************************************************
  2. C       D X M 2 D 
  3. C
  4. C       Joachim Falk , Uni Hamburg , (0049)/(0) - 40 / 4123 5047
  5. C       Institute for Geophysics
  6. C       11.07.1996
  7. C       email   falk@dkrz.de
  8. C
  9. C       ZA, ZE = first and last grid-row    of differentiation
  10. C       SA, SE = first and last grid-column of dIFferentiation
  11. C**********************************************************************
  12.       SUBROUTINE DXM2D(IN,OUT, nx, nz, dx, ZA, ZE, SA, SE)  
  13.       IMPLICIT NONE  
  15.       INCLUDE 'param.f'
  17.       REAL IN(MAXCOLS,MAXROWS), OUT(MAXCOLS,MAXROWS),
  18.      &     G1, G2, G3, F1, G4, G5, G6, dx
  19.       INTEGER nx, nz, r, c, INIT,
  20.      &        ZA, ZE, SA, SE
  21.       
  22.       SAVE F1, G1, G2, G3, G4, G5, G6, INIT
  23.       DATA INIT /0/  
  24.            
  25.       IF (INIT .eq. 0) THEN  
  26.          INIT =   1  
  27.          F1   =   1. / dx          
  28.          G1   =   1.22133636474609     / dx
  29.          G2   =  -0.09693145751953116 / dx
  30.          G3   =   0.01744766235351561 / dx
  31.          G4   =  -2.967289515904012E-3 /dx
  32.          G5   =   3.590053982204854E-4 /dx
  33.          G6   =  -2.184781161221585E-5 /dx
  34.       ENDIF  
  36.        DO r=1,nz
  37.           DO c = 1 , SA-1
  38.              OUT(c,r) = 0.  
  39.           ENDDO  
  40.           DO c = SE+1 , nx
  41.              OUT(c,r) = 0.  
  42.           ENDDO  
  43.        ENDDO
  45.        DO r=1, ZA-1
  46.           DO c=SA,SE
  47.              OUT(c,r) = 0.  
  48.           ENDDO
  49.        ENDDO  
  50.       
  51.       DO r = ZA , ZE  
  52.          DO c = SA , SA+5
  53.             OUT(c,r) = F1 * (IN(c,r) - IN(c-1,r))  
  54.          ENDDO  
  55.          DO c = SA+6 , SE-5
  56.             OUT(c,r) =        
  57.      &          G1 * (IN(c  ,r) - IN(c-1,r)) +  
  58.      &          G2 * (IN(c+1,r) - IN(c-2,r)) +  
  59.      &          G3 * (IN(c+2,r) - IN(c-3,r)) +
  60.      &          G4 * (IN(c+3,r) - IN(c-4,r)) +
  61.      &          G5 * (IN(c+4,r) - IN(c-5,r)) +
  62.      &          G6 * (IN(c+5,r) - IN(c-6,r))
  63.          ENDDO
  64.          DO c = SE-4 , SE
  65.             OUT(c,r) = F1 * (IN(c,r) - IN(c-1,r))  
  66.          ENDDO  
  67.       ENDDO 
  68.       
  69.        DO r = ZE+1 , nz
  70.           DO c = SA , SE
  71.              OUT(c,r) = 0.0
  72.           ENDDO
  73.        ENDDO
  75.       RETURN
  76.       END