Ignore:
Timestamp:
2013-11-28T12:41:27+01:00 (7 years ago)
Author:
flavoni
Message:

update IDL_scripts for ORCA2_LIM and ORCA2_LIM3, add plot of max_mld output 5days, see ticket: #724

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceThick.pro

    r2751 r4320  
    1 pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    1212  IF keyword_set(arc) THEN var = var+'_Arc_' 
    1313  IF keyword_set(ant) THEN var = var+'_Ant_' 
    14   IF keyword_set(march) THEN var = var+'March' 
     14  IF keyword_set(april) THEN var = var+'April' 
    1515  IF keyword_set(sept) THEN var = var+'Sept' 
    1616 
    1717  filename = cdti3 + '_'+var+'_'+std_file1_I 
    18   if std_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I 
    19   if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
     18  IF std_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I 
     19  IF KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
    2020   
    2121  IF keyword_set(arc) THEN BEGIN  
    2222    domdef, 20, 380, 50, 90 
    2323    vmin = 0. 
    24     vmax = 10. 
     24    vmax = 6. 
    2525    vint = 0.5 
    26     fmt = '(i2)' 
     26    fmtd = '(f4.1)' 
    2727    vmind = -3. 
    2828    vmaxd = 3. 
    2929    vintd = 0.5 
    30     fmtd = '(i2)' 
     30    fmtd = '(f4.1)' 
     31    div = 6 
    3132    map = [90, 0, 0] 
    3233  ENDIF 
     
    4142    vintd = 0.2 
    4243    fmtd = '(f4.1)' 
     44    div = 5 
    4345    map = [-90, 0, 0] 
    4446  ENDIF 
     
    4850  title = var+'!C'+std_file1_I 
    4951  plt, (Ithi1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $ 
    50        , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ 
     52       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, divisions = div $ 
    5153       , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex 
    5254;                                ; 
     
    5456    title = title + std_file2_I 
    5557    plt, Ithi1.arr - Ithi2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 
    56          , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 
     58         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title  $ 
    5759         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 
    5860  endif                         
Note: See TracChangeset for help on using the changeset viewer.