Changeset 310 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Timestamp:

Sep 29, 2011, 3:59:01 PM (14 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

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars

Files:

• STORM_JULIEN_LAST/aeropacity_tachemobile_z.F (modified) (1 diff)
• STORM_JULIEN_LAST/initracer.F (added)
• initracer.F (modified) (1 diff)
• initracer.F.changealphalift (deleted)
• initracer.F.link (added)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/STORM_JULIEN_LAST/aeropacity_tachemobile_z.F

@@ -497 +497 @@
       ztoploc = 10      ! target pseudo-altitude of local storm (km)
       radloc = 0.5      ! radius of dust storm (degree)
-      lonloc = 80       ! center longitude of storm (deg)
-      latloc = -25      ! center latitude of storm (deg)
+      lonloc = 25.      ! center longitude of storm (deg)
+      latloc = -3.      ! center latitude of storm (deg)
 
       DO ig=1,ngrid

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/initracer.F

@@ -1 @@
-link ../../../mars_lmd_new/libf/phymars/initracer.F
+link STORM_JULIEN_LAST/initracer.F