羅敬

銀蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 81.3
散金: 19
紅花: 2
帖子: 43
在線: 18.3小時
蟲號: 2318513
注冊: 2013-03-04
性別: GG
專業(yè): 光學和光電子材料

[交流] 誰有美國標準大氣模型 1976 跪求呀已有3人參與

找了好久沒找到，希望小木蟲的大神可以幫幫忙！！感激不盡

回復此樓

» 猜你喜歡

地球物理學和空間物理學論文潤色/翻譯怎么收費? 已經(jīng)有203人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

請問美國哪個學校大氣污染這塊專業(yè)最好？已經(jīng)有18人回復
【求助】6S模型大氣校正求助已經(jīng)有13人回復

1樓 2014-07-13 21:03:43

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

mingkeng

至尊木蟲 (正式寫手)

應助: 11 (小學生)
金幣: 9386.4
帖子: 620
在線: 213.2小時
蟲號: 194003
注冊: 2006-02-22
專業(yè): 大氣科學

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

c  1976年美國標準大氣
   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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

mingkeng

至尊木蟲 (正式寫手)

應助: 11 (小學生)
金幣: 9386.4
帖子: 620
在線: 213.2小時
蟲號: 194003
注冊: 2006-02-22
專業(yè): 大氣科學

引用回帖:

2樓: Originally posted by mingkeng at 2014-07-26 19:26:23
c  1976年美國標準大氣
   PROGRAM STD
   DIMENSION HBASE(9)
   DATA HBASE/ 0.,1.1E4,2.E4,3.2E4,4.7E4,5.1E4,7.1E4,2*84852.0/
...

執(zhí)行這個程序，就可以算出美國標準大氣的參數(shù)，然后隨便找個軟件繪圖即可。

贊一下

回復此樓

3樓2014-07-26 19:26:55

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

水晶杯喝茶

新蟲 (正式寫手)

應助: 85 (初中生)
金幣: 762.3
散金: 454
紅花: 8
帖子: 301
在線: 460.4小時
蟲號: 2512282
注冊: 2013-06-18
專業(yè): 爆炸與沖擊動力學

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

2樓: Originally posted by mingkeng at 2014-07-26 19:26:23
c  1976年美國標準大氣
   PROGRAM STD
   DIMENSION HBASE(9)
   DATA HBASE/ 0.,1.1E4,2.E4,3.2E4,4.7E4,5.1E4,7.1E4,2*84852.0/
...

這用什么程序執(zhí)行啊

贊一下

回復此樓

4樓2014-07-26 20:38:17

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

mingkeng

至尊木蟲 (正式寫手)

應助: 11 (小學生)
金幣: 9386.4
帖子: 620
在線: 213.2小時
蟲號: 194003
注冊: 2006-02-22
專業(yè): 大氣科學

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

4樓: Originally posted by 水晶杯喝茶 at 2014-07-26 20:38:17
這用什么程序執(zhí)行啊...

這是Fortran程序

贊一下

回復此樓

5樓2014-07-27 08:18:38

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

愛笑的螃蟹

銅蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 337.7
帖子: 104
在線: 21.7小時
蟲號: 2020637
注冊: 2012-09-22
性別: GG
專業(yè): 計算機科學的基礎理論

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

請問是浙大光電帥男羅敬？

贊一下

回復此樓

你要盡心、盡性、盡意愛主你的神。

6樓2015-02-12 11:27:36

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 318求調(diào)劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 求調(diào)劑一志愿海大，0703化學學碩304分，有大創(chuàng)項目，四級已過 +3	幸運哩哩 2026-03-22	6/300	2026-03-22 14:12 by 30660438
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +5	生物工程調(diào)劑 2026-03-17	9/450	2026-03-21 23:32 by zhujy1982
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 材料學學碩080502 337求調(diào)劑-一志愿華中科技大學 +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 一志愿山大07化學 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 機械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 求調(diào)劑 +3	Ma_xt 2026-03-17	3/150	2026-03-21 02:05 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 考研調(diào)劑求學校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 藥學383 求調(diào)劑 +3	藥學chy 2026-03-15	5/250	2026-03-20 22:11 by 云游重陽
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 295復試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3	123阿標 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111