【正文】 t=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 !氣候場(chǎng)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),ix=1,N),iy=1,M) end do 1000 FORMAT() end subroutine qhc !求距平場(chǎng) subroutine jpc(h,s,qh,jp) implicit none parameter N=37,M=17,nyr=4,nmon=12real h(N,M,nmon,nyr),s(N,M,nmon),qh(N,M,nmon),jp(N,M,nmon,nyr)integer ix,iy,it,ittcall qhc(h,qh) do it=1,nmon do iy=1,M do ix=1,N do itt=1,nyr jp(ix,iy,it,itt)=h(ix,iy,it,itt)qh(ix,iy,it) end do end do end do end do do itt=1,nyr do it=1,nmon write(4) ((jp(ix,iy,itt,itt),ix=1,N),iy=1,M) end do end do 1000 FORMAT()end subroutine jpc !求標(biāo)準(zhǔn)差場(chǎng) subroutine bzcc(h,s,qh,jp,bzc) implicit none param