Changeset 117 for trunk/mesoscale/PLOT/SPEC

Timestamp:

May 10, 2011, 6:33:29 PM (14 years ago)

Author:

aslmd

Message:

M 116 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
* CORRECTION MAJEURE DE BUG : PSFC et TRACEURS PAS BIEN MIS A JOUR *

• pdq au lieu de pdq*dt (dt = pas de temps dynamique)
• idem pour pdpsrf

IMPACT SURTOUT SUR LES SIMULATIONS LONGUES ET LES SIMULATIONS TACHE DE POUSSIERE (stage Julien)
TEST A EFFECTUER POUR LES SIMULATIONS CYCLE DE L'EAU

M 116 mars/libf/phymars/newsedim.F
M 116 mars/README
Correction de deux bugs dans newsedim.F: 1-e(-x) trop faible et endif mal place

M 116 mesoscale/LMDZ.MARS.new/myGCM/start.nc
A 0 mesoscale/LMDZ.MARS.new/myGCM/DEFS_JB/SP51_HR_dq_r3.0n0.5a0.5_MY26_TM
A 0 mesoscale/LMDZ.MARS.new/myGCM/DEFS_JB/SP51_HR_dq_r3.0n0.5a0.5_MY26_TM/run.def
A 0 mesoscale/LMDZ.MARS.new/myGCM/DEFS_JB/SP51_HR_dq_r3.0n0.5a0.5_MY26_TM/traceur.def
A 0 mesoscale/LMDZ.MARS.new/myGCM/DEFS_JB/SP51_HR_dq_r3.0n0.5a0.5_MY26_TM/startfi72.nc
A 0 mesoscale/LMDZ.MARS.new/myGCM/DEFS_JB/SP51_HR_dq_r3.0n0.5a0.5_MY26_TM/start72.nc
A 0 mesoscale/LMDZ.MARS.new/myGCM/DEFS_JB/SP51_HR_dq_r3.0n0.5a0.5_MY26_TM/callphys.def
M 116 mesoscale/LMDZ.MARS.new/myGCM/traceur.def
M 116 mesoscale/LMDZ.MARS.new/myGCM/startfi.nc
M 116 mesoscale/LMDZ.MARS.new/myGCM/callphys.def
M 116 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/compile
M 116 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/launch_gcm
M 116 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/run_mcd_3days
Nouvelle base de donnees d'etats initiaux sans les nuages radiativement actifs

A 0 mesoscale/PLOT/SPEC/LES/turb_period_psfc.pro
A 0 mesoscale/PLOT/SPEC/LES/turb_inc.pro.new
M 116 mesoscale/PLOT/SPEC/GW/gravitwave2.pro
M 116 mesoscale/PLOT/SPEC/GW/gravitwaveprof.pro
Petites MAJ routines graphiques

Location:

trunk/mesoscale/PLOT/SPEC

Files:

• GW/gravitwave2.pro (modified) (1 diff)
• GW/gravitwaveprof.pro (modified) (4 diffs)
• LES/turb_inc.pro.new (added)
• LES/turb_period_psfc.pro (added)

trunk/mesoscale/PLOT/SPEC/GW/gravitwave2.pro

@@ -13 +13 @@
 folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_narrowmountain_2D.morepoints.160818/'
 folder='/d5/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_LT15.161544/'
+folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_LTcst.168440/'
+folder='/home/aslmd/GRAVITWAVE/GW_MARS_highwind_3D_loweropacity_morepoints_widehill.168468/'
 file=folder+'wrfout_d01_9999-09-09_09:00:00_z'
 charvar='W' & charvarc='W'

trunk/mesoscale/PLOT/SPEC/GW/gravitwaveprof.pro

@@ -9 +9 @@
         path='../TESTGW/'
         path='/tmp7/aslmd/'
-path='/d5/aslmd/GRAVITWAVE/'
+path='/home/aslmd/GRAVITWAVE/'
         experiment='mcd_lt16_rcpcst'
         ;experiment='mcd_lt16'
@@ -29 +29 @@
 experiment='GW_MARS_highwind_3D_loweropacity_morepoints.161542'
 experiment='GW_MARS_highwind_3D_loweropacity_morepoints_LT15.161544'
+experiment='GW_MARS_highwind_3D_loweropacity_morepoints_LTcst.168440'
+experiment='GW_MARS_highwind_3D_loweropacity_morepoints_widehill.168468'
         file=path+experiment+'_wrfout_d01_9999-09-09_09:00:00_z'
 file=path+experiment+'/wrfout_d01_9999-09-09_09:00:00_z'
@@ -88 +90 @@
         ;mean_columnp = TOTAL(columnp(*,*,*,nts:nte),4) / float(nte-nts+1)
         anomal_invar = invar & for i=0,s[4]-1 do anomal_invar(*,*,*,i) = anomal_invar(*,*,*,i) - mean_invar
-;
-;
-;
-        ;goto, perturb
+;;
+;;
+;;
+;       goto, perturb
 for nt=nts,nte do begin
 for nx=nxs,nxe do begin
@@ -257 +259 @@
         oplot, overplot, overplot_column
      endfor
+
+        ;temppp = reform(mean_invar(nxs,middle,*))
+        ;hache = 191.*temppp/3.72/1000.
+        ;ok = column / ( 2. * hache )
+        ;print, exp(ok)
+        ;yeah = exp(ok)
+        ;oplot, yeah, column
+
 PS_End, /PNG