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

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

?? fdtd3d_cpml.f90

?? 基于matlab的三維時域有限差分程序（pec邊界條件）
?? 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一区二区三区免费野_久草精品视频
国产福利精品导航| 欧美精选一区二区| 欧美日韩免费电影| 国产欧美视频一区二区| 亚洲成a人片在线观看中文| 国产精品性做久久久久久| 欧美日韩亚洲丝袜制服| 中文字幕在线一区| 精品一区二区在线视频| 国产99久久久久久免费看农村| 欧美色中文字幕| 欧美专区亚洲专区| 911精品产国品一二三产区| 亚洲视频在线观看三级| 国产+成+人+亚洲欧洲自线| 欧美成人免费网站| 天天色天天操综合| 欧美午夜不卡视频| 亚洲综合激情网| 色婷婷综合五月| 国产精品毛片久久久久久| 国产精华液一区二区三区| 日韩精品在线一区| 日产欧产美韩系列久久99| 欧美四级电影网| 亚洲国产视频直播| 欧美日韩一区在线观看| 亚洲摸摸操操av| 91黄色免费版| 亚洲国产精品尤物yw在线观看| 97精品电影院| 一区二区三区蜜桃网| 欧洲精品视频在线观看| 亚洲福利一二三区| 欧美理论在线播放| 蜜桃视频在线观看一区二区| 欧美三区在线观看| 国产成人精品亚洲777人妖| 国产美女在线精品| 国产乱码精品一区二区三区忘忧草| 男女视频一区二区| 国产成人一区在线| 蜜桃精品视频在线观看| 久久99精品久久久久久久久久久久| 日日摸夜夜添夜夜添亚洲女人| 婷婷成人激情在线网| 国产在线国偷精品产拍免费yy| 成人激情视频网站| 欧美经典一区二区| 成人免费毛片嘿嘿连载视频| 最新热久久免费视频| 欧美羞羞免费网站| 狠狠色丁香久久婷婷综合丁香| 久久久久久久一区| 色av一区二区| 日韩黄色小视频| 国产日韩欧美制服另类| 色伊人久久综合中文字幕| 亚洲a一区二区| 久久久亚洲欧洲日产国码αv| 成人一区二区三区视频在线观看| 一区二区日韩av| 欧美刺激午夜性久久久久久久 | 99视频在线精品| 一区二区三区在线播| 精品国产3级a| 91久久精品一区二区二区| 久久精品国产77777蜜臀| 国产精品视频免费| 欧美精品乱码久久久久久| 成人午夜看片网址| 蜜桃一区二区三区四区| 日韩一区日韩二区| 2020日本不卡一区二区视频| 色香色香欲天天天影视综合网| 麻豆国产91在线播放| 中文字幕一区二区三区精华液 | 成人免费高清在线观看| 亚洲福利一二三区| 最近日韩中文字幕| www久久精品| 欧美人与性动xxxx| 99riav一区二区三区| 久久国产欧美日韩精品| 亚洲国产精品一区二区www在线| 国产亚洲精品精华液| 日韩一区二区免费高清| 日本精品免费观看高清观看| 久久国产免费看| 亚洲色图清纯唯美| 久久久久久久综合日本| 国产suv一区二区三区88区| 日韩毛片在线免费观看| 国产**成人网毛片九色| 成人欧美一区二区三区| 欧美三级日韩三级国产三级| 日韩av一级电影| 日韩精品一区二区三区中文不卡 | 欧美系列在线观看| 国内精品免费在线观看| 五月婷婷欧美视频| 亚洲高清免费一级二级三级| 亚洲婷婷综合色高清在线| 国产精品网站在线观看| 久久精品亚洲国产奇米99| 精品久久一区二区| 欧美一区二区美女| 91精品国产91热久久久做人人| 777久久久精品| 欧美影院午夜播放| 在线免费不卡视频| 日本韩国视频一区二区| 91丨porny丨国产| 99精品久久99久久久久| 91蝌蚪porny九色| 91免费国产视频网站| 色拍拍在线精品视频8848| 91在线精品一区二区| 91美女在线观看| 在线观看精品一区| 欧美吻胸吃奶大尺度电影| 欧美日韩一区二区三区视频 | 亚洲bt欧美bt精品777| 亚洲1区2区3区4区| 日本不卡一区二区三区 | 日韩电影一二三区| 美女视频黄a大片欧美| 精品一区二区三区在线观看国产 | 欧美国产精品专区| 综合电影一区二区三区| 亚洲一卡二卡三卡四卡五卡| 日本欧洲一区二区| 国产成人免费视频网站| 91视视频在线观看入口直接观看www| 欧美久久久久免费| 久久久国产精华| 久久精品二区亚洲w码| 久久精品国产亚洲a| 麻豆91免费观看| 国产成人在线网站| 国产精品99久久久久久似苏梦涵| 亚洲综合免费观看高清在线观看| 亚洲一二三四在线| 首页国产欧美日韩丝袜| 国产成人在线看| 欧美精品vⅰdeose4hd| 2024国产精品视频| 久久久精品免费观看| 亚洲男女一区二区三区| 美女视频黄频大全不卡视频在线播放| 国产成人精品免费| 欧美日韩国产片| 欧美国产国产综合| 日韩经典中文字幕一区| 波波电影院一区二区三区| 欧美人体做爰大胆视频| 亚洲国产精品成人综合色在线婷婷| 依依成人综合视频| 国产成人av电影在线观看| 欧美久久久久久久久| 成人欧美一区二区三区白人| 蜜桃免费网站一区二区三区| 91麻豆高清视频| 久久久久久夜精品精品免费| 亚洲丰满少妇videoshd| 99在线热播精品免费| 精品成人在线观看| 日韩成人一级片| 日本精品免费观看高清观看| 国产人成亚洲第一网站在线播放| 午夜精品在线看| 色婷婷久久久综合中文字幕 | 99视频精品在线| 2欧美一区二区三区在线观看视频 337p粉嫩大胆噜噜噜噜噜91av | 久久99精品久久久久| 色哟哟国产精品| 中文字幕成人在线观看| 国产揄拍国内精品对白| 日韩午夜在线播放| 天堂影院一区二区| 在线亚洲一区观看| 亚洲日本欧美天堂| 不卡视频在线看| 国产精品美女www爽爽爽| 国产成人免费网站| 国产日韩精品一区二区浪潮av| 日本亚洲电影天堂| 欧美日韩大陆一区二区| 亚洲午夜影视影院在线观看| 97久久久精品综合88久久| 国产精品久久久99| 不卡在线观看av| 国产精品色一区二区三区| 福利视频网站一区二区三区| 久久精品夜夜夜夜久久| 精品一区二区三区在线视频| 亚洲精品一区在线观看| 国内精品国产成人国产三级粉色| 精品日韩欧美在线| 国产精一品亚洲二区在线视频|