Ignore:
Timestamp:
2011-04-26T10:45:52+02:00 (10 years ago)
Author:
smasson
Message:

large update of IDL_scripts for Is-ENES

File:
1 moved

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_ArcSal.pro

    r2733 r2751  
    1 pro xxx_EXP1_S_ARC, S1, SLev, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_ArcSal, S1, S2in, SLevin, Z100 = z100, POSTSCRIPT = postscript, _extra = ex 
     2 
     3  compile_opt idl2, strictarrsubs 
     4 
     5@common   
     6@std_common 
     7 
     8  CASE n_params() OF 
     9    2:BEGIN 
     10      Slev = S2in 
     11    END 
     12    3:BEGIN 
     13      IF S2in.arr[0] EQ -1 THEN return 
     14      S2 = S2in 
     15      Slev = Slevin 
     16    END 
     17  ENDCASE 
     18; 
     19  cdti3 = string(cnt, format = '(i3.3)') 
     20  print, cdti3 + ') ' + blabla 
     21; 
     22  IF keyword_set(z100) THEN tmp = min(abs(100 - gdept), ind) ELSE ind = 0 
     23  sdepref = strtrim(round(gdept[ind]), 1)+'m' 
     24  filename = cdti3 + '_Arctic_Sal'+sdepref+'_'+std_file1_T 
     25  IF keyword_set(S2) THEN filename = filename + '_' + std_file2_T 
     26  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
     27;  
     28  domdef, 20, 380, 60, 90 
     29; 
     30  varunit = S1.unit 
     31  titleorg = 'Salinity ('+sdepref+')!C' 
     32; 
     33  IF keyword_set(S2) THEN BEGIN 
     34    title = titleorg+std_file1_T+ ' - '+std_file2_T 
     35    plt, S1.arr - S2.arr, MIN = -4., MAX = 4., INTER = 0.2, CELL_FILL = 2,  STYLE = 'so0so', format = '(f4.1)'  $      
     36         , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $       
     37         , /ORTHO, MAP = [90, 0, 0], LATDEL = 5, boxzoom = [ind, ind], /zindex, /portrait, _extra = ex 
     38  ENDIF ELSE BEGIN  
     39    title = titleorg+std_file1_T 
     40    plt, S1, MIN = 30.0, MAX = 36., INTER = 0.2, CELL_FILL = 2, format = '(f4.1)' $ 
     41         , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $ 
     42         , /ORTHO, MAP = [90, 0, 0], LATDEL = 5, boxzoom = [ind, ind], /zindex, /PORTRAIT, _extra = ex 
     43  ENDELSE  
     44;  
     45  IF keyword_set(S2) THEN BEGIN 
     46    title = titleorg+std_file2_T+ ' - Levitus' 
     47    tmp = S2.arr - SLev.arr 
     48  ENDIF ELSE BEGIN  
     49    title = titleorg+std_file1_T+ ' - Levitus' 
     50    tmp = S1.arr - SLev.arr 
     51  ENDELSE 
     52  plt, temporary(tmp), MIN = -4., MAX = 4., INTER = 0.2, CELL_FILL = 2,  STYLE = 'so0so', format = '(f4.1)'  $      
     53       , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $       
     54       , /ORTHO, MAP = [90, 0, 0], LATDEL = 5, boxzoom = [ind, ind], /zindex, /NOERASE, _extra = ex 
     55; 
     56  domdef 
     57; 
     58  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
     59  if KEYWORD_SET(postscript) then closeps 
    260   
    3  @common  
    4  compile_opt idl2, strictarrsubs 
    5  
    6  
    7  if KEYWORD_SET(POSTSCRIPT) then begin 
    8     openps, '008_S_ARC_'+sEXP1+'_y'+start_end+'.ps', portrait = 1 
    9  endif 
    10  
    11                                 ;  
    12  domdef,20,380,60,90 
    13                                 ; 
    14  title=sEXP1+' Salinity SSS' 
    15  plt, S1, MIN = 30.0, MAX = 36., INTER = 0.2, CELL_FILL = 2  $ 
    16       , small = [1,2,1], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $ 
    17       , /ORTHO, MAP = [90,0,0], LATDEL = 5, boxzoom = [4., 6.], /PORTRAIT, _extra = ex 
    18                                 ;  
    19  
    20  title=sEXP1+' - Levitus Salinity SSS m / year '+start_end 
    21  plt, S1.arr-SLev.arr, MIN = -4., MAX = 4., INTER = 0.2, CELL_FILL = 2,  STYLE = 'so0so'  $      
    22       , small = [1,2,2], COAST_THICK = 2, TITLE = title, CHARSIZE = -0.55, GLINETHICK = 2. $       
    23       , /ORTHO, MAP = [90,0,0], LATDEL = 5, boxzoom = [4., 6.], /NOERASE, _extra = ex 
    24                                 ; 
    25  domdef 
    26  
    27  title = 'Artic salinity at SSS meters' 
    28  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/'+sEXP1+'_'+sEXP2+'/ps/008_S_ARC_'+sEXP1+'_y'+start_end+'.png  />  ' ] 
    29  putfile, '/tmp/prova.txt', text 
    30  
    31  
    32  if KEYWORD_SET(POSTSCRIPT) then begin 
    33     closeps 
    34  endif 
    35  
    36   
    37  return 
     61  return 
    3862end 
Note: See TracChangeset for help on using the changeset viewer.