AminFazlKazemi

Hi Dr. Huang This is Amin Fazl kazemi from I.R.Iran. I've just downloaded your description of SOR solving method to find UREF. It is published February 2020 and refers to...

Thanks for description. My PHD thesis topic is about studying diabatic/disssipative processes effects using your LWA conservation law. By the way, I have tried to solve equation 12 Nakamura/Solomon 2010...

> Hi Dr. Huang I used this code to solve UREF based on Nakamura_Solomon 2010: ` @numba.guvectorize(['(float64[:],int64,float64,float64[:,:],float64[:,:],float64[:],float64[:],float64[:],float64[:,:],int64,int64,float64,float64[:,:])'] , '(z),(),(),(m,z),(m,z),(m),(z),(z),(m,z),(),(),()->(m,z)',target='cpu', nopython=True) def density_weighted_elliptic_solver(A,max_iter,acceptable_precision,first_guess,epsilon,dmu2,dz2,ru0,forcing,ny,nz,alfa,result): dz=dz2**.5 dmu=dmu2**.5 precision=1000. iter_=0 old=first_guess.copy() for iter_ in...

Hi You are correct the units are squared meters. `Area=(numpy.radians(resolution)*earth_radius)**2 *(numpy.cos(numpy.radians(data.latitude)))+0*data.longitude` well , yes I used earth-radius as not complete sphere but the effect of this is partial. I checked...

Hi. for Latitude=75 degrees North negative_intrusion_area= 9.08332754e+12 squared meters positive_intrusion_area = 4.88416361e+11 squared meters

This is the code I used. Itried but didn't manage to add indentations correctly. The data can be downloaded from era interim dataset with 0.75 degrees resolution for April 2011...

[LWA.txt](https://github.com/miniufo/xcontour/files/8640338/LWA.txt) This is the code

In the finite amplitude limit, we see enclosed or cut-off contours forming. I suggest test definition of area within contours visually in such a case. i also suggest to read...