大家来找茬，一个微小的语法错误

方自之

各位老铁、Crossin先生好：
        py新手，由于之前都是用fortran做数值模拟较多，有些习惯延续到python.我的问题如下：
1、如代码，我遍历数组p[3000,2000]依然是fortran的习惯是：两个嵌套循环。python中是否建议这样使用，是否有更规范便捷的方法？
2、如代码，我在两个for循环内的if语句，报语法错误，很捉急，没有看出来，请大家帮忙找茬,第35行。
3、由于数据较大，我起初想创建的数组维度是p[240,3000,2000],但报错memory error,自己查了（其实根据字面意思也能猜出来）是数组太大，我无奈改为p[3000,2000]。请问python中该如何解决？
4、代码很简单，都是一元函数运算，如果大家觉得在书写方法上有可以优化的地方，或者有些坏习惯，也请告诉我。
先谢谢各位了

1. # -*- coding: utf-8 -*-
   
2. #本程序用于实现理想台风模拟，参考防灾减灾手册。
   
3. import numpy as np
   
4. import matplotlib.pyplot as plt
   
5. import matplotlib.mlab as mlab
   
6. from mpl_toolkits.basemap import Basemap, cm
   
7. import math
   
8. 
   
9. #参数设定
   
10. R_max = 39.0                                 #最大风速半径（km）
    
11. pi    = 3.1415926535897932384626             #pai
    
12. P_max = 1010.0                               #无穷远处最大气压（hPa）
    
13. P_0   = 965.0                                #台风中心气压
    
14. V_0   = 25.0                                 #台风移速（km/h）
    
15. x_def = 2000
    
16. y_def = 3000                                 #区域范围（km）
    
17. t_def = 120                                  #模拟时间区间(h)
    
18. theta_1 = pi/18.0
    
19. theta_2 = pi/7.2                             #流入角度
    
20. 
    
21. W_r     = 3.029*(P_max-P_0)**0.644
    
22. #Surface Pressure(hPa)
    
23. P    = np.zeros([y_def,x_def])
    
24. U    = np.zeros([y_def,x_def])
    
25. V    = np.zeros([y_def,x_def])
    
26. 
    
27. pit_x = 1000
    
28. pit_y = 1500
    
29. P[pit_y,pit_x] = P_0                           #设置台风中心
    
30. 
    
31. for i in range(x_def):
    
32.          for j in range(y_def):
    
33.                   r=math.sqrt(abs(i-pit_x)**2+abs(j-pit_y)**2)
    
34.                   theta_3 = pi/12.0*((r-R_max)/(0.2*R_max)+(pi/18)
    
35.                   if r>0 and r<= R_max:                                      
    
36.                            A      = -((i-pit_x)*math.sin(theta_1)+(j-pit_y)*math.cos(theta_1))
    
37.                            B      =(i-pit_x)*math.cos(theta_1)-(j-pit_y)*math.sin(theta_1)
    
38.                            P[j,i] = P_0+(0.25*(P_max-P_0)*(r/R_max)**3)
    
39.                            U[j,i] = (R_max/(r-R_max)0+W_r*(r/R_max)**1.5*A/r
    
40.                            V[j,i] = (R_max/(r-R_max)*V_0)+W_r*(r/R_max)**0.5*B/r
    
41.                   elif(r>R_max and r<R_max*1.2):                         #theta_3                                             
    
42.                            P[j,i] =P_max-(0.75*(P_max-P_0)*(R_max/r))
    
43.                            A      = -((i-pit_x)*math.sin(theta_3)+(j-pit_y)*math.cos(theta_3))
    
44.                            B      =(i-pit_x)*math.cos(theta_3)-(j-pit_y)*math.sin(theta_3)
    
45.                            U[j,i] = (R_max/(r-R_max)0+W_r*(r/R_max)**1.5*A/r
    
46.                            V[j,i] = (R_max/(r-R_max)*V_0)+W_r*(r/R_max)**0.5*B/r
    
47.                   elif(r>=1.2*R_max):                                #theta_2                                    
    
48.                            P[j,i] =P_max-(0.75*(P_max-P_0)*(R_max/r))
    
49.                            A      = -((i-pit_x)*math.sin(theta_2)+(j-pit_y)*math.cos(theta_2))
    
50.                            B      =(i-pit_x)*math.cos(theta_2)-(j-pit_y)*math.sin(theta_2)
    
51.                            U[j,i] = (R_max/(r-R_max)0+W_r*(r/R_max)**1.5*A/r
    
52.                            V[j,i] = (R_max/(r-R_max)*V_0)+W_r*(r/R_max)**0.5*B/r
    
53. 
    
54. #Print Projection

复制代码

crossin先生

你问代码错误的时候请把报错信息也附上。看报错是调试的第一步。

这里是因为34行你少了半个括号。

两层循环遍历没问题。

数组过大建议分段处理，或者配合循环动态处理，一般都可以通过代码来解决。

方自之

crossin先生 发表于 2017-6-28 22:45
你问代码错误的时候请把报错信息也附上。看报错是调试的第一步。

这里是因为34行你少了半个括号。

谢谢Crossin先生！以后我会注意提问的完整性，方便老司机带路。这次错误如图，不是那种准确的报错，所以直接口述了。现问题已解决。绝大多数数值计算问题都可以用numpy解决，我知道这只是python的冰山一角，也不是你的授课重点，但你课程中的遍历文件夹，搜索内容等教程对我十分有用，使我在处理大量数据时更加自动化。PS. 最近闲时在玩你推荐的jieba中文分词器。谢谢你提供的所有平台和资源，我很受用。

crossin先生

方自之 发表于 2017-6-29 08:13
谢谢Crossin先生！以后我会注意提问的完整性，方便老司机带路。这次错误如图，不是那种准确的报错，所以 ...