- Timestamp:
- 09/28/05 14:23:30 (3 years ago)
- Files:
-
- 1 modified
-
MarsGCM/+source/compile/MG8M0F1.f (modified) (25 diffs)
Legend:
- Unmodified
- Added
- Removed
-
MarsGCM/+source/compile/MG8M0F1.f
r150 r158 1 c ********************************************************************* 2 c ********************************************************************* 3 c ** 4 c ** Model II Mars 5 c ** Based on GCMII code for Mars from Noam on Athena /u/noam/+NEWMARS 6 c ** which was in turn from David's 8x10x23 model II. 7 c ** 8 c ** CHANGE HISTORY: 9 c ** 10 c ** 09/28/05 fixed includes, reordered a little code 11 c ** 12 c ** NOTES: 13 c ** 14 c ********************************************************************* 15 c ********************************************************************* 16 1 17 C**** MG8M0F1 BM08M0M23 MG8M0T noam 3/12/99 0.1 2 18 C**** 0.2 … … 22 38 C**** CALL FINITE DIFFERENCE EPFLUX DIAGNOSTICS 0.44 23 39 C**** UMA,VMA,PAA,TAA,PAVG SAVED ON UNIT 20 EVERY USET HOURS 0.5 24 INCLUDE ' ./BM08M0M23.COM'1.40 INCLUDE 'BM08M0M23.COM' 1. 25 41 COMMON U,V,T,SX,SY,SZ,P,Q 2. 26 42 COMMON/WORK1/WORKX(IM0,JM0,69),UA(IM0,JM0,LM0),VA(IM0,JM0,LM0), 3. … … 412 428 close (30) 336.61 413 429 c**** try using UNIX copy ('cp') to save the uset file (fort.20) 336.62 414 cmm d = 'cp fort.20 '//336.63430 cmmnd = 'cp fort.20 '// 336.63 415 431 * jmnth0(1:3)//cyear//'.use'//label1(1:llab1) 336.64 416 432 IF(USET.GT.0.) CALL SYSTEM(CMMND) 336.87 … … 470 486 C**** INITIAL CONDITIONS, AND CALCULATES THE DISTANCE PROJECTION ARRAYS 504. 471 487 C**** 505. 472 INCLUDE ' ./BM08M0M23.COM'506.488 INCLUDE 'BM08M0M23.COM' 506. 473 489 COMMON U,V,T,SX,SY,SZ,P,Q 507. 474 490 CHARACTER*4 C,C1,NAMD60,DISK,RUNID 508. … … 1080 1096 C**** DETERMINED BY U, V AND P 1005. 1081 1097 C**** 1006. 1082 INCLUDE ' ./BM08M0M23.COM'1007.1098 INCLUDE 'BM08M0M23.COM' 1007. 1083 1099 COMMON/WORK1/PIT(IM0,JM0),SD(IM0,JM0,LM0-1),PU(IM0,JM0,LM0), 1008. 1084 1100 * PV(IM0,JM0,LM0) 1008.1 … … 1196 1212 C**** DETERMINED BY DT1 AND THE CURRENT AIR MASS FLUXES 1104. 1197 1213 C**** 1105. 1198 INCLUDE ' ./BM08M0M23.COM'1106.1214 INCLUDE 'BM08M0M23.COM' 1106. 1199 1215 COMMON/WORK1/PIT(IM0,JM0),SD(IM0,JM0,LM0-1),PU(IM0,JM0,LM0), 1107. 1200 1216 * PV(IM0,JM0,LM0) 1107.01 … … 1242 1258 C**** AS DETERMINED BY DT1 AND THE CURRENT AIR MASS FLUXES 1204. 1243 1259 C**** 1205. 1244 INCLUDE ' ./BM08M0M23.COM'1206.1260 INCLUDE 'BM08M0M23.COM' 1206. 1245 1261 COMMON/WORK1/PIT(IM0,JM0),SD(IM0,JM0,LM0-1),PU(IM0,JM0,LM0), 1207. 1246 1262 * PV(IM0,JM0,LM0) 1207.1 … … 1403 1419 C**** THE PRESSURE GRADIENT FORCE 1404. 1404 1420 C**** 1405. 1405 REAL*8 KAPAP1,KAPAP2 1405.5 1406 INCLUDE './BM08M0M23.COM' 1406. 1421 C REAL*8 KAPAP1,KAPAP2 1405.5 1422 INCLUDE 'BM08M0M23.COM' 1406. 1423 REAL*8 KAPAP1,KAPAP2 1406.5 1407 1424 COMMON/WORK3/PHI(IM0,JM0,LM0),SPA(IM0,JM0,LM0) 1407. 1408 1425 COMMON/WORKR/PU(IM0,JM0) 1408. … … 1575 1592 C**** CURRENT AIR MASS FLUXES 1804. 1576 1593 C**** 1805. 1577 INCLUDE ' ./BM08M0M23.COM'1806.1594 INCLUDE 'BM08M0M23.COM' 1806. 1578 1595 COMMON/WORK1/PIT(IM0,JM0),SD(IM0,JM0,LM0-1),PU(IM0,JM0,LM0), 1807. 1579 1596 * PV(IM0,JM0,LM0) 1807.1 … … 1692 1709 C**** THIS VERSION OF AVRX DOES SO BY TRUNCATING THE FOURIER SERIES. 1905. 1693 1710 C**** 1906. 1694 INCLUDE ' ./BM08M0M23.COM'1907.1711 INCLUDE 'BM08M0M23.COM' 1907. 1695 1712 DIMENSION PU(IM0,JM0) 1908. 1696 1713 COMMON/AVRXX/AVCOS(72,36),AVSIN(72,36),BYSN(36),CN(2,37) 1909. … … 1734 1751 C**** HUMIDITY CAUSED BY CONDENSATION. 2004. 1735 1752 C**** 2005. 1736 INCLUDE ' ./BM08M0M23.COM'2006.1753 INCLUDE 'BM08M0M23.COM' 2006. 1737 1754 COMMON U,V,T,SX,SY,SZ,P,Q 2007. 1738 1755 COMMON/WORK1/CONV(IM0,JM0,LM0),PK(IM0,JM0,LM0),PREC(IM0,JM0), 2008. … … 2303 2320 C**** RUN1 IS NOT ACUMULATED IN ADAILY FOR DIAG6 3006. 2304 2321 C**** 3007. 2305 INCLUDE ' ./BM08M0M23.COM'3008.2322 INCLUDE 'BM08M0M23.COM' 3008. 2306 2323 COMMON U,V,T,SX,SY,SZ,P,Q 3009. 2307 2324 COMMON/WORK1/CONV(IM0,JM0,LM0),PK(IM0,JM0,LM0),PREC(IM0,JM0), 3010. … … 3043 3060 C**** THIS SUBROUTINES ADDS THE RADIATION HEATING TO THE TEMPERATURES 4003. 3044 3061 C**** 4004. 3045 INCLUDE ' ./BM08M0M23.COM'4005.3062 INCLUDE 'BM08M0M23.COM' 4005. 3046 3063 COMMON U,V,T,SX,SY,SZ,P,Q 4006. 3047 3064 COMMON/WORK1/CONV(IM0,JM0,LM0),PK(IM0,JM0,LM0),PREC(IM0,JM0), 4007. … … 3628 3645 C**** SURFACE SPECIFIC HUMIDITY, AND SURFACE WIND COMPONENTS. 4806. 3629 3646 C**** 4807. 3630 INCLUDE ' ./BM08M0M23.COM'4808.3647 INCLUDE 'BM08M0M23.COM' 4808. 3631 3648 COMMON U,V,T,SX,SY,SZ,P,Q 4809. 3632 3649 COMMON/WORK1/CONV(IM0,JM0,LM0),PK(IM0,JM0,LM0),PREC(IM0,JM0), 4810. … … 4378 4395 C**** GROUND TEMPERATURE, GROUND WATER AND ICE, AND SNOW MELTING. 6004. 4379 4396 C**** 6005. 4380 INCLUDE ' ./BM08M0M23.COM'6006.4397 INCLUDE 'BM08M0M23.COM' 6006. 4381 4398 COMMON U,V,T,SX,SY,SZ,P,Q 6007. 4382 4399 COMMON/WORK3/E0(IM0,JM0,4),E1(IM0,JM0,4),EVAPOR(IM0,JM0,4) 6008. … … 5039 5056 C**** THIS ROUTINE ONLY CHECKS LAYERS 2 TO LM. 6805. 5040 5057 C**** 6806. 5041 INCLUDE ' ./BM08M0M23.COM'6807.5058 INCLUDE 'BM08M0M23.COM' 6807. 5042 5059 COMMON U,V,T,SX,SY,SZ,P,Q 6808. 5043 5060 COMMON/WORK1/CONV(IM0,JM0,LM0),PK(IM0,JM0,LM0) 6809. … … 5187 5204 C**** THE MESOSPHERE. 7004. 5188 5205 C**** 7005. 5189 INCLUDE ' ./BM08M0M23.COM'7006.5206 INCLUDE 'BM08M0M23.COM' 7006. 5190 5207 LOGICAL QWRITE,QTELLE 7007. 5191 5208 LOGICAL INFINT, INFIN4 7007.5 … … 5814 5831 C**** 3 SMOOTH P AND T 7508. 5815 5832 C**** 7509. 5816 INCLUDE ' ./BM08M0M23.COM'7510.5833 INCLUDE 'BM08M0M23.COM' 7510. 5817 5834 COMMON U,V,T,SX,SY,SZ,P,Q 7511. 5818 5835 REAL*8 X,XS,Y 7512. … … 5885 5902 C**** TAKE PLACE AT THE BEGINNING OF A NEW DAY. 8004. 5886 5903 C**** 8005. 5887 INCLUDE ' ./BM08M0M23.COM'8006.5904 INCLUDE 'BM08M0M23.COM' 8006. 5888 5905 COMMON U,V,T,SX,SY,SZ,P,Q 8007. 5889 5906 COMMON/WORK2/Z1OOLD(IM0,JM0),XO(IM0,JM0,3),XZO(IM0,JM0) 8008. … … 6170 6187 C**** ELIMINATES SMALL ICE BERGS WHEN ROICE < .01 AND TGO > 1 (C). 8506. 6171 6188 C**** 8507. 6172 INCLUDE ' ./BM08M0M23.COM'8508.6189 INCLUDE 'BM08M0M23.COM' 8508. 6173 6190 COMMON U,V,T,SX,SY,SZ,P,Q 8509. 6174 6191 COMMON/WORK2/Z1OOLD(IM0,JM0) 8510. … … 6269 6286 C**** THE ERRORS ARE CORRECTED. 9006. 6270 6287 C**** 9007. 6271 INCLUDE ' ./BM08M0M23.COM'9008.6288 INCLUDE 'BM08M0M23.COM' 9008. 6272 6289 COMMON U,V,T,SX,SY,SZ,P,Q 9008.001 6273 6290 COMMON/WORK1/WORKX(IM0,JM0,69),UA(IM0,JM0,LM0),VA(IM0,JM0,LM0), 9008.004 … … 6311 6328 C**** AND TO USE THE ALTERNATE TEMPERATURES 9105. 6312 6329 C**** 9106. 6313 INCLUDE ' ./BM08M0M23.COM'9107.6330 INCLUDE 'BM08M0M23.COM' 9107. 6314 6331 DIMENSION PA(IM0,JM0),SM(IM0,JM0,LM0),SXM(IM0,JM0,LM0), 9108. 6315 6332 * SYM(IM0,JM0,LM0),SZM(IM0,JM0,LM0),PB(IM0,JM0),S(IM0,JM0,LM0) 9109. … … 6491 6508 C**** ITS SLOPES IN THE EAST-WEST DIRECTION USING THE SLOPES SCHEME. 9304. 6492 6509 C**** 9305. 6493 INCLUDE ' ./BM08M0M23.COM'9306.6510 INCLUDE 'BM08M0M23.COM' 9306. 6494 6511 DIMENSION PA(IM0,JM0),SM(IM0,JM0,1),SXM(IM0,JM0,1),SYM(IM0,JM0,1),9307. 6495 6512 * SZM(IM0,JM0,LM0),PB(IM0,JM0) 9308. … … 6562 6579 C**** ITS SLOPES IN THE NORTH-SOUTH DIRECTION USING THE SLOPES SCHEME. 9504. 6563 6580 C**** 9505. 6564 INCLUDE ' ./BM08M0M23.COM'9506.6581 INCLUDE 'BM08M0M23.COM' 9506. 6565 6582 DIMENSION PA(IM0,JM0),SM(IM0,JM0,LM0),SXM(IM0,JM0,LM0), 9507. 6566 6583 * SYM(IM0,JM0,LM0),SZM(IM0,JM0,LM0),PB(IM0,JM0) 9508. … … 6676 6693 C**** IN THE Z DIRECTION USING THE SLOPES SCHEME. 9704. 6677 6694 C**** 9705. 6678 INCLUDE ' ./BM08M0M23.COM'9706.6695 INCLUDE 'BM08M0M23.COM' 9706. 6679 6696 DIMENSION PA(IM0,JM0),SM(IM0,JM0,LM0),SXM(IM0,JM0,LM0), 9707. 6680 6697 * SYM(IM0,JM0,LM0),SZM(IM0,JM0,LM0) 9708.
