source: trunk/mesoscale/LMD_MM_MARS/SRC/ARWpost/atlas.gw.def.ref @ 13

;;---------------------------------------
;
;  fill here plot parameters then type:
;
;  @trace.idl
;  -or- @report.idl
;  -or- @movie.idl
;
;  this file will be saved !
;
;;---------------------------------------


;-------------------------------------;
; DATA LOCATION and PLOT DESTINATION  ;
; default is ./                       ;
;-------------------------------------;
datafolder='/tmp7/aslmd/MERIDIANI_FLUSH'
plotfolder=datafolder
;; for web anim and tex report (default is 'plot')
textitle='atlas_2_1416_FLUSH'
texcomments='sol'

;;; load/save ARWpost-generated text files ? default is no
;backup_data='yes'
;already_data='yes'

;-------------------------------
;-------------------------------
; namelist.ARWPOST parameters
;-------------------------------
;-------------------------------

nam1=" fields = 'HGT,U,V,TSURF' "
;; ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^
;; CHANGE HERE FIELDS TO BE EXTRACTED 
;; ex: tk,HGT,PSFC,U,V,W,T,tk,ICETOT,RAVE,QH2O_ICE,QH2O,PTOT,TSURF

nam2=" start_date = '2024-09-07_12:00:00'
nam3=" end_date = '2024-09-08_18:00:00'
;nam2=" start_date = '2024-09-13_12:00:00'
;nam3=" end_date = '2024-09-13_18:00:00'
;; ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^
;; CHANGE HERE START/END DATE

;nam4=" interval_seconds = 3700 "
nam4=" interval_seconds = 7400 "
;nam4=" interval_seconds = 88800 "              
;; ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 
;; SET INTERVAL   2h: 7400s, 4h: 14800s, 6h: 22200s, 12h: 44400s        

nam5=" interp_method = 1 "      
nam5=" interp_method = 0 "
;; ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^              
;; SET VERT. COORD. 0: eta levels, 1: height (inc) or pressure (dec)

tabnam=[0.0, 5.0, 1.0]
tabnam=[999.,0.,1.,4.,6.,10.,20.]
;; ^ ^ ^ ^ ^ ^ ^ ^ ^
;; SET HEIGHT/PRESSURE LEVELS WITH [first alt, last alt, step]

;nest='d03'  
;; ^ ^ ^ ^
;; SET NEST



;----------------------------------------------------
;----------------------------------------------------
; what do you want to plot ?
; THIS MUST BE SET ACCORDINGLY TO PREVIOUS SECTION
;   - do not plot a field not extracted
;   - do not ask for more plot than data available
;----------------------------------------------------
;----------------------------------------------------

num=23                  ;; how many plots ?
;num=7
;num=1

field1='U'              ;; FILLED
;field2='TSURF'         ;; CONTOURED
winds=['U','V']         ;; VECTORS
topo=1

;;; VORTICITY
;field1='HGT' & winds=['U','V'] & topo=1

;----------------------------------
;----------------------------------
; what kind of plot do you want ? 
;----------------------------------
;----------------------------------

;; first is 1
plot='map' & level=2 ;1 
;plot='meridional' & nlon=50
;plot='zonal' & nlat=66

;---------;
; special ;
;---------;
;;; SEA-LEVEL PRESSURE
;field1='PSFC' & winds=['tk','HGT'] & topo=1 & plot='map' & level=2 ;6 ;2 ;7
;;; HORIZONTAL WIND VELOCITY
;field1='U' & winds=['U','V'] & topo=1 & plot='map' & level=2 ;6 ;2 ;9


;;;useidlfile='true' ;; if you are just changing what is below...


;----------------------------
;----------------------------
; change plot settings 
;----------------------------
;----------------------------

;; uncomment to modify default plot settings
;; | | | | | | | | | | | | | | | | | | | | |
;; v v v v v v v v v v v v v v v v v v v v v

;;;
;;; PLOT TITLES
;;;
;subtitle_user='-- LMD Martian Mesoscale Model --'            ;; grepALL
;title_user='Vertical Wind (!Nm!N.s!U-1!N)'     ;; grepALL
;title_user='Dust opacity' ;; grepALL
;title_axis=['x coord (km)','y coord (km)']     ;; grepALL

;;;
;;; COLOR TABLES
;;;
;flag_cb='false'                 ;; grepALL
format='(F5.1)'                  ;; grepALL
;format='(I0)'                   ;; grepALL
;colors=128                      ;; grepALL
pal=22                           ;; grepALL
;pal=33                           ;; grepALL
poscb=0.85                       ;; grepALL
poscb=0.68                       ;; grepALL
;; GOOD COLOR TABLES: 4, 18, 22, 16, 37, 33, 39

;;;
;;; FILL LIMITS
;;;
minfield_init= 0.    ;; grepALL
maxfield_init= 25.    ;; grepALL

lim_max = maxfield_init & w=where((what_I_plot ge lim_max) and (what_I_plot le 1e9)) & if (w[0] ne -1) then what_I_plot[w]=lim_max ;; grepALL
lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min ;; grepALL
;;lim_blank = 0.2 & w=where(abs(what_I_plot) le lim_blank) & if (w[0] ne -1) then what_I_plot[w]=missing_value ;; grepALL

;;;
;;; WINDS
;;;
windex=25.      ;; grepALL
stride=3.       ;; grepALL

;;;
;;; CONTOUR
;;;
;overcontour=overcontour-mean(overcontour)      ;; grepALL
lev=-10. + 0.5*findgen(100) ;; grepALL
;lev=120.+5.*findgen(50) ;; grepALL

;;;
;;; AXIS
;;;
;isotropic='false'       ;; grepALL
intervalx=2.            ;; grepALL
intervaly=2.            ;; grepALL

;;;
;;; MAP LIMITS 
;;;
windowx=[-20.,08.]       ;; grepMAP
windowy=[-16.,12.]       ;; grepMAP


;;;
;;; VERTICAL SECTION LIMITS
;;;
;minalt=-4.             ;; grepSEC
;maxalt=16.             ;; grepSEC
;minspace=-30.          ;; grepSEC
;maxspace=0.            ;; grepSEC

;;
;; METRIC UNITS FOR VERTICAL SECTION
;;
;factor=20.        ;; grepSEC
;space=space*60. & spacekm='true'       ;; grepSEC
;minspace=minspace*60./factor           ;; grepSEC
;maxspace=maxspace*60./factor           ;; grepSEC
;intervalx=round(intervalx*60./factor)   ;; grepSEC

;;;
;;; TRICKS
;;;

;; pour tracer juste les vecteurs sur un fond uni ou vide
;pal=0                                          ;; 1/4 grepMAP
;what_I_plot(*,*)=what_I_plot(*,*)*0.+0.2       ;; 2/4 grepMAP
;what_I_plot(0,0)=0.                            ;; 3/4 grepMAP
;flag_cb='false'                                        ;; 4/4 grepMAP

;;; truc pour tracer juste les contours et une zone grisee de topo 
;pal=0                                           ;; 1/9 grepALL
;w=where(abs(what_I_plot) lt missing_value)     ;; 2/9 grepALL
;what_I_plot[w]=0.                              ;; 3/9 grepALL
;w=where(abs(what_I_plot) gt missing_value)     ;; 4/9 grepALL
;what_I_plot[w]=0.75                            ;; 5/9 grepALL
;w=where(what_I_plot eq 0.)                     ;; 6/9 grepALL
;what_I_plot[w]=missing_value                           ;; 7/9 grepALL
;what_I_plot(0,0)=1. & what_I_plot(1,0)=0.      ;; 8/9 grepALL
;flag_cb='false'                                ;; 9/9 grepALL



;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;


;;; POUBELLE


;; DEFAULT: ARWPOST.EXE ALREADY RUN (i.e. input.ctl and input.dat already there)
;; UNCOMMENT FOLLOWING LINE TO EXTRACT FILES 
;extract='yes'   ;;   EXECUTE ARWPOST.EXE
                ;;   WITH THE FOLLOWING PARAMETERS :


;@map_sea_level.idl
;@map_ertel_vorticity.idl
;field1='HGT'
;field2='U'
;winds=['U','V']
;topo=1


;winds=['U','V'] & field1='SWDOWN' & topo=1
;winds=['V','W'] & field1='V' 
;field1='ICETOT' & topo=1
;field1='RAVE' & topo=1
;field1='QH2O_ICE' & topo=1
;field1='QH2O' & winds=['U','W']
;field1='QH2O' & field2='W'
;field1='QH2O' & winds=['U','W']
;field1='QH2O_ICE' & field2='tk'
;field1='tk' & topo=1
;field1='V' & field2='V'
;field1='W' & winds=['V','W']
;field1='tk' & winds=['V','W']
;field1='tk'
;field1='W' & topo=1
;field1='TSK' & topo=1
;field1='U' & winds=['U','V'] & topo=1
;winds=['V','W'] & field1='V'
;winds=['V','W'] & field1='W'
;field1='U' & winds=['U','V'] & topo=1
;field1='HGT' & field2='HGT' 
;field1='U' & winds=['U','V'] ;& topo=1
;field1='PSFC' & winds=['U','V'] & field2='W'
;field1='W' & winds=['U','V'] & field2='PSFC'
;field1='TSURF' & topo=1 ;& winds=['U','V']
;field1='T' & topo=1 & winds=['U','V']
;field1='PSFC' & winds=['tk','HGT'] & topo=1
;field1='tk' & topo=1
;field1='tk' & winds=['U','V'] & topo=1
;field1='T' & field2='T'
;field1='T' & winds=['U','V'] 
;field1='HGT' & winds=['HGT','HGT']