资源描述：
《短期气候预测实习程序总结.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、本人为南京信息工程大学大气科学系学生，在我大四上时，不幸选修课选了短期气候预测实习课程。其烦人程度超乎我的想象，一边在准备考研和找工作，一边还要花心思完成实习内容，当时根本喘不过气来。为让学弟学妹们不重蹈覆辙，在此共享出我所用的程序，给你们一些帮助。不过需要明确一点，这次实习对于个人编程水平的提高帮助很大，学有余力的同学，应仅把本文当做参考，理解基础上使用，而非不劳而获的资本。以下为正文：具体实习要求参见课本，在此不赘述实习一：大气环流状况的表征programEX1reala(144,73,12,65),ave1(144,73),ave7(144,73),asum(144,73
2、)realdev(144,73,65),latave(73,12,65),latsum(73,12,65)reallatdev(144,73,12,65)open(2,file='d:1hgt500.grd',form='binary')!补充正确路径open(4,file='d:1ave7.grd',form='binary')!补充数据输出路径open(5,file='d:1dev.grd',form='binary')open(6,file='d:1latdev.grd',form='binary')doit=1,65doimo=1,12doj=1,73d
3、oi=1,144read(2)a(i,j,imo,it)enddo;enddo;enddo;enddocccccccccccccc请完成以下的程序!月时间平均ave7imo=7doi=1,144doj=1,73doit=1,65asum(i,j)=asum(i,j)+a(i,j,imo,it)enddoave7(i,j)=asum(i,j)/65.0enddoenddo!7月距平deviationdoi=1,144doj=1,73doit=1,65dev(i,j,it)=a(i,j,imo,it)-ave7(i,j)enddoenddoenddo!纬圈平均latitudeave
4、ragedoit=1,65doimo=1,12doj=1,73doi=1,144latsum(j,imo,it)=latsum(j,imo,it)+a(i,j,imo,it)enddolatave(j,imo,it)=latsum(j,imo,it)/144.0enddoenddoenddo!纬向偏差latitudedeviationdoit=1,65doimo=1,12doj=1,73doi=1,144latdev(i,j,imo,it)=a(i,j,imo,it)-latave(j,imo,it)enddoenddoenddoenddo!写数据write(4)((ave7(
5、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实习二：大气环流分型PROGRAMEOFCTHISPROGRAMUSESEOFFORANALYSINGTIMESERIESCOFMETEOROLOGICALFIELDCM:LENTHOFTIMESERIES!!!!!!!!!!m:时间序列长度CN:NUMBER
6、OFGRID-POINTS!!!!!!!!!!n:格点数CKS=-1:SELF;KS=0:DEPATURE;KS=1:STANDERDLIZEDDEPATURECKV:NUMBEROFEIGENVALUESWILLBEOUTPUTCKVT:NUMBEROFEIGENVECTORSANDTIMESERIESWILLBEOUTPUTCMNH=MIN(M,N)CEGVT=EIGENVACTORS,ECOF=TIMECOEFFICIENTSFOREGVT.CER(KV,1)=LAMDA,LAMDAEIGENVALUECER(KV,2)=ACCUMULATELAMDACER(KV,3)=
7、THESUMOFCOMPONENTSVECTORSPROJECTEDONTOcEIGENVACTOR.CER(KV,4)=ACCUMULATEER(KV,3)CPARAMETER(M=61,N=41*21,MNH=61,KS=1,KV=8,KVT=8,pi=3.1415926)CDIMENSIONF(N,M),A(MNH,MNH),S(MNH,MNH),ER(MNH,4),*DF(N),V(MNH),AVF(N),EGVT(N,KVT),ECOF(M,KVT)dimensionhh(144,
