ocssw V2020
oci_kernal Namespace Reference

Functions

def normalized (a, order=2, axis=-1)
 

Variables

int color_lvl = 8
 
 rgb = list(permutations(range(0,256,color_lvl),3))
 
 colors = np.array(sample(rgb,120))
 
list awl
 
int start = 200
 
int stop = 20000
 
 wl = np.copy(awl[8:-7])
 
 nwl = len(wl)
 
 x = np.arange(start,stop,1, )
 
 nx = len(x)
 
 lsz = np.linspace(np.log(start),np.log(stop),nwl)
 
 sz = np.exp(lsz)
 
 nsz = len(sz)
 
 I = np.identity(nsz, dtype=np.float64)
 
 H0 = -np.identity(nsz, dtype=np.float64)
 
 H1 = np.dot(H0.T,H0)
 
 H2 = H1.copy()
 
 H3 = np.zeros([nsz,nwl])
 
 sig = np.zeros(nsz)
 
 pdf = np.zeros((nx, nsz))
 
 spdf = np.sum(pdf,axis=1)
 
float m = 1.5+0.00001j
 
int nang = 1
 
int ng = 7
 
int std = 1
 
 of = np.linspace(-(ng-1)/2,(ng-1)/2, ng)
 
 gsz = np.zeros(ng)
 
 kbh = np.zeros((nwl, nsz), dtype=np.float64)
 
 kbtz = np.zeros((nwl, nsz), dtype=np.float64)
 
 kbhg = np.ones((ng, nwl, nsz), dtype=np.float64)
 
 kbht = np.ones((nwl, nsz, ng), dtype=np.float64)
 
 wnd = np.array(signal.gaussian(7, std=std))/(std * np.sqrt(2*np.pi))
 
 qb
 
 gs
 
float R = 0.5
 
int spect = np.ones(nwl)/5 + R*kbh[:,0]
 
int guess = np.ones(nsz)/10
 
 sd = np.zeros(nsz)
 
 aot = np.dot(kbh,sd)
 
 ax1
 
 ax2
 
 figsize
 
 xlabel
 
 ylabel
 
 title
 
 b
 
 True
 
 which
 
 axis
 
int rp0 = 1e-6
 
float rp1 = 0.001
 
float rp2 = 0.001
 
float rp3 = 0.001
 
float rp1n = rp1
 
float rp2n = rp2
 
float rp3n = rp3
 
 rikern = la.inv(np.dot(kbh.T,kbh) + rp0*I + rp1n*H1 + rp2n*H2 + rp3n*H3)
 
 psd = np.tensordot(np.dot(rikern,kbh.T),spect,(((1),(0))))
 
 rguess = np.tensordot(kbh,psd,(((1),(0))))
 
 color
 
 fin = tables.open_file('OCI2020084005000.L1B_PACE.nc', 'r')
 
 oci
 
string opath = "figures"
 
 roci = np.tensordot(kbh,psd,(((1),(0))))
 
string fname = opath + "/SZ_" + "{: .0e}".format(rp0) + "_" + "{: .0e}".format(rp1) + \
 

Function Documentation

◆ normalized()

def oci_kernal.normalized (   a,
  order = 2,
  axis = -1 
)

Definition at line 16 of file oci_kernal.py.

Variable Documentation

◆ aot

aot = np.dot(kbh,sd)

Definition at line 141 of file oci_kernal.py.

◆ awl

list awl
Initial value:
1 = [310, 315, 320, 325, 330, 335, 340, 345, 350, 355, 360, 365, 370, 375, 380, 385,
2  390, 395, 400, 405, 410, 415, 420, 425, 430, 435, 440, 445, 450, 455, 460, 465,
3  470, 475, 480, 485, 490, 495, 500, 505, 510, 515, 520, 525, 530, 535, 540, 545,
4  550, 555, 560, 565, 570, 575, 580, 585, 590, 595, 600, 605, 610, 615, 620, 625,
5  630, 635, 640, 645, 650, 655, 660, 665, 670, 675, 680, 685, 690, 695, 700, 705,
6  710, 715, 720, 725, 730, 735, 740, 745, 750, 755, 760, 765, 770, 775, 780, 785,
7  790, 795, 800, 805, 810, 815, 820, 825, 830, 835, 840, 845, 850, 855, 860, 865,
8  870, 875, 880, 885, 890, 895, 940, 1040, 1250, 1378, 1615, 2130, 2260 ]

Definition at line 26 of file oci_kernal.py.

◆ ax1

ax1

Definition at line 146 of file oci_kernal.py.

◆ ax2

ax2

Definition at line 146 of file oci_kernal.py.

◆ axis

axis

Definition at line 150 of file oci_kernal.py.

◆ b

b

Definition at line 150 of file oci_kernal.py.

◆ color

color

Definition at line 183 of file oci_kernal.py.

◆ color_lvl

int color_lvl = 8

Definition at line 21 of file oci_kernal.py.

◆ colors

float colors = np.array(sample(rgb,120))

Definition at line 23 of file oci_kernal.py.

◆ figsize

figsize

Definition at line 146 of file oci_kernal.py.

◆ fin

fin = tables.open_file('OCI2020084005000.L1B_PACE.nc', 'r')

Definition at line 191 of file oci_kernal.py.

◆ fname

string fname = opath + "/SZ_" + "{: .0e}".format(rp0) + "_" + "{: .0e}".format(rp1) + \

Definition at line 238 of file oci_kernal.py.

◆ gs

gs

Definition at line 115 of file oci_kernal.py.

◆ gsz

gsz = np.zeros(ng)

Definition at line 104 of file oci_kernal.py.

◆ guess

int guess = np.ones(nsz)/10

Definition at line 131 of file oci_kernal.py.

◆ H0

H0 = -np.identity(nsz, dtype=np.float64)

Definition at line 47 of file oci_kernal.py.

◆ H1

H1 = np.dot(H0.T,H0)

Definition at line 53 of file oci_kernal.py.

◆ H2

H2 = H1.copy()

Definition at line 55 of file oci_kernal.py.

◆ H3

H3 = np.zeros([nsz,nwl])

Definition at line 59 of file oci_kernal.py.

◆ I

I = np.identity(nsz, dtype=np.float64)

Definition at line 45 of file oci_kernal.py.

◆ kbh

def kbh = np.zeros((nwl, nsz), dtype=np.float64)

Definition at line 105 of file oci_kernal.py.

◆ kbhg

kbhg = np.ones((ng, nwl, nsz), dtype=np.float64)

Definition at line 107 of file oci_kernal.py.

◆ kbht

kbht = np.ones((nwl, nsz, ng), dtype=np.float64)

Definition at line 108 of file oci_kernal.py.

◆ kbtz

kbtz = np.zeros((nwl, nsz), dtype=np.float64)

Definition at line 106 of file oci_kernal.py.

◆ lsz

lsz = np.linspace(np.log(start),np.log(stop),nwl)

Definition at line 41 of file oci_kernal.py.

◆ m

float m = 1.5+0.00001j

Definition at line 90 of file oci_kernal.py.

◆ nang

int nang = 1

Definition at line 91 of file oci_kernal.py.

◆ ng

int ng = 7

Definition at line 101 of file oci_kernal.py.

◆ nsz

nsz = len(sz)

Definition at line 43 of file oci_kernal.py.

◆ nwl

nwl = len(wl)

Definition at line 38 of file oci_kernal.py.

◆ nx

nx = len(x)

Definition at line 40 of file oci_kernal.py.

◆ oci

oci
Initial value:
1 = np.append(fin.get_node("/", "observation_data/Lt_blue")[8:59], \
2  fin.get_node("/", "observation_data/Lt_red")[1:,:,:],axis=0)

Definition at line 192 of file oci_kernal.py.

◆ of

of = np.linspace(-(ng-1)/2,(ng-1)/2, ng)

Definition at line 103 of file oci_kernal.py.

◆ opath

string opath = "figures"

Definition at line 197 of file oci_kernal.py.

◆ pdf

pdf = np.zeros((nx, nsz))

Definition at line 80 of file oci_kernal.py.

◆ psd

psd = np.tensordot(np.dot(rikern,kbh.T),spect,(((1),(0))))

Definition at line 179 of file oci_kernal.py.

◆ qb

qb

Definition at line 115 of file oci_kernal.py.

◆ R

float R = 0.5

Definition at line 126 of file oci_kernal.py.

◆ rgb

rgb = list(permutations(range(0,256,color_lvl),3))

Definition at line 22 of file oci_kernal.py.

◆ rguess

rguess = np.tensordot(kbh,psd,(((1),(0))))

Definition at line 180 of file oci_kernal.py.

◆ rikern

rikern = la.inv(np.dot(kbh.T,kbh) + rp0*I + rp1n*H1 + rp2n*H2 + rp3n*H3)

Definition at line 176 of file oci_kernal.py.

◆ roci

roci = np.tensordot(kbh,psd,(((1),(0))))

Definition at line 208 of file oci_kernal.py.

◆ rp0

int rp0 = 1e-6

Definition at line 154 of file oci_kernal.py.

◆ rp1

int rp1 = 0.001

Definition at line 155 of file oci_kernal.py.

◆ rp1n

float rp1n = rp1

Definition at line 160 of file oci_kernal.py.

◆ rp2

int rp2 = 0.001

Definition at line 156 of file oci_kernal.py.

◆ rp2n

float rp2n = rp2

Definition at line 161 of file oci_kernal.py.

◆ rp3

int rp3 = 0.001

Definition at line 157 of file oci_kernal.py.

◆ rp3n

float rp3n = rp3

Definition at line 162 of file oci_kernal.py.

◆ sd

sd = np.zeros(nsz)

Definition at line 134 of file oci_kernal.py.

◆ sig

sig = np.zeros(nsz)

Definition at line 74 of file oci_kernal.py.

◆ spdf

spdf = np.sum(pdf,axis=1)

Definition at line 84 of file oci_kernal.py.

◆ spect

spect = np.ones(nwl)/5 + R*kbh[:,0]

Definition at line 127 of file oci_kernal.py.

◆ start

int start = 200

Definition at line 35 of file oci_kernal.py.

◆ std

int std = 1

Definition at line 102 of file oci_kernal.py.

◆ stop

int stop = 20000

Definition at line 36 of file oci_kernal.py.

◆ sz

sz = np.exp(lsz)

Definition at line 42 of file oci_kernal.py.

◆ title

title

Definition at line 149 of file oci_kernal.py.

◆ True

True

Definition at line 150 of file oci_kernal.py.

◆ which

which

Definition at line 150 of file oci_kernal.py.

◆ wl

wl = np.copy(awl[8:-7])

Definition at line 37 of file oci_kernal.py.

◆ wnd

wnd = np.array(signal.gaussian(7, std=std))/(std * np.sqrt(2*np.pi))

Definition at line 109 of file oci_kernal.py.

◆ x

x = np.arange(start,stop,1, )

Definition at line 39 of file oci_kernal.py.

◆ xlabel

xlabel

Definition at line 148 of file oci_kernal.py.

◆ ylabel

ylabel

Definition at line 148 of file oci_kernal.py.