亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

? 歡迎來到蟲蟲下載站！ | ?? 資源下載 ?? 資源專輯 ?? 關于我們
? 蟲蟲下載站

?? fdtd3d_cpml.f90

?? FDTD的fortran程序
?? F90
?? 第 1 頁 / 共 3 頁
字號: 
!.....................................................................
         DO i = 1,nxPML_1-1
	      psi_Hyx_1(i,j,k) = bh_x_1(i)*psi_Hyx_1(i,j,k)                   &
				+ ch_x_1(i)*(Ez(i+1,j,k) - Ez(i,j,k))/dx
	      Hy(i,j,k) = Hy(i,j,k) + DB*psi_Hyx_1(i,j,k)
         ENDDO
!.....................................................................
!  PML for top Hy, i-direction
!.....................................................................
         ii = nxPML_2-1
         DO i = Imax+1-nxPML_2,Imax-1
	      psi_Hyx_2(ii,j,k) = bh_x_2(ii)*psi_Hyx_2(ii,j,k)     &
				+ ch_x_2(ii)*(Ez(i+1,j,k) -    &
                                Ez(i,j,k))/dx
	      Hy(i,j,k) = Hy(i,j,k) + DB*psi_Hyx_2(ii,j,k)
            ii = ii-1
	   ENDDO
     ENDDO
   ENDDO
   DO i = 1,Imax-1
      DO j = 1,Jmax-1
!.....................................................................
!  PML for bottom Hy, k-direction
!.....................................................................
         DO k = 2,nzPML_1
	      psi_Hyz_1(i,j,k-1) = bh_z_1(k-1)*psi_Hyz_1(i,j,k-1)          &
				+ ch_z_1(k-1)*(Ex(i,j,k-1) - Ex(i,j,k))/dz
	      Hy(i,j,k) = Hy(i,j,k) + DB*psi_Hyz_1(i,j,k-1)
         ENDDO
!.....................................................................
!  PML for top Hy, k-direction
!.....................................................................
         kk = nzPML_2-1
         DO k = Kmax+1-nzPML_2,Kmax-1
	    psi_Hyz_2(i,j,kk) = bh_z_2(kk)*psi_Hyz_2(i,j,kk)               &
				+ ch_z_2(kk)*(Ex(i,j,k-1) -    &
                                Ex(i,j,k))/dz
	    Hy(i,j,k) = Hy(i,j,k) + DB*psi_Hyz_2(i,j,kk)
            kk = kk-1
         ENDDO
     ENDDO
   ENDDO

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  UPDATE Hz
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   DO k = 1,Kmax-1
      DO i = 1,Imax-1
        DO j = 1,Jmax-1
            Hz(i,j,k) = DA * Hz(i,j,k) + DB       &
                  * ((Ey(i,j,k) - Ey(i+1,j,k))*den_hx(i) +        &
			    (Ex(i,j+1,k) - Ex(i,j,k))*den_hy(j))
	   ENDDO
      ENDDO
      DO j = 1,Jmax-1
!.....................................................................
!  PML for bottom Hz, x-direction
!.....................................................................
         DO i = 1,nxPML_1-1
   	      psi_Hzx_1(i,j,k) = bh_x_1(i)*psi_Hzx_1(i,j,k)                 &
	 			+ ch_x_1(i) *(Ey(i,j,k) - Ey(i+1,j,k))/dx
	      Hz(i,j,k) = Hz(i,j,k) + DB*psi_Hzx_1(i,j,k)
         ENDDO
!.....................................................................
!  PML for top Hz, x-direction
!.....................................................................
         ii = nxPML_2-1
         DO i = Imax+1-nxPML_2,Imax-1
   	      psi_Hzx_2(ii,j,k) = bh_x_2(ii)*psi_Hzx_2(ii,j,k)            &
	 			+ ch_x_2(ii) *(Ey(i,j,k) -       &
                                Ey(i+1,j,k))/dx
	      Hz(i,j,k) = Hz(i,j,k) + DB*psi_Hzx_2(ii,j,k)
            ii = ii-1
	   ENDDO
      ENDDO
      DO i = 1,Imax-1
!.....................................................................
!  PML for bottom Hz, y-direction
!.....................................................................
         DO j = 1,nyPML_1-1
            psi_Hzy_1(i,j,k) = bh_y_1(j)*psi_Hzy_1(i,j,k)                   &
				+ ch_y_1(j)*(Ex(i,j+1,k) - Ex(i,j,k))/dy
	      Hz(i,j,k) = Hz(i,j,k) + DB*psi_Hzy_1(i,j,k)
         ENDDO
!.....................................................................
!  PML for top Hz, y-direction
!.....................................................................
         jj = nyPML_2-1
         DO j = Jmax+1-nyPML_2,Jmax-1
            psi_Hzy_2(i,jj,k) = bh_y_2(jj)*psi_Hzy_2(i,jj,k)               &
				+ ch_y_2(jj)*(Ex(i,j+1,k) -    &
                                Ex(i,j,k))/dy
	      Hz(i,j,k) = Hz(i,j,k) + DB*psi_Hzy_2(i,jj,k)
            jj = jj-1
	   ENDDO
      ENDDO
   ENDDO

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  UPDATE Ex
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   DO k = 1,Kmax-1
      DO i = 1,Imax-1
	   DO j = 2,Jmax-1
              IF (i >= istart-1 .and. i <= iend .and. j >= jstart .and.  &
                   j <= jend .and. k >= kstart .and. k <= kend) THEN
                 Ex(i,j,k) = 0.0
              ELSE
	         Ex(i,j,k) = CA(i,j,k) * Ex(i,j,k) + CB(i,j,k) *       &
  			( (Hz(i,j,k) - Hz(i,j-1,k))*den_ey(j)  +    &
			  (Hy(i,j,k) - Hy(i,j,k+1))*den_ez(k) )
              ENDIF
	   ENDDO
	ENDDO 
      DO i = 1,Imax-1
!.....................................................................
!  PML for bottom Ex, j-direction
!.....................................................................
         DO j = 2,nyPML_1
  	      psi_Exy_1(i,j,k) = be_y_1(j)*psi_Exy_1(i,j,k)                 &
	 			+ ce_y_1(j) *(Hz(i,j,k) - Hz(i,j-1,k))/dy
            Ex(i,j,k) = Ex(i,j,k) + CB(i,j,k)*psi_Exy_1(i,j,k)
         ENDDO
!.....................................................................
!  PML for top Ex, j-direction
!.....................................................................
         jj = nyPML_2
         DO j = Jmax+1-nyPML_2,Jmax-1
  	      psi_Exy_2(i,jj,k) = be_y_2(jj)*psi_Exy_2(i,jj,k)       &
	 			+ ce_y_2(jj) *(Hz(i,j,k) -       &
                                Hz(i,(j-1),k))/dy
            Ex(i,j,k) = Ex(i,j,k) + CB(i,j,k)*psi_Exy_2(i,jj,k)
            jj = jj-1
         ENDDO
	ENDDO
   ENDDO
   DO i = 1,Imax-1
      DO j = 2,Jmax-1
!.....................................................................
!  PML for bottom Ex, k-direction
!.....................................................................
         DO k = 1,nzPML_1
  	      psi_Exz_1(i,j,k) = be_z_1(k)*psi_Exz_1(i,j,k)                 &
	 			+ ce_z_1(k) *(Hy(i,j,k) - Hy(i,j,k+1))/dz
            Ex(i,j,k) = Ex(i,j,k) + CB(i,j,k)*psi_Exz_1(i,j,k)
         ENDDO
!.....................................................................
!  PML for top Ex, k-direction
!.....................................................................
         kk = nzPML_2
         DO k = Kmax-nzPML_2,Kmax-1
  	      psi_Exz_2(i,j,kk) = be_z_2(kk)*psi_Exz_2(i,j,kk)             &
	 			+ ce_z_2(kk) *(Hy(i,j,k) -       &
                                Hy(i,j,k+1))/dz
            Ex(i,j,k) = Ex(i,j,k) + CB(i,j,k)*psi_Exz_2(i,j,kk)
            kk = kk-1
         ENDDO
	ENDDO
   ENDDO

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  UPDATE Ey
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   DO k = 1,Kmax-1
      DO i = 2,Imax-1
	   DO j = 1,Jmax-1
              IF (i >= istart .and. i <= iend .and. j >= jstart-1 .and. &
                   j <= jend .and. k >= kstart .and. k <= kend) THEN
                 Ey(i,j,k) = 0.0
              ELSE
                 Ey(i,j,k) = CA(i,j,k) * Ey(i,j,k) + CB(i,j,k) *    &
			( (Hz(i-1,j,k) - Hz(i,j,k))*den_ex(i) +         &
			  (Hx(i,j,k+1) - Hx(i,j,k))*den_ez(k) )
              ENDIF
	   ENDDO 
      ENDDO
      DO j = 1,Jmax-1
!.....................................................................
!  PML for bottom Ey, i-direction
!.....................................................................
         DO i = 2,nxPML_1
	      psi_Eyx_1(i,j,k) = be_x_1(i)*psi_Eyx_1(i,j,k)                   &
				+ ce_x_1(i)*(Hz(i-1,j,k) - Hz(i,j,k))/dx
	      Ey(i,j,k) = Ey(i,j,k) + CB(i,j,k)*psi_Eyx_1(i,j,k)
         ENDDO
!.....................................................................
!  PML for top Ey, i-direction
!.....................................................................
         ii = nxPML_2
         DO i = Imax+1-nxPML_2,Imax-1
	      psi_Eyx_2(ii,j,k) = be_x_2(ii)*psi_Eyx_2(ii,j,k)              &
				+ ce_x_2(ii)*(Hz(i-1,j,k) -    &
                               Hz(i,j,k))/dx
	      Ey(i,j,k) = Ey(i,j,k) + CB(i,j,k)*psi_Eyx_2(ii,j,k)
            ii = ii-1
	   ENDDO
     ENDDO
   ENDDO
   DO i = 2,Imax-1
      DO j = 1,Jmax-1
!.....................................................................
!  PML for bottom Ey, k-direction
!.....................................................................
         DO k = 1,nzPML_1
	      psi_Eyz_1(i,j,k) = be_z_1(k)*psi_Eyz_1(i,j,k)                   &
				+ ce_z_1(k)*(Hx(i,j,k+1) - Hx(i,j,k))/dz
	      Ey(i,j,k) = Ey(i,j,k) + CB(i,j,k)*psi_Eyz_1(i,j,k)
         ENDDO
!.....................................................................
!  PML for top Ey, k-direction
!.....................................................................
         kk = nzPML_2
         DO k = Kmax-nzPML_2,Kmax-1
	    psi_Eyz_2(i,j,kk) = be_z_2(kk)*psi_Eyz_2(i,j,kk)               &
				+ ce_z_2(kk)*(Hx(i,j,k+1) -    &
                                Hx(i,j,k))/dz
	    Ey(i,j,k) = Ey(i,j,k) + CB(i,j,k)*psi_Eyz_2(i,j,kk)
            kk = kk-1
         ENDDO
     ENDDO
   ENDDO
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  UPDATE Ez
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   DO k = 2,Kmax-1
      DO i = 2,Imax-1
         DO j = 2,Jmax-1
            Ez(i,j,k) = CA(i,j,k) * Ez(i,j,k) + CB(i,j,k)       &
                  * ((Hy(i,j,k) - Hy(i-1,j,k))*den_ex(i) +        &
			    (Hx(i,j-1,k) - Hx(i,j,k))*den_ey(j))
	   ENDDO
      ENDDO
      DO j = 2,Jmax-1
!.....................................................................
!  PML for bottom Ez, x-direction
!.....................................................................
         DO i = 2,nxPML_1
   	      psi_Ezx_1(i,j,k) = be_x_1(i)*psi_Ezx_1(i,j,k)             &
	 			+ ce_x_1(i) *(Hy(i,j,k) - Hy(i-1,j,k))/dx
	      Ez(i,j,k) = Ez(i,j,k) + CB(i,j,k)*psi_Ezx_1(i,j,k)
         ENDDO
!.....................................................................
!  PML for top Ez, x-direction
!.....................................................................
         ii = nxPML_2
         DO i = Imax+1-nxPML_2,Imax-1
   	      psi_Ezx_2(ii,j,k) = be_x_2(ii)*psi_Ezx_2(ii,j,k)       &
	 			+ ce_x_2(ii) *(Hy(i,j,k) -       &
                                Hy(i-1,j,k))/dx
	      Ez(i,j,k) = Ez(i,j,k) + CB(i,j,k)*psi_Ezx_2(ii,j,k)
            ii = ii-1
	   ENDDO
      ENDDO
      DO i = 2,Imax-1
!.....................................................................
!  PML for bottom Ez, y-direction
!.....................................................................
         DO j = 2,nyPML_1
            psi_Ezy_1(i,j,k) = be_y_1(j)*psi_Ezy_1(i,j,k)            &
				+ ce_y_1(j)*(Hx(i,j-1,k) - Hx(i,j,k))/dy
	      Ez(i,j,k) = Ez(i,j,k) + CB(i,j,k)*psi_Ezy_1(i,j,k)
         ENDDO
!.....................................................................
!  PML for top Ez, y-direction
!.....................................................................
         jj = nyPML_2
         DO j = Jmax+1-nyPML_2,Jmax-1
            psi_Ezy_2(i,jj,k) = be_y_2(jj)*psi_Ezy_2(i,jj,k)         &
				+ ce_y_2(jj)*(Hx(i,j-1,k) -    &
                                Hx(i,j,k))/dy
	      Ez(i,j,k) = Ez(i,j,k) + CB(i,j,k)*psi_Ezy_2(i,jj,k)
            jj = jj-1
	   ENDDO
      ENDDO
   ENDDO

!-----------------------------------------------------------------------
!   SOURCE
!-----------------------------------------------------------------------
   i = isource
   j = jsource
   k = ksource
   source = -2.0*((n*dt-tO)/tw) * exp(-((n*dt-tO)/tw)**2.0)
   Ez(i,j,k) = Ez(i,j,k) - CB(i,j,k)*source

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  RECORD GRID FOR VISUALIZATION
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   IF (n == record_grid) then
      DO j = 1,Jmax
         DO i = 1,Imax      
            write(33,*)Ez(i,j,ksource+1)
	   ENDDO
	ENDDO
   CLOSE(UNIT = 33)
   endif   
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  WRITE TO OUTPUT FILES
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  
   P1 = Ez(isource,jsource,ksource)
   write(30,*)P1      
   P2 = Ey(irecv1,jrecv1,krecv1)
   write(31,*)P2      

   IF (mod(n,10) == 0) then
    WRITE(*,*)n, " of ", nmax
   ENDIF

   ENDDO
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  END TIME STEP
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!.:. .:. .:. .:. .:. .:. .:. .:. .:. .:. .:. .:. .:. .:. .:. .:. .:. .:.
    WRITE(*,*)"done time-stepping"

!-----------------------------------------------------------------------
! CLOSE OUTPUT FILES
!-----------------------------------------------------------------------   
   CLOSE(UNIT = 30)
   CLOSE(UNIT = 31)

   END PROGRAM fdtd3D_CPML
!cccccccc1ccccccccc2ccccccccc3ccccccccc4ccccccccc5cccgtgccc6ccccccccc7cc

?? 快捷鍵說明

復制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
亚洲一区免费在线观看| 欧美精品一区二区久久婷婷| 中文久久乱码一区二区| 国产精品一区二区x88av| 国产日韩欧美精品综合| 成人午夜免费av| 亚洲啪啪综合av一区二区三区| 色噜噜狠狠色综合中国| 亚洲午夜精品在线| 欧美一区二区免费| 国产成人免费网站| 亚洲精品日韩综合观看成人91| 色999日韩国产欧美一区二区| 亚洲成人免费电影| 欧美精品一区二区三区在线 | 免费成人av资源网| 国产亚洲va综合人人澡精品| 91女人视频在线观看| 香蕉乱码成人久久天堂爱免费| 日韩女优av电影| 成人激情电影免费在线观看| 亚洲一区欧美一区| 一区二区三区在线看| 91官网在线免费观看| 精品在线播放免费| 亚洲图片另类小说| 制服丝袜在线91| 成人毛片老司机大片| 亚洲午夜久久久久| 国产欧美精品在线观看| 欧美三区免费完整视频在线观看| 精品在线一区二区三区| 亚洲视频一区二区在线| 日韩视频在线你懂得| 97久久超碰国产精品| 捆绑紧缚一区二区三区视频| 国产精品视频一二三区| 91精品国产91综合久久蜜臀| 本田岬高潮一区二区三区| 午夜精品视频一区| 国产精品久久久久久户外露出| 欧美年轻男男videosbes| 国产成人日日夜夜| 免费美女久久99| 亚洲图片欧美综合| 国产精品看片你懂得 | 欧美一区二区三区不卡| 成人99免费视频| 国产在线精品视频| 天堂蜜桃一区二区三区| 亚洲欧美日韩人成在线播放| 久久精品一区二区三区四区| 欧美巨大另类极品videosbest| 99国产精品一区| 国产成人8x视频一区二区| 久久国产人妖系列| 午夜久久久久久电影| 亚洲精品欧美专区| 国产精品国产三级国产aⅴ无密码| 精品国产乱码久久久久久夜甘婷婷 | 国产日韩欧美一区二区三区综合 | 国产自产高清不卡| 亚洲h在线观看| 一区二区三区欧美视频| 中文字幕亚洲在| 国产精品天美传媒| 日本一二三不卡| 国产婷婷精品av在线| 久久亚洲精精品中文字幕早川悠里| 欧美日韩你懂得| 欧美精品一卡二卡| 欧美日韩精品欧美日韩精品 | 秋霞电影网一区二区| 午夜欧美电影在线观看| 国产精品不卡在线| 国产精品欧美久久久久一区二区 | 中文字幕欧美日韩一区| 国产色产综合产在线视频| 久久综合久久久久88| 久久综合久色欧美综合狠狠| 久久免费电影网| 2023国产一二三区日本精品2022| 精品88久久久久88久久久| www激情久久| 中文欧美字幕免费| 日韩一区中文字幕| 一区二区三区在线免费播放| 一区二区三区日韩| 日韩国产高清在线| 精品一区二区三区免费| 国产精品996| 91亚洲精华国产精华精华液| 一本到一区二区三区| 欧美私模裸体表演在线观看| 91麻豆精品国产91久久久久久 | 亚洲综合图片区| 五月综合激情婷婷六月色窝| 美女任你摸久久| 国产91在线|亚洲| 91视频免费观看| 欧美喷水一区二区| 精品区一区二区| 欧美激情一区在线| 亚洲色图视频免费播放| 一级特黄大欧美久久久| 午夜成人免费电影| 激情六月婷婷久久| 91影院在线免费观看| 欧美日韩精品综合在线| 2024国产精品| 亚洲男人的天堂在线aⅴ视频| 天堂在线亚洲视频| 国产伦精品一区二区三区免费| 91视视频在线直接观看在线看网页在线看| 在线中文字幕一区| 久久综合丝袜日本网| 亚洲卡通欧美制服中文| 麻豆精品在线视频| 色哟哟欧美精品| 久久先锋资源网| 中文字幕一区在线观看视频| 日韩经典中文字幕一区| 岛国一区二区三区| 日韩一区二区精品| 亚洲女人小视频在线观看| 九九九久久久精品| 欧日韩精品视频| 亚洲国产成人一区二区三区| 全部av―极品视觉盛宴亚洲| 成+人+亚洲+综合天堂| 欧美videossexotv100| 一区二区在线电影| 国产电影一区二区三区| 欧美精品日韩一区| 亚洲精品中文字幕在线观看| 久久99精品久久只有精品| 99久久久久久99| 久久夜色精品国产噜噜av| 亚洲国产欧美日韩另类综合| 国产suv精品一区二区6| 欧美一区二区视频在线观看| 亚洲免费观看高清| 成人美女视频在线观看18| 欧美r级电影在线观看| 一区二区三区影院| 成人av资源在线| 国产亚洲欧洲997久久综合| 免费观看久久久4p| 欧美精品一卡两卡| 亚洲国产日韩在线一区模特| 色婷婷综合五月| 国产精品久久久久久久浪潮网站| 日本亚洲三级在线| 欧美色大人视频| 亚洲国产日韩一区二区| 色视频欧美一区二区三区| 中文字幕一区二区三| 国产不卡高清在线观看视频| 26uuu另类欧美| 国产精品一区二区黑丝 | 欧美激情在线一区二区| 国产在线不卡一卡二卡三卡四卡| 欧美一级一级性生活免费录像| 亚洲成人激情综合网| 欧美午夜精品久久久久久孕妇| 日韩毛片视频在线看| 成人国产在线观看| 国产精品青草久久| 97成人超碰视| 一区二区三区在线影院| 欧美综合一区二区| 亚洲成av人**亚洲成av**| 欧美日韩精品一区视频| 日本一区中文字幕| 欧美成人伊人久久综合网| 伦理电影国产精品| 久久久久久久久久看片| 粉嫩aⅴ一区二区三区四区| 欧美国产激情一区二区三区蜜月| 丁香亚洲综合激情啪啪综合| 国产精品国产三级国产aⅴ入口| 91影视在线播放| 亚洲国产精品天堂| 日韩欧美国产系列| 国产99久久久久| 伊人婷婷欧美激情| 91麻豆精品国产91久久久久久| 久久国产福利国产秒拍| 国产日韩欧美在线一区| 97精品久久久午夜一区二区三区 | 精品88久久久久88久久久| 丁香网亚洲国际| 亚洲精品亚洲人成人网在线播放| 欧美日精品一区视频| 久久国产精品无码网站| 亚洲国产成人午夜在线一区| 在线亚洲人成电影网站色www| 视频一区二区三区在线| 久久亚洲一区二区三区明星换脸 | 久久一日本道色综合|