分类： 学习

https://en.wikipedia.org/wiki/Colatitude

In spherical coordinates, colatitude is the complementary angle of the latitude, i.e. the difference between 90° and the latitude,^[1] where southern latitudes are denoted with a minus sign.

Examples

Latitude and colatitude sum up to 90°.

place                latitude               colatitude

north pole​        90°                           0°

equator              ​0°                           90° 

south pole        -90°                       180°

定义了一个结构体数组，每个成员都想从另一个相同的结构体中拷贝过来，遂用了memcpy，比如

typedef struct cata *CATA, single_cata;

int size=​sizeof(single_cata);

memcpy(CATA+i*size,&single_cata,size);​

运行时发现除了结构体数组CATA的第一个值是正确的，后面的全错了，后来搜了一下，发现要如是改：

​memcpy((char*)CATA+i*size,&single_cata,size);​

perl在进行汉字操作之前要先转码，不然无法正确操作。

​use Encode;

​$s=’汉字’;

$str=​decode ‘utf8’,$s;   #######或者gb2312之类的

$len=length($str); ######这样就可以正确得到汉字字符串的长度了​

参考http://www.cnblogs.com/tsingke/p/6194737.html

采用box-muller方法:

​

#include

#include ​

double gaussrand2(double mu,double sigma){​

  static double u,v;

  static int phase=0;

  double z;

  if(phase==0){

    u=rand()/(RAND_MAX+1.0);

    v=rand()/(RAND_MAX+1.0);

    z=sqrt(-2.0*log(u))*sin(2.0*pi*v);

  }

  else{

    z=sqrt(-2.0*log(u))*cos(2.0*pi*v);

  }

  phase=1-phase;

  return(mu+z*sigma);

}

1. fortran中的函数名全部用小写，记得带下划线，如func_

2. 参数必须是地址

3. 编译时加上-lgfortran​

import numpy as np

​import matplotlib.pyplot as plt 

 from mpl_toolkits.basemap import Basemap 

 from obspy import read_events 

from obspy.imaging.beachball import beach 

 event = read_events( 'https://earthquake.usgs.gov/archive/product/moment-tensor/' 'us_20005ysu_mww/us/1470868224040/quakeml.xml', format='QUAKEML')[0] 

origin = event.preferred_origin() or event.origins[0] 

ocmec = event.preferred_focal_mechanism() or event.focal_mechanisms[0] 

tensor = focmec.moment_tensor.tensor 

import time
a='1987-06-05 12:31:28'
tm=time.strptime(a,'%Y-%m-%d %H:%M:%S')
epo_time=time.mktime(tm)

import matplotlib.pyplot as plt
import matplotlib.dates as dates

ax.xaxis.set_minor_locator(dates.YearLocator())
ax.xaxis.set_major_locator(dates.YearLocator(5))

可以自己写个简单函数，也可以使用math模块下的modf函数
import math
a,b=math.modf(32.005)
print a,b

a=0.005000000000002558
b=32.0

*****注意a，b都是浮点型