Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/README
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/README	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/README	(revision 2199)
@@ -0,0 +1,2 @@
+To use these idl tools you need to download some cliatogies and mask files;
+that you can find here: http://dodsp.idris.fr/reee512/NEMO_OUT/ORCA2_LIM/
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/fig_ht_fromQ.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/fig_ht_fromQ.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/fig_ht_fromQ.pro	(revision 2199)
@@ -0,0 +1,86 @@
+pro fig_ht_fromQ,  Q1, Q2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, MASK_FILENAME = mask_filename, _extra = ex
+  
+  compile_opt idl2, strictarrsubs
+
+@common
+
+  tmask_save = tmask  
+                                ; suppress 1/2 north ligne dans TMASK
+  tmask[jpi/2:jpi-1,jpj-1,0] = 0b ; 0 byte
+
+  Qave1 = moyenne(Q1,'xy') 
+  Qave2 = moyenne(Q2,'xy')
+
+  Qnet1 = Q1 - Qave1
+  Qnet2 = Q2 - Qave2
+
+  index=where(gphit eq max(gphit))
+  nx=index[0] mod jpi
+  gphi_save=gphit
+  gphit[0,*]=gphit[nx,*]
+                                ; 
+  msk = read_ncdf( 'atlmsk_nomed', filename = mask_filename, _extra = ex)
+  
+                                ; from 2D array to 3D array 
+  msk = msk.arr[*]#replicate(1., nzt)
+  
+                                ; x average (zonal mean)
+                                ; *1.E-15 to have PetaWatt
+  Q1x = moyenne(Qnet1*e2t,'x',/integration)*1.E-15
+  Q2x = moyenne(Qnet2*e2t,'x',/integration)*1.E-15
+  Q1x_atl = moyenne(Qnet1*e2t*msk,'x',/integration)*1.E-15
+  Q2x_atl = moyenne(Qnet2*e2t*msk,'x',/integration)*1.E-15
+                                ; northward heat flux transport from antartic  
+  htr1 = total(Q1x, /cumulative)
+  htr2 = total(Q2x, /cumulative)
+  htr1_atl = total(reverse(Q1x_atl), /cumulative)
+  htr1_atl = -reverse(htr1_atl)
+  htr2_atl = total(reverse(Q2x_atl), /cumulative)
+  htr2_atl = -reverse(htr2_atl)
+
+  bad_gphit = where(gphit[0,*] lt -30.)
+  htr1_atl[bad_gphit] = !Values.F_NaN   ; we take from values north till 30° South
+  htr2_atl[bad_gphit] = !Values.F_NaN   ; we take from values north till 30° South
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '004_MHT_'+sEXP1+'-'+sEXP2+'_y'+start_end+ '.ps', /landscape
+  endif
+  
+; update data informations
+  varname = 'MHT'
+  varunit = 'Pw'
+  vargrid = 'T' 
+
+  ibce1 = STRTRIM(Qave1, 1)
+  ibce2 = STRTRIM(Qave2, 1)
+
+  title = sEXP1+' (Red) MHT (PetaW) - (Blue) Atlantic MHT (PetaW) / year '+start_end
+  subtitle = sEXP1+' (Qnet='+ibce1+' W/m2) -  (Qnet='+ibce2+' W/m2)'
+
+  plt1d, htr1,'y',min = -2., max = 2.5, TITLE=title, SUBTITLE = subtitle $ 
+        , SMALL = [2,2,1], YTITLE = 'PW',  XGRIDSTYLE = 2, TICKLEN = 1, YGRIDSTYLE = 2,  _extra = ex   
+  plt1d, htr1_atl,'y', COLOR = 50, /ov1d,  _extra = ex
+
+ if sEXP1 ne sEXP2 then begin
+  
+  title = sEXP1+' - '+sEXP2+' MHT (TeraW) differences / year '+start_end
+  plt1d, (htr1-htr2)*1.e3,'y',min = 0., max = 500., TITLE = title $
+         , SMALL = [2,2,3], TICKLEN = 1, XGRIDSTYLE = 2 $
+         , YTITLE = 'PW', YGRIDSTYLE = 2,   /NOERASE, _extra = ex
+  plt1d, (htr1_atl - htr2_atl),'y', COLOR = 50, /ov1d,  _extra = ex
+ 
+ endif 
+
+  title = 'Qnet'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/004_MHT_'+sEXP1+'-'+sEXP2+'_y'+start_end+ '.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+  
+  tmask=tmask_save
+  gphit=gphi_save
+  
+  return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/getname.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/getname.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/getname.pro	(revision 2199)
@@ -0,0 +1,26 @@
+function getname,dir,exp,freq,suff,date
+  
+ file_list = file_search (dir+'/'+exp+'_*_*_'+freq+'_'+suff+'.nc', count = nb_file )
+  
+ file_list = file_basename( file_list )
+
+ filename = ''
+ i = 0
+
+ while  i lt nb_file  and filename eq '' do begin
+  
+    split_name = STRSPLIT(file_list[i],'_',/EXTRACT ) 
+    date1 = long(split_name[1])
+    date2 = long(split_name[2])
+
+
+    if date1 - date le 0 and  date - date2 le 0 then filename = dir+'/'+file_list[i]
+    i = i + 1
+    
+ endwhile 
+
+ if filename eq '' then print,'filename does NOT FOUND' 
+
+ return,filename
+
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/plt_bsf.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/plt_bsf.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/plt_bsf.pro	(revision 2199)
@@ -0,0 +1,41 @@
+pro plt_bsf, U1, U2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+
+  @common
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '005_bsf_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps',  portrait = 1
+  endif
+  domdef, 0, 6000
+
+
+  bb1 = bsf(U1.arr, refvalue = 0., refpoint = [25, 0] )
+  bb2 = bsf(U2.arr, refvalue = 0., refpoint = [25, 0])
+ 
+
+  title = sEXP1+' Barotropic Stream Function'+' / year '+start_end
+  plt, bb1, min =-200., max = 200., int = 10.,/portrait, FORMAT = '(I4)', STYLE= 'so0so' $
+       , small = [1,2,1], COAST_THICK = 2, TITLE = title,  _extra = ex
+
+
+  if sEXP1 ne sEXP2 then begin
+
+     title = sEXP1+' - '+sEXP2+' Barotropic Stream Function'+' / year '+start_end
+     plt, bb1.arr - bb2.arr , min = -20., max = 20., int = 2.,  FORMAT = '(I2)', STYLE= 'so0so' $
+          , small = [1,2,2], COAST_THICK = 2, TITLE = title, /NOERASE, _extra = ex
+   
+  endif
+  
+  title = 'Barotropic Stream Function'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/005_bsf_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+
+  return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/plt_msf.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/plt_msf.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/plt_msf.pro	(revision 2199)
@@ -0,0 +1,66 @@
+pro plt_msf, V1, V2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, SUBBASIN = subbasin, MASK_FILENAME = mask_filename,_extra = ex
+
+  compile_opt idl2, strictarrsubs
+
+  @common
+
+  if KEYWORD_SET(SUBBASIN) then subname = subbasin 
+
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '019_msf_'+subname+'_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps',  portrait = 1
+  endif
+
+  if KEYWORD_SET(SUBBASIN) then begin
+     
+     CASE 1 of 
+        subbasin eq 'Atl' : var = 'atlmsk_nomed'
+        subbasin eq 'Ind' : var = 'indmsk'
+        subbasin eq 'Pac' : var = 'indpacmsk'
+        subbasin eq 'Glo' : var = 'glomsk_nomed'
+     ENDCASE
+                                ;
+     msk = read_ncdf( var, filename = mask_filename, /nostruct, _extra = ex )
+
+  endif else msk = tmask[*,*,0]
+;
+
+  
+  CASE subname OF
+     'Glo':lat_ext=[-80,90]
+     'Atl':lat_ext=[-30,90]
+     'Ind':lat_ext=[-30,30]
+     'Pac':lat_ext=[-30,70]
+     
+  ENDCASE
+
+ if subbasin eq 'Pac' then subname = "IndoPac"
+  
+  domdef, 0, 6000
+
+  mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma)
+  mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma)
+  title=sEXP1+' Meridional Stream Function '+subname+'/ year '+start_end
+  pltz, mm1, -20., 20., int = 1., boxzoom = [ind[0:1],lat_ext, 0, 5500], /xindex  , FORMAT = '(I3)' $
+        , small=[1,2,1], COAST_THICK=2 , zoom = 5500, maskdta = ma, /no_partial, TITLE = title, style = 'so0so'
+
+  if sEXP1 ne sEXP2 then begin
+     
+     title=sEXP1+' - '+sEXP2+subname+'/ year '+start_end
+     pltz, mm1.arr-mm2.arr, -10., 10., int = 1., boxzoom = [ind[0:1],lat_ext, 0, 5500], /xindex,  FORMAT = '(I3)'  $
+           , small=[1,2,2], COAST_THICK=2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, /NOERASE, style = 'so0so'
+     
+  endif
+  
+  title = + subname+ 'Meridional Stream Function'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/019_msf_'+subname+'_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+
+  return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_com.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_com.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_com.pro	(revision 2199)
@@ -0,0 +1,5 @@
+;COMMON STD_VAR, std_EXP1, std_EXP2, std_YEAR, std_iodir_data
+;COMMON STD_FIL1, std_file1_T, std_file1_U, std_file1_V, std_file1_W, std_file1_I
+;COMMON STD_FIL2, std_file2_T, std_file2_U, std_file2_V, std_file2_W, std_file2_I
+;COMMON STD_LEVITUS, std_file_Levitus_T, std_file_Levitus_S
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_plots.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_plots.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_plots.pro	(revision 2199)
@@ -0,0 +1,375 @@
+pro std_plots, std_EXP1, std_EXP2, ystart, _extra=ex
+
+
+PRINT, 'usage: postscript '
+PRINT, '========================================================================'
+PRINT, 'std_plots, ''std_EXP1'', ''std_EXP2'', ystart, /postscript'
+PRINT, '========================================================================'
+PRINT, 'std_plots, ''core2000'', ''core2000'', 1991, /postscript'
+PRINT, '========================================================================'
+
+
+  compile_opt idl2, strictarrsubs
+
+  PRINT, ''
+  PRINT, ' 	############################################'
+  PRINT, ''
+  PRINT, '                    LAUNCH of std_plots'
+  PRINT, ''
+  PRINT, ' 	############################################'
+  PRINT, ''
+
+
+@common
+@initorca2
+@std_com
+
+  year_start = ystart * 10000L + 101L                                                 ; year of start (corresponding to filename)
+  year_end = (ystart + 9L) * 10000L + 1231L                                           ; year of end (corresponding to filename)
+  start_end = string(year_start,format= '(I08)')+'_'+string(year_end,format ='(I08)') ; variable string for filename
+;;;;; ======================CHANGE PATH================================================================= ;;;;;;;;;;
+  std_iodir_data    = '/path_of_data_in_NetCdf_format/'        			      ; path of data in NetCdf format
+  std_iodir_climato = '/path_of_climatological_data/' 	                              ; path of climatological data
+  std_iodir_mask    = '/path_of_mask_files/'	                                      ; path of mask files (ex: subbasins)
+;;;;; ======================END CHANGE PATH================================================================= ;;;;;;;;;;
+
+  
+;;;;; ==================================================== WARNING ================================================================= ;;;;;;;;;;
+  d1 = (ystart + 9L) * 10000L + 101L ;  date to select last record of file with 10 years of record
+  d2 = year_end  ; date to select last record of file with 10 years of record
+  id1 = (ystart + 9L) * 10000L + 1201L ; date to select starting date for (monthly) ice output 
+  id2 = d2 ; date to select last record for (monthly) ice output
+  id_mars_1 = (ystart + 9L) * 10000L + 301L ; date to select mars ice output 
+  id_mars_2 = (ystart + 9L) * 10000L + 331L  ; date to select mars ice output
+  id_sept_1 = (ystart + 9L) * 10000L + 901L ; date to select september ice output 
+  id_sept_2 = (ystart + 9L) * 10000L + 930L ; date to select september ice output
+;
+; Levitus 98
+  std_file_Levitus_T =  'data_1y_potential_temperature_nomask.nc'
+  std_file_Levitus_S =  'data_1y_salinity_nomask.nc'
+
+; Output run experience1
+  std_file1_T     = std_EXP1+'_'+start_end+'_1Y_grid_T.nc' ; grid T input file
+  std_file1_U     = std_EXP1+'_'+start_end+'_1Y_grid_U.nc' ; grid U input file
+  std_file1_V     = std_EXP1+'_'+start_end+'_1Y_grid_V.nc' ; grid V input file
+;  std_file1_W    = std_EXP1+'_'+start_end+'_1Y_grid_W.nc' ; grid W input file
+  std_file1_I     = std_EXP1+'_'+start_end+'_1M_icemod.nc' ; ice    input file
+  
+; Output run experience2
+  std_file2_T     = std_EXP2+'_'+start_end+'_1Y_grid_T.nc' ; grid T input file
+  std_file2_U     = std_EXP2+'_'+start_end+'_1Y_grid_U.nc' ; grid U input file
+  std_file2_V     = std_EXP2+'_'+start_end+'_1Y_grid_V.nc' ; grid V input file
+;  std_file2_W    = std_EXP2+'_'+start_end+'_1Y_grid_W.nc' ; grid W
+  std_file2_I     = std_EXP2+'_'+start_end+'_1M_icemod.nc' ; ice    input file
+
+
+  PRINT, ''
+  PRINT, '	std_EXP1       : ', std_EXP1
+  PRINT, '	std_EXP2       : ', std_EXP2
+  PRINT, '	std_iodir_data : ', std_iodir_data 
+  PRINT, '	std_file1T : ', std_file1_T
+  PRINT, '	std_file1U : ', std_file1_U
+  PRINT, '	std_file1V : ', std_file1_V
+;  PRINT, '	std_file1W : ', std_file1_W
+  PRINT, '	std_file2I : ', std_file1_I
+  PRINT, '	std_file2T : ', std_file2_T
+  PRINT, '	std_file2U : ', std_file2_U
+  PRINT, '	std_file2V : ', std_file2_V
+;  PRINT, '	std_file2W : ', std_file2_W
+  PRINT, '	std_file2I : ', std_file2_I
+  PRINT, ''
+
+;#########################################################################
+; reading variables
+
+;;; 3D ;;;
+; temperature
+  T1   = read_ncdf('thetao', d1, d2, filename = std_iodir_data+std_file1_T                          ) 
+  T2   = read_ncdf('thetao', d1, d2, filename = std_iodir_data+std_file2_T                          ) 
+  TLev = read_ncdf('votemper', filename = std_iodir_climato+std_file_Levitus_T                ) 
+  TRey = read_ncdf('sst'     , filename = std_iodir_climato+'NewREY_ORCA2_1991_2000_1y.nc'    ) 
+
+;help,T1,T2,TLev,TRey
+
+; salinity
+  S1   = read_ncdf('so', d1, d2, filename = std_iodir_data+std_file1_T             ) 
+  S2   = read_ncdf('so', d1, d2, filename = std_iodir_data+std_file2_T             ) 
+  SLev = read_ncdf('vosaline', filename = std_iodir_climato+std_file_Levitus_S     ) 
+
+; Net Downward heat flux
+  Q1   = read_ncdf('qns+qsr', d1, d2, filename = std_iodir_data+std_file1_T           ) 
+  Q2   = read_ncdf('qns+qsr', d1, d2, filename = std_iodir_data+std_file2_T           )
+  geo  = read_ncdf('heatflow', filename=  std_iodir_climato+'geothermal_heating.nc'   )
+  geo  = geo.arr*1.e-3      ; convert into W/m2
+;climatology 
+  QNET  = read_ncdf('qnet', filename = std_iodir_climato+'OAFlux_1my_01_12_1984_2004_orca2_qnet.nc'   )
+
+; erp (evaporation damping)
+  ERP1   = read_ncdf('wfcorr', d1, d2, filename = std_iodir_data+std_file1_T  )
+  ERP2   = read_ncdf('wfcorr', d1, d2, filename = std_iodir_data+std_file2_T  )
+  ERP1   = ERP1.arr * 86400
+  ERP2   = ERP2.arr * 86400
+
+; emp (evaporation minus precipitation)
+  EMP1   = read_ncdf('wfo', d1, d2, filename = std_iodir_data+std_file1_T     )
+  EMP2   = read_ncdf('wfo', d1, d2, filename = std_iodir_data+std_file2_T     )
+  EMP1   = EMP1.arr * 86400
+  EMP2   = EMP2.arr * 86400
+
+ 
+;mixed layer depth
+  MXL1 = read_ncdf('mldr10_1', d1, d2, filename = std_iodir_data+std_file1_T         )  ; 10 m
+  MXL2 = read_ncdf('mldr10_1', d1, d2, filename = std_iodir_data+std_file2_T         )  ; 10 m
+;climatology 
+  MLD  = read_ncdf('mld', filename = std_iodir_climato+'mld_DR003_c1m_ORCA2_1y.nc'   )
+
+; current
+  U1   = read_ncdf('uoce_eff', d1, d2, filename = std_iodir_data+std_file1_U         )
+  U2   = read_ncdf('uoce_eff', d1, d2, filename = std_iodir_data+std_file2_U         )
+
+; Take into account the Gent & McWilliams "Eddy Induced velocity"
+  V1   = read_ncdf('voce_eff', d1, d2, filename = std_iodir_data+std_file1_V         )
+  V2   = read_ncdf('voce_eff', d1, d2, filename = std_iodir_data+std_file2_V         )
+
+; ice
+  Ithi_mars_1   = read_ncdf('sit', id_mars_1, id_mars_2, filename = std_iodir_data+std_file1_I ) 
+  Ithi_mars_2   = read_ncdf('sit', id_mars_1, id_mars_2, filename = std_iodir_data+std_file2_I )
+
+  Ithi_sept_1   = read_ncdf('sit', id_sept_1, id_sept_2, filename = std_iodir_data+std_file1_I ) 
+  Ithi_sept_2   = read_ncdf('sit', id_sept_1, id_sept_2, filename = std_iodir_data+std_file2_I )
+
+  Ifra_mars_1   = read_ncdf('sic', id_mars_1, id_mars_2, filename = std_iodir_data+std_file1_I )
+  Ifra_mars_2   = read_ncdf('sic', id_mars_1, id_mars_2, filename = std_iodir_data+std_file2_I )
+
+  Ifra_sept_1   = read_ncdf('sic', id_sept_1, id_sept_2, filename = std_iodir_data+std_file1_I )
+  Ifra_sept_2   = read_ncdf('sic', id_sept_1, id_sept_2, filename = std_iodir_data+std_file2_I )
+
+ 
+;#########################################################################
+
+;#########################################################################
+;######################  STANDARD PLOTS   ################################
+;#########################################################################
+
+; fixed color tabled
+  lct,64
+
+  PRINT, '  ' 
+  PRINT, '				Experience treated:', std_EXP1
+  PRINT, '  ' 
+  PRINT, '  ' 
+  PRINT, '  ' 
+
+       ;;;;;;;; ErP & Emp salinity damping term
+       ;;;;;;;; =============================================
+  PRINT, '============================================= '
+  PRINT, ' PLOTS 2D '
+  PRINT, '============================================= '
+  PRINT, '                                ErP salinity damping term'
+  PRINT, '1'
+  xxx_Erp, ERP1, ERP2, std_EXP1, std_EXP2, start_end, _extra = ex
+  PRINT, '2'
+  xxx_Emp, EMP1, EMP2, std_EXP1, std_EXP2, start_end, _extra = ex
+ 
+       ;;;;; Net heat flux
+        ;;;;; =============
+  PRINT, ' '
+  PRINT, '                               Net heat flux'
+  PRINT, '3 '
+  xxx_Qnet, Q1, Q2, QNET, std_EXP1, std_EXP2, start_end, _extra = ex
+
+       ;;; Meridionnal Heat Transport
+        ;;; ===========================
+  PRINT, ' '
+  PRINT, '                               Meridionnal Heat Transport'
+  PRINT, '4 '
+  fig_ht_fromQ, Q1.arr+geo, Q2.arr+geo, std_EXP1, std_EXP2, start_end, MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+
+        ;;; Global Barotropic Function
+        ;;; ===========================
+  PRINT, ' '
+  PRINT, ' 				Global Barotropic Function ' 
+  PRINT, '5 '
+  plt_bsf, U1, U2, std_EXP1, std_EXP2, start_end, _extra=ex
+
+  ;; mean Temperature diff with New Reynolds
+  ;; =======================================
+  PRINT, ' '
+  PRINT, ' 				mean Temperature diff with New Reynolds'
+  PRINT, '6 '
+  xxx_SST, T1[*,*,0], T2[*,*,0], TRey, std_EXP1, std_EXP2, start_end, _extra=ex
+
+
+        ;;;; mean Salinity diff with Levitus
+        ;;;; ===============================
+  PRINT, ' '
+  PRINT, '                               mean Salinity diff with Levitus'
+  PRINT, '7 '
+  xxx_SSS, S1[*,*,0], S2[*,*,0], SLev[*,*,0], std_EXP1, std_EXP2, start_end, _extra=ex
+
+      ;;;;;;; Arctic mean Salinity diff with Levitus at z=105 meters
+        ;;;;;;; ======================================================
+  PRINT, ' '
+  PRINT, '                               Arctic mean Salinity diff with Levitus at z=105 meters'
+  PRINT, '8 exp1-levitus'
+  xxx_EXP1_S_ARC_z105, S1, SLev, std_EXP1, std_EXP2, start_end, _extra = ex       
+  PRINT, '9 exp2-levitus'
+  xxx_EXP2_S_ARC_z105, S1, S2, SLev, std_EXP1, std_EXP2, start_end, _extra = ex
+
+
+        ;;;;;;;; mean Temperature diff with Levitus at z=105 meters
+        ;;;;;;;; ==================================================
+  PRINT, ' '
+  PRINT, '                               mean Temperature diff with Levitus at z=105 meters'
+  PRINT, '10'
+  xxx_T_z105, T1, T2, TRey, std_EXP1, std_EXP2, start_end, _extra=ex
+                              
+
+        ;;;;;;;; mean Salinity diff with Levitus at z=105 meters
+        ;;;;;;;; ==================================================
+  PRINT, ' '
+  PRINT, '                               mean Salinity diff with Levitus at z=105 meters'
+  PRINT, '11'
+  xxx_S_z105, S1, S2, SLev, std_EXP1, std_EXP2, start_end, _extra=ex
+                            
+
+  PRINT, '============================================= '
+  PRINT, ' PLOTS 3D '
+  PRINT, '============================================= '
+ 
+  PRINT, '                         Mixed layer depth'
+  PRINT, '12 '
+  xxx_mxl10, MXL1, MXL2, MLD, std_EXP1, std_EXP2, start_end,  _extra = ex
+  PRINT, '13'
+  xxx_ZonMxl, MXL1, MXL2, MLD, std_EXP1, std_EXP2, start_end, _extra = ex
+ 
+
+        ;;;;;;;;; Zonal mean Temperature diff with Levitus
+        ;;;;;;;;; ========================================
+  PRINT, ' '
+  PRINT, '                               Zonal mean Temperature diff with Levitus'
+  PRINT, '14 '
+  xxx_zonal_T, T1, T2, TLev, std_EXP1, std_EXP2, start_end, MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+
+  PRINT, '14 '                      
+  xxx_zonal_T, T1, T2, TLev, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Atl', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+  
+  PRINT, '14 '
+  xxx_zonal_T, T1, T2, TLev, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Ind', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+  
+  PRINT, '14 '
+  xxx_zonal_T, T1, T2, TLev, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Pac', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+ 
+
+        ;;;;;;;;;; Zonal mean Salinity diff with Levitus
+        ;;;;;;;;;; ========================================
+  PRINT, ' '
+  PRINT, '                               Zonal mean Salinity diff with Levitus'
+  PRINT, '15 '
+  xxx_zonal_S, S1, S2, SLev, std_EXP1, std_EXP2, start_end, MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+                               
+  PRINT, '15 '
+  xxx_zonal_S, S1, S2, SLev, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Atl', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+                              
+  PRINT, '15'
+  xxx_zonal_S, S1, S2, SLev, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Ind', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex 
+
+  PRINT, '15 '
+  xxx_zonal_S, S1, S2, SLev, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Pac', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex 
+  
+        ;;;;;;; Seasonal & Annual mean Arctic/Antarctic Ice Thickness & Ice Concentration
+        ;;;;;; ==========================================================================
+  PRINT, ' '
+  PRINT, '                               Arctic/Antarctic Ice Thickness'
+  PRINT, ' 16'
+  xxx_ARC_Icethick_MARS, Ithi_mars_1, Ithi_mars_2, std_EXP1, std_EXP2, start_end, _extra = ex
+ 
+  PRINT, ' 16'
+  xxx_ARC_Icethick_SEPT, Ithi_sept_1, Ithi_sept_2, std_EXP1, std_EXP2, start_end, _extra = ex
+ 
+  PRINT, ' 17'
+  xxx_ANT_Icethick_MARS, Ithi_mars_1, Ithi_mars_2, std_EXP1, std_EXP2, start_end, _extra = ex
+
+  PRINT, ' 17'
+  xxx_ANT_Icethick_SEPT, Ithi_sept_1, Ithi_sept_2, std_EXP1, std_EXP2, start_end, _extra = ex
+
+  PRINT, '                               Arctic/Antarctic Ice Concentration (area of sea-ice per grid cell area'
+  PRINT, ' 18'
+  xxx_Iceleadfrac_MARS, Ifra_mars_1, Ifra_mars_2, std_EXP1, std_EXP2, start_end, _extra = ex
+
+  PRINT, ' 18'
+  xxx_Iceleadfrac_SEPT, Ifra_sept_1, Ifra_sept_2, std_EXP1, std_EXP2, start_end, _extra = ex
+
+                       
+        ;;; Global and Atlantic Meridional Function
+        ;;; =======================================
+  PRINT, ' '
+  PRINT, ' 				Global and Atlantic Meridional Function'
+  PRINT, '19 '
+  plt_msf, V1, V2, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Glo', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+  PRINT, '19 '
+  plt_msf, V1, V2, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Atl', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+  PRINT, '19 '
+  plt_msf, V1, V2, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Ind', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+  PRINT, '19 '
+  plt_msf, V1, V2, std_EXP1, std_EXP2, start_end, SUBBASIN = 'Pac', MASK_FILENAME = std_iodir_mask+'subbasins_orca21_nored.nc', _extra = ex
+
+
+  PRINT, '============================================= '
+  PRINT, ' Local PLOTS '
+  PRINT, '============================================= '
+
+        ;;;;;;;;; Equatorial Temperature 
+        ;;;;;;;;; ======================
+  PRINT, ' '
+  PRINT, '                               Equatorial Temperature'
+  PRINT, ' '
+  PRINT, '20 '
+  xxx_Eq_T, T1, T2, Tlev, std_EXP1, std_EXP2, start_end, _extra = ex
+
+        ;;;;;;;; Equatorial Salinity 
+        ;;;;;;;; ===================
+  PRINT, ' '
+  PRINT, '                               Equatorial Salinity'
+  PRINT, '21'
+  xxx_Eq_S, S1, S2, SLev, std_EXP1, std_EXP2, start_end, _extra = ex
+
+        ;;; Equatorial zonal velocity
+        ;;; =========================
+  PRINT, ' '
+  PRINT, '                               Equatorial zonal velocity'
+  PRINT, '22 '
+  xxx_Eq_U, U1, U2, std_EXP1, std_EXP2, start_end, _extra = ex
+  
+       ;;;;;;;;;;;; Mediterranean salt tongue 
+       ;;;;;;;;;;;; =============================================== 
+  
+  PRINT, ' '
+  PRINT, '                               Mediterranean salt tongue at depth=700 and lat 40°N'
+  PRINT, '23'
+  xxx_Med_S_tongue, S1, S2, SLev, std_EXP1, std_EXP2, start_end, DEPTH = 700,  LAT = 40, _extra = ex
+  ;
+  PRINT, '                               Mediterranean salt tongue at depth=1000 and lat 38°N'
+  PRINT, '23'
+  xxx_Med_S_tongue, S1, S2, SLev, std_EXP1, std_EXP2, start_end, DEPTH = 1000, LAT = 38, _extra = ex
+
+        ;;;;;;;;; Mediterranean 
+        ;;;;;;;;; ===============================================
+  PRINT, ' '
+  PRINT, '                               Mediterranean water at 38°N'
+  PRINT, '24'
+  xxx_Med_S_depth, S1, S2, SLev, std_EXP1, std_EXP2, start_end, DEPTH = 1000, LAT = 38, _extra = ex
+  PRINT, '                               Mediterranean water at 38°N'
+  PRINT, '24 '
+  xxx_Med_S_depth, S1, S2, SLev, std_EXP1, std_EXP2, start_end, DEPTH = 700,  LAT = 40, _extra = ex
+
+        ;;;;;;;;;;; Vertical Global mean T & S
+        ;;;;;;;;;;; ===========================
+  PRINT, ' '
+  PRINT, '                                Vertical Global mean T & S'
+  PRINT, '25 '
+  xxx_TS_ZGlobm, T1, T2, TLev, S1, S2, SLev, std_EXP1, std_EXP2, start_end, _extra = ex
+
+ 
+  return
+END
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_plots.sh
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_plots.sh	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_plots.sh	(revision 2199)
@@ -0,0 +1,198 @@
+#!/bin/sh 
+#
+# AUTHOR - date
+# ===========
+# Françoise Pinsard - 01/2010 - LOCEAN
+# Simona Flavoni  - 01/2010 - LOCEAN
+#
+# DESCRIPTION
+# ===========
+# Graphics outputs files are written in ./exp1_exp2/ps/.
+#
+# A global document with all plots in Postscript, PDF or HTML file can be produced.
+# 
+# EXAMPLES
+# ========
+# $ ./std_plots.sh -exp1 v32colib -exp2 v32colib -y 19800101_19801231 -format html
+# 
+# $ ./std_plots.sh -exp1 core100 -exp2 core100 -y 00910101_01001231 -format pdf
+#
+# PROBLEMS
+# ========
+# no batch processing possible because of SAXO interactiviy Seb is working on it ...
+#
+# TODO
+# ====
+# portability on linux
+# define overriding policy
+# finalize ps and pdf production 
+#
+usage=" Usage : ${command} -exp1 exp1 -exp2 exp2 -y start -format html
+ example:  ./std_plots.sh -exp1 core2000 -exp2 core2000 -y 1991 -format pdf 
+ exp1 is your experience name
+ exp2 is the reference experience (if you don't have one exp2 = exp1)
+ ystart is the first year of your 10years file output; example: 361 if you want to take file core2000_03610101_03701231_1Y_grid_T.nc"
+#
+minargcount=8
+if [ ${#} -lt ${minargcount} ]; then
+ echo "Not enought arguments"
+ echo "${usage}"
+ exit 1
+fi
+#
+while [ ! -z "${1}" ] # ++ pb bash
+do
+ case ${1} in
+     -exp1)	 exp1=${2}        shift  ;;
+     -exp2)	 exp2=${2}        shift  ;;
+     -y)	 ystart=${2}	  shift	 ;;
+     -format)	 format=${2}	  shift	 ;;
+     -h)	
+	 echo "${usage}"
+	 exit 0	
+	 ;;
+     *) # other choice
+	 echo "${usage}"
+	 exit 1
+	 ;;
+ esac
+ shift # next flag
+done
+#
+# ++ check parameters
+set -u
+#
+# define output directory for POSTCRIPT files
+#++ change directory
+#===============CHANGE DIRECTORIES==========================
+psdir=~/PS/${exp1}_${exp2}/ps/
+pdfdir=~/PDF/${exp1}_${exp2}/pdf
+htmldir=~/HTML/${exp1}_${exp2}/pdf
+#===============END CHANGE DIRECTORIES=======================
+# test if ${psdir} already exists
+if [ -d ${psdir} ]; then
+ echo "${psdir} already exist"
+ #+++exit 1
+else
+ mkdir -p ${psdir}
+fi 
+# test if ${pdfdir} already exists
+if [ -d ${pdfdir} ]; then
+ echo "${psdir} already exist"
+ #+++exit 1
+else
+ mkdir -p ${pdfdir}
+fi 
+
+# memorize current IDL_STARTUP file
+memo_idl_startup=${IDL_STARTUP}
+
+# generation of IDL commands file to launch several std_plots commands
+# NOTE:  $(( 10#${ystart} )) is used to have decimal, if not the number that start with 0 (for example if you put 0361) would be interpeted like alphanumerical number 
+rm /tmp/std_plots_multi.pro 2> /dev/null
+cat << EOF >> /tmp/std_plots_multi.pro
+ psdir='${psdir}'
+ std_plots,'${exp1}','${exp2}', $(( 10#${ystart} )), /postscript
+ exit
+EOF
+#
+idl_command='/Applications/itt/idl64/bin/idl'  # !! only on my mac
+#
+# check if IDL available
+type ${idl_command}
+status_type=${?}
+if [ ${status_type} -ne 0 ]
+then
+ echo "eee : idl not found"
+ exit 1
+fi
+# debug
+more /tmp/std_plots_multi.pro
+#
+# launch IDL
+${idl_command} /tmp/std_plots_multi.pro
+status_idl=${?}
+# test if idl returns ok
+#echo "iii : ${status_idl}"
+#
+# restore IDL_STARTUP file
+IDL_STARTUP=${memo_idl_startup}
+export IDL_STARTUP
+
+# produce the final document
+case ${format} in
+    pdf)
+	rm ${pdfdir}/all_plots.pdf 2> /dev/null
+	for file in ${psdir}/*.ps ; do
+	    echo "ps2pdf"
+#           echo "nome file: ${file}"
+		ps2pdf -sPAPERSIZE=a4 ${file}
+	done
+#	echo "dentro pdf"
+	ls ${pdfdir}
+	mv *.pdf ${pdfdir}/.
+	texexec --pdfarrange --result=/tmp/all_plots.pdf ${pdfdir}/*.pdf
+	rm -rf *.aux *.log /tmp/*.log
+	mv /tmp/all_plots.pdf ${pdfdir}
+#commented because convert gives images of bad quality
+#	convert ${psdir}/*.ps /tmp/all_plots.pdf
+#       convert -resize 800x600 ${psdir}/${exp1}_${exp2}/ps/*.ps /tmp/all_plots.pdf
+	;;
+    html)
+	rm /tmp/std_plots.html 2> /dev/null
+# check if convert available
+	convert_command='/sw/bin/convert'  # !! only on my mac
+	type ${convert_command}
+	status_type=${?}
+	if [ ${status_type} -ne 0 ]
+	then
+	    echo "eee : convert not found"
+	    exit 1
+	fi
+# debug
+	
+	for file in `ls ${psdir}/*.ps | sed -e "s/\.ps//"` ; do
+	    ${convert_command} -antialias  ${file}.ps  ${file}.png
+	done
+	 # creation of a HTML file
+	cat << EOF > /tmp/std_plots.html
+	<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
+	<html xmlns="http://www.w3.org/1999/xhtml">
+	<head>
+	<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
+        <hr>
+       	</head>
+	<body>
+EOF
+	cat /tmp/std_plots.txt >>  /tmp/std_plots.html
+	cat << EOF >> /tmp/std_plots.html
+<hr/>
+EOF
+	
+	cat << EOF >> /tmp/std_plots.html
+
+  </body>
+  </html>
+EOF
+	more /tmp/std_plots.html
+	;;
+    *)	
+	echo " format ${format} not implemented"
+	exit 1
+	;;
+esac
+
+    #
+#=====================
+# prepare to put images on dods
+# +++ to finish it
+#LOGIN=xxxx
+#Tag_Name=ORCA2_LIM2
+#rsh ${LOGIN}@gaya.idris.fr exec /bin/ksh <<EOF
+#       cd ${homegaya}/IGCM_OUT/${Tag_Name}/CORE2/INTERAN/plots_pdf/${exp1}
+#       /usr/local/bin/dods_rm DODS/pub/${LOGIN}/${Tag_Name}/CORE2/INTERAN/plots_pdf/${exp1} > /dev/null 2>&1
+#       /usr/local/bin/dods_cp ${exp1} DODS/pub/${LOGIN}/${Tag_Name}/CORE2/INTERAN/plots_pdf/${exp1} > /dev/null 2>&1
+#EOF
+#=====================
+# end
+exit 0
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_ts.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_ts.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/std_ts.pro	(revision 2199)
@@ -0,0 +1,99 @@
+pro std_ts, std_EXP1, std_EXP2, year1, year2, _extra=ex
+
+
+PRINT, 'usage: postscript '
+PRINT, '========================================================================'
+PRINT, 'std_ts, ''std_EXP1'', ''std_EXP2'', ''year_start'', ''year_end'', /postscript'
+PRINT, '========================================================================'
+PRINT, 'std_ts, ''core2000'', ''core2000'', ''1'', ''2000'', /postscript'
+PRINT, '========================================================================'
+
+
+  compile_opt idl2, strictarrsubs
+
+  PRINT, ''
+  PRINT, ' 	############################################'
+  PRINT, ''
+  PRINT, '                   plot of all TIME SERIES'
+  PRINT, ''
+  PRINT, ' 	############################################'
+  PRINT, ''
+
+
+@common
+@initorca2
+@std_com
+
+
+;  std_iodir_data    = '/Users/sflod/idl_PLOTS/DATA_STORE/RUN_INTERAN/mercure/' ; path of data in NetCdf format
+  std_iodir_data    = '/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/brodie/' ; path of data in NetCdf format
+
+
+;#########################################################################
+;######################  TIME SERIES   ################################
+;#########################################################################
+
+; fixed color tabled
+  lct,64
+
+ 
+
+  ;; Temperature Mean
+  ;; =======================================
+  PRINT, ' '
+  PRINT, ' 			 Temperature Mean'
+  PRINT, '1 '
+  ts_mean_T, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex
+
+  ;; Salinity Mean
+  ;; =======================================
+  PRINT, ' '
+  PRINT, ' 			 Salinity Mean'
+  PRINT, '2 '
+  ts_mean_S, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex
+  
+  ;; Sea Surface Height
+  ;; =======================================
+  PRINT, ' '
+  PRINT, ' 			 SSH'
+  PRINT, '3 '
+  ts_mean_SSH, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex
+
+   ;; Heat flux Mean
+  ;; =======================================
+  PRINT, ' '
+  PRINT, ' 			 Heat Mean'
+  PRINT, '4 '
+  ts_mean_Q, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex
+
+  ;; EMP Mean
+  ;; =======================================
+  PRINT, ' '
+  PRINT, ' 			 EMP Mean'
+  PRINT, '5 '
+  ts_mean_EMP, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex
+
+
+  ;; Drake transport
+  ;; =======================================
+  PRINT, ' '
+  PRINT, ' 		         Drake Transport'
+  PRINT, '6 '
+  ts_mean_U, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex
+
+  ;; Max AMOC
+  ;; =======================================
+  PRINT, ' '
+  PRINT, ' 		         Max AMOC'
+  PRINT, '7 '
+  ts_mean_V, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex
+
+  ;; Sea-Ice Area 
+  ;; =======================================
+  PRINT, ' '
+  PRINT, ' 		         Sea-Ice Area'
+  PRINT, '8 '
+  ts_mean_ICE, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex
+
+  return
+END
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/time_series.sh
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/time_series.sh	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/time_series.sh	(revision 2199)
@@ -0,0 +1,183 @@
+#!/bin/sh 
+#
+# AUTHOR - date
+# ===========
+# Simona Flavoni  - 02/2010 - LOCEAN
+#
+# DESCRIPTION
+# ===========
+# Graphics outputs files are written in my_IDL/PS/ps_ts/.
+#
+# A global document with all plots in Postscript, PDF or HTML file can be produced.
+# 
+# EXAMPLES
+# ========
+# $ ./time_series.sh -exp1 sEXP1 -exp2 sEXP2 -ystart year1 -yend year2 -format pdf
+# 
+# PROBLEMS
+# ========
+# no batch processing possible because of SAXO interactiviy Seb is working on it ...
+#
+# TODO
+# ====
+# portability on linux
+# define overriding policy
+# finalize ps and pdf production 
+#
+usage=" Usage : ${command} -exp1 exp1 -exp2 exp2 -ystart start -yend end -format pdf
+ exp1 is your experience name
+ exp2 is name of experience of reference 
+ ystart is the first year from which you want to compute Time Series 
+ yend is the last year for which you want to compute Time Series 
+ example:  ./time_series.sh -exp1 core2000 -exp2 core2000 -ystart 1 -yend 2000 -format pdf "
+#
+minargcount=10
+if [ ${#} -lt ${minargcount} ]; then
+ echo "Not enought arguments"
+ echo "${usage}"
+ exit 1
+fi
+#
+while [ ! -z "${1}" ] # ++ pb bash
+do
+ case ${1} in
+     -exp1)	 exp1=${2}        shift  ;;
+     -exp2)	 exp2=${2}        shift  ;;
+     -ystart)	 start=${2}	  shift	 ;;
+     -yend)	 end=${2}	  shift	 ;;
+     -format)	 format=${2}	  shift	 ;;
+     -h)	
+	 echo "${usage}"
+	 exit 0	
+	 ;;
+     *) # other choice
+	 echo "${usage}"
+	 exit 1
+	 ;;
+ esac
+ shift # next flag
+done
+#
+# ++ check parameters
+set -u
+#
+# define output directory for POSTCRIPT files
+#++ change directory
+#================CHANGE DIRECTORIES=====================
+psdir=/path_output_PS_time_series
+pdfdir=/path_output_PDF_time_series
+htmldir=/path_output_HTML_time_series
+#================END CHANGE DIRECTORIES=================
+# test if ${psdir} already exists
+if [ -d ${psdir} ]; then
+ echo "${psdir} already exist"
+ #+++exit 1
+else
+ mkdir -p ${psdir}
+fi 
+# test if ${pdfdir} already exists
+if [ -d ${pdfdir} ]; then
+ echo "${psdir} already exist"
+ #+++exit 1
+else
+ mkdir -p ${pdfdir}
+fi 
+
+# memorize current IDL_STARTUP file
+memo_idl_startup=${IDL_STARTUP}
+
+# generation of IDL commands file to launch several std_plots commands
+rm /tmp/time_series_multi.pro 2> /dev/null
+cat << EOF >> /tmp/time_series_multi.pro
+ psdir='${psdir}'
+ std_ts,'${exp1}','${exp2}',${start},${end},/postscript
+ exit
+EOF
+#
+idl_command='/Applications/itt/idl64/bin/idl'  # !! only on my mac
+#
+# check if IDL available
+type ${idl_command}
+status_type=${?}
+if [ ${status_type} -ne 0 ]
+then
+ echo "eee : idl not found"
+ exit 1
+fi
+# debug
+more /tmp/time_series_multi.pro
+#
+# launch IDL
+${idl_command} /tmp/time_series_multi.pro
+status_idl=${?}
+# test if idl returns ok
+# echo "iii : ${status_idl}"
+#
+# restore IDL_STARTUP file
+IDL_STARTUP=${memo_idl_startup}
+export IDL_STARTUP
+
+# produce the final document
+case ${format} in
+    pdf)
+	rm ${pdfdir}/all_ts.pdf 2> /dev/null
+	echo "psdir: ${psdir}"
+	for file in ${psdir}/*.ps ; do
+	    echo "ps2pdf"
+		ps2pdf -sPAPERSIZE=a4 ${file}
+	done
+	ls ${pdfdir}
+	mv *.pdf ${pdfdir}/.
+	texexec --pdfarrange --result=/tmp/all_ts.pdf ${pdfdir}/*.pdf
+	rm -rf *.aux *.log /tmp/*.log
+	mv /tmp/all_ts.pdf ${pdfdir}/all_ts.pdf
+#commented because convert gives images of bad quality
+#	convert ${psdir}/*.ps /tmp/all_ts.pdf
+#       convert -resize 800x600 ${psdir}/${exp1}_${exp2}/ps/*.ps /tmp/all_ts.pdf
+	;;
+    html)
+	rm /tmp/time_series.html 2> /dev/null
+# check if convert available
+	convert_command='/sw/bin/convert'  # !! only on my mac
+	type ${convert_command}
+	status_type=${?}
+	if [ ${status_type} -ne 0 ]
+	then
+	    echo "eee : convert not found"
+	    exit 1
+	fi
+# debug
+	
+	for file in `ls ${psdir}/*.ps | sed -e "s/\.ps//"` ; do
+	    ${convert_command} -antialias  ${file}.ps  ${file}.png
+	done
+	 # creation of a HTML file
+	cat << EOF > /tmp/time_series.html
+	<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
+	<html xmlns="http://www.w3.org/1999/xhtml">
+	<head>
+	<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
+        <hr>
+       	</head>
+	<body>
+EOF
+	cat /tmp/time_series.txt >>  /tmp/time_series.html
+	cat << EOF >> /tmp/time_series.html
+<hr/>
+EOF
+	
+	cat << EOF >> /tmp/time_series.html
+
+  </body>
+  </html>
+EOF
+	more /tmp/time_series.html
+	;;
+    *)	
+	echo " format ${format} not implemented"
+	exit 1
+	;;
+esac
+
+    #
+exit 0
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_EMP.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_EMP.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_EMP.pro	(revision 2199)
@@ -0,0 +1,66 @@
+pro ts_mean_EMP, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+  
+  @common
+
+ @initorca2
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, sEXP1+'_EMP_mean.ps', portrait = 1
+ endif
+
+ ts_EMP = 0.
+ ts_EMP2 = 0.
+ ts_Time = 0.
+
+ for year=year1,year2 do begin
+
+ date1 = year * 10000L + 101
+ date2 = year * 10000L + 1231
+ filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1)
+ filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1)
+
+ if filename eq '' then stop 
+ EMP = read_ncdf('wfo',date1 ,date2 , filename = filename, /nostruct)
+ EMP2 = read_ncdf('wfo',date1 ,date2 , filename = filename2, /nostruct)
+
+ ts_EMP = [ ts_EMP, moyenne(EMP, 'xy')]
+ ts_EMP2 = [ ts_EMP2, moyenne(EMP2, 'xy')]
+ ts_Time = [ ts_Time, Time] ; Time = date
+
+ endfor
+
+  surf_oce =  e1t * e2t  * tmask[*,*,0]
+  surf_oce = total(surf_oce)
+
+  ts_EMP = ts_EMP[1:*] ; remove first record of 0
+  ts_EMP = ts_EMP * 1.E-09 * surf_oce ; to have Sv
+
+  ts_EMP2 = ts_EMP2[1:*]                ; remove first record of 0
+  ts_EMP2 = ts_EMP2 * 1.E-09 * surf_oce ; to have Sv
+
+  varunit= 'Sv'
+  vol_oce = moyenne( replicate(1., jpi, jpj, jpk), /integration, 'xyz')
+
+  Time = ts_Time[1:*] ; remove first record of 0
+  jpt = n_elements(Time)
+
+  title = sEXP1+' Global Annual EMP Budget (Sv) --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt,ts_EMP,'t', -.01, 0.00, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
+       , small = [1,2,1], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex
+
+  if sEXP2 ne sEXP1 then begin
+     
+     title = sEXP1+' - '+sEXP2+' Global Annual EMP Budget (Sv) --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+     pltt,ts_EMP - ts_EMP2,'t', -.002, .002, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI, /NOERASE $
+          , COLOR = 250, small = [1,2,2], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex
+  endif
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+
+  return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_ICE.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_ICE.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_ICE.pro	(revision 2199)
@@ -0,0 +1,100 @@
+pro ts_mean_ICE, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+  
+  @common
+  
+  @initorca2
+  
+  std_iodir_climato = '/Users/sflod/idl_PLOTS/CLIMATOLOGIES/'
+  std_iodir_interan = '/Users/sflod/idl_PLOTS/INTERAN/'
+
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, sEXP1+'_ice_area_mean.ps', portrait = 1
+  endif
+  
+  
+
+;NORTHERN Hemisphere
+  ts_ICE_N = 0.
+  ts_ICE_N_2 = 0.
+
+;SOUTHERN Hemisphere
+  ts_ICE_S = 0.
+  ts_ICE_S_2 = 0.
+
+  ts_Time = 0.
+  
+  for year=year1,year2 do begin
+     date1 = year * 10000L + 101
+     date2 = year * 10000L + 1231
+     filename = getname(s_iodir_data,sEXP1,'1M','icemod',date1)
+     filename2 = getname(s_iodir_data,sEXP2,'1M','icemod',date1)
+
+     if filename eq '' then stop
+     
+     domdef, 0, jpi-1, 30, 90, /xindex
+
+; computation of annual mean
+     ICE_N = read_ncdf('sic',date1 ,date2 , filename = filename, /nostruct, dir = 't')
+     ICE_N_2 = read_ncdf('sic',date1 ,date2 , filename = filename2, /nostruct, dir = 't')
+     
+     ts_ICE_N   = [ ts_ICE_N, moyenne(ICE_N, 'xy', /integration) ]
+     ts_ICE_N_2 = [ ts_ICE_N_2, moyenne(ICE_N_2, 'xy', /integration) ]
+    
+     domdef, 0, jpi-1, -30, -90, /xindex
+     
+     ICE_S = read_ncdf('sic',date1 ,date2 , filename = filename, /nostruct, dir = 't')
+     ICE_S_2 = read_ncdf('sic',date1 ,date2 , filename = filename2, /nostruct, dir = 't')
+     
+     ts_ICE_S   = [ ts_ICE_S, moyenne(ICE_S, 'xy', /integration) ]
+     ts_ICE_S_2 = [ ts_ICE_S_2, moyenne(ICE_s_2, 'xy', /integration) ]
+ 
+    ts_Time  = [ ts_Time, Time] ; Time = date
+ 
+  endfor
+
+;  surf_oce =  e1t[firstxt:lastxt,firstyt:lastyt] * e2t[firstxt:lastxt,firstyt:lastyt]  * tmask[firstxt:lastxt,firstyt:lastyt,0]
+;  surf_oce = total(surf_oce)
+
+  ts_ICE_N   = ts_ICE_N[1:*]        ; remove first record of 0 
+  ts_ICE_N_2 = ts_ICE_N_2[1:*]      ; remove first record of 0
+  ts_ICE_N   = ts_ICE_N * 1.E-12    ; remove first record of 0
+  ts_ICE_N_2 = ts_ICE_N_2 * 1.E-12  ; remove first record of 0
+  
+  ts_ICE_S   = ts_ICE_S[1:*]        ; remove first record of 0 
+  ts_ICE_S_2 = ts_ICE_S_2[1:*]      ; remove first record of 0
+  ts_ICE_S   = ts_ICE_S * 1.E-12    ; remove first record of 0
+  ts_ICE_S_2 = ts_ICE_S_2 * 1.E-12  ; remove first record of 0
+  
+  Time = ts_Time[1:*]           ; remove first record of 0
+  jpt = n_elements(Time)
+  
+  title = sEXP1+' Global Annual Mean Ice Area (North. Hemisp.) --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt, ts_ICE_N, 't', 0., 15., year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
+       , small = [1,2,1], YTITLE='10^12 m^2 ', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+  
+  title = sEXP1+' Global Annual Mean Ice Area (South. Hemisp.) --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt, ts_ICE_S, 't', 0., 15., year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI, /NOERASE $
+        , small = [1,2,2], YTITLE='10^12 m^2 ', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+  
+  if sEXP2 ne sEXP1 then begin
+     title = sEXP1+'  - '+sEXP2+' Global Annual Mean Ice Area (North. Hemisp.) --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+     pltt, ts_ICE_N - ts_ICE_N_2 , 't', -2., 2., year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
+           , COLOR = 250, small = [1,2,1], YTITLE='10^12 m^2 ', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+     
+     title = sEXP1+' - '+sEXP2+' Global Annual Mean Ice Area (South. Hemisp.) --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+     pltt, ts_ICE_S - ts_ICE_S_2, 't', -2., 2., year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI, /NOERASE $
+           , COLOR = 250, small = [1,2,2], YTITLE='10^12 m^2 ', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+
+  endif
+
+  
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+  
+  
+  return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_Q.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_Q.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_Q.pro	(revision 2199)
@@ -0,0 +1,59 @@
+pro ts_mean_Q, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+  
+  @common
+  
+ @initorca2
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, sEXP1+'_Heat_mean.ps', portrait = 1
+ endif
+ 
+ ts_Heat = 0.
+ ts_Heat2 = 0.
+ ts_Time = 0.
+
+ for year=year1,year2 do begin
+    
+    date1 = year * 10000L + 101
+    date2 = year * 10000L + 1231
+    filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1)
+    filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1)
+
+    if filename eq '' then stop 
+    Heat = read_ncdf('qns+qsr',date1 ,date2 , filename = filename, /nostruct)
+    Heat2 = read_ncdf('qns+qsr',date1 ,date2 , filename = filename2, /nostruct)
+
+    ts_Heat = [ ts_Heat, moyenne(Heat, 'xy')]
+    ts_Heat2 = [ ts_Heat2, moyenne(Heat2, 'xy')]
+    ts_Time = [ ts_Time, Time]  ; Time = date
+
+ endfor
+ 
+ ts_Heat = ts_Heat[1:*]         ; remove first record of 0
+ ts_Heat2 = ts_Heat2[1:*]         ; remove first record of 0
+ Time = ts_Time[1:*]            ; remove first record of 0
+ 
+ jpt = n_elements(Time)
+ 
+ title = sEXP1+' Global Annual Mean Surface Heat Flux --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+ pltt,ts_Heat,'t', -4., 4., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
+      , small = [1,2,1], YTITLE = 'W/m2',  TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+
+ if sEXP2 ne sEXP1 then begin
+
+    title = sEXP1+' - '+sEXP2+' Global Annual Mean Surface Heat Flux --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+    pltt,ts_Heat - ts_Heat2,'t', -4., 4., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI, /NOERASE $
+         ,  COLOR = 250, small = [1,2,2], YTITLE = 'W/m2', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+    
+ endif
+
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+ 
+ 
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_S.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_S.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_S.pro	(revision 2199)
@@ -0,0 +1,103 @@
+pro ts_mean_S, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+  
+  @common
+  
+  @initorca2
+  
+  std_iodir_climato = '/Users/sflod/idl_PLOTS/CLIMATOLOGIES/'
+  std_iodir_interan = '/Users/sflod/idl_PLOTS/INTERAN/'
+
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, sEXP1+'_Salinity_mean.ps', portrait = 1
+  endif
+  
+  filename = std_iodir_climato+'data_1y_salinity_nomask.nc'
+  Lev = read_ncdf('vosaline', filename = filename, /nostruct )
+  Lev_xyz = moyenne(Lev, 'xyz')
+  Lev_z = moyenne(Lev, 'xy',/KEEPBOTTOM)
+
+  ts_Sal = 0.
+  ts_Sal_2 = 0.
+  ts_Sal_Lev = 0.
+  ts_Sal_Lev_2 = 0.
+  ts_Time = 0.
+  ts_z = fltarr(jpk)
+  ts_z_2 = fltarr(jpk)
+
+  for year=year1,year2 do begin
+     
+     date1 = year * 10000L + 101
+     date2 = year * 10000L + 1231
+     filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1)
+     filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1)
+
+     if filename eq '' then stop 
+     Sal   = read_ncdf('so',date1 ,date2 , filename = filename, /nostruct)
+     Sal_2 = read_ncdf('so',date1 ,date2 , filename = filename2, /nostruct)
+     SSH   =  read_ncdf('zos',date1 ,date2 , filename = filename, /nostruct)
+     SSH_2 =  read_ncdf('zos',date1 ,date2 , filename = filename2, /nostruct)
+
+     ts_Sal = [ ts_Sal, moyenne(Sal, 'xyz', ssh = SSH) ]
+     ts_Sal_2 = [ ts_Sal_2, moyenne(Sal_2, 'xyz', ssh = SSH_2) ]
+     ts_Sal_Lev = [ ts_Sal_Lev, moyenne(Sal, 'xyz', ssh = SSH) - Lev_xyz ]
+     ts_Sal_Lev_2 = [ ts_Sal_Lev_2, moyenne(Sal_2, 'xyz', ssh = SSH_2) - Lev_xyz ]
+     ts_Time = [ ts_Time, Time] ; Time = date
+
+     ts_z = [ [ts_z], [moyenne(Sal, 'xy',/KEEPBOTTOM, ssh = SSH)- Lev_z] ]
+     ts_z_2 = [ [ts_z_2], [moyenne(Sal_2, 'xy',/KEEPBOTTOM, ssh = SSH_2)- Lev_z] ]
+
+  endfor
+  
+  ts_Sal = ts_Sal[1:*]             ; remove first record of 0
+  ts_Sal_2 = ts_Sal_2[1:*]         ; remove first record of 0
+  ts_Sal_Lev = ts_Sal_Lev[1:*]     ; remove first record of 0
+  ts_Sal_Lev_2 = ts_Sal_Lev_2[1:*] ; remove first record of 0
+  
+  ts_z = ts_z[*,1:*]            ; remove first record of 0
+  ts_z_2 = ts_z_2[*,1:*]            ; remove first record of 0
+  
+  Time = ts_Time[1:*]           ; remove first record of 0
+  jpt = n_elements(Time)
+
+
+  title = sEXP1+' Global Annual Mean Salinity --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt, ts_Sal, 't', 34.54, 34.76, INTER = .2, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
+       , small = [1,2,1], YTITLE='psu', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+
+  title = sEXP1+' Global Annual Mean Salinity - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt, ts_Sal_Lev, 't', -.05, .05, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
+       , small = [1,2,2], YTITLE='psu',  TITLE = title, FORMAT = '(I3)', /NOERASE, _extra = ex
+
+;;;;;;;;;;;;;;;;; new page ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+ERASE
+ 
+  title = sEXP1+' ZT-plot Salinity Mean (0-1500m) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt,ts_z ,'zt', -.5, .5, INTER = .025, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $ 
+       , small = [1,2,1], YTITLE='psu', STYLE = 'so0so',  TITLE = title, boxzoom = 1500., /PORTRAIT, _extra = ex
+
+  title = sEXP1+' ZT-plot Salinity Mean (0-6000) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt,ts_z,'zt', -.5, .5, INTER = .025, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $ 
+       , small = [1,2,2], YTITLE='psu', STYLE = 'so0so',  TITLE = title, boxzoom = 6000., zoom = 1500.,  /NOERASE, _extra = ex
+  
+
+if sEXP2 ne sEXP1 then begin
+
+ERASE
+
+   title = sEXP1+' - '+sEXP2+' Global Annual Mean Salinity --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+   pltt, ts_Sal - ts_Sal_2, 't', -1., 1., INTER = .2, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $
+         , COLOR = 250, small = [1,2,1], YTITLE='psu', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+   
+endif
+
+  
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+  
+  
+  return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_SSH.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_SSH.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_SSH.pro	(revision 2199)
@@ -0,0 +1,60 @@
+pro ts_mean_SSH, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+  
+  @common
+  
+  @initorca2
+  
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, sEXP1+'_SSH_mean.ps', portrait = 1
+  endif
+  
+  ts_SSH = 0.
+  ts_SSH2 = 0.
+  ts_Time = 0.
+  
+ 
+  for year=year1,year2 do begin
+     
+     date1 = year * 10000L + 101
+     date2 = year * 10000L + 1231
+     filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1)
+     filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1)
+     
+     SSH = read_ncdf('zos',date1 ,date2 , filename = filename, /nostruct)
+     SSH2 = read_ncdf('zos',date1 ,date2 , filename = filename2, /nostruct)
+
+     ts_SSH = [ ts_SSH, moyenne(SSH, 'xy')]
+     ts_SSH2 = [ ts_SSH2, moyenne(SSH2, 'xy')]
+     ts_Time = [ ts_Time, Time] ; Time = date
+     
+  endfor
+  
+  ts_SSH = ts_SSH[1:*]          ; remove first record of 0
+  ts_SSH2 = ts_SSH2[1:*]          ; remove first record of 0
+  Time = ts_Time[1:*]           ; remove first record of 0
+  
+  jpt = n_elements(Time)
+  
+  title = sEXP1+' Global Annual Mean SSH --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt,ts_SSH,'t', -.001, .001, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
+       , small = [1,2,1], YTITLE = 'm', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+  
+  if sEXP2 ne sEXP1 then begin
+     
+     title = sEXP1+' - '+sEXP2+' Global Annual Mean SSH --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+     pltt, ts_SSH - ts_SSH2,'t', -.0001, .0001, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI, /NOERASE $
+           , COLOR = 250, small = [1,2,2], YTITLE = 'm', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+     
+  endif
+  
+
+ 
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+
+  return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_T.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_T.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_T.pro	(revision 2199)
@@ -0,0 +1,105 @@
+pro ts_mean_T, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+  
+  @common
+  
+  @initorca2
+  
+  std_iodir_climato = '/Users/sflod/idl_PLOTS/CLIMATOLOGIES/'
+;  std_iodir_interan = '/Users/sflod/idl_PLOTS/INTERAN/'
+  std_iodir_data = '/Users/sflod/idl_PLOTS/INTERAN/mercure/'
+
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, sEXP1+'_Temperature_mean.ps', portrait = 1
+  endif
+  
+  filename = std_iodir_climato+'data_1y_potential_temperature_nomask.nc'
+  Lev = read_ncdf('votemper', filename = filename, /nostruct )
+  Lev_xyz = moyenne(Lev, 'xyz')
+  Lev_z = moyenne(Lev, 'xy',/KEEPBOTTOM)
+
+  ts_Temp = 0.
+  ts_Temp_2 = 0.
+  ts_Temp_Lev = 0.
+  ts_Temp_Lev_2 = 0.
+  ts_Time = 0.
+  ts_z = fltarr(jpk)
+  ts_z_2 = fltarr(jpk)
+  
+  for year=year1,year2 do begin
+     
+     date1 = year * 10000L + 101L
+     date2 = year * 10000L + 1231L
+ 
+print,s_iodir_data
+
+     filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1)
+     filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1)
+
+     if filename eq '' then stop 
+     Temp   = read_ncdf('thetao',date1 ,date2 , filename = filename, /nostruct)
+     Temp_2 = read_ncdf('thetao',date1 ,date2 , filename = filename2, /nostruct)
+     SSH    =  read_ncdf('zos',date1 ,date2 , filename = filename, /nostruct)
+     SSH_2  =  read_ncdf('zos',date1 ,date2 , filename = filename2, /nostruct)
+
+     ts_Temp = [ ts_Temp, moyenne(Temp, 'xyz',ssh = SSH) ]
+     ts_Temp_2 = [ ts_Temp_2, moyenne(Temp_2, 'xyz',ssh = SSH_2) ]
+     ts_Temp_Lev = [ ts_Temp_Lev, moyenne(Temp, 'xyz',ssh = SSH) - Lev_xyz ]
+     ts_Temp_Lev_2 = [ ts_Temp_Lev_2, moyenne(Temp_2, 'xyz',ssh = SSH_2) - Lev_xyz ]
+     ts_Time = [ ts_Time, Time] ; Time = date
+
+     ts_z = [ [ts_z], [moyenne(Temp, 'xy',/KEEPBOTTOM, ssh = SSH)- Lev_z] ]
+     ts_z_2 = [ [ts_z_2], [moyenne(Temp_2, 'xy',/KEEPBOTTOM, ssh = SSH_2)- Lev_z] ]
+
+  endfor
+
+  ts_Temp = ts_Temp[1:*]        ; remove first record of 0 
+  ts_Temp_2 = ts_Temp_2[1:*]        ; remove first record of 0
+  ts_Temp_Lev = ts_Temp_Lev[1:*] ; remove first record of 0
+  ts_Temp_Lev_2 = ts_Temp_Lev_2[1:*] ; remove first record of 0
+  ts_z = ts_z[*,1:*]             ; remove first record of 0
+  ts_z_2 = ts_z_2[*,1:*]             ; remove first record of 0
+  
+  Time = ts_Time[1:*]           ; remove first record of 0
+  jpt = n_elements(Time)
+  
+  title = sEXP1+' Global Annual Mean Temperature --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt, ts_Temp, 't', 1., 4.5, INTER = .5, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
+       , small = [1,2,1], YTITLE='Deg C', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+
+  title = sEXP1+' Global Annual Mean Temperature - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt,ts_Temp_Lev, 't', -1., 1., INTER = .25,  year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
+       , small = [1,2,2], YTITLE='Deg C', TITLE = title, FORMAT = '(I3)', /NOERASE, _extra = ex
+
+  ;;;;;;;;;;;;;;;;; new page  ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+
+ERASE
+
+  title = sEXP1+' ZT-plot Temperature Mean (0-1500m) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt,ts_z ,'zt', -2., 2., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $ 
+       , small = [1,2,1], YTITLE='Deg C', TITLE = title, boxzoom = 1500., /NOERASE, _extra = ex
+  
+  title = sEXP1+' ZT-plot Temperature Mean (0-6000m) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+  pltt,ts_z ,'zt', -2., 2., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $ 
+       , small = [1,2,2], YTITLE='Deg C', TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex
+  
+if sEXP2 ne sEXP1 then begin
+  
+ERASE
+ 
+   title = sEXP1+' - '+sEXP2+' Global Annual Mean Temperature --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+   pltt, ts_Temp - ts_Temp_2, 't', -1., 1., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
+         , COLOR = 250, small = [1,2,1], YTITLE='Deg C', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+   
+endif
+
+  
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+  
+  
+  return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_U.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_U.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_U.pro	(revision 2199)
@@ -0,0 +1,77 @@
+pro ts_mean_U, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
+
+
+  compile_opt idl2, strictarrsubs
+  
+  @common
+  
+ @initorca2
+ 
+ if max(glamt) gt 300 then testlam = abs(glamt - 300) else testlam = abs(glamt + 60)
+
+ testlat = abs(gphit + 30)
+
+ index=where(abs(testlam - min(testlam)) lt 1 and abs(testlat - min(testlat)) lt 1 )
+ xindex=index[0] mod jpi
+ yindex=index[0]/jpi
+
+ domdef, xindex, xindex+1, 0, yindex, 0, jpk-1, /index
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, sEXP1+'_U_mean.ps', portrait = 1
+ endif
+ 
+ ts_bbD = 0.
+ ts_bbD2 = 0.
+ ts_Time = 0.
+ 
+ for year=year1,year2 do begin
+    
+    date1 = year * 10000L + 101
+    date2 = year * 10000L + 1231
+    filename = getname(s_iodir_data,sEXP1,'1Y','grid_U',date1)
+    filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_U',date1)
+
+    if filename eq '' then stop 
+    U = read_ncdf('uoce_eff',date1 ,date2 , filename = filename, /nostruct)
+    U2 = read_ncdf('uoce_eff',date1 ,date2 , filename = filename2, /nostruct)
+     
+    bb1 = bsf(U)
+    bb2 = bsf(U2)
+    bb1=bb1.arr[0,*]
+    bb2=bb2.arr[0,*]
+
+    bbDRAKE = bb1[nyt-1]-bb1[0]
+    bb2DRAKE = bb2[nyt-1]-bb2[0]
+
+    ts_bbD = [ ts_bbD, bbDRAKE ]
+    ts_bbD2 = [ ts_bbD2, bb2DRAKE ]
+    ts_Time = [ ts_Time, Time]  ; Time = date
+       
+ endfor
+ 
+ ts_bbD = ts_bbD[1:*]           ; remove first record of 0
+ ts_bbD2 = ts_bbD2[1:*]         ; remove first record of 0
+ Time = ts_Time[1:*]            ; remove first record of 0
+ 
+ jpt = n_elements(Time)
+
+ title = sEXP1+' DRAKE transport --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+ pltt,ts_bbD,'t', 0., 200., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
+      , small = [1,2,1],  TITLE = title, YTITLE='Sv', FORMAT = '(I3)', /PORTRAIT, _extra = ex
+ 
+ if sEXP2 ne sEXP1 then begin
+    
+    title = sEXP1+' - '+sEXP2+' DRAKE transport --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+    pltt,ts_bbD - ts_bbD2,'t', -30., 30., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI, /NOERASE $
+        , COLOR = 250, small = [1,2,2],  TITLE = title, YTITLE='Sv', FORMAT = '(I3)', /PORTRAIT, _extra = ex
+    
+ endif
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+ 
+ 
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_V.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_V.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/ts_mean_V.pro	(revision 2199)
@@ -0,0 +1,132 @@
+pro ts_mean_V, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex
+
+  
+  compile_opt idl2, strictarrsubs
+  
+  @common
+  
+ @initorca2
+ 
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, sEXP1+'_V_mean.ps', portrait = 1
+ endif
+ 
+; inizialization of vectors (for experience 1 and 2) 
+ ts_mAMOC45 = 0.
+ ts_mAABW70 = 0.
+ ts_mAABC50 = 0.
+ ts_mAMOC45_2 = 0.
+ ts_mAABW70_2 = 0.
+ ts_mAABC50_2 = 0.
+
+ ts_Time = 0.
+ 
+ for year=year1,year2 do begin
+    
+    date1 = year * 10000L + 101
+    date2 = year * 10000L + 1231
+    filename = getname(s_iodir_data,sEXP1,'1Y','grid_V',date1)
+    filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_V',date1)
+    
+    if filename eq '' then stop 
+    
+    domdef, 0, jpk-1,/zindex
+
+    V1 = read_ncdf('voce_eff',date1 ,date2 , filename = filename, /nostruct)
+    V2 = read_ncdf('voce_eff',date1 ,date2 , filename = filename2, /nostruct)
+    
+    msf_tot = msf(V1, indexboxzoom = ind)
+    msf2_tot = msf(V2, indexboxzoom = ind)
+    
+    yaxis = gphit[ind[0], ind[2]:ind[3]]
+ 
+; computation of max Atlatic Meridional Overturninc Circulation at 40°N and 50°N
+;;  see plot msf global
+    indy = where(yaxis gt 40 and yaxis le 50)
+    domdef, 0, 3500
+    
+    mAMOC45 = max(msf_tot.arr[indy, firstzw:lastzw],/NaN)
+    mAMOC45_2 = max(msf2_tot.arr[indy, firstzw:lastzw],/NaN)
+
+; computation of max atlantic Antarctic Bottom Water between 80°S and 65°S
+;; see plot msf global
+    indy = where(yaxis gt -80 and yaxis le -65)
+    domdef, 300, 3500
+    
+    mAABW70 = min(msf_tot.arr[indy, firstzw:lastzw],/NaN)
+    mAABW70_2 = min(msf2_tot.arr[indy, firstzw:lastzw],/NaN)
+
+;  computation of max Antarctic Abyssal Bottom Cell between 65°S and 30°N , 
+;; see plot msf global
+    indy = where(yaxis gt -65 and yaxis le 30)
+    domdef, 2500, 5000
+    
+    mAABC50 = min(msf_tot.arr[indy, firstzw:lastzw],/NaN)
+    mAABC50_2 = min(msf2_tot.arr[indy, firstzw:lastzw],/NaN)
+
+; creation of Time Series
+; for EXPERIENCE 1
+    ts_mAMOC45 = [ ts_mAMOC45, mAMOC45 ]
+    ts_mAABW70 = [ ts_mAABW70, mAABW70 ]
+    ts_mAABC50 = [ ts_mAABC50, mAABC50 ]
+
+; for EXPERIENCE 2
+    ts_mAMOC45_2 = [ ts_mAMOC45_2, mAMOC45_2 ]
+    ts_mAABW70_2 = [ ts_mAABW70_2, mAABW70 ]
+    ts_mAABC50_2 = [ ts_mAABC50_2, mAABC50 ]
+   
+    ts_Time = [ ts_Time, Time]  ; Time = date
+    
+ endfor
+
+ ts_mAMOC45 = ts_mAMOC45[1:*]   ; remove first record of 0
+ ts_mAABW70 = ts_mAABW70[1:*]   ; remove first record of 0
+ ts_mAABC50 = ts_mAABC50[1:*]   ; remove first record of 0
+ 
+ ts_mAMOC45_2 = ts_mAMOC45_2[1:*]   ; remove first record of 0
+ ts_mAABW70_2 = ts_mAABW70_2[1:*]   ; remove first record of 0
+ ts_mAABC50_2 = ts_mAABC50_2[1:*]   ; remove first record of 0
+
+ Time = ts_Time[1:*]            ; remove first record of 0
+ 
+ jpt = n_elements(Time)
+ 
+ title = sEXP1+' Max Atlantic MOC between 40N and 50N --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+ pltt, ts_mAMOC45,'t', 0., 30., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
+      ,  small = [1,2,1],  TITLE = title, YTITLE='Sv', FORMAT = '(I3)', /PORTRAIT, _extra = ex
+ 
+ if sEXP2 ne sEXP1 then begin
+    title = sEXP1+' - '+sEXP2+' Max Atlantic MOC between 40N and 50N --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+    pltt, ts_mAMOC45 - ts_mAMOC45_2,'t', -9., 9., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI, /NOERASE $
+          ,  COLOR = 250, small = [1,2,2],  TITLE = title, YTITLE='Sv', FORMAT = '(I3)', /PORTRAIT, _extra = ex
+ endif
+ 
+ title = sEXP1+' Max AntArctic Bottom Water between 80S and 65S --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+ pltt,-ts_mAABW70,'t', 0., 20., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
+      , small = [1,2,1],  TITLE = title, YTITLE='Sv', FORMAT = '(I3)', /PORTRAIT, _extra = ex
+ 
+ if sEXP2 ne sEXP1 then begin 
+    title = sEXP1+' - '+sEXP2+' Max AntArctic Bottom Water between 80S and 65S --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+    pltt,-ts_mAABW70 + ts_mAABW70_2 ,'t', 0., 20., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI, /NOERASE $
+         ,  COLOR = 250, small = [1,2,2],  TITLE = title, YTITLE='Sv', FORMAT = '(I3)', /PORTRAIT, _extra = ex
+ endif
+ 
+ title = sEXP1+' Max AntArctic Bottom Cell between 65S and 30N --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+ pltt,-ts_mAABC50,'t', 5., 30., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $
+      ,  small = [1,2,1],  TITLE = title, YTITLE='Sv', FORMAT = '(I3)', /PORTRAIT, _extra = ex
+
+ if sEXP2 ne sEXP1 then begin
+
+    title = sEXP1+' - '+sEXP2+' Max AntArctic Bottom Cell between 65S and 30N --- year '+strtrim(year1,1)+' - '+strtrim(year2,1)
+    pltt,-ts_mAABC50 + ts_mAABC50_2,'t', 5., 30., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI,/NOERASE  $
+         ,  COLOR = 250, small = [1,2,2],  TITLE = title, YTITLE='Sv', FORMAT = '(I3)', /PORTRAIT, _extra = ex
+ endif
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+ 
+ 
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ANT_Icethick_MARS.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ANT_Icethick_MARS.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ANT_Icethick_MARS.pro	(revision 2199)
@@ -0,0 +1,44 @@
+pro xxx_ANT_Icethick_MARS,  Ithi_mars_1, Ithi_mars_2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+
+ @common
+
+ compile_opt idl2, strictarrsubs
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '017_ANT_IceThick_MARS_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+ endif
+ 
+                                ;
+ domdef,20,380,-50,-90
+                                ; 
+ title=sEXP1+' Antartic Ice Thickness MARS / year '+start_end
+
+ Ithi_mars_1.arr = Ithi_mars_1.arr < 1.
+ plt, Ithi_mars_1.arr - 1.E-04, MIN = 0., MAX = 1., INTER = 0.1 , STYLE = 'so0so',CELL_FILL = 2, /STRICTFILL $ 
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title  $
+      , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [-90,0,0], /PORTRAIT, _extra = ex
+
+ if sEXP1 ne sEXP2 then begin
+    
+    title=sEXP1+' - '+sEXP2+' Antartic Ice Thickness MARS / year '+start_end
+    plt, Ithi_mars_1.arr - Ithi_mars_2.arr, MIN = -1., MAX = 1., INTER = 0.1, /STRICTFILL , STYLE = 'so0so' $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title $
+         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [-90,0,0], /NOERASE, _extra = ex
+    
+ endif
+ 
+ domdef
+
+ title = 'Antarctic Sea-ice Thickness'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/017_ANT_IceThick_MARS_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+
+ return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ANT_Icethick_SEPT.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ANT_Icethick_SEPT.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ANT_Icethick_SEPT.pro	(revision 2199)
@@ -0,0 +1,44 @@
+pro xxx_ANT_Icethick_SEPT,  Ithi_sept_1, Ithi_sept_2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+
+ @common
+
+ compile_opt idl2, strictarrsubs
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '017_ANT_IceThick_SEPT_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+ endif
+ 
+                                ;
+ domdef,20,380,-50,-90
+                                ; 
+ title=sEXP1+' Antartic Ice Thickness SEPT / year '+start_end
+
+ Ithi_sept_1.arr = Ithi_sept_1.arr < 1.
+ plt, Ithi_sept_1.arr - 1.E-04, MIN = 0., MAX = 1., INTER = 0.1 , STYLE = 'so0so',CELL_FILL = 2, /STRICTFILL $ 
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title  $
+      , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [-90,0,0], /PORTRAIT, _extra = ex
+
+ if sEXP1 ne sEXP2 then begin
+    
+    title=sEXP1+' - '+sEXP2+' Antartic Ice Thickness SEPT / year '+start_end
+    plt, Ithi_sept_1.arr - Ithi_sept_2.arr, MIN = -1., MAX = 1., INTER = 0.1, /STRICTFILL , STYLE = 'so0so' $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title $
+         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [-90,0,0], /NOERASE, _extra = ex
+    
+ endif
+ 
+ domdef
+
+ title = 'Antarctic Sea-ice Thickness'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/017_ANT_IceThick_SEPT_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+
+ return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ARC_Icethick_MARS.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ARC_Icethick_MARS.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ARC_Icethick_MARS.pro	(revision 2199)
@@ -0,0 +1,41 @@
+pro xxx_ARC_Icethick_MARS, Ithi_mars_1, Ithi_mars_2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ @common
+
+ compile_opt idl2, strictarrsubs
+
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '016_ARC_IceThick_MARS_'+sEXP1+'_'+sEXP2+'y'+start_end+'.ps', portrait = 1
+ endif
+ 
+ domdef,20,380,55,90
+           
+ title=sEXP1+' Arctic Ice Thickness MARS/ year '+start_end
+
+ Ithi_mars_1.arr = Ithi_mars_1.arr < 6.
+ plt, Ithi_mars_1.arr - 1.E-04, MIN = 0., MAX = 6., INTER = 0.25 ,/STRICTFILL , CELL_FILL = 2, STYLE = 'so0so'  $
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title $
+      , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [90,0,0], /PORTRAIT, _extra = ex
+ ;
+ if sEXP1 ne sEXP2 then begin            
+    title = sEXP1+' - '+sEXP2+' Arctic Ice Thickness / year '+start_end
+    plt, Ithi_mars_1.arr - Ithi_mars_2.arr, MIN = -3., MAX = 3., INTER = .5, /STRICTFILL , STYLE = 'so0so' $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title $
+         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [90,0,0], /NOERASE, _extra = ex
+    
+ endif                          ; end of 2
+ 
+ domdef
+
+ title = 'Arctic Sea-ice Thickness MARS'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/016_ARC_IceThick_MARS_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ARC_Icethick_SEPT.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ARC_Icethick_SEPT.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ARC_Icethick_SEPT.pro	(revision 2199)
@@ -0,0 +1,41 @@
+pro xxx_ARC_Icethick_SEPT, Ithi_sept_1, Ithi_sept_2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ @common
+
+ compile_opt idl2, strictarrsubs
+
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '016_ARC_IceThick_SEPT_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+ endif
+ 
+ domdef,20,380,55,90
+           
+ title=sEXP1+' Arctic Ice Thickness SEPT / year '+start_end
+
+ Ithi_sept_1.arr = Ithi_sept_1.arr < 6.
+ plt, Ithi_sept_1.arr - 1.E-04, MIN = 0., MAX = 6., INTER = 0.25 ,/STRICTFILL , CELL_FILL = 2, STYLE = 'so0so'  $
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title $
+      , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [90,0,0], /PORTRAIT, _extra = ex
+ ;
+ if sEXP1 ne sEXP2 then begin            
+    title = sEXP1+' - '+sEXP2+' Arctic Ice Thickness SEPT / year '+start_end
+    plt, Ithi_sept_1.arr - Ithi_sept_2.arr, MIN = -3., MAX = 3., INTER = .5, /STRICTFILL , STYLE = 'so0so' $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title $
+         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [90,0,0], /NOERASE, _extra = ex
+    
+ endif                          ; end of 2
+ 
+ domdef
+
+ title = 'Arctic Sea-ice Thickness'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/016_ARC_IceThick_SEPT_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_EXP1_S_ARC_z105.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_EXP1_S_ARC_z105.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_EXP1_S_ARC_z105.pro	(revision 2199)
@@ -0,0 +1,38 @@
+pro xxx_EXP1_S_ARC_z105, S1, SLev, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ @common 
+ compile_opt idl2, strictarrsubs
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '008_S_ARC_z105_'+sEXP1+'_y'+start_end+'.ps', portrait = 1
+ endif
+
+                                ; 
+ domdef,20,380,60,90
+                                ;
+ title=sEXP1+' Salinity 105 m'
+ plt, S1, MIN = 30.0, MAX = 35., INTER = 0.25, CELL_FILL = 2  $
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $
+      , /ORTHO, MAP = [90,0,0], LATDEL = 5, boxzoom = [105., 106.], /PORTRAIT, _extra = ex
+                                ; 
+
+ title=sEXP1+' - Levitus Salinity 105 m / year '+start_end
+ plt, S1.arr-SLev.arr, MIN = -2., MAX = 2., INTER = 0.2, CELL_FILL = 2,  STYLE = 'so0so'  $     
+      , small = [1,2,2], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $      
+      , /ORTHO, MAP = [90,0,0], LATDEL = 5, boxzoom = [105., 106.], /NOERASE, _extra = ex
+                                ;
+ domdef
+
+ title = 'Artic salinity at 105 meters'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/008_S_ARC_z105_'+sEXP1+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+ 
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_EXP2_S_ARC_z105.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_EXP2_S_ARC_z105.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_EXP2_S_ARC_z105.pro	(revision 2199)
@@ -0,0 +1,47 @@
+pro xxx_EXP2_S_ARC_z105, S1, S2, SLev, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ @common 
+ compile_opt idl2, strictarrsubs
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '009_S_ARC_z105_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+ endif
+
+                                ; 
+ domdef,20,380,60,90
+                                ;
+ if sEXP1 ne sEXP2 then begin
+
+    title=sEXP1+' - '+sEXP2+' Salinity 105 m / year '+start_end
+    plt, S1.arr-S2.arr, MIN = -2., MAX = 2., INTER = 0.2, CELL_FILL = 2,  STYLE = 'so0so'  $     
+         , small = [1,2,1], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $      
+         , /ORTHO, MAP = [90,0,0], LATDEL = 5, boxzoom = [105., 106.], /NOERASE, _extra = ex
+    
+ endif else begin
+    
+    title='Levitus Salinity 105 m'
+    plt, SLev.arr, MIN = 33., MAX = 35., INTER = 0.25 , CELL_FILL = 2 $ 
+         , small = [1,2,1], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $ 
+         , /ORTHO, MAP = [90,0,0], LATDEL = 5, boxzoom = [105., 106.], _extra = ex
+ endelse
+                                ;
+ title=sEXP2+' - Levitus Salinity 105 m / year '+start_end
+ plt, S2.arr-SLev.arr, MIN = -2., MAX = 2., INTER = 0.2, CELL_FILL = 2,  STYLE = 'so0so'  $     
+      , small = [1,2,2], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $      
+      , /ORTHO, MAP = [90,0,0], LATDEL = 5, boxzoom = [105., 106.], /NOERASE, _extra = ex
+                                ;
+ domdef
+;
+ title = 'Artic salinity at 105 meters'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/009_S_ARC_z105_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+ 
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Emp.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Emp.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Emp.pro	(revision 2199)
@@ -0,0 +1,32 @@
+pro xxx_Emp, EMP1, EMP2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ 
+  compile_opt idl2, strictarrsubs
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '002_Emp_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+  endif
+  
+  title=sEXP1+' Emp (mm/day) / year '+start_end
+  plt, EMP1,  MIN = -15., MAX = 15., INTER = 2.5 $
+       , small = [1,2,1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra=ex
+  
+if sEXP1 ne sEXP2 then begin
+
+  title=sEXP1+' - '+sEXP2+' / Emp year '+start_end
+  plt, EMP1-EMP2, MIN = -5., MAX = 5., INTER = 1. $
+       , small = [1,2,2], COAST_THICK = 2, TITLE = title,  /NOERASE, _extra=ex
+endif
+
+  title = 'Emp field'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/002_Emp_'+sEXP1+'_'+sEXP2+'_y '+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+  return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Eq_S.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Eq_S.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Eq_S.pro	(revision 2199)
@@ -0,0 +1,41 @@
+pro xxx_Eq_S, S1, S2, SLev, sEXP1, sEXP2, start_end,  POSTSCRIPT = postscript, _extra = ex
+  
+  compile_opt idl2, strictarrsubs
+
+  S = S1.arr-S2.arr
+  sEXP = sEXP1+' - '+sEXP2
+  min =  -1.
+  max = -min
+  inter = .1
+  
+  if sEXP2 eq sEXP1 then begin
+     sEXP = sEXP1+' - Levitus' 
+     S = S1.arr-Slev.arr
+     min = -1.
+     max = -min
+     inter = .1
+  endif
+  
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '021_Eq_S_'+sEXP+'_y'+start_end+'.ps', portrait = 1
+  endif
+                                ;
+  title = sEXP1+' Equatorial S  / year '+start_end
+  pltz, S1, MININ = 33., MAXIN = 37., INTER = .2, typein='xz'  $
+        , small=[1,2,1], COAST_THICK = 2, TITLE = title, boxzoom = [20., 380., -1., 1., 0., 500.], /PORTRAIT, _extra = ex
+                                ; 
+  title = sEXP+' Eq. S / year '+start_end
+  pltz, S, MININ = min, MAXIN = max, INTER = inter, typein = 'xz', STYLE = 'so0so' $
+        , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = [20., 380., -1., 1., 0., 500.], /PORTRAIT, _extra = ex
+;                              
+  title = 'Equatorial differences with Levitus and '+sEXP+' of salinity'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/021_Eq_S_'+sEXP+'_y'+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+  return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Eq_T.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Eq_T.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Eq_T.pro	(revision 2199)
@@ -0,0 +1,42 @@
+pro xxx_Eq_T, T1, T2, TLev, sEXP1, sEXP2, start_end,  POSTSCRIPT = postscript, _extra = ex
+  
+  compile_opt idl2, strictarrsubs
+
+  T = T1.arr-T2.arr
+  sEXP = sEXP1+' - '+sEXP2
+  min =  -2.
+  max = -min
+  inter = .25
+  
+  if sEXP2 eq sEXP1 then begin
+     sEXP = sEXP1+' - Levitus'
+     T = T1.arr-TLev.arr
+     min = -4.
+     max = -min
+     inter = .5
+   endif
+  
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '020_Eq_T_'+sEXP+'_y'+start_end+'.ps', portrait = 1
+  endif
+                                ;
+  title = sEXP1+' Equatorial T  / year '+start_end
+  pltz, T1, MININ = 2., MAXIN = 30., INTER = 1., typein = 'xz'  $
+        , small = [1,2,1], COAST_THICK = 2, TITLE = title, boxzoom = [20., 380., -1., 1., 0., 500.], /PORTRAIT, _extra = ex
+                                ; 
+ 
+  title = sEXP+' Eq. T / year '+start_end
+  pltz, T, MININ = min, MAXIN = max, INTER = inter, typein = 'xz', STYLE = 'so0so' $
+        , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = [20., 380., -1., 1., 0., 500.], /NOERASE, _extra = ex
+                                ; 
+  title ='Equatorial differences with Levitus and '+sEXP+'of temperature'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/020_Eq_T_'+sEXP+'_y'+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+  
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+  return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Eq_U.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Eq_U.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Eq_U.pro	(revision 2199)
@@ -0,0 +1,30 @@
+pro xxx_Eq_U, U1, U2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+  compile_opt idl2, strictarrsubs
+  
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '022_Eq_U_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+  endif
+                                ;
+  title = sEXP1+' Equatorial U_zonal / year '+start_end
+  pltz, U1, MININ = -1., MAXIN = 1., INTER = .1, typein = 'xz', STYLE = 'so0so', UNIT = varunit $
+        , small = [1,2,1], COAST_THICK = 2, TITLE = title, boxzoom = [20., 380., -1., 1., 0., 500.], ZOOM = 500, /PORTRAIT, _extra = ex
+
+  if sEXP1 ne sEXP2 then begin
+                                ;
+     title = sEXP1+' - '+sEXP2+' U_zonal / year '+start_end
+     pltz, U1.arr-U2.arr, MININ = -.5, MAXIN = .5, INTER = .1, typein = 'xz', STYLE = 'so0so', UNIT = varunit $
+           , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = [20., 380., -1., 1., 0., 500.], ZOOM = 500, /PORTRAIT, _extra = ex
+                                ;
+  endif
+  
+  title ='Equatorial differences with Levitus and '+sEXP2+'of zonal velocity U'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/021_Eq_U_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+  
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+  return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Erp.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Erp.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Erp.pro	(revision 2199)
@@ -0,0 +1,33 @@
+pro xxx_Erp, ERP1, ERP2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ 
+  compile_opt idl2, strictarrsubs
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '001_Erp_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+  endif
+  
+  title=sEXP1+' Erp (mm/day) / year '+start_end
+  plt, ERP1,  MIN = -5., MAX = 5., INTER = .5 $
+       , small = [1,2,1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra=ex
+  
+  if sEXP1 ne sEXP2 then begin
+
+     title=sEXP1+' - '+sEXP2+' / Erp year '+start_end
+     plt, ERP1-ERP2, MIN = -2., MAX = 2., INTER = .25 $
+          , small = [1,2,2], COAST_THICK = 2, TITLE = title,  /NOERASE, _extra=ex
+     
+  endif
+
+  title = 'Erp field'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/001_Erp_'+sEXP1+'_'+sEXP2+'_y '+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+  return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Iceleadfrac_MARS.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Iceleadfrac_MARS.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Iceleadfrac_MARS.pro	(revision 2199)
@@ -0,0 +1,31 @@
+pro xxx_Iceleadfrac_MARS, Ifra_mars_1, Ifra_mars_2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+ 
+  compile_opt idl2, strictarrsubs
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '018_Iceleadfrac_MARS_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+  endif
+  
+  title=sEXP1+' Iceleadfrac MARS / year '+start_end
+  plt, Ifra_mars_1,  MIN = 0., MAX = 1., INTER = .1 $
+       , small = [1,2,1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra=ex
+  
+if sEXP1 ne sEXP2 then begin
+
+  title=sEXP1+' - '+sEXP2+' / Iceleadfrac MARS year '+start_end
+  plt, Ifra_mars_1.arr - Ifra_mars_2.arr, MIN = -1., MAX = 1., INTER = .1 $
+       , small = [1,2,2], COAST_THICK = 2, TITLE = title,  /NOERASE, _extra=ex
+endif
+
+  title = 'Iceldeafrac  field'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/018_Iceleadfrac_MARS_'+sEXP1+'_'+sEXP2+'_y '+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+  return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Iceleadfrac_SEPT.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Iceleadfrac_SEPT.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Iceleadfrac_SEPT.pro	(revision 2199)
@@ -0,0 +1,31 @@
+pro xxx_Iceleadfrac_SEPT, Ifra_sept_1, Ifra_sept_2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+  compile_opt idl2, strictarrsubs
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     openps, '018_Iceleadfrac_SEPT_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+  endif
+  
+  title=sEXP1+' Iceleadfrac SEPT / year '+start_end
+  plt, Ifra_sept_1,  MIN = 0., MAX = 1., INTER = .1 $
+       , small = [1,2,1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra=ex
+  
+if sEXP1 ne sEXP2 then begin
+
+  title=sEXP1+' - '+sEXP2+' / Iceleadfrac SEPT year '+start_end
+  plt, Ifra_sept_1.arr - Ifra_sept_2.arr, MIN = -1., MAX = 1., INTER = .1 $
+       , small = [1,2,2], COAST_THICK = 2, TITLE = title,  /NOERASE, _extra=ex
+endif
+
+  title = 'Iceldeafrac field'
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/018_Iceleadfrac_SEPT_'+sEXP1+'_'+sEXP2+'_y '+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+
+
+  if KEYWORD_SET(POSTSCRIPT) then begin
+     closeps
+  endif
+
+  return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Med_S_depth.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Med_S_depth.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Med_S_depth.pro	(revision 2199)
@@ -0,0 +1,45 @@
+pro xxx_Med_S_depth, S1, S2, SLev, sEXP1, sEXP2, start_end, DEPTH=depth, LAT=lat, POSTSCRIPT = postscript,  _extra = ex
+  
+  compile_opt idl2, strictarrsubs
+
+   @common
+
+
+   if KEYWORD_SET(POSTSCRIPT) then begin
+      openps, '024_Med_S_depth_'+sEXP1+'_'+strtrim(lat,1)+'_d'+strtrim(depth,1)+'_S_y'+start_end+'.ps'   
+   endif
+;domdef, 0, 6000
+
+;domdef,270,365,5,70
+   minval = min(abs(gdept -  depth),index)
+   
+   title = sEXP1+'  Mediterranean water at ' +strtrim(lat,1)+'°N / year '+start_end 
+   pltz, S1.arr, MININ = 35., MAXIN = 37., INTER = 0.1,  UNIT = varunit $
+         , small = [1,2,1], COAST_THICK = 2, endpoints = [300., lat, 357., lat], TITLE = title, boxzoom = [2000.], ZOOM = 2000., /PORTRAIT, _extra = ex 
+;
+
+   if sEXP1 ne sEXP2 then begin
+      
+      title = sEXP1+' - '+sEXP2+' Mediterranean  water at '+strtrim(lat,1)+'°N / year '+start_end
+      pltz, S1.arr - S2.arr, MININ = -.4, MAXIN = .4, INTER = .05, UNIT = varunit $ 
+            , small = [1,2,2], COAST_THICK = 2, endpoints = [300., lat, 357., lat], TITLE = title, boxzoom = [2000.], ZOOM = 2000.,  _extra = ex 
+      
+   endif else begin
+   
+      title = sEXP1+' - Levitus Mediterranean  water at '+strtrim(lat,1)+'°N / year '+start_end
+      pltz, S1.arr-SLev.arr, MININ = -1., MAXIN = 1., INTER=0.1, UNIT = varunit $ 
+            , small = [1,2,2], COAST_THICK = 2, endpoints = [300., lat, 357., lat], TITLE = title, boxzoom = [2000.], ZOOM = 2000., _extra = ex 
+   endelse
+
+
+   title ='Mediterranean water Lat: '+strtrim(lat,1)+', '+strtrim(depth,1)
+   text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/024_Med_Levitus_S_depth_'+sEXP1+'_'+strtrim(lat,1)+'_d'+strtrim(depth,1)+'_S_y'+start_end+'.png  />  ' ]
+   putfile, '/tmp/prova.txt', text
+
+   if KEYWORD_SET(POSTSCRIPT) then begin
+      closeps
+   endif
+
+
+   return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Med_S_tongue.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Med_S_tongue.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Med_S_tongue.pro	(revision 2199)
@@ -0,0 +1,44 @@
+pro xxx_Med_S_tongue, S1, S2, SLev, sEXP1, sEXP2, start_end, DEPTH=depth, LAT=lat, POSTSCRIPT = postscript,  _extra = ex
+  
+  compile_opt idl2, strictarrsubs
+
+   @common
+
+
+   if KEYWORD_SET(POSTSCRIPT) then begin
+      openps, '023_Med_S_tongue_'+sEXP1+'_'+strtrim(lat,1)+'_d'+strtrim(depth,1)+'_S_y'+start_end+'.ps'   
+   endif
+;domdef, 0, 6000
+
+;domdef,270,365,5,70
+minval = min(abs(gdept -  depth),index)
+  
+   title = sEXP1+' Mediterranean salt tongue at depth=' + strtrim(gdept[index],1)+' m / year '+start_end
+   plt, S1.arr, MININ=35., MAXIN=37., INTER=0.1,  UNIT= varunit $
+         , small=[1,2,1], boxzoom = [270, 365, 5,70, index,index],/zindex, COAST_THICK=2, TITLE=title,  /PORTRAIT, _extra=ex
+;
+  if sEXP1 ne sEXP2 then begin
+     
+     title = sEXP1+' - '+sEXP2+' Mediterranean salt tongue at depth=' + strtrim(gdept[index],1)+' m / year '+start_end
+     plt, S1.arr - S2.arr, MININ = -.4, MAXIN = .4, INTER = .05,  UNIT = varunit $
+          , small=[1,2,2], boxzoom = [270, 365, 5,70, index,index],/zindex, COAST_THICK=2, TITLE=title,  /PORTRAIT, _extra=ex
+  
+  endif else begin
+     
+     title = sEXP1+' - Levitus Mediterranean salt tongue at depth=' + strtrim(gdept[index],1)+' m / year '+start_end
+     plt, S1.arr - SLev.arr, MININ=-1., MAXIN=1., INTER=0.1,  UNIT = varunit $
+          , small=[1,2,2], boxzoom = [270, 365, 5,70, index,index],/zindex, COAST_THICK=2, TITLE=title,  /PORTRAIT, _extra=ex
+     
+  endelse
+
+  title ='Mediterranean water Lat: '+strtrim(lat,1)+', '+strtrim(depth,1)
+  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/023_Med_S_tongue_'+sEXP1+'_'+strtrim(lat,1)+'_d'+strtrim(depth,1)+'_S_y'+start_end+'.png  />  ' ]
+  putfile, '/tmp/prova.txt', text
+   
+  if KEYWORD_SET(POSTSCRIPT) then begin
+      closeps
+   endif
+
+
+   return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_PHCXX_PHCYY_ANT_Icethick.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_PHCXX_PHCYY_ANT_Icethick.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_PHCXX_PHCYY_ANT_Icethick.pro	(revision 2199)
@@ -0,0 +1,44 @@
+pro xxx_PHCXX_PHCYY_ANT_Icethick,  I1, I2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+
+ @common
+
+ compile_opt idl2, strictarrsubs
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '0003_'+sEXP1+'_'+sEXP2+'_ANT_IceThick_y'+start_end+'.ps', portrait = 1
+ endif
+ 
+                                ;
+ domdef,20,380,-50,-90
+                                ; 
+ title=sEXP1+' Antartic Ice Thickness / year '+start_end
+
+ I1.arr = I1.arr < 1.
+ plt, I1.arr - 1.E-04, MIN = 0., MAX = 1., INTER = 0.1 , STYLE = 'so0so',CELL_FILL = 2, /STRICTFILL $ 
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title  $
+      , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [-90,0,0], /PORTRAIT, _extra = ex
+
+ if sEXP1 ne sEXP2 then begin
+    
+    title=sEXP1+' - '+sEXP2+' Antartic Ice Thickness / year '+start_end
+    plt, I1.arr - I2.arr, MIN = -1., MAX = 1., INTER = 0.1, /STRICTFILL , STYLE = 'so0so' $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title $
+         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [-90,0,0], /NOERASE, _extra = ex
+    
+ endif
+ 
+ domdef
+
+ title = 'Antarctic Sea-ice Thickness'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/0003_'+sEXP1+'_'+sEXP2+'_ANT_IceThick_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+
+ return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_PHCXX_PHCYY_ARC_Icethick.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_PHCXX_PHCYY_ARC_Icethick.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_PHCXX_PHCYY_ARC_Icethick.pro	(revision 2199)
@@ -0,0 +1,41 @@
+pro xxx_PHCXX_PHCYY_ARC_Icethick,  I1, I2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ @common
+
+ compile_opt idl2, strictarrsubs
+
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '0003_'+sEXP1+'_'+sEXP2+'_ARC_IceThick_y'+start_end+'.ps', portrait = 1
+ endif
+ 
+ domdef,20,380,55,90
+           
+ title=sEXP1+' Arctic Ice Thickness / year '+start_end
+
+ I1.arr = I1.arr < 6.
+ plt, I1.arr - 1.E-04, MIN = 0., MAX = 6., INTER = 0.25 ,/STRICTFILL , CELL_FILL = 2, STYLE = 'so0so'  $
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title $
+      , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [90,0,0], /PORTRAIT, _extra = ex
+ ;
+ if sEXP1 ne sEXP2 then begin            
+    title = sEXP1+' - '+sEXP2+' Arctic Ice Thickness / year '+start_end
+    plt, I1.arr - I2.arr, MIN = -3., MAX = 3., INTER = .5, /STRICTFILL , STYLE = 'so0so' $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title $
+         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = [90,0,0], /NOERASE, _extra = ex
+    
+ endif                          ; end of 2
+ 
+ domdef
+
+ title = 'Arctic Sea-ice Thickness'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/0003_'+sEXP1+'_'+sEXP2+'_ARC_IceThick_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Qnet.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Qnet.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_Qnet.pro	(revision 2199)
@@ -0,0 +1,43 @@
+pro xxx_Qnet,  Q1, Q2, QNET, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ @common
+
+ compile_opt idl2, strictarrsubs
+
+ Q = Q1.arr-Q2.arr
+ sEXP = sEXP1+' - '+sEXP2
+ min =  -200.
+ max = -min
+ inter = 20.
+ 
+ if sEXP2 eq sEXP1 then begin
+    sEXP = sEXP1+' - OAFlux'
+    Q = Q1.arr-QNET.arr
+    min = -200.
+    max = -min
+    inter = 20. 
+ endif
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '003_Qnet_'+sEXP+'_y'+start_end+'.ps', portrait = 1
+ endif
+ 
+ title=sEXP1+' / Qnet year '+start_end
+ plt,  Q1,  MIN = -200., MAX = 200., INTER = 20., STYLE = 'so0so' $
+       , small = [1,2,1], COAST_THICK = 2, TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex
+                                ;
+ title=sEXP+' / Qnet year '+start_end
+ plt, Q, MIN = min, MAX = max, INTER = inter, STYLE = 'so0so'  $
+      , small = [1,2,2], COAST_THICK = 2, TITLE = title, FORMAT = '(I3)', /NOERASE, _extra = ex
+                                ;
+ title = 'Qnet differences with '+sEXP
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/003_Qnet_'+sEXP+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+ return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_SSS.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_SSS.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_SSS.pro	(revision 2199)
@@ -0,0 +1,40 @@
+pro xxx_SSS, S1, S2, SLev, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ @common
+ compile_opt idl2, strictarrsubs
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '007_SSS_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+ endif 
+      
+ varunit = S1.unit  
+ 
+ title = sEXP1+' / SSS year '+start_end
+ plt, S1.arr[*,*,0] , MIN = 33., MAX = 41., INTER = .25, STYLE = 'so0so' $
+      , small = [1,2,1], COAST_THICK = 2, TITLE=title, /NOERASE, _extra = ex 
+                                ;
+ if sEXP1 ne sEXP2 then begin
+    
+    title = sEXP1+' - '+sEXP2+' / SSS year '+start_end
+    plt, S1.arr[*,*,0] - S2.arr[*,*,0], MIN = -2., MAX = 2., INTER = .2 $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title, _extra=ex
+
+ endif else begin
+
+    title = sEXP1+' - Levitus / SSS year '+start_end
+    plt, S1.arr[*,*,0]-SLev.arr[*,*,0], MIN = -3., MAX = 3., INTER = 0.2, STYLE = 'so0so' $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex
+
+ endelse
+ 
+                                ;
+   title = 'SSS differences with Levitus and '+sEXP1+'-'+sEXP2
+   text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/007_SSS_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.png  />  ' ]
+   putfile, '/tmp/prova.txt', text
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif 
+
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_SST.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_SST.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_SST.pro	(revision 2199)
@@ -0,0 +1,44 @@
+pro xxx_SST, T1, T2, TRey, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+
+ @common 
+ compile_opt idl2, strictarrsubs
+
+ T = T1.arr[*,*,0]-T2.arr[*,*,0]
+ sEXP = sEXP1+' - '+sEXP2
+ min =  -1.
+ max = -min
+ inter = 0.1
+
+ if sEXP2 eq sEXP1 then begin
+    sEXP = sEXP1
+    T = T1.arr[*,*,0]
+    min = -2.
+    max = 32.
+    inter = 1. 
+ endif
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '006_SST_'+sEXP+'_y'+start_end+'.ps', portrait = 1
+ endif
+   
+ varunit = T1.unit
+
+                               ; 
+ title = sEXP+' / SST year '+start_end
+ plt, T , MIN = min, MAX = max, INTER = inter ,  STYLE = 'so0so' $
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title, FORMAT = '(I2)', /NOERASE, _extra = ex
+                                ; 
+ title=sEXP1+' - NewReynolds  / SST year '+start_end
+ plt, T1.arr[*,*,0]-TRey.arr[*,*,0], MIN = -8., MAX = 8., INTER = 0.5 , STYLE = 'so0so' $
+      , small = [1,2,2], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex
+                 
+ title = 'SST differences with Reynolds and  '+sEXP
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/006_SST_'+sEXP+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_S_z105.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_S_z105.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_S_z105.pro	(revision 2199)
@@ -0,0 +1,44 @@
+pro xxx_S_z105, S1, S2, SLev, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+
+ @common
+ compile_opt idl2, strictarrsubs
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '011_S_z105_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+ endif
+ 
+;
+ varunit = S1.unit
+;
+
+ title=sEXP1+' S (PSU) depth=105m / Sal year '+start_end
+ plt, S1.arr, MIN = 33., MAX = 38.2, INTER = .2, STYLE = 'so0so' $
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title, boxzoom = [105., 106.], /NOERASE, _extra = ex
+;
+ if sEXP1 ne sEXP2 then begin
+    
+    title = sEXP1+' - '+sEXP2+' S (PSU) depth=105m / Sal year '+start_end
+    plt, S1.arr[*,*,0] - S2.arr[*,*,0], MIN = -1., MAX = 1., INTER = .1 $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title, _extra=ex
+    
+ endif else begin
+    
+    title=sEXP1+' - Levitus S (PSU) depth=105m / Sal year '+start_end
+    plt, S1.arr-SLev.arr, MIN = -1., MAX = 1., INTER=0.1, STYLE='so0so' $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = [105., 106.], /PORTRAIT, _extra = ex                           
+    
+ endelse
+ 
+
+ 
+ title = 'Salinity differences at z=105 m with Levitus and '+sEXP2
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/011_S_z105_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text ;
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_TS_ZGlobm.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_TS_ZGlobm.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_TS_ZGlobm.pro	(revision 2199)
@@ -0,0 +1,36 @@
+pro xxx_TS_ZGlobm, T1, T2, TLev, S1, S2, SLev, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+
+ @common
+ compile_opt idl2, strictarrsubs
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '025_TS_ZGlobmean_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', /portrait
+ endif
+
+                                ;
+ title=sEXP1+' (Black) & '+sEXP2+' (Red) - Levitus Temp. Global mean / year '+start_end
+ plt1d, T1.arr - TLev.arr, typein = 'z', ticklen = 1 , MIN = -2., MAX = 2., boxzoom = [4.,5300.], /KEEPBOTTOM $
+        , small = [1,2,1], XGRIDSTYLE = 2, YGRIDSTYLE = 2, TITLE = title, _extra = ex
+ plt1d, T2.arr - TLev.arr, typein = 'z',  ticklen = 1 , MIN = -2., MAX = 2., boxzoom = [4.,5300.], /KEEPBOTTOM  $
+        , /ov1d, COLOR = 250,TITLE = title, _extra = ex 
+                                ;
+
+ title=sEXP1+'(Black) & '+sEXP2+' (Red) - Levitus Sal. Global mean / year '+start_end
+ plt1d, S1.arr - SLev.arr, typein = 'z', ticklen = 1 , MIN = -.2, MAX = .2, boxzoom = [4.,5300.], /KEEPBOTTOM $
+        , small = [1,2,2], XGRIDSTYLE = 2, YGRIDSTYLE = 2, TITLE = title, /NOERASE, _extra = ex
+ plt1d, S2.arr - SLev.arr, typein = 'z',  ticklen = 1 , MIN = -2., MAX = 2., boxzoom = [4.,5300.], /KEEPBOTTOM $
+        ,  /ov1d, COLOR = 250, TITLE = title, _extra = ex 
+
+
+ title = 'Temperature & Salinity horizontal mean'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/025_TS_ZGlobmean_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+ return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_T_z105.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_T_z105.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_T_z105.pro	(revision 2199)
@@ -0,0 +1,41 @@
+pro xxx_T_z105,  T1, T2, TLev, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ @common 
+ compile_opt idl2, strictarrsubs
+
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '010_T_z105_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+ endif
+;
+ varunit = T1.unit     
+           ;                              
+ title = sEXP1+' T (°C) depth=105m / Temp year '+start_end
+ plt, T1.arr, MIN = -2., MAX = 32., INTER = 1., STYLE = 'so0so'  $
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title, FORMAT = '(I2)', boxzoom = [105., 106.], /NOERASE, _extra = ex
+                                ;
+if sEXP1 ne sEXP2 then begin
+
+  title = sEXP1+' - '+sEXP2+' T (°C) depth=105m / Temp year '+start_end
+  plt, T1.arr-T2.arr, MIN = -2., MAX = 2., INTER = .2 $
+       , small = [1,2,2], COAST_THICK = 2, TITLE = title,  _extra=ex
+
+endif else begin
+
+ title = sEXP1+' - Levitus T (°C) depth=105m / Temp year '+start_end
+ plt, T1.arr-TLev.arr, MIN = -4., MAX = 4., INTER = .5, STYLE = 'so0so'  $
+      , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = [105., 106.], /PORTRAIT, _extra = ex   
+endelse
+                    
+
+ title = 'Temperature differences at z=105 m with Levitus and '+sEXP1+'-'+sEXP2
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/010_T_z105_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text ;
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+ 
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ZonMxl.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ZonMxl.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_ZonMxl.pro	(revision 2199)
@@ -0,0 +1,34 @@
+pro xxx_ZonMxl, MXL1, MXL2, MLD, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+  
+ @common
+
+ compile_opt idl2, strictarrsubs
+
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '013_ZonMxl_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+ endif
+ 
+
+                                ;
+ title=sEXP1+' (Black) , '+sEXP2+' (Red) & Clim (Blue) Mxl10 / year '+start_end
+ plt1d, -MXL1.arr, MIN = -300., MAX = 0., INTER = 10., typein = 'y' $
+        , small=[1,1,1], boxzoom=5500, ZOOM=500, CHARSIZE=.8, TITLE = title, /PORTRAIT, _extra=ex  ; color by default = black
+                                ;
+ plt1d, -MXL2.arr, MIN = -300., MAX = 0., INTER = 10., typein = 'y' $
+        , /ov1d, COLOR = 250, TITLE = title, /NOERASE, _extra = ex  ; color 250 = red
+                                ;
+ plt1d, -MLD.arr, MIN = -300., MAX = 0., INTER = 10., typein = 'y' $
+        ,  /ov1d, COLOR = 100, THICK = 3, TITLE = title, /NOERASE, _extra = ex ; color 100 = blue
+
+ title = 'Mxl field'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/013_ZonMxl_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+ return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_mxl10.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_mxl10.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_mxl10.pro	(revision 2199)
@@ -0,0 +1,46 @@
+pro xxx_mxl10, MXL1, MXL2, MLD, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex
+
+ @common
+
+ compile_opt idl2, strictarrsubs
+
+ MXL = MXL1.arr-MXL2.arr
+ sEXP = sEXP1+' - '+sEXP2
+ min =  -100.
+ max = -min
+ inter = 10.
+ 
+ if sEXP2 eq sEXP1 then begin
+    sEXP = sEXP1+' - MLD '
+    MXL = MXL1.arr-MLD.arr
+    min = -80.
+    max = -min
+    inter = 10. 
+ endif
+
+ 
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    openps, '012_Mxl10_'+sEXP+'_y'+start_end+'.ps', portrait = 1
+ endif
+ 
+ title=sEXP1+' Mxl10 / year '+start_end
+ plt, MXL1, MIN = 0., MAX = 500., INTER = 25., /NOCONTOUR ,FORMAT = '(I3)' $
+      , small = [1,2,1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex
+ ; 
+
+ title=sEXP+' Mxl10 / year '+start_end
+ plt, MXL, MIN = min, MAX = max, INTER = inter  $
+      , small = [1,2,2], COAST_THICK = 2, TITLE = title, /NOCONTOUR, /NOERASE, _extra = ex
+                                ; end of 2
+ title = 'Mxl10 field'
+ text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/012_Mxl10_'+sEXP+'_y'+start_end+'.png  />  ' ]
+ putfile, '/tmp/prova.txt', text
+
+
+ if KEYWORD_SET(POSTSCRIPT) then begin
+    closeps
+ endif
+
+ return
+end
+
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_zonal_S.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_zonal_S.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_zonal_S.pro	(revision 2199)
@@ -0,0 +1,70 @@
+pro xxx_zonal_S, S1, S2, SLev, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, SUBBASIN = subbasin, MASK_FILENAME = mask_filename, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+
+   @common
+
+   if KEYWORD_SET(SUBBASIN) then  subname = subbasin else  subname = "Glo"
+   
+
+   if KEYWORD_SET(POSTSCRIPT) then begin
+      openps, '015_Zonal_S_'+subname+'_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+   endif
+
+                                ; looking for longitudinal index corresponding to the highest latitude (closest to the North Pole)
+   index=where(gphit eq max(gphit))
+   nx=index[0] mod jpi
+   gphi_save=gphit
+   gphit[0,*]=gphit[nx,*]
+
+   if KEYWORD_SET(SUBBASIN) then begin
+                                ; save tmask
+      tmasksv = tmask
+                                ; 
+      CASE 1 of 
+         subbasin eq 'Atl' : var = 'atlmsk_nomed'
+         subbasin eq 'Ind' : var = 'indmsk_nored'
+         subbasin eq 'Pac' : var = 'pacmsk'
+      ENDCASE
+                                ;
+      msk = read_ncdf( var, filename = mask_filename, _extra =ex)
+
+; from 2D array to 3D array
+      msk = msk.arr[*]#replicate(1., nzt)
+      tmask[firstxt:lastxt,firstyt:lastyt,firstzt:lastzt] = tmask[firstxt:lastxt,firstyt:lastyt,firstzt:lastzt] * msk
+   endif 
+
+   varunit = S1.unit
+;
+   title=sEXP1+' '+subname+' / S year '+start_end
+   pltz, S1.arr, MININ = 32., MAXIN = 37., INTER = .2, typein = 'yz' $ 
+         , small = [1,2,1], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /PORTRAIT,  _extra=ex 
+; 
+   if sEXP1 ne sEXP2 then begin
+      
+      title = sEXP1+' - '+sEXP2+' '+subname+' / S year '+start_end
+      pltz, S1.arr-S2.arr, MIN = -1., MAX = 1., INTER = .1, typein = 'yz' $
+           , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /NOERASE,  _extra=ex
+  
+   endif else begin
+      
+      title=sEXP1+' - Levitus '+subname+'/ S year'+start_end
+      pltz, S1.arr-SLev.arr, MININ = -1, MAXIN = 1., INTER = .1, typein = 'yz'  $
+            , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /NOERASE, _extra = ex
+   endelse
+;
+   title = 'Salinity differences in zonal average ' + subname +' with Levitus and '+sEXP1+'-'+sEXP2
+   text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/015_Zonal_S_'+subname+'_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.png  />  ' ]
+   putfile, '/tmp/prova.txt', text
+
+   if KEYWORD_SET(POSTSCRIPT) then begin
+      closeps
+   endif
+
+
+   if KEYWORD_SET(SUBBASIN) then tmask = tmasksv
+
+   gphit=gphi_save
+
+   return
+end
Index: branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_zonal_T.pro
===================================================================
--- branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_zonal_T.pro	(revision 2199)
+++ branches/devukmo2010/CONFIG/ORCA2_LIM/IDL_scripts/xxx_zonal_T.pro	(revision 2199)
@@ -0,0 +1,77 @@
+pro xxx_zonal_T, T1, T2, TLev, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, SUBBASIN = subbasin, MASK_FILENAME = mask_filename, _extra = ex
+  
+  compile_opt idl2, strictarrsubs
+
+   @common
+
+   if KEYWORD_SET(SUBBASIN) then subname = subbasin else subname = "Glo"
+
+
+   if KEYWORD_SET(POSTSCRIPT) then begin
+      openps, '014_Zonal_T_'+subname+'_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.ps', portrait = 1
+   endif
+
+                                ; looking for longitudinal index corresponding to the highest latitude (closest to the North Pole)
+   index=where(gphit eq max(gphit))
+   nx=index[0] mod jpi
+   gphi_save=gphit
+   gphit[0,*]=gphit[nx,*]
+
+   if KEYWORD_SET(SUBBASIN) then begin
+ ; save tmask
+      tmasksv = tmask
+ ; 
+      CASE 1 of 
+         subbasin eq 'Atl' : var = 'atlmsk_nomed'
+         subbasin eq 'Ind' : var = 'indmsk_nored'
+         subbasin eq 'Pac' : var = 'pacmsk'
+      ENDCASE
+                                ;
+      msk = read_ncdf( var, filename = mask_filename, _extra =ex)
+
+; from 2D array to 3D array
+      msk = msk.arr[*]#replicate(1., nzt)
+      tmask[firstxt:lastxt,firstyt:lastyt,firstzt:lastzt] = tmask[firstxt:lastxt,firstyt:lastyt,firstzt:lastzt] * msk
+
+   endif
+
+varunit = T1.unit
+;
+ 
+   title=sEXP1+' '+subname+' / T year '+start_end
+   pltz, T1.arr, MININ = -2., MAXIN = 30., INTER = 1., typein = 'yz', FORMAT = '(I2)' $ 
+         , small =[ 1,2,1], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /PORTRAIT, _extra = ex 
+; 
+if sEXP1 ne sEXP2 then begin
+
+  title = sEXP1+' - '+sEXP2+' T (°C) depth=105m / Temp year '+start_end
+  pltz, T1.arr-T2.arr, MIN = -2., MAX = 2., INTER = .2, typein = 'yz'  $
+       , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /PORTRAIT, _extra=ex
+
+endif else begin
+
+   title=sEXP1+' - Levitus '+subname+'/ T year'+start_end
+   pltz, T1.arr-TLev.arr, MININ = -4., MAXIN = 4., INTER = .5, typein='yz', FORMAT = '(I2)' $
+         , small = [1,2,2], COAST_THICK = 2, TITLE = title, boxzoom = 5500, ZOOM = 1000, /NOERASE, _extra = ex
+   
+endelse
+         
+   
+;
+   title = 'Temperature differences in zonal average ' + subname+' with Levitus and '+sEXP1+'-'+sEXP2
+   text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/v32colib_v32colib/ps/014_Zonal_T_'+subname+'_'+sEXP1+'-'+sEXP2+'_y'+start_end+'.png  />  ' ]
+   putfile, '/tmp/prova.txt', text
+
+
+   if KEYWORD_SET(POSTSCRIPT) then begin
+      closeps
+   endif
+
+   if KEYWORD_SET(SUBBASIN) then tmask = tmasksv
+   
+  
+   gphit=gphi_save
+
+
+   return
+end
Index: branches/devukmo2010/CONFIG/POMME/EXP00/namelist
===================================================================
--- branches/devukmo2010/CONFIG/POMME/EXP00/namelist	(revision 2189)
+++ branches/devukmo2010/CONFIG/POMME/EXP00/namelist	(revision 2199)
@@ -178,4 +178,5 @@
    sn_prec     = 'precip_1988.nc' ,       -1          , 'precip'   ,    .true.      , .true. ,   'yearly' , ''
    sn_snow     = 'precip_1988.nc' ,       -1          , 'precip'   ,    .true.      , .true. ,   'yearly' , ''
+   sn_tdif     = 'taudif_core'    ,       24          , 'taudif'   ,    .true.      , .true.  ,   'yearly'  ,''
 !
    cn_dir      = './'      !  root directory for the location of the bulk files
@@ -318,7 +319,7 @@
     rn_dpsin   =    1.     !     -           -         -       south   -      -
     rn_dpeob   = 1500.     !  time relaxation (days) for the east  open boundary
-    rn_dpwob   =   15.     !     -           -         -     west    -      -
+    rn_dpwob   =  150.     !     -           -         -     west    -      -
     rn_dpnob   =  150.     !     -           -         -     north   -      -
-    rn_dpsob   =   15.     !     -           -         -     south   -      -
+    rn_dpsob   =  150.     !     -           -         -     south   -      -
     rn_volemp  =  1.       !  = 0 the total volume change with the surface flux (E-P-R)
                            !  = 1 the total volume remains constant
