source: trunk/MESOSCALE_DEV/PLOT/RESERVE/obsolete/trace.idl @ 207

Last change on this file since 207 was 86, checked in by aslmd, 14 years ago

*
mars + LMD_MM_MARS
* Precompilation flag MESOSCALE for better transparency

* in shared phymars between GCM and mesoscale model

*

M 85 mars/libf/phymars/meso_physiq.F
M 85 mars/libf/phymars/meso_inifis.F
Added a pre-compilation flag MESOSCALE so that the LMDZ.MARS GCM
will compile without stating errors because of mesoscale routines.

M 85 mars/libf/phymars/newcondens.F
M 85 mars/libf/phymars/testphys1d.F
M 85 mars/libf/phymars/dustlift.F
D 85 mars/libf/phymars/meso_testphys1d.F
D 85 mars/libf/phymars/meso_dustlift.F
D 85 mars/libf/phymars/meso_newcondens.F
Now, this MESOSCALE precompilation flag can be used to lower
the number of meso_* routines when adaptations for mesoscale
applications are not very extended.
--> Three meso_* routines were deleted and changes are
now impacted under the MESOSCALE flag in the original GCM routines
--> Completely transparent for GCM compilation since it is devoid
of the -DMESOSCALE option
--> Very good for syncing because changes in dustlift, newcondens
will be directly available in the mesoscale model

M 84 mesoscale/LMD_MM_MARS/makemeso
Changed meso_testphys1d in testphys1d

M 84 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new/makegcm_pgf
M 84 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new/makegcm_mpifort
M 84 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new/makegcm_ifort
M 84 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new/makegcm_g95
M 84 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new/makegcm_mpi
Added the option -DMESOSCALE in these scripts

*
LMD_MM_MARS
* Various minor changes related to water cycle and plotting routines

* Also included the GW test case

*

A 0 mesoscale/LMDZ.MARS.new/myGCM/DEFS_JB/callphys.def.orig
M 84 mesoscale/NOTES.txt
D 84 mesoscale/LMD_MM_MARS/SRC/ARWpost/idl
M 84 mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
M 84 mesoscale/LMD_MM_MARS/SIMU/gnome_launch.meso
M 85 mesoscale/PLOT/MINIMAL/map_latlon.pro
D 85 mesoscale/PLOT/SPEC/LES/getget.pro
M 85 mesoscale/PLOT/SPEC/MAP/map_uvt.pro
A + - mesoscale/PLOT/SPEC/getget.pro
A 0 mesoscale/PLOT/RESERVE/obsolete
A 0 mesoscale/TESTS/TESTGW.tar.gz
M 84 000-USERS

File size: 2.6 KB
Line 
1;;---------------------------------------
2;  TRACE
3;
4;  **** DO NOT MODIFY ****
5;
6;  set parameters in plot.def then type:
7;
8;  @trace.idl
9;
10;;---------------------------------------
11
12field1='' & field2='' & winds=''
13backup_data='no' & already_data='no'   
14datafolder='./' & plotfolder='./'
15topo=0
16extract='no'
17
18SPAWN, '\rm user.idl'
19SPAWN, 'grep -v grep plot.def > user.idl'
20
21;; to avoid dumb and useless warnings from IDL
22space=0. & altitude=0. & lon=0. & lat=0
23what_I_plot=0. & title_user='' & missing_value=0.
24
25@user.idl
26if (nam5 eq " interp_method = 0 ") then coord='model_level'
27if ((nam5 eq " interp_method = 1 ") and (tabnam(0) lt tabnam(1))) then coord='height'
28
29.compile call_arwpost
30if (extract eq 'yes') then call_arwpost, nam1, nam2, nam3, nam4, nam5, tabnam
31
32denom=plot+'_'+coord+'_'+field1
33if ((field2 ne '') and (topo eq 0)) then denom=denom+'_'+field2
34if (topo eq 1) then denom=denom+'_HGT'
35if (winds(0) ne '') then denom=denom+'_'+winds(0)+winds(1)
36
37if (already_data eq 'yes') then SPAWN, 'ln -sf  '+datafolder+'/'+coord+'.input.ctl input.ctl'
38if (already_data eq 'yes') then SPAWN, 'ln -sf  '+datafolder+'/'+coord+'.input.dat input.dat'
39
40SPAWN, '\rm param_plot.idl'
41SPAWN, 'grep grepALL plot.def > param_plot.idl'
42if ((plot eq 'meridional') or (plot eq 'zonal')) then SPAWN, 'grep grepSEC plot.def >> param_plot.idl'
43if (plot eq 'map') then SPAWN, 'grep grepMAP plot.def >> param_plot.idl'
44.compile map_latlon
45.compile section
46
47!p.charthick = 2.0
48!p.thick = 3.0
49!x.thick = 2.0
50!y.thick = 2.0
51
52out_wrf, $
53        plot=plot, $
54        field1=field1, $
55        field2=field2, $
56        when=1, $
57        winds=winds, $
58        level=level,$
59        save_data=yeah,$        ;; data is loaded only once
60        topo=topo,$
61;        range=range,$
62        colors=colors, $
63        nlon=nlon, $
64        nlat=nlat, $
65        save_ps=denom+string(1001,'(I0)')
66;
67; display the first one as a sample
68;
69SPAWN, 'display ./'+denom+'1001.ps &'
70
71if (num ge 2) then for i=2,num do $
72        out_wrf, $
73        plot=plot, $
74        field1=field1, $
75        field2=field2, $
76        when=i, $
77        winds=winds, $
78        level=level,$
79        save_data=yeah,$        ;; data is loaded only once
80        topo=topo,$
81        range=range,$
82        colors=colors, $
83        nlon=nlon, $
84        nlat=nlat, $
85        save_ps=denom+string(1000+i,'(I0)')
86
87SPAWN, 'mv -f ./'+denom+'????.ps '+plotfolder+'/'
88if (backup_data eq 'yes') then SPAWN, 'cp -f namelist.ARWpost '+datafolder+'/'+coord+'.namelist.ARWpost'
89if (backup_data eq 'yes') then SPAWN, 'cp -f input.ctl '+datafolder+'/'+coord+'.input.ctl'
90if (backup_data eq 'yes') then SPAWN, 'cp -f input.dat '+datafolder+'/'+coord+'.input.dat'
91
92SPAWN, 'cp -f plot.def '+plotfolder+'/'+denom+'.def'
93
Note: See TracBrowser for help on using the repository browser.