Crossin的编程教室

标题: 大家来找茬，一个微小的语法错误 [打印本页]
作者: 方自之    时间: 2017-6-28 18:55
标题: 大家来找茬，一个微小的语法错误
本帖最后由 方自之 于 2017-6-28 19:00 编辑

各位老铁、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

  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                             #流入角度

  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])

  27. pit_x = 1000
  28. pit_y = 1500
  29. P[pit_y,pit_x] = P_0                           #设置台风中心

  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

  54. #Print Projection
复制代码
作者: crossin先生    时间: 2017-6-28 22:45
你问代码错误的时候请把报错信息也附上。看报错是调试的第一步。

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

两层循环遍历没问题。

数组过大建议分段处理，或者配合循环动态处理，一般都可以通过代码来解决。
作者: 方自之    时间: 2017-6-29 08:13
crossin先生 发表于 2017-6-28 22:45
你问代码错误的时候请把报错信息也附上。看报错是调试的第一步。

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

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


微信截图_20170629080344.png (28.21 KB, 下载次数: 417)

本次报错（现已解决）

本次报错（现已解决）

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





欢迎光临 Crossin的编程教室 (https://bbs.crossincode.com/) Powered by Discuz! X2.5