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

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

并行计算

开发平台：

Matlab

dzp2d_12.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, G4, G5, G6
  19.       INTEGER nx, nz, r, c, INIT, ZA, ZE, SA, SE, fs
  20.       SAVE F1, G1, G2, G3, G4, G5, G6,  INIT
  22.       DATA INIT /0/  
  23.   
  24.       IF (INIT .eq. 0) THEN  
  25.          INIT =   1  
  26.          F1   =   1. / dz           
  27.          G1   =   1.22133636474609     / dz
  28.          G2   =  -0.09693145751953116 / dz
  29.          G3   =   0.01744766235351561 / dz
  30.          G4   =  -2.967289515904012E-3 /dz
  31.          G5   =   3.590053982204854E-4 /dz
  32.          G6   =  -2.184781161221585E-5 /dz
  33.       ENDIF  
  35.       DO c = 1 , SA-1
  36.          DO r = 1 , nz
  37.             OUT(c,r) = 0.
  38.          ENDDO
  39.       ENDDO
  40.       DO c = SE+1 , nx
  41.          DO r = 1 , nz
  42.             OUT(c,r) = 0.
  43.          ENDDO
  44.       ENDDO
  45.       DO r = 1 , ZA-1
  46.          DO c = SA , SE
  47.             OUT(c,r) = 0.
  48.          ENDDO
  49.       ENDDO
  51.       DO r =  ZA , ZA+4
  52.          DO c = SA , SE
  53.             OUT(c,r) = F1 * ( IN(c,r+1)- IN(c,r) ) 
  54.          ENDDO
  55.       ENDDO  
  58.       DO r = ZA+5 , ZE-6
  59.          DO c = SA , SE
  60.             OUT(c,r) =        
  61.      &          G1 * (IN(c,r+1) - IN(c,r  )) + 
  62.      &          G2 * (IN(c,r+2) - IN(c,r-1)) +  
  63.      &          G3 * (IN(c,r+3) - IN(c,r-2)) +
  64.      &          G4 * (IN(c,r+4) - IN(c,r-3)) +
  65.      &          G5 * (IN(c,r+5) - IN(c,r-4)) +
  66.      &          G6 * (IN(c,r+6) - IN(c,r-5))
  67.          ENDDO
  68.       ENDDO  
  70.       DO r = ZE-5 , ZE
  71.          DO c = SA , SE
  72.             OUT(c,r) =  F1 * ( IN(c,r+1) - IN(c,r) ) 
  73.          ENDDO
  74.       ENDDO  
  75.          
  76. C      DO r = ZE+1 , nz
  77. C         DO c = SA , SE
  78. C            OUT(c,r) = 0.
  79. C         ENDDO
  80. C      ENDDO
  82.       RETURN
  83.       END