Changeset 310 for trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/mymath.py

Timestamp:

Sep 29, 2011, 3:59:01 PM (13 years ago)

Author:

aslmd

Message:

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:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/mymath.py (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SRC/PYTHON/mymath.py

@@ -14 +14 @@
         return u'\u00b0'
 
+def writeascii ( tab, filename ):
+    mydata = tab
+    myfile = open(filename, 'w')
+    for line in mydata:
+        myfile.write(str(line) + '\n')
+    myfile.close()
+    return
+