我想用底图创建一个密度图.我拥有的数据是无人值守的,重复的或彼此非常接近的.我已尝试对数据进行网格化,然后使用pcolor绘制bin的数量,但即使所有数据集的长度相等,我仍然会收到缓冲区大小错误.我最初的想法是使用下面的底图脚本,但我只能分散工作,虽然这不会给我一个密度图.
m = Basemap(resolution='f',projection='merc',
lon_0=160,
llcrnrlat=-30.0,
urcrnrlat=30.0,
llcrnrlon=100.,
urcrnrlon=270.0,
lat_ts=0.0)
m.drawmapboundary(fill_color='white')
m.fillcontinents(color='#F5DEB3',lake_color='#85A6D9')
m.drawcoastlines(color='black', linewidth=.4)
m.drawcountries(color='#6D5F47', linewidth=.4)
m.drawmeridians(np.arange(0, 360, 20),
labels=[0,0,0,1],
color='black',
dashes=[1,0],
labelstyle='+/-',
linewidth=0.2)
m.drawparallels(np.arange(-90, 90, 10),
labels=[1,0,0,0],
color='black',
dashes=[1,0],
labelstyle='+/-',
linewidth=0.2)
lats =
array([ 14.375, 14.125, 14.125, 9.375, 13.625, 14.375, 5.625,
13.875, 14.625, 5.875, 8.875, 5.625, 8.875, 13.375,
11.125, 8.375, 12.375, 6.125, 5.375, 8.375, 7.375,
7.875, 14.375, 14.875, 9.875, 11.125, 14.875, 7.875,
9.125, 11.625, 5.125, 10.875, 5.125, 12.125, 12.625,
8.625, 5.125, 8.375, 11.625, 11.375, 12.875, 14.375,
8.875, 8.375, 6.375, …我正在将一些F77文件转换为F90.但是,有一些常见的块未转换.我是F90的新手,没有F77的经验.有人可以告诉我如何将以下示例代码转换为F90吗？
BLOCK DATA SETUP
INTEGER A,B,C
REAL I,J,K,L
COMMON /AREA1/ A,B,C
COMMON /AREA2/ I,J,K,L
DATA A,B,C,I,J,K,L/0,1,2,10.0,-20.0,30.0,-40.0/
END
我的想法是将数组A,B和C放在一个模块中.我没有得到的主要是AREA1和AREA2.如何在F77中使用它们以及如何翻译它们？我的第一个猜测是丢弃它们,只需在模块中定义A,B和C. 但是,它们是一种包含A,B和C的派生类型吗？
我已将一个 Maven 项目导入 Intellij 并向该src目录添加一个模块。但现在我发现我无法像其他项目一样将 java 包添加到新目录中，即每当我右键单击新添加的目录时，Intellij 都不会给我添加 java 包的可能性。这似乎是一个设置问题，但我找不到控制此问题的设置。我使用的是 Intellij 的社区版：2018.3
我有一个问题,我的工作时间和时间差异很大.到目前为止,我已经使用许多if语句解决了这个问题,但这些都容易出错.在寻找更好的解决方案而不重新发明轮子时,我遇到了时间,日期和时间.但这些对我来说似乎太不灵活,所以我正在寻找有关如何使用python的内置函数解决这个问题的想法和提示 - 如果可能的话.这就是问题:将从条带到规则网格的卫星数据分箱涉及大量简单的时间计算.如果我有日期整数:200702100700我想从这样的时间格式中添加和减去小时和分钟,这样日期 - 18会给我一个考虑闰年的新日期,可以跨越,年,月,日向后和转发.这是我对python的datetime包的理解:
from datetime import timedelta, date, time
t1 = date(2007,01,01); t = time(12,0); t2 = datetime.combine(t1,t)
t3 = timedelta(hours=18); t4 = t2 - t3; t5 = t2 + t3
print t4
>>>datetime.timedelta(2006, 23401, 3001)
这个时间在用于进行数千次计算的脚本中不是很有用.我的问题是t4没有年,月或日对象,而2006似乎被视为小时.这一切都非常令人困惑.这可能在python/numpy/scipy中吗？如果是这样,有人可以给我一个提示吗？
我对 Web 端点进行了一个简单的测试，并且正在使用restAssured 测试 json 架构。我不断得到：
org.junit.jupiter.api.extension.ParameterResolutionException: No ParameterResolver registered for parameter [java.lang.String arg0] in method [void People.init(java.lang.String,java.lang.String)].
我读过JUnit 5手册，我真的不明白这个错误。据我所知，该init方法的输入字符串没有任何问题。我真的不明白这个错误。有人可以向我解释一下导致此错误的原因以及如何在这种情况下解决它吗？
import Utils;
import org.junit.jupiter.api.*;
import org.springframework.test.context.TestPropertySource;
import static io.restassured.RestAssured.get;
import static io.restassured.module.jsv.JsonSchemaValidator.matchesJsonSchemaInClasspath;
@DisplayName(value = "Tests endpoint")
@TestInstance(TestInstance.Lifecycle.PER_CLASS)
@TestMethodOrder(MethodOrderer.OrderAnnotation.class)
@TestPropertySource(locations = "classpath:resources")
class People {
private String PEOPLE_QUERY = "PeopleQStar.json";
private String host = "http://dev-dev/";
private String endPoint = "people";
Utils utils = new Utils();
@BeforeAll
void init(String host, String endPoint)
{
utils.setHostAndPath(host, endPoint);
utils.setCommonSettings();
}
@AfterAll
void cleanUp() …我正在将一些函数从Excel电子表格翻译成python.根据numpy docs:arctan2,我需要使用aranan2的atan2 .问题是两个结果甚至没有关闭:
oc = 23.4384863405
sal = 89.7814630647
sra = np.arctan2(np.cos(np.deg2rad(sal)),
np.cos(np.deg2rad(oc))*np.sin(np.deg2rad(sal)))
results: Excel = 1.566714757 Numpy = 0.00415720646 ??
我相信Excel结果是正确的.这是错误的numpy.
现在要么我没有正确使用arctan2,要么atan2在numpy中不是arctan2,或者numpy中有一个bug,或者我刚刚完全迷失在这里.
我正在使用python版本2.7.2和numpy 1.6.2
请问有什么想法吗？谢谢
在 Fortran 中，可以为指针分配内存，或者不能：
real(kind=jp), target :: bt(100,100)
real(kind=jp), pointer :: pt(:,:)
但是你可以为指针分配内存pt：
allocate(pt(100,100))
我的问题是：利弊是什么？就我所见，为指针分配内存违背了指针的目的并占用了更多内存。当然，我对指针的了解是有限的，所以如果有人能向我解释这里发生了什么，我将不胜感激。
我正在使用混合 FORTRAN 77 和 Fortran 90 代码的模型，此外，我正在使用英特尔编译器编译代码。