|
Last change
on this file since 315 was
310,
checked in by aslmd, 14 years ago
|
|
MESOSCALE:
-- Code in storm scenario corresponds to 'OMEGA' reference case [Julien Faure]
-- Easier settings for dust lifting without the need to recompile [see below]
-- A few modifications to plot and dust-devil detection PYTHON routines
29/09/11 == AS
--> To easily explore sensitivity to lifting thresholds: in dustlift.F, ustar_seuil=sqrt(stress/rho)
and alpha_lift[dust_mass] can be prescribed through an external stress.def parameter file.
--- alpha_lift[dust_number] is computed from alpha_lift[dust_mass] as in initracer.F
--- ustar_seuil is more user-friendly than stress because direct comparison with ustar from model
--> For the moment this is MESOSCALE only, but potentially useful to everyone
|
|
File size:
522 bytes
|
| Rev | Line | |
|---|
| [180] | 1 | def min (field,axis=None): |
|---|
| 2 | import numpy as np |
|---|
| 3 | return np.array(field).min(axis=axis) |
|---|
| 4 | |
|---|
| 5 | def max (field,axis=None): |
|---|
| 6 | import numpy as np |
|---|
| 7 | return np.array(field).max(axis=axis) |
|---|
| 8 | |
|---|
| 9 | def mean (field,axis=None): |
|---|
| 10 | import numpy as np |
|---|
| 11 | return np.array(field).mean(axis=axis) |
|---|
| 12 | |
|---|
| 13 | def deg (): |
|---|
| 14 | return u'\u00b0' |
|---|
| 15 | |
|---|
| [310] | 16 | def writeascii ( tab, filename ): |
|---|
| 17 | mydata = tab |
|---|
| 18 | myfile = open(filename, 'w') |
|---|
| 19 | for line in mydata: |
|---|
| 20 | myfile.write(str(line) + '\n') |
|---|
| 21 | myfile.close() |
|---|
| 22 | return |
|---|
| 23 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
