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


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

;; 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 :

;-------------------------------
;-------------------------------

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

;nam2=" start_date = '2024-05-03_06:00:00' "	;; CHANGE HERE START DATE
;nam3=" end_date   = '2024-05-04_06:00:00' "	;; CHANGE HERE END DATE
;nam2=" start_date = '2024-06-53_06:00:00' "	;; CHANGE HERE START DATE
;nam3=" end_date = '2024-06-53_12:00:00' "	;; CHANGE HERE END DATE

nam2=" start_date = '2024-07-17_00:00:00' "	;; CHANGE HERE START DATE
nam3=" end_date = '2024-07-18_00:00:00' "	;; CHANGE HERE END DATE


;nam4=" interval_seconds = 14800 "		;; SET INTERVAL   2h: 7400s, 4h: 14800s, 6h: 22200s	
;nam4=" interval_seconds = 7400 "               ;; SET INTERVAL   2h: 7400s, 4h: 14800s, 6h: 22200s
;nam4=" interval_seconds = 3700 "
;nam4=" interval_seconds = 88800 "
nam4=" interval_seconds = 22200 "
;nam4=" interval_seconds = 44400 "
nam5=" interp_method = 0 "			;; SET VERT. COORD. 0: eta levels, 1: height (inc) or pressure (dec)
;nam5=" interp_method = 1 "
tabnam=[0.,6.,0.1]   	;; SET HEIGHT/PRESSURE LEVELS WITH [first alt, last alt, step]
;tabnam=[-3.,5.,0.1]
tabnam=[-5.,10.,0.2]

nest='d03'  ;; SET NEST
nest='d01'  ;; 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=1		   	;; how many plots ?
num=6                   ;; how many plots ?
;num=2                   ;; how many plots ?
;num=3			;; how many plots ?
;num=12
;num=13
;num=4
;num=7

;@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='W' ;& 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='W' & field2='U'
;field2='tk' & field1='W'
field1='tk'

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

plot='map' 	
level=1		; first level is 1
;level=20
;;level=41
;level=10
;level=11
;;level=1
;level=9
;level=8
;level=3
;level=49
;level=5
;level=15
;level=17
;level=5
;level=10
;level=7  ;;SEA LEVEL
;level=30
level=7

;plot='meridional'  	
;nlon=75	; first is 1
;;;;nlon=87
;nlon=72
;nlon=70
;nlon=69
;nlon=60  
;nlon=162
;nlon=100
nlon=60
nlon=150
nlon=250
nlon=30

;plot='zonal' 		
;;;nlat=145	; first is 1
;;nlat=24
;;nlat=25
;nlat=100.
;nlat=165

;@map_wind_stress.idl


;----------------------------
;----------------------------
; 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

minfield_init=200.              ;; grepALL
maxfield_init=260.              ;; grepALL
;minfield_init=0.              ;; grepALL
;maxfield_init=48.               ;; grepALL
;intervalx=5.                   ;; grepALL
;intervaly=5.                    ;; grepALL
;title_user='Vertical Wind (!Nm!N.s!U-1!N), LT 13, 980m ABG'	;; grepALL
;title_axis=['x coord (km)','y coord (km)']	;; grepALL
;subtitle_user='LT 13, 980m ABG, Contours: surface pressure anomaly each 0.5 Pa'	;; grepALL
;colors=64			;; grepALL
;pal=33				;; grepALL
;pal=0				;; grepALL
;pal=22				;; grepALL
;format='(I0)'			;; grepALL
;windex=30.			;; grepALL
;stride=3			;; grepALL
;flag_cb='false'		;; grepALL
;what_I_plot=abs(what_I_plot)    ;; grepALL ??! minfield
;what_I_plot=what_I_plot+220.
;isotropic='false'		;; grepALL
;overcontour=overcontour-mean(overcontour)	;; grepALL
;lev=[-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.]		;; grepALL
;lev=[-5.,0.,5.,10.,15.,20.,25.] ;; grepALL
;lev=[-1.,5.,10.,15.,20.,25.] ;; grepALL


;minalt=-3.			;; grepSEC
;maxalt=5.			;; grepSEC
;;;altitude=altitude*1000.  	;; grepSEC
;minspace=-18.0			;; grepSEC
;maxspace=-9.0			;; grepSEC
;minspace=-11.0                 ;; grepSEC
;maxspace=-7.0                  ;; grepSEC
;factor=20.  ;10.			   ;; grepSEC
;space=space*59. & spacekm='true'   ;; grepSEC
;minspace=minspace*59./factor       ;; grepSEC
;maxspace=maxspace*59./factor       ;; grepSEC
;intervalx=round(intervalx*59./factor)       ;; grepSEC


;windowy=[10.,13.]        ;; grepMAP
;windowx=[14.,17.]	;; grepMAP
;windowy=[-12.,-2.]        ;; grepMAP
;windowx=[-78.,-68.]     ;; grepMAP
;windowy=[45.,55.]        ;; grepMAP
;windowx=[-45.,-35.]     ;; grepMAP
;lon=lon*59.-lon(0)*59.		;; grepMAP ;;attention plots successifs 
;lat=lat*59.-lat(0)*59.		;; grepMAP
;windowy=[0.,14.]	;; grepMAP
;windowx=[0.,14.]   ;; grepMAP


;lim_max = maxfield_init				;; 1/3 grepALL
;w=where(what_I_plot ge lim_max) 		;; 2/3 grepALL
;if (w[0] ne -1) then what_I_plot[w]=lim_max	;; 3/3 grepALL

;lim_min = minfield_init                         ;; 1/3 grepALL
;w=where(what_I_plot le lim_min)                 ;; 2/3 grepALL
;if (w[0] ne -1) then what_I_plot[w]=lim_min     ;; 3/3 grepALL


;lim_blank = 10 					;; 1/3 grepALL
;w=where(abs(what_I_plot) le lim_blank)			;; 2/3 grepALL
;if (w[0] ne -1) then what_I_plot[w]=missing_value	;; 3/3 grepALL

;title_user='Sea-level surface pressure'+title_user	;; grepALL
subtitle_user='LMD mesoscale model'     		;; grepALL


;; 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


;------------------------------
;------------------------------
; some last possible settings 
;------------------------------
;------------------------------

backup_data='no'   	;; save ARWpost-generated text files ? default is no
already_data='no'  	;; load ARWpost-generated text files ? default is no
;datafolder='./WHERE_RESULTS_ARE/'	;; default is ./
;plotfolder='./WHERE_RESULTS_ARE/'       ;; default is ./
;plotfolder='./WHERE_PLOTS_ARE/'	        ;; default is ./
;datafolder='/d5/aslmd/WRF_Mars/EN_COURS/newTH20WC/month5/'
;datafolder='/d5/aslmd/WRF_Mars/ARTICLE/VM/'
datafolder='/d5/aslmd/WRF_Mars/DATA/VM/'
plotfolder=datafolder

textitle='generic'		;; one word.  default is 'plot'
texcomments='sol'	;; one word.  default is 'plot'