资源描述：
《短期气候预测实习程序总结材料》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实用文档本人为南京信息工程大学大气科学系学生，在我大四上时，不幸选修课选了短期气候预测实习课程。其烦人程度超乎我的想象，一边在准备考研和找工作，一边还要花心思完成实习内容，当时根本喘不过气来。为让学弟学妹们不重蹈覆辙，在此共享出我所用的程序，给你们一些帮助。不过需要明确一点，这次实习对于个人编程水平的提高帮助很大，学有余力的同学，应仅把本文当做参考，理解基础上使用，而非不劳而获的资本。以下为正文：具体实习要求参见课本，在此不赘述实习一：大气环流状况的表征programEX1reala(144,73,12,65),ave1(144,73),ave7(144,73
2、),asum(144,73)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
3、,65doimo=1,12doj=1,73doi=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.0enddo文案大全实用文档enddo!7月距平deviationdoi=1,144doj=1,73doit=1,65dev(i,j,it)=a(i,j,imo,it)-a
4、ve7(i,j)enddoenddoenddo!纬圈平均latitudeaveragedoit=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(
5、j,imo,it)enddoenddoenddoenddo!写数据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实习二：大气环流分型PROGRAMEOFCTHISPROGRAMUSESEOFFORANALYSINGTIMESERIESCOFM
6、ETEOROLOGICALFIELDCM:LENTHOFTIMESERIES!!!!!!!!!!m:时间序列长度CN:NUMBEROFGRID-POINTS!!!!!!!!!!n:格点数CKS=-1:SELF;KS=0:DEPATURE;KS=1:STANDERDLIZEDDEPATURECKV:NUMBEROFEIGENVALUESWILLBEOUTPUTCKVT:NUMBEROFEIGENVECTORSANDTIMESERIESWILLBEOUTPUTCMNH=MIN(M,N)CEGVT=EIGENVACTORS,ECOF=TIMECOEFFICIENTS
7、FOREGVT.CER(KV,1)=LAMDA,LAMDAEIGENVALUECER(KV,2)=ACCUMULATELAMDACER(KV,3)=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),
8、EGVT(N,KVT),ECOF(M,
