 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2254 | 回復(fù): 5 | |||
[交流]
誰有美國標(biāo)準(zhǔn)大氣模型 1976 跪求呀 已有3人參與
|
| 找了好久沒找到，希望小木蟲的大神可以幫幫忙！！感激不盡 |
回復(fù)此樓» 猜你喜歡
招聘博士研究生
已經(jīng)有45人回復(fù)
-
地球物理學(xué)和空間物理學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有164人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 請問美國哪個(gè)學(xué)校大氣污染這塊專業(yè)最好？
已經(jīng)有18人回復(fù)
- 【求助】6S模型大氣校正求助
已經(jīng)有13人回復(fù)
mingkeng
至尊木蟲 (正式寫手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 9386.4
- 帖子: 620
- 在線: 213.2小時(shí)
- 蟲號: 194003
- 注冊: 2006-02-22
- 專業(yè): 大氣科學(xué)
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
c 1976年美國標(biāo)準(zhǔn)大氣 PROGRAM STD DIMENSION HBASE(9) DATA HBASE/ 0.,1.1E4,2.E4,3.2E4,4.7E4,5.1E4,7.1E4,2*84852.0/ PRINT 1002 1002 FORMAT(1H1,1X,'STANDARD ATMOSPHERE LAYERS - ', 2 2X,' FROM US STANDARD ATMOSPHERE, 1976 ',/) PRINT 1008 PRINT 1009 1008 FORMAT(10X,' P Z T TGRAD 2 DENSITY') 1009 FORMAT(10X,' MB METERS FEET DEG C DEG/M 2 KG/M3') HMT=HBASE(1) CALL STDZ2P(P,HMT, 0, TCENT, DENS) DO1I= 2,8 PP=P PHMT=HMT PTCENT=TCENT PDENS=DENS HMT=HBASE(I) CALL STDZ2P(P,HMT, 0, TCENT, DENS) TGRAD=(TCENT-PTCENT)/(HMT-PHMT) HFT=PHMT/.3048 WRITE (6,1007)PP,PHMT,HFT,PTCENT,TGRAD,PDENS,PP 1007 FORMAT(10X,F9.4,2F12.2,F10.2,F10.5,F20.8,E16.8) 1 CONTINUE TGRAD=0. HFT=HMT/.3048 WRITE (6,1007)P,HMT,HFT,TCENT,TGRAD,DENS,P C C DO P2Z SUBROUTINE PRINT 1003 1003 FORMAT(1H1,1X,'PRINT FOR PRESSURE LEVELS') PRINT 1005 PRINT 1006 PINC=50. P=1000. 10 CALL STDP2Z(P,Z,ZFT,TCENT,DENS) PRINT 1004,P,Z,ZFT,TCENT,DENS 1004 FORMAT(10X,F7.2,2F10.0,F10.2,F16.6) 1005 FORMAT(10X,' P Z T DENSITY') 1006 FORMAT(10X,' MB METERS FEET DEG C KG/M3') IF(P.LE.100.) PINC=10. IF(P.LE.10.)PINC=1. IF(P.LE.1.)PINC=.1 P=P-PINC IF(P.GE..09) GO TO 10 END SUBROUTINE STDZ2P(P,H1,L,TCENT,RHO) C C COMPUTES PRESSURE, TEMPERATURE, AND DENSITY VALUES FROM INPUT HEIGHTS C BASED ON US STANDARD ATMOSPHERE, 1976 C VALUES NOT VALID ABOVE 84852 M. C C INPUT C H1 - HEIGHT IN FEET OR METERS C L - UNITS FLAG - 0=METERS, 1=FEET C OUTPUT C P - PRESSURE IN MB C TCENT - TEMPERATURE IN DEG C C RHO - DENSITY IN KG/M3 C DIMENSION HBASE(10),HTOP(10),TB(10),GRAD(10),ABS(10) DIMENSION PB(11),RB(10) DATA HBASE/ 0.,1.1E4,2.E4,3.2E4,4.7E4,5.1E4,7.1E4,3*84852.0/ DATA HTOP /1.1E4,2.E4,3.2E4,4.7E4,5.1E4,7.1E4,4*84852.0/ DATA TB / 15.,2*-56.5,-44.5,2*-2.5,-58.5,3*-86.2/ DATA GRAD/-.0065,0.,.001,.0028,0.,-.0028,-.002,3*0./ DATA CHECK/0./ C H1 IS ALTITUDE, IF IN METERS SET L .LE. 0. IF IN FEET SET L .GE. 1 C P IS PRESSURE IN MB. TCENT IS TEMP CELSIUS. RHO IS DENSITY IN KG/M**3 C IF(CHECK.NE.0.) GO TO 16 ABSZ = 273.15 PZERO=1013.250 G=980.665 R = 83143200./28.9644 PB(1) = PZERO DO 3 I=1,8 ABS(I) = TB(I) + ABSZ 3 CONTINUE DO 15 I=1,8 IF (GRAD(I) .EQ. 0) GO TO 10 C HERE FOR TEMP GRADIENT NOT ZERO 5 PB(I+1) = PB(I)*((ABS(I+1)/ABS(I))**(100.*G/(-GRAD(I)*R))) GO TO 12 C HERE FOR TEMP GRADIENT ZERO 10 PB(I+1) = PB(I)*2.7182818**((-100.*G/(R*ABS(I)))*(HTOP(I)-HBASE(I) X )) 12 RB(I) = 1000.*PB(I)/(R*ABS(I)) * 1000. 15 CONTINUE CHECK=10. 16 CONTINUE C IF(L. LT.1) GO TO 30 C H = H1*.3048 GO TO 4 30 H = H1 4 CONTINUE C FIND THE ATMOSPHERIC LAYER WE ARE IN 35 DO 55 I = 1,8 IF (H .GT. HBASE(I+1)) GO TO 55 37 IF (GRAD(I)) 39,43,39 C TEMP GRADIENT NOT ZERO 39 P= PB(I)*((H-HBASE(I))*GRAD(I)/ABS(I)+1)**(-100.*G/(GRAD(I)*R)) GO TO 45 C ISOTHERMAL LAYER 43 P= PB(I)*EXP((H-HBASE(I))*(-G)/(.01*R*ABS(I))) 45 T = ABS(I) + GRAD(I)*(H-HBASE(I)) TCENT = T - ABSZ RHO = 1.E6* P/(R*T) RETURN 55 CONTINUE 56 WRITE(6,98)H 98 FORMAT(' STDATM - ALTITUDE OUT OF RANGE ',F12.2,//) END SUBROUTINE STDP2Z(P,H,HFT,TCENT,DENS) C C COMPUTES HEIGHT, TEMPERATURE, AND DENSITY VALUES FROM INPUT PRESSURES C BASED ON US STANDARD ATMOSPHERE, 1976 C VALUES NOT VALID ABOVE 84852 KM. C C INPUT C P - PRESSURE IN MB C OUTPUT C H - HEIGHT IN METERS C HFT - HEIGHT IN FEET C TCENT - TEMPERATURE IN DEG C C DENS - DENSITY IN KG/M3 C DIMENSION HBASE(10),HTOP(10),TB(10),GRAD(10),ABS(10),PB(10),RB(10) DATA HBASE/ 0.,1.1E4,2.E4,3.2E4,4.7E4,5.1E4,7.1E4,3*84852.0/ DATA HTOP /1.1E4,2.E4,3.2E4,4.7E4,5.1E4,7.1E4,4*84852.0/ DATA GRAD/-.0065,0.,.001,.0028,0.,-.0028,-.002,3*0./ DATA ABS/288.15,216.65,216.65,228.65,270.65,270.65, 2 214.65,186.95,2*0./ DATA PB/1013.25,0.22632064E+03,0.54748887E+02,0.86801870E+01, 2 0.11090631E+01,0.66938874E+00,0.39564205E-01,0.37338360E-02, 3 2*0./ ABSZ = 273.15 G=980.665 R = 83143200./28.9644 IF(P .GT. PB(1)) GO TO 56 35 DO 55 I = 1,8 IF (P .GE. PB(I+1)) GO TO 37 55 CONTINUE 56 PRINT 1001,P 1001 FORMAT (' STDATM- PRESSURE OUT OF EXCEEDED ',F8.2) 37 IF (GRAD(I)) 39,43,39 C TEMP GRADIENT NOT ZERO 39 H = ABS(I)/(-GRAD(I)) *(1. -(P/PB(I))**(-GRAD(I)*R/(100.*G))) X + HBASE(I) GO TO 45 C ISOTHERMAL LAYER 43 H = HBASE(I) + .01*R*ABS(I)/G *(ALOG(PB(I)/P)) 45 T = ABS(I) + GRAD(I)*(H-HBASE(I)) TCENT = T - ABSZ DENS = 1000*1000.*P/(R*T) HFT = H/.3048 RETURN END |
2樓2014-07-26 19:26:23
mingkeng
至尊木蟲 (正式寫手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 9386.4
- 帖子: 620
- 在線: 213.2小時(shí)
- 蟲號: 194003
- 注冊: 2006-02-22
- 專業(yè): 大氣科學(xué)
3樓2014-07-26 19:26:55
4樓2014-07-26 20:38:17
mingkeng
至尊木蟲 (正式寫手)
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 9386.4
- 帖子: 620
- 在線: 213.2小時(shí)
- 蟲號: 194003
- 注冊: 2006-02-22
- 專業(yè): 大氣科學(xué)
5樓2014-07-27 08:18:38
愛笑的螃蟹
銅蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 337.7
- 帖子: 104
- 在線: 21.7小時(shí)
- 蟲號: 2020637
- 注冊: 2012-09-22
- 性別: GG
- 專業(yè): 計(jì)算機(jī)科學(xué)的基礎(chǔ)理論
6樓2015-02-12 11:27:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 羅敬 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +3 | NIFFFfff 2026-03-20 | 3/150 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 289求調(diào)劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|