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

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

并行计算

开发平台：

Matlab

dzp2d.f：源码内容
  1. C**********************************************************************
  2. C     D Z P 2 D
  3. C
  4. C     Joachim Falk , Uni Hamburg , (0049)/(0) - 40 / 4123 5050
  5. C     Institute for Geophysics
  6. C     email falk@dkrz.de
  7. C
  8. C       fs = jumper for free surface (1=on, 2=off)
  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 DZP2D(IN, OUT, nx, nz, dz, ZA, ZE, SA, SE, fs)  
  13.       IMPLICIT NONE
  15.       INCLUDE 'param.f'
  17.       REAL IN(MAXCOLS,MAXROWS), OUT(MAXCOLS,MAXROWS), 
  18.      &     dz, G1, G2, G3, F1
  19.       INTEGER nx, nz, r, c, INIT, ZA, ZE, SA, SE, fs
  20.       SAVE F1, G1, G2, G3, INIT
  22.       DATA INIT /0/  
  23.   
  24.       IF (INIT .eq. 0) THEN  
  25.          INIT =   1  
  26.          F1   =   1. / dz           
  27.          G1   =   75. / (64. * dz)  
  28.          G2   =   -25. / (384. * dz)    
  29.          G3   =   3. / (640. * dz)      
  30.       ENDIF  
  32.       DO c = 1 , SA-1
  33.          DO r = 1 , nz
  34.             OUT(c,r) = 0.
  35.          ENDDO
  36.       ENDDO
  37.       DO c = SE+1 , nx
  38.          DO r = 1 , nz
  39.             OUT(c,r) = 0.
  40.          ENDDO
  41.       ENDDO
  42.       DO r = 1 , ZA-1
  43.          DO c = SA , SE
  44.             OUT(c,r) = 0.
  45.          ENDDO
  46.       ENDDO
  48.       DO r =  ZA , ZA+1
  49.          DO c = SA , SE
  50.             OUT(c,r) = F1 * ( IN(c,r+1)- IN(c,r) ) 
  51.          ENDDO
  52.       ENDDO  
  55.       DO r = ZA+2 , ZE-3
  56.          DO c = SA , SE
  57. C            OUT(c,r) =        
  58. C     &          F1 * (IN(c,r+1) - IN(c,r  ))  ! second order
  59.             OUT(c,r) =        
  60.      &          G1 * (IN(c,r+1) - IN(c,r  )) + 
  61.      &          G2 * (IN(c,r+2) - IN(c,r-1)) +  
  62.      &          G3 * (IN(c,r+3) - IN(c,r-2))
  63.          ENDDO
  64.       ENDDO  
  66.       DO r = ZE-2 , ZE
  67.          DO c = SA , SE
  68.             OUT(c,r) =  F1 * ( IN(c,r+1) - IN(c,r) ) 
  69.          ENDDO
  70.       ENDDO  
  71.          
  72. C      DO r = ZE+1 , nz
  73. C         DO c = SA , SE
  74. C            OUT(c,r) = 0.
  75. C         ENDDO
  76. C      ENDDO
  78.       RETURN
  79.       END