New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2751 for trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_sst.pro – NEMO

Ignore:
Timestamp:
2011-04-26T10:45:52+02:00 (13 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_sst.pro

    r2733 r2751  
    1 pro xxx_SST, T1, T2, TRey, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_sst, T1, T2, TRey, POSTSCRIPT = postscript, _extra = ex 
    22 
    3  @common  
    4  compile_opt idl2, strictarrsubs 
     3  compile_opt idl2, strictarrsubs 
    54 
    6  T = T1.arr[*,*,0]-T2.arr[*,*,0] 
    7  sEXP = sEXP1+' - '+sEXP2 
    8  min =  -1. 
    9  max = -min 
    10  inter = 0.1 
     5@common   
     6@std_common 
    117 
    12  if sEXP2 eq sEXP1 then begin 
    13     sEXP = sEXP1 
    14     T = T1.arr[*,*,0] 
     8  cdti3 = string(cnt, format = '(i3.3)') 
     9  print, cdti3 + ') ' + blabla 
     10  filename = cdti3 + '_SST_'+std_file1_T 
     11  if std_file1_T NE std_file2_T then filename = filename + '_'+std_file2_T 
     12  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
     13 
     14  if std_file1_T EQ std_file2_T then begin 
     15    T = T1.arr[*, *, 0] 
    1516    min = -2. 
    1617    max = 32. 
    1718    inter = 1.  
    18  endif 
     19    STYLE = 0 
     20    title = 'SST!C'+std_file1_T 
     21  ENDIF ELSE BEGIN  
     22    T = T1.arr[*, *, 0] - T2.arr[*, *, 0] 
     23    min =  -1. 
     24    max = -min 
     25    inter = 0.1 
     26    STYLE = 'so0so' 
     27    title = 'SST!C'+std_file1_T+' - '+std_file2_T 
     28  ENDELSE  
     29; 
     30  varunit = T1.unit 
     31;  
     32  plt, T, MIN = min, MAX = max, INTER = inter, STYLE = STYLE $ 
     33       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, FORMAT = '(I2)', /PORTRAIT, _extra = ex 
     34;  
     35  title = 'SST!C'+std_file1_T+' - NewReynolds' 
     36  plt, T1.arr[*, *, 0] - TRey.arr, MIN = -8., MAX = 8., INTER = 0.5, STYLE = 'so0so', FORMAT = '(I2)' $ 
     37       , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, /NOERASE, _extra = ex 
     38   
     39  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
     40  if KEYWORD_SET(postscript) then closeps 
    1941 
    20  if KEYWORD_SET(POSTSCRIPT) then begin 
    21     openps, '006_SST_'+sEXP+'_y'+start_end+'.ps', portrait = 1 
    22  endif 
    23     
    24  varunit = T1.unit 
    25  
    26                                ;  
    27  title = sEXP+' / SST year '+start_end 
    28  plt, T , MIN = min, MAX = max, INTER = inter ,  STYLE = 'so0so' $ 
    29       , small = [1,2,1], COAST_THICK = 2, TITLE = title, FORMAT = '(I2)', /NOERASE, _extra = ex 
    30                                 ;  
    31  title=sEXP1+' - NewReynolds  / SST year '+start_end 
    32  plt, T1.arr[*,*,0]-TRey.arr[*,*,0], MIN = -8., MAX = 8., INTER = 0.5 , STYLE = 'so0so' $ 
    33       , small = [1,2,2], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 
    34                   
    35  title = 'SST differences with Reynolds and  '+sEXP 
    36  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/'+sEXP1+'_'+sEXP2+'/ps/006_SST_'+sEXP+'_y'+start_end+'.png  />  ' ] 
    37  putfile, '/tmp/prova.txt', text 
    38  
    39  if KEYWORD_SET(POSTSCRIPT) then begin 
    40     closeps 
    41  endif 
    42  
    43  return 
     42  return 
    4443end 
Note: See TracChangeset for help on using the changeset viewer.