天氣學(xué)診斷實習(xí)報告-全文預(yù)覽
【正文】 i=1,nx),j=1,ny),iz=1,nz),it=1,nt)open(7,file=39。)write(5)((((sd(i,j,iz,it),i=1,nx),j=1,ny),iz=1,nz),it=1,nt)!計算渦度平流 do it=1,nt do iz=1,nz do j=2,ny1 do i=2,nx1 wopl(i,j,iz,it)=1./(2*a)*(u(i,j,iz,it)*(wd(i+1,j,iz,it)wd(i1,j,iz,it))/cos((10+d*(j1))*)/dd+v(i,j,iz,it)*(wd(i,j+1,iz,it)wd(i,j1,iz,it))/dd)enddo enddo enddo enddo!計算散度平流 do it=1,nt do iz=1,nz do j=2,ny1 do i=2,nx1 wnpl(i,j,iz,it)=1./(2*a)*(u(i,j,iz,it)*(t(i+1,j,iz,it)t(i1,j,iz,it))/cos((10+d*(j1))*)/dd+v(i,j,iz,it)*(t(i,j+1,iz,it)t(i,j1,iz,it))/dd)enddo enddo enddo enddo!保存渦度平流和散度平流open(6,file=39。)write(4)((((wd(i,j,iz,it),i=1,nx),j=1,ny),iz=1,nz),it=1,nt)open(5,file=39。)do it=1,nt do iz=1,nz do j=1,ny do i=1,nxread(1)u(i,j,iz,it)read(2)v(i,j,iz,it)read(3)t(i,j,iz,it)enddo enddo enddo enddo!計算渦度 do it=1,nt do iz=1,nz do j=2,ny1 do i=2,nx1 wd(i,j,iz,it)=1./(2*a)*((v(i+1,j,iz,it)v(i1,j,iz,it))/cos((10+d*(j1))*)/dd(u(i,j+1,iz,it)u(i,j1,iz,it))/dd+2*u(i,j,iz,it)*tan((10+d*(j1))*))enddo enddo enddo enddo!計算散度 do it=1,nt do iz=1,nz do j=2,ny1 do i=2,nx1 sd(i,j,iz,it)=1./(2*a)*((u(i+1,j,iz,it)u(i1,j,iz,it))/cos((10+d*(j1))*)/dd+(v(i,j+1,iz,it)v(i,j1,iz,it))/dd2*v(i,j,iz,it)*tan((10+d*(j1))*))enddo enddo enddo enddo!保存渦度和散度open(4,file=39。)open(3,file=39。)open(2,file=39。)write(9)((((v(i,j,iz,it),i=1,nx),j=1,ny),iz=1,nz),it=1,nt)end程序2 渦度、散度、渦度平流和溫度平流program woduwendupingliu parameter(a=6371000.,nx=53,ny=29,nz=11,nt=8,d=,dd=*)real wd(nx,ny,nz,nt),sd(nx,ny,nz,nt),wopl(nx,ny,nz,nt),wnpl(nx,ny,nz,nt),u(nx,ny,nz,nt),v(nx,ny,nz,nt),t(nx,ny,nz,nt)open(1,file=39。)write(8)((((u(i,j,iz,it),i=1,nx),j=1,ny),iz=1,nz),it=1,nt)open(9,file=39。)write(7)((((ttd(i,j,iz,it),i=1,nx),j=1,ny),iz=1,nz),it=1,nt)open(8,file=39。)write(6)((((t(i,j,iz,it),i=1,nx),j=1,ny),iz=1,nz),it=1,nt)open(7,file=39。)write(5)((((h(i,j,iz,it),i=1,nx),j=1,ny),iz=1,nz),it=1,nt)open(6,file=39。//sj(it))do i=1,3 read(4,*)enddo do j=ny,1,1 read(4,*)(u(i,j,iz,it),i=1,nx)enddo do j=ny,1,1 read(4,*)(v(i,j,iz,it),i=1,nx)enddo enddo enddo!存放數(shù)據(jù)open(5,file=39。//sj(it))do i=1,4 read(3,*)enddo do j=ny,1,1 read(3,*)(ttd(i,j,iz,it),i=1,nx)enddo open(4,file=39。//sj(it))do i=1,4 read(2,*)enddo do j=ny,1,1 read(2,*)(t(i,j,iz,it),i=1,nx)enddo open(3,file=39。//sj(it))do i=1,4 read(1,*)enddo do j=ny,1,1 read(1,*)(h(i,j,iz,it),i=1,nx)enddo open(2,file=39。 do iz=1,nz do it=1,ntopen(1,file=39。 z(10)=39。 z(8)=39。 z(6)=39。 z(4)=39。 z(2)=39。 sj(8)=39。 sj(6)=39。 sj(4)=39。 sj(2)=39。用第二種修正方案，其中大氣層頂?shù)拇怪彼俣瓤梢灾苯硬捎?，也可以用絕熱法。實習(xí)內(nèi)容：、等壓線、等溫線，了解當(dāng)時環(huán)流形勢（主要是500和850hPa）。所給要素為高度場，風(fēng)場，溫度以及溫度露點差，共11層。實習(xí)過程中遇到一些困難，但經(jīng)過翻閱fortran書籍和向老師請教，我明白了我遇到的困難大都出于同一些問題并解決了它們。26日上升運動達到最強，強烈的垂直上升運動為此次降水的形成提供了充足的動力條件，此時降雨強度也是最強的。⑥500hPa垂直速度分析：在w方程中，負值代表上升，正值代表下沉。5月25日20時，水汽通量輻合，之后最大輻合區(qū)隨著氣旋東移。分析高低空的散度：在我國東部高空有輻散，地空有輻合，對應(yīng)著上升運動，為降水提供有利的條件。②500hPa渦度平流分析：500hpa的8個時次渦度平流差異較大，高低值中心有些散亂，但從正負值所處的區(qū)域來分析，可大致得出，5月25日東經(jīng)95度地區(qū)有正渦度平流中心，對應(yīng)有利于氣旋的形成，隨后正渦度平流中心向東移動，氣旋也隨之向東移動。①環(huán)流形勢分析從850hpa的高度場可以看出，5月25日20時大約在北緯30度、東經(jīng)100度貴州附近有一個氣旋生成，氣旋強度較弱，之后開始加強并向東移動。用第二種修正方案，其中大氣層頂?shù)拇怪彼俣瓤梢灾苯硬捎?，也可以用絕熱法。四、實習(xí)內(nèi)容：、等壓線、等溫線，了解當(dāng)時環(huán)流形勢（主要是500和850hPa）。所給要素為高度場，風(fēng)場，溫度以及溫度露點差，共11層。掌握如何分析得到的各物理量的圖形。第一篇：天氣學(xué)診斷實習(xí)報告南京信息工程大學(xué)天氣學(xué)診斷分析實習(xí)報告實習(xí)名稱天氣學(xué)診斷分析日期2016年12月17日得
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1