氣候變化及其診斷實習報告(已修改)
2025-08-26 23:01 本頁面
【正文】 氣候變化及其診斷實習報告姓名： 學號：20091349026班級：2009級氣候一班實習一（1）實習目的：編fortran程序，求500hPa高度場的氣候場、距平場、標準化場，標準差場，并能用Grads做出圖形。（2）程序：program dlimplicit none parameter N=37,M=17,nyr=4,nmon=12real h(N,M,nmon,nyr),s(N,M,nmon),qh(N,M,nmon),bzc(N,M,nmon,nyr),jp(N,M,nmon,nyr)integer ix,iy,it,ittopen(1,file=39。f:/sx/h500/39。,form=39。binary39。)open(2,file=39。f:/sx/h500/39。)open(3,file=39。f:/sx/h500/39。,form=39。binary39。)open(4,file=39。f:/sx/h500/39。,form=39。binary39。)open(5,file=39。f:/sx/h500/39。,form=39。binary39。)open(6,file=39。f:/sx/h500/39。,form=39。binary39。) do itt=1,nyr do it=1,nmon read(1) ((h(ix,iy,it,itt),ix=1,N),iy=1,M) write(2,1000) ((h(ix,iy,itt,itt),ix=1,N),iy=1,M) end do end docall qhc(h,qh)call jpc(h,s,qh,jp)call bzcc(h,s,qh,jp,bzc)call bzhh(h,qh,s,jp,bzc)1000 FORMAT()end program !氣候場subroutine qhc(h,qh)implicit noneparameter N=37,M=17,nyr=4,nmon=12real h(N,M,nmon,nyr),s(N,M,nmon),qh(N,M,nmon)integer ix,iy,it,itts=0do it=1,nmon do ix=1,N do iy=1,M do itt=1,nyr s(ix,iy,it)=s(ix,iy,it)+h(ix,iy,it,itt) end do qh(ix,iy,it)=s(ix,iy,it)/nyr end do end do end do do it=1,nmon write(3) ((qh(ix,iy,it),i
環(huán)評公示相關推薦
鄂ICP備17016276號-1