氣象程序設(shè)計(jì)及繪圖實(shí)習(xí)報告

............wd....wd..實(shí)習(xí)要求：依據(jù)所供給的實(shí)習(xí)說明和資料，按要求完成實(shí)習(xí)內(nèi)容；將完成的實(shí)習(xí)內(nèi)容按步驟寫入實(shí)習(xí)報告中，包括資料說明、Fortran程序、CTL文件、gs文件、繪制的圖形〔要求在圖形上方用drawtitle的方式添加標(biāo)題，注明姓名拼音及學(xué)號〕依據(jù)所繪圖形及查找相關(guān)文獻(xiàn)，說明本次實(shí)習(xí)所求物理量場的意義。實(shí)習(xí)內(nèi)容：Nino34海溫指數(shù)與海平面氣壓相關(guān)圖Nino34160站降水相關(guān)圖1Nino341月份海平面氣壓相關(guān)圖資料說明Nino34指數(shù)：Nino34.txt19511202312Nino3.4區(qū)〔熱帶太平170度~1205度~5度〕區(qū)域平均海溫指數(shù)資料，共63121-12月的海溫指數(shù)；全球海平面氣壓月平均數(shù)據(jù)：slp.jan.grd195120231月全球海平面氣63144*732.5?*2.5?；步驟：Fortrancorr.grid.f901951-2023年1Nino34海溫指數(shù)與海平面氣corr.slp.grd文件中；corr.grid.f90文件：programcx open(7,file=”D:\expriments\grads\corr.slp.implicitnone txt”)realt(63),p(73,144,63),r(73,144),relation, open(8,file=”D:\expriments\grads\corr.slp.z(63) grd”,form=”binary”)character*12c(63) write(7,”(144F5.2)”)character*7a(63) ((r(i,j),j=1,144),i=1,73)integeri,j,k write(8)((r(i,j),j=1,144),i=1,73)open(5,file=”D:\expriments\grads\slp.jan. endgrd”,form=”binary”) functionrelation(a,b,n)open(6,file=”D:\expriments\grads\nino34. implicitnonetxt”) integer,intent(in)::ndoi=1,63 real,intent(in)::a(n),b(n)read(6,”(A12)”)c(i) real::relation!返回的相關(guān)系數(shù)enddo integer::i,j !循環(huán)掌握變量doi=1,63 real::sfenzi,sfenmu1,sfenmu2,s!加法器read(c(i),”(A7,F5.2)”)a(i),t(i) real::amean,bmean!a,b向量的平均值enddo !計(jì)算平均值dok=1,63 s=0.read(5)((p(i,j,k),j=1,144),i=1,73) doi=1,nenddo s=s+a(i)!print*,p enddodoi=1,73 amean=s/ndoj=1,144 s=0.dok=1,63 doi=1,nz(k)=p(i,j,k) s=s+b(i)enddo enddoenddor(i,j)=relation(z,t,63) enddoenddo!print*,r

sfenzi=0.sfenmu1=0.sfenmu2=0. sfenmu1=sfenmu1+(a(i)-amean)**2doi=1,n sfenmu2=sfenmu2+(b(i)-bmean)**2sfenzi=sfenzi+(a(i)-amean)*(b(i)- enddobmean)relation=sfenzi/sqrt(sfenmu1*sfenmu2)endfunctionrelationcorr.slp.grd文件中。corr.slp.grd1.ctl：dsete:\grads\20230706qimo\test1\corr.slp.grdtitleslpundef-9.99e33xdef144linear-1802.5ydef73linear-902.5zdef1levels1000tdef63linearjun19511yrvars1slp00sealevelpressureendvars1.gs文件：”reinit””opene:\grads\20230706qimo\test1\1.ctl””setlev1000””setlon-180180””setlat-9090””sett1””setgxoutshaded””setgradsoff””setblack-0.20.2””dslp””setgxoutcontour””dslp””drawtitleLuYu””cbar1.00””printime:\grads\20230706qimo\test1\1.pngwhite”;1951-2023年1Nino34海溫指數(shù)與海平面氣壓相關(guān)系數(shù)圖：1951-2023年1Nino34海溫指數(shù)與海平面氣壓相關(guān)系數(shù)圖1Nino347160站降水相關(guān)圖：資料說明Nino34指數(shù)：Nino34.txt19511202312Nino3.4區(qū)〔熱帶太平170度~1205度~5度〕區(qū)域平均海溫指數(shù)資料，共63121-12月的海溫指數(shù)；160站降水資料：r1607.txt1601951~2023637月降水資料，數(shù)19517id.txt160站都寫入后再19527月各站降水量，...................20237月各站降水量依據(jù)id.txt中各站挨次寫入；160站站點(diǎn)資料：id.txt160160行。每行第一8位字符型的站號，后面兩個數(shù)字分別為該站點(diǎn)緯度和經(jīng)度；步驟：Fortrancorr.station.f901951-20231Nino347月160站降水相關(guān)系數(shù)，計(jì)算結(jié)果保存于corr.7.txt和corr.7.grd文件中；programcx doi=1,160implicitnone read(10,*)stid(i),lat(i),lon(i)integerp(160,63),i,j,k enddoreal print*,stids(63),r(160),relation,t(63),lon(160),lat(16 doj=1,630),tim read(5,*)(p(i,j),i=1,160)character*12c(63) enddocharacter*7a(63) doi=1,63character*8stid(160) read(6,”(A12)”)c(i)integernlev,nflag enddoopen(6,file=”E:\grads\20230706qimo\test2 doi=1,63\nino34.txt”) read(c(i),”(A7,F5.2)”)a(i),t(i)open(5,file=”E:\grads\20230706qimo\test2 enddo\r1607.txt”) doi=1,160open(10,file=E:\grads\20230706qimo\test doj=1,631\id.txt”) s(j)=p(i,j)enddo real::relation!返回的相關(guān)系數(shù)r(i)=relation(s,t,63) integer::i,j !循環(huán)掌握變量enddo real::sfenzi,sfenmu1,sfenmu2,s!加法器!計(jì)算平均值!print*,r real::amean,bmean!a,b向量的平均值!計(jì)算平均值doi=1,ns=0.nlev=1doi=1,ns=0.s=s+a(i)open(7,file=”E:\grads\20230706qimo\test2s=s+a(i)amean=s/nenddo\corr.7.txt”)amean=s/nenddodoi=1,ns=0.\corr.7.grd”,form=”binary”)dodoi=1,ns=0.bmean=s/nenddos=s+b(i)write(7,*)stid(i),lon(i),lat(i),tim,nlev,nflag,r(i)endbmean=s/nenddos=s+b(i)doi=1,160r(i)write(8)stid(i),lat(i),lon(i),tim,nlev,nflag, !計(jì)算相關(guān)系數(shù)r(i)enddonlev=0write(8)stid(160-1),lat(160-1),lon(160-1),tim,nlev,nflag,r(160-1)endfunctionrelation(a,b,n)!本程序計(jì)算兩列向量的相關(guān)系數(shù)!a,b分別是待計(jì)算的向量!n是向量的長度，要求兩列向量等長implicitnoneinteger,intent(in)::nreal,intent(in)::a(n),b(n)corr.7.grd的數(shù)據(jù)描述文件2.ctl；dsetE:\grads\20230706qimo\test2\corr.7.datdtypestation

sfenzi=0.sfenmu1=0.sfenmu2=0.doi=1,nsfenzi=sfenzi+(a(i)-amean)*(b(i)-bmean)sfenmu1=sfenmu1+(a(i)-amean)**2sfenmu2=sfenmu2+(b(i)-bmean)**2enddorelation=sfenzi/sqrt(sfenmu1*sfenmu2)endfunctionrelationstnmapE:\grads\20230706qimo\test2\corr7.mapUNDEF2.5E33TITLEcorrTDEF1linearjan19511moVARS1r099relationendvarscorr.7.grdcorr7.map；預(yù)備GrADSGrid.grdGrid161.ctl〔已供給〕；2.gs1951-2023年1Nino347160站降水相關(guān)系數(shù)corr.7.gmf文件中。”reinit” E:\grads\20230706qimo\test2\grid161.ctl”‘open ”openE:\grads\20230706qimo\test2\2.ctl””setlat354””setlon70135””definea=oacres(g,r.2)””definea1=maskout(a,g-0.5)””defineaa=smth9(a1)””setgradsoff””setmpdsetworld””setmap111””enable printE:\grads\20230706qimo\test1\corr7.gmf”

”setgxoutshaded””setcmin0””daa””setgxoutcontour””setccolor2””daa””drawtitleLuYu””print””disableprint”分析：

1951-2023年1Nino347160站降水相關(guān)系數(shù)圖11951-2023年1Nino34海溫指數(shù)與海平面氣壓相關(guān)系數(shù)圖，從圖中可以看出，赤道中東太平洋地區(qū)海溫指數(shù)和海平面氣壓成明顯正相關(guān)，最大相關(guān)系數(shù)已超過0.8；而赤道錫太平洋地區(qū)海溫指數(shù)和海平面氣壓則成負(fù)相關(guān)，最小相關(guān)系數(shù)低于-0.6。由此可以解釋當(dāng)發(fā)生厄爾尼諾大事時，由于赤道中東太平洋和秘魯沿岸等地區(qū)海溫上升，使沃克環(huán)流減弱，對中東太平洋地區(qū)的下沉氣流產(chǎn)生負(fù)的影響，使中東太平洋地區(qū)海平面氣壓降低；與之相反，西太平洋地區(qū)上升