短期气候预测实习程序总结.doc

本人为南京信息工程大学大气科学系学生，在我大四上时，不幸选修课选了短期气候预测实习课程。其烦人程度超乎我的想象，一边在准备考研和找工作，一边还要花心思完成实习内容，当时根本喘不过气来。 为让学弟学妹们不重蹈覆辙，在此共享出我所用的程序，给你们一些帮助。 不过需要明确一点，这次实习对于个人编程水平的提高帮助很大，学有余力的同学，应仅把本文当做参考，理解基础上使用，而非不劳而获的资本。 以下为正文： 具体实习要求参见课本，在此不赘述 实习一：大气环流状况的表征 program EX1 real a(144,73,12,65),ave1(144,73),ave7(144,73),asum(144,73) real dev(144,73,65),latave(73,12,65),latsum(73,12,65) real latdev(144,73,12,65) open(2,file=d:\1\hgt500.grd,form=binary) !补充正确路径 open(4,file=d:\1\ave7.grd,form=binary) !补充数据输出路径 open(5,file=d:\1\dev.grd,form=binary) open(6,file=d:\1\latdev.grd,form=binary) do it=1,65 do imo=1,12 do j=1,73 do i=1,144 read(2)a(i,j,imo,it) enddo;enddo;enddo;enddo cccccccccccccc 请完成以下的程序 !月时间平均 ave7 imo=7 do i=1,144 do j=1,73 do it=1,65 asum(i,j)=asum(i,j)+a(i,j,imo,it) enddo ave7(i,j)=asum(i,j)/65.0 enddo enddo !7月距平 deviation do i=1,144 do j=1,73 do it=1,65 dev(i,j,it)=a(i,j,imo,it)-ave7(i,j) enddo enddo enddo !纬圈平均 latitude average do it=1,65 do imo=1,12 do j=1,73 do i=1,144 latsum(j,imo,it)=latsum(j,imo,it)+a(i,j,imo,it) enddo latave(j,imo,it)=latsum(j,imo,it)/144.0 enddo enddo enddo !纬向偏差 latitude deviation do it=1,65 do imo=1,12 do j=1,73 do i=1,144 latdev(i,j,imo,it)=a(i,j,imo,it)-latave(j,imo,it) enddo enddo enddo enddo !写数据 write(4) ((ave7(i,j),i=1,144),j=1,73) write(5) (((dev(i,j,it),i=1,144),j=1,73),it=1,65) write(6) ((((latdev(i,j,imo,it),i=1,144),j=1,73),imo=1,12),it=1,65) close(4) close(5) close(6) End 实习二：大气环流分型 PROGRAM EOF C THIS PROGRAM USES EOF FOR ANALYSING TIME SERIES C OF METEOROLOGICAL FIELD C M:LENTH OF TIME SERIES !!!!!!!!!! m:时间序列长度 C N:NUMBER OF GRID-POINTS !!!!!!!!!! n:格点数 C KS=-1:SELF; KS=0:DEPATURE; KS=1:STANDERDLIZED DEPATURE C KV:NUMBER OF EIGENVALUES WILL BE OUTPUT C KVT:NUMBER OF EIGENVECTORS AND TIME SERIES WILL BE OUTPUT C MNH=MIN(M,N) C EGVT=EIGENVACTORS, ECOF=TIME COEFFICIENTS FOR EGVT. C ER(KV,1)=LAMDA,LAMDA EIGENVALUE C