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 4559 for trunk/NEMOGCM/CONFIG/ORCA2_LIM – NEMO

Ignore:
Timestamp:
2014-03-21T10:22:46+01:00 (10 years ago)
Author:
flavoni
Message:

add plot & ts of ice velocity for orca2 idl plots, see ticket #724

Location:
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceAge.pro

    r4519 r4559  
    1 pro std_plot_IceAge, Iage1, Iage2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_IceAge, Iage1, Iage2, Ifra1, Ifra2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
    … …  
    5050; 
    5151  title = var+'!C'+std_file1_I 
     52   
     53  Iage1.arr = Iage1.arr * ( Ifra1.arr gt 0.15 ) 
     54  
    5255  plt, Iage1.arr - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $ 
    5356       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, divisions = div  $ 
    … …  
    5659  if std_file1_I NE std_file2_I then begin             
    5760    title = title + std_file2_I 
     61  
     62    Iage2.arr = Iage2.arr * ( Ifra2.arr gt 0.15 ) 
     63   
    5864    plt, Iage1.arr - Iage2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 
    5965         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title  $ 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro

    r4518 r4559  
    3232; 
    3333  title = var+'!C'+std_file1_I 
    34   plt, (Ifra1.arr < 1. ) - 1.E-04, MIN = 0., MAX = 1., INTER = .1, /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 
     34  Ifra1.arr = Ifra1.arr * ( Ifra1.arr gt 0.15 ) 
     35  plt, Ifra1.arr - 1.E-04, MIN = 0., MAX = 1., INTER = .1, /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 
    3536       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ 
    3637       , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex 
    … …  
    3839  if std_file1_I NE std_file2_I then begin             
    3940    title = title + std_file2_I 
     41    Ifra2.arr = Ifra2.arr * ( Ifra2.arr gt 0.15 ) 
    4042    plt, Ifra1.arr - Ifra2.arr, MIN = -1, MAX = 1, INTER = .1, STYLE = 'so0so', format = '(f4.1)' $ 
    4143         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceSal.pro

    r4519 r4559  
    1 pro std_plot_IceSal, Isal1, Isal2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_IceSal, Isal1, Isal2, Ifra1, Ifra2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
    … …  
    3232  ; 
    3333   title = var+'!C'+std_file1_I 
    34    plt, Isal1.arr - 1.E-04, MIN = 0., MAX = 12., INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 
     34 
     35   Isal1.arr = Isal1.arr * ( Ifra1.arr gt 0.15 ) 
     36   plt, Isal1.arr - 1.E-04, MIN = 0., MAX = 10., INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 
    3537        , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ 
    3638        , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex 
    … …  
    3840   if std_file1_I NE std_file2_I then begin 
    3941     title = title + std_file2_I 
    40      plt, Isal1.arr - Isal2.arr, MIN = -4, MAX = 4, INTER = 1., STYLE = 'so0so', format = '(f4.1)' $ 
     42      
     43     Isal2.arr = Isal2.arr * ( Ifra2.arr gt 0.15 ) 
     44      
     45     plt, Isal1.arr - Isal2.arr, MIN = -2., MAX = 2., INTER = 1., STYLE = 'so0so', format = '(f4.1)' $ 
    4146          , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 
    4247          , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceThick.pro

    r4518 r4559  
    1 pro std_plot_IceThick, Ithi1, Ithi2, ARC = arc, ANT = ant, APRIL = april, JAN = jan, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_IceThick, Ithi1, Ithi2, Ifra1, Ifra2, ARC = arc, ANT = ant, APRIL = april, JAN = jan, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
    … …  
    5050; 
    5151  title = var+'!C'+std_file1_I 
     52 
     53  Ithi1.arr = Ithi1.arr * ( Ifra1.arr gt 0.15 ) 
     54  
    5255  plt, (Ithi1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $ 
    5356       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, divisions = div  $ 
    … …  
    5659  if std_file1_I NE std_file2_I then begin             
    5760    title = title + std_file2_I 
     61 
     62    Ithi2.arr = Ithi2.arr * ( Ifra2.arr gt 0.15 ) 
     63 
    5864    plt, Ithi1.arr - Ithi2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 
    5965         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title  $ 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro

    r4518 r4559  
    1 pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_SnowThick, Isnow1, Isnow2, Ifra1, Ifra2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
    … …  
    5050; 
    5151  title = var+'!C'+std_file1_I 
     52   
     53  Isnow1.arr = Isnow1.arr * ( Ifra1.arr gt 0.15 ) 
     54   
    5255  plt, (Isnow1.arr < 10. ) - 1.E-04, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $ 
    5356       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ 
    … …  
    5659  if std_file1_I NE std_file2_I then begin             
    5760    title = title + std_file2_I 
     61    
     62    Isnow2.arr = Isnow2.arr * ( Ifra2.arr gt 0.15 ) 
     63      
    5864    plt, Isnow1.arr - Isnow2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 
    5965         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r4518 r4559  
    3636; 
    3737  allrec =  1 - keyword_set(long(getenv('READ_ONLY_FIRST_RECORD'))) 
     38   
    3839; Output run experience1 
    3940  std_file1_T     = isafile(getenv('FILE1_T'), title = 'exp1 grid T input file', iodir = std_iodir_data) 
    … …  
    192193  ENDELSE 
    193194; 
    194   Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
    195   caldat, time, mm 
    196   febr = where(mm EQ 2, cnt) 
    197   Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
    198   march = where(mm EQ 3, cnt) 
    199   Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
    200   sept = where(mm EQ 9, cnt) 
    201   Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit}  
    202   undefine, Iage_1 
    203 ; 
    204   IF std_file2_I NE std_file1_I THEN BEGIN 
    205      Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
    206      caldat, time, mm 
    207      febr = where(mm EQ 2, cnt) 
    208      Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
    209      march = where(mm EQ 3, cnt) 
    210      Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
    211      sept = where(mm EQ 9, cnt) 
    212      Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
    213      undefine, Iage_2 
    214   ENDIF ELSE BEGIN  
    215      Iage_febr_2 = {arr:-1} 
    216      Iage_march_2 = {arr:-1} 
    217      Iage_sept_2 = {arr:-1} 
    218   ENDELSE 
    219 ; 
    220195  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    221196  caldat, time, mm 
     197  jan = where(mm EQ 1, cnt) 
     198  Ifra_jan_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
    222199  febr = where(mm EQ 2, cnt) 
    223200  Ifra_febr_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
    224201  march = where(mm EQ 3, cnt) 
    225202  Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
     203  april = where(mm EQ 4, cnt) 
     204  Ifra_april_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ifra_1.unit}  
    226205  sept = where(mm EQ 9, cnt) 
    227206  Ifra_sept_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 
    … …  
    231210    Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )  
    232211    caldat, time, mm 
     212    jan = where(mm EQ 1, cnt) 
     213    Ifra_jan_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
    233214    febr = where(mm EQ 2, cnt) 
    234215    Ifra_febr_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
    235216    march = where(mm EQ 3, cnt) 
    236217    Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
     218    april = where(mm EQ 4, cnt) 
     219    Ifra_april_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
    237220    sept = where(mm EQ 9, cnt) 
    238221    Ifra_sept_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 
    239222    undefine, Ifra_2 
    240223  ENDIF ELSE BEGIN  
     224    Ifra_jan_2 = {arr:-1}  
    241225    Ifra_febr_2 = {arr:-1} 
    242226    Ifra_march_2 = {arr:-1} 
     227    Ifra_april_2 = {arr:-1} 
    243228    Ifra_sept_2 = {arr:-1} 
    244229  ENDELSE 
    … …  
    271256; 
    272257   Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I )  
    273    Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
    274    msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations  
    275    help, Ifra_1 
    276    help, mask 
     258;SF  ready for mask : remove 0.15% for observations 
     259;SF 
     260;SF   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
     261;SF   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations  
    277262   caldat, time, mm 
    278263   march = where(mm EQ 3, cnt) 
    279264   febr = where(mm EQ 2, cnt) 
    280    ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    281    Isal_1.arr = Isal_1.arr * msk 
    282    help, Isal_1 
    283265   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    284    help, Isal_march_1 
     266;SF 
     267;SF   Isal_1.arr = Isal_1.arr * msk 
     268;SF  
     269;SF   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    285270   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    286271   sept = where(mm EQ 9, cnt) 
    … …  
    289274 ; 
    290275   Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
    291    Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 
    292    msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations  
    293276   IF std_file2_I NE std_file1_I THEN BEGIN 
    294      Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
    295      Isal_2.arr = Isal_2.arr * msk 
    296277     caldat, time, mm 
    297278     march = where(mm EQ 3, cnt) 
    … …  
    307288     Isal_sept_2 = {arr:-1} 
    308289   ENDELSE 
     290; 
     291  IvelU_1 = read_ncdf(getenv('VAR1_IvelU'), allrecords = allrec, filename = std_file1_I )  
     292  IvelV_1 = read_ncdf(getenv('VAR1_IvelV'), allrecords = allrec, filename = std_file1_I )  
     293  Ivelo_1 = read_ncdf(getenv('VAR1_Ivelo'), allrecords = allrec, filename = std_file1_I )  
     294  caldat, time, mm 
     295  febr = where(mm EQ 2, cnt) 
     296  IvelU_febr_1 = {arr:1./float(cnt) * total(reform(IvelU_1.arr[*, *, febr],nxt,nyt,cnt), 3), unit:IvelU_1.unit, g: 'T'} 
     297  IvelV_febr_1 = {arr:1./float(cnt) * total(reform(IvelV_1.arr[*, *, febr],nxt,nyt,cnt), 3), unit:IvelV_1.unit, g: 'T'} 
     298  Ivelo_febr_1 = {arr:1./float(cnt) * total(reform(Ivelo_1.arr[*, *, febr],nxt,nyt,cnt), 3), unit:Ivelo_1.unit, g: 'T'}  
     299  march = where(mm EQ 3, cnt) 
     300  IvelU_march_1 = {arr:1./float(cnt) * total(reform(IvelU_1.arr[*, *, march],nxt,nyt,cnt), 3), unit:IvelU_1.unit, g: 'T'} 
     301  IvelV_march_1 = {arr:1./float(cnt) * total(reform(IvelV_1.arr[*, *, march],nxt,nyt,cnt), 3), unit:IvelV_1.unit, g: 'T'} 
     302  Ivelo_march_1 = {arr:1./float(cnt) * total(reform(Ivelo_1.arr[*, *, march],nxt,nyt,cnt), 3), unit:Ivelo_1.unit, g: 'T'}  
     303  sept = where(mm EQ 9, cnt) 
     304  IvelU_sept_1 = {arr:1./float(cnt) * total(reform(IvelU_1.arr[*, *, sept],nxt,nyt,cnt), 3), unit:IvelU_1.unit, g: 'T'} 
     305  IvelV_sept_1 = {arr:1./float(cnt) * total(reform(IvelV_1.arr[*, *, sept],nxt,nyt,cnt), 3), unit:IvelV_1.unit, g: 'T'} 
     306  Ivelo_sept_1 = {arr:1./float(cnt) * total(reform(Ivelo_1.arr[*, *, sept],nxt,nyt,cnt), 3), unit:Ivelo_1.unit, g: 'T'} 
     307; 
     308  IF std_file2_I NE std_file1_I THEN BEGIN 
     309    IvelU_2 = read_ncdf(getenv('VAR2_IvelU'), allrecords = allrec, filename = std_file2_I )  
     310    IvelV_2 = read_ncdf(getenv('VAR2_IvelV'), allrecords = allrec, filename = std_file2_I )  
     311    Ivelo_2 = read_ncdf(getenv('VAR2_Ivelo'), allrecords = allrec, filename = std_file2_I )  
     312    caldat, time, mm 
     313    febr = where(mm EQ 2, cnt) 
     314    IvelU_febr_2 = {arr:1./float(cnt) * total(reform(IvelU_2.arr[*, *, febr],nxt,nyt,cnt), 3), unit:IvelU_2.unit} 
     315    IvelV_febr_2 = {arr:1./float(cnt) * total(reform(IvelV_2.arr[*, *, febr],nxt,nyt,cnt), 3), unit:IvelV_2.unit} 
     316    Ivelo_febr_2 = {arr:1./float(cnt) * total(reform(Ivelo_2.arr[*, *, febr],nxt,nyt,cnt), 3), unit:Ivelo_2.unit} 
     317    march = where(mm EQ 3, cnt) 
     318    IvelU_march_2 = {arr:1./float(cnt) * total(reform(IvelU_2.arr[*, *, march],nxt,nyt,cnt), 3), unit:IvelU_2.unit} 
     319    IvelV_march_2 = {arr:1./float(cnt) * total(reform(IvelV_2.arr[*, *, march],nxt,nyt,cnt), 3), unit:IvelV_2.unit} 
     320    Ivelo_march_2 = {arr:1./float(cnt) * total(reform(Ivelo_2.arr[*, *, march],nxt,nyt,cnt), 3), unit:Ivelo_2.unit} 
     321    sept = where(mm EQ 9, cnt) 
     322    IvelU_sept_2 = {arr:1./float(cnt) * total(reform(IvelU_2.arr[*, *, sept],nxt,nyt,cnt), 3), unit:IvelU_2.unit} 
     323    IvelV_sept_2 = {arr:1./float(cnt) * total(reform(IvelV_2.arr[*, *, sept],nxt,nyt,cnt), 3), unit:IvelV_2.unit}    
     324    Ivelo_sept_2 = {arr:1./float(cnt) * total(reform(Ivelo_2.arr[*, *, sept],nxt,nyt,cnt), 3), unit:Ivelo_2.unit} 
     325    undefine, Ifra_2 
     326  ENDIF ELSE BEGIN  
     327    IvelU_febr_2 = {arr:-1} 
     328    IvelV_febr_2 = {arr:-1} 
     329    Ivelo_febr_2 = {arr:-1} 
     330    IvelU_march_2 = {arr:-1} 
     331    IvelV_march_2 = {arr:-1} 
     332    Ivelo_march_2 = {arr:-1} 
     333    IvelU_sept_2 = {arr:-1} 
     334    IvelV_sept_2 = {arr:-1} 
     335    Ivelo_sept_2 = {arr:-1} 
     336  ENDELSE 
     337; 
     338  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
     339  caldat, time, mm 
     340  febr = where(mm EQ 2, cnt) 
     341  Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     342  march = where(mm EQ 3, cnt) 
     343  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     344  sept = where(mm EQ 9, cnt) 
     345  Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit}  
     346  undefine, Iage_1 
     347; 
     348  IF std_file2_I NE std_file1_I THEN BEGIN 
     349     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
     350     caldat, time, mm 
     351     febr = where(mm EQ 2, cnt) 
     352     Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     353     march = where(mm EQ 3, cnt) 
     354     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     355     sept = where(mm EQ 9, cnt) 
     356     Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     357     undefine, Iage_2 
     358  ENDIF ELSE BEGIN  
     359     Iage_febr_2 = {arr:-1} 
     360     Iage_march_2 = {arr:-1} 
     361     Iage_sept_2 = {arr:-1} 
     362  ENDELSE 
    309363; 
    310364  jpt = 1 
    … …  
    442496  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex 
    443497; 
    444   cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40°N' 
     498  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40¡ N' 
    445499  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex 
    446500; 
    447   cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38°N' 
     501  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38¡ N' 
    448502  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 
    449503; 
    450   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: JAN' 
    451   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, /ARC, /JAN, _extra = ex 
    452 ; 
    453   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
    454   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 
    455 ;   
    456   cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
    457   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 
    458 ; 
    459   cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT' 
    460   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 
    461 ; 
    462   cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
    463   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 
    464 ;   
    465   cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT' 
    466   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 
    467 ;   
    468   cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH' 
    469   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, /ANT, /FEBR, _extra = ex 
    470 ; 
    471   cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
    472   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 
    473 ; 
     504;; 
     505; all plot are done for ice fraction > 0.15% 
     506;; 
    474507  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH' 
    475508  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
    … …  
    484517  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
    485518; 
     519  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: JAN' 
     520  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, Ifra_jan_1, Ifra_jan_2, /ARC, /JAN, _extra = ex 
     521; 
     522  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
     523  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, Ifra_april_1, Ifra_april_2, /ARC, /APRIL, _extra = ex 
     524;   
     525  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
     526  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, Ifra_april_1, Ifra_april_2, /ANT, /APRIL, _extra = ex 
     527; 
     528  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT' 
     529  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
     530; 
    486531  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: MARCH' 
    487   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 
     532  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
    488533;  
    489534  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
    490   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
     535  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 
    491536;   
    492537  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH' 
    493   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, /ANT, /FEBR, _extra = ex 
     538  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 
    494539;  
    495540  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
    496   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
     541  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
    497542;  
    498543  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: MARCH' 
    499   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 
     544  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
    500545;  
    501546  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT' 
    502   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 
     547  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 
    503548;   
    504549  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: FEBRUARY' 
    505   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 
     550  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 
    506551;  
    507552  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT' 
    508   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex 
    509 ;  
     553  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
     554;  
     555  cnt = cnt+1   &   blabla = 'Arctic Ice Velocity: FEBRUARY' 
     556  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_febr_1, IvelU_febr_2, IvelV_febr_1, IvelV_febr_2, Ivelo_febr_1, Ivelo_febr_2, /ARC, /FEBR, _extra = ex 
     557; 
     558  cnt = cnt+1   &   blabla = 'Arctic Ice Velocity: SEPT' 
     559  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_sept_1, IvelU_sept_2, IvelV_sept_1, IvelV_sept_2, Ivelo_sept_1, Ivelo_sept_2, /ARC,  /SEPT, _extra = ex 
     560; 
     561  cnt = cnt+1   &   blabla = 'Antartic Ice Velocity: MARCH' 
     562  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_march_1, IvelU_march_2, IvelV_march_1, IvelV_march_2, Ivelo_march_1, Ivelo_march_2, /ANT,  /MARCH, _extra = ex 
     563; 
     564  cnt = cnt+1   &   blabla = 'Antartic Ice Velocity: SEPT' 
     565  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_sept_1, IvelU_sept_2, IvelV_sept_1, IvelV_sept_2, Ivelo_sept_1, Ivelo_sept_2, /ANT,  /SEPT, _extra = ex 
     566; 
     567  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
     568  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
     569;   
     570  cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT' 
     571  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2,  Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 
     572;   
     573  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: FEBR' 
     574  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 
     575; 
     576  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
     577  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2,  Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
     578; 
    510579  
    511580  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*] 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_1Y.pro

    r4518 r4559  
    273273   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
    274274   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations  
    275    help, Ifra_1 
    276    help, mask 
    277275   caldat, time, mm 
    278276   march = where(mm EQ 3, cnt) 
    … …  
    280278   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    281279   Isal_1.arr = Isal_1.arr * msk 
    282    help, Isal_1 
    283280   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    284    help, Isal_march_1 
    285281   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    286282   sept = where(mm EQ 9, cnt) 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro

    r4518 r4559  
    289289   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
    290290   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations  
    291    help, Ifra_1 
    292    help, mask 
    293291   caldat, time, mm 
    294292   march = where(mm EQ 3, cnt) 
    … …  
    296294   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    297295   Isal_1.arr = Isal_1.arr * msk 
    298    help, Isal_1 
    299296   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    300    help, Isal_march_1 
    301297   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    302298   sept = where(mm EQ 9, cnt) 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh

    r4516 r4559  
    7272VAR1_Isal=sisali 
    7373VAR1_Iage=siages 
     74VAR1_IvelU=sivelu 
     75VAR1_IvelV=sivelv 
     76VAR1_Ivelo=sivelo 
    7477# 
    7578#===================== EXP2 ===================== 
    … …  
    9295VAR2_Isal=sisali 
    9396VAR2_Iage=siages 
     97VAR2_IvelU=sivelu 
     98VAR2_IvelV=sivelv 
     99VAR2_Ivelo=sivelo 
    94100# 
    95101######################### Export Variables ############################### 
    … …  
    115121export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    116122export VAR1_U      VAR1_V 
    117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
     123export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 
    118124#===================== EXP2 ===================== 
    119125export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    120126export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    121127export VAR2_U      VAR2_V 
    122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
     128export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 
    123129# 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.last_year_5D

    r4514 r4559  
    7272VAR1_Isal=iicesali 
    7373VAR1_Iage=iiceages 
     74VAR1_IvelU=iicevelu 
     75VAR1_IvelV=iicevelv 
     76VAR1_Ivelo=iicevelo 
    7477# 
    7578#===================== EXP2 ===================== 
    … …  
    9295VAR2_Isal=iicesali 
    9396VAR2_Iage=iiceages 
     97VAR2_IvelU=iicevelu 
     98VAR2_IvelV=iicevelv 
     99VAR2_Ivelo=iicevelo 
    94100# 
    95101######################### Export Variables ############################### 
    … …  
    115121export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    116122export VAR1_U      VAR1_V 
    117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
     123export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 
    118124#===================== EXP2 ===================== 
    119125export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    120126export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    121127export VAR2_U      VAR2_V 
    122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
     128export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 
    123129# 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.new_names

    r4514 r4559  
    7272VAR1_Isal=sisali 
    7373VAR1_Iage=siages 
     74VAR1_IvelU=sivelu 
     75VAR1_IvelV=sivelv 
     76VAR1_Ivelo=sivelo 
    7477# 
    7578#===================== EXP2 ===================== 
    … …  
    9295VAR2_Isal=sisali 
    9396VAR2_Iage=siages 
     97VAR2_IvelU=sivelu 
     98VAR2_IvelV=sivelv 
     99VAR2_Ivelo=sivelo 
    94100# 
    95101######################### Export Variables ############################### 
    … …  
    115121export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    116122export VAR1_U      VAR1_V 
    117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
     123export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 
    118124#===================== EXP2 ===================== 
    119125export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    120126export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    121127export VAR2_U      VAR2_V 
    122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
     128export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 
    123129# 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names

    r4514 r4559  
    7272VAR1_Isal=iicesali 
    7373VAR1_Iage=iiceages 
     74VAR1_IvelU=iicevelu 
     75VAR1_IvelV=iicevelv 
     76VAR1_Ivelo=iicevelo 
    7477# 
    7578#===================== EXP2 ===================== 
    … …  
    9497VAR2_Isal=iicesali 
    9598VAR2_Iage=iiceages 
     99VAR2_IvelU=iicevelu 
     100VAR2_IvelV=iicevelv 
     101VAR2_Ivelo=iicevelo 
    96102# 
    97103######################### Export Variables ############################### 
    … …  
    117123export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    118124export VAR1_U      VAR1_V 
    119 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
     125export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 
    120126#===================== EXP2 ===================== 
    121127export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    122128export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    123129export VAR2_U      VAR2_V 
    124 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
     130export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 
    125131# 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1

    r4514 r4559  
    7272VAR1_Isal=iicesali 
    7373VAR1_Iage=iiceages 
     74VAR1_IvelU=iicevelu 
     75VAR1_IvelV=iicevelv 
     76VAR1_Ivelo=iicevelo 
    7477# 
    7578#===================== EXP2 ===================== 
    … …  
    9295VAR2_Isal=iicesali 
    9396VAR2_Iage=iiceages 
     97VAR2_IvelU=iicevelu 
     98VAR2_IvelV=iicevelv 
     99VAR2_Ivelo=iicevelo 
    94100# 
    95101######################### Export Variables ############################### 
    … …  
    115121export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    116122export VAR1_U      VAR1_V 
    117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
     123export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 
    118124#===================== EXP2 ===================== 
    119125export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    120126export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    121127export VAR2_U      VAR2_V 
    122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
     128export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 
    123129# 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2

    r4514 r4559  
    7272VAR1_Isal=iicesali 
    7373VAR1_Iage=iiceages 
     74VAR1_IvelU=iicevelu 
     75VAR1_IvelV=iicevelv 
     76VAR1_Ivelo=iicevelo 
    7477# 
    7578#===================== EXP2 ===================== 
    … …  
    9295VAR1_Isal=iicesali 
    9396VAR1_Iage=iiceages 
     97VAR2_IvelU=iicevelu 
     98VAR2_IvelV=iicevelv 
     99VAR2_Ivelo=iicevelo 
    94100# 
    95101######################### Export Variables ############################### 
    … …  
    115121export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    116122export VAR1_U      VAR1_V 
    117 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
     123export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage VAR1_IvelU VAR1_IvelV VAR1_Ivelo 
    118124#===================== EXP2 ===================== 
    119125export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    120126export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    121127export VAR2_U      VAR2_V 
    122 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
     128export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage VAR2_IvelU VAR2_IvelV VAR2_Ivelo 
    123129# 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vel.pro

    r4514 r4559  
    1717@common 
    1818@std_common 
    19  
     19; get exp1 info 
     20  vICE1 = getenv('VAR1_ICE')   &   prefix = getenv('V1ICE_PREF')    &   suffix = getenv('V1ICE_SUFF') 
     21; get exp2 info 
     22  vICE2 = getenv('VAR2_ICE')   &   prefix2 = getenv('V2ICE_PREF')   &   suffix2 = getenv('V2ICE_SUFF') 
    2023; get exp1 info 
    2124  vICE_vel_1 = getenv('VAR1_Ivel')   &   prefix = getenv('V1Iv_PREF')  &   suffix = getenv('V1Iv_SUFF') 
    … …  
    3437  ; ICE Velocity in NORTH Hemisphere 
    3538  domdef, 0, jpi-1, 30, 90, /xindex 
    36   Velo_N = rseries_ncdf(vICE_vel_1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec,/nostruct) ;!! warning positive northward 
     39  Velo_N = rseries_ncdf(vICE_vel_1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, /nostruct) ;!! warning positive northward 
     40  ICE_N_15 = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec,  /nostruct)  
     41  print, 'N15', max(ICE_N_15) 
    3742   
    38   ;;;;;; if grossempoyenne is not present dimension of array is not ok: 54   Array[180, 54, 120]  
    39   ;;;;;; with the domain dimensions [jpi/nx, jpj/ny, jpk/nz, jpt] = [180/180, 149/50, 31/31, 12] 
    4043   
    41   ICE_velo_N = grossemoyenne(Velo_N, 'xy', /integration, mask2d = masknp) 
    42   help, ICE_velo_N 
    43   print, 'max ice velo', max(ICE_velo_N) 
     44  ICE_N_15[where(ICE_N_15 lt 0.15)] = 0. 
     45   
     46  ICE_velo_N = grossemoyenne( (Velo_N < 1.e10) * (ICE_N_15 < 1.e10), 'xy',/integration, mask2d = masknp) 
     47  ICE_N_15 = grossemoyenne(ICE_N_15, 'xy',/integration, mask2d = masknp) 
     48   
     49  ICE_velo_N = ICE_velo_N / ICE_N_15 
    4450   
    4551  if jpt mod 12 ne 0 then stop 
    … …  
    4753  ICE_velo_N = reform(ICE_velo_N, 12, nyr) 
    4854  ICE_velo_N = total(ICE_velo_N,2)/nyr 
    49   ICE_velo_N = {arr:ICE_velo_N * 1.e-19 * 86400 * 365, unit : 'm/year'}   
    50   help, ICE_velo_N 
    51   print, 'max ice velo', max(ICE_velo_N.arr) 
    52      
     55  ICE_velo_N = {arr:ICE_velo_N , unit : 'm/s'}   
     56    
     57  
    5358  ;ICE Velocity in SOUTH Hemisphere 
    5459  domdef, 0, jpi-1, -90, -30, /xindex 
    5560  Velo_S = rseries_ncdf(vICE_vel_1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec,/nostruct) ;!! warning positive northward 
     61 
    5662    
    5763     
    … …  
    5965  jpt=12 
    6066  time=julday(1,15,1900)+30*lindgen(12) 
    61   pltt, ICE_velo_N, 't', /REMPLI, /PORTRAIT, MIN = 0., MAX = 60000000. , XGRIDSTYLE = 1 $ 
     67 
     68 
     69  pltt, ICE_velo_N, 't', /REMPLI, /PORTRAIT, XGRIDSTYLE = 1 $ 
    6270        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, _extra = ex     
    6371  ; 
    6472  title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Ice Velocity (Black SOLID simulation)' 
    65   pltt, Velo_S, 't',  /REMPLI, MIN = 0., MAX =  60000000., /NOERASE, XGRIDSTYLE = 1 $ 
     73  pltt, Velo_S, 't',  /REMPLI, /NOERASE, XGRIDSTYLE = 1 $ 
    6674       , small = [1, 2, 2], YTITLE = varunit, TITLE = title, _extra = ex 
    6775  ; 
    … …  
    7785    domdef, 0, jpi-1, 30, 90, /xindex 
    7886    Velo_N_2 = rseries_ncdf(vICE_vel_2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec,/nostruct) 
    79     help, Velo_N_2  
    8087     
    8188    ;;;;;; if grossempoyenne is not present dimension of array is not ok: 54   Array[180, 54, 120]  
    … …  
    8390   
    8491    ICE_velo_N_2 = grossemoyenne(Velo_N_2, 'xy', /integration, mask2d = masknp) 
    85     help, ICE_velo_N_2 
    8692    print, 'max ice velo', max(ICE_velo_N_2) 
    8793   
    … …  
    9197    ICE_velo_N_2 = total(ICE_velo_N_2,2)/nyr 
    9298    ICE_velo_N_2 = {arr:ICE_velo_N_2 * 86400 * 365, unit : 'm/year'}   
    93     help, ICE_velo_N_2 
    9499    print, 'max ice velo', max(ICE_velo_N_2.arr) 
    95100     
    … …  
    98103    domdef, 0, jpi-1, -90, -30, /xindex 
    99104    Velo_S_2 = rseries_ncdf(vICE_vel_2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec,/nostruct) 
    100 ;SF    help, Velo_S_2 
    101105       
    102106    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
    103107 
     108    jpt=12 
     109    time=julday(1,15,1900)+30*lindgen(12) 
     110 
    104111    title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 +'!C'+'Ice Velocity (BLACK) ' 
    105     pltt, ICE_velo_N , 't', MIN = 0., MAX =  60000000. , date1, date2, /REMPLI, XGRIDSTYLE = 1, window = 2  $ 
     112    pltt, ICE_velo_N , 't', /REMPLI, XGRIDSTYLE = 1, window = 2  $ 
    106113          , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 
    107     pltt, ICE_velo_N_2, 't', date1, date2, /REMPLI $ 
     114    pltt, ICE_velo_N_2, 't', /REMPLI $ 
    108115          , /ov1d, color = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 
    109116  
    110117    title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 +'!C'+'Ice Velocity (BLACK) ' 
    111     pltt, Velo_S , 't', MIN = -.04, MAX = .04, date1, date2, /REMPLI, XGRIDSTYLE = 1 $ 
     118    pltt, Velo_S , 't', /REMPLI, XGRIDSTYLE = 1 $ 
    112119          , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 
    113     pltt, Velo_S_2, 't', date1, date2, /REMPLI $ 
     120    pltt, Velo_S_2, 't', /REMPLI $ 
    114121         , /ov1d, color = 250, small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex 
    115122 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro

    r4514 r4559  
    7878  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 
    7979; 
    80 ;SF    cnt = cnt+1   &   blabla = 'ICE Velocity ' 
    81 ;SF    IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 
     80  cnt = cnt+1   &   blabla = 'ICE Velocity ' 
     81  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 
    8282; 
    8383  cnt = cnt+1   &   blabla = 'Sea-Ice Volume Export at Fram Strait' 

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh

    r4516 r4559  
    1818idl_command=/Applications/itt/idl64/bin/idl 
    1919# 
    20 PS_DIR=$( pwd )/ts_heatbug-sbcmodMV_ps 
    21 PDF_DIR=$( pwd )/ts_heatbug-sbcmodMV_pdf 
     20PS_DIR=$( pwd )/ts_heatbug_ps 
     21PDF_DIR=$( pwd )/ts_heatbug_pdf 
    2222HTML_DIR=$( pwd )/html 
    2323SAXO_DIR=/Users/sflod/SAXO_DIR 
    … …  
    5353# 
    5454#DATE1=20010101   ;   DATE2=25001231 
    55 DATE1=20010101   ;   DATE2=20501231 
     55DATE1=20310101   ;   DATE2=20401231 
    5656#DATE1=00010101      ;   DATE2=00101231 
    5757# 
    … …  
    7272# 
    7373#DATE1_2=20010101   ;   DATE2_2=20101231 
    74 DATE1_2=20010101   ;   DATE2_2=20501231 
     74DATE1_2=20310101   ;   DATE2_2=20401231 
    7575#DATE1_2=20010101   ;   DATE2_2=25001231 
    7676# 
    77 VAR2_T=thetao      ;   V2T_PREF=sbcmodMV      ;   V2T_SUFF=_1Y_grid_T.nc 
    78 VAR2_S=so          ;   V2S_PREF=sbcmodMV      ;   V2S_SUFF=_1Y_grid_T.nc 
    79 VAR2_SSH=zos       ;   V2SSH_PREF=sbcmodMV    ;   V2SSH_SUFF=_1Y_grid_T.nc 
    80 VAR2_Q=tohfls      ;   V2Q_PREF=sbcmodMV      ;   V2Q_SUFF=_1Y_grid_T.nc 
    81 VAR2_EMP=wfo       ;   V2EMP_PREF=sbcmodMV    ;   V2EMP_SUFF=_1Y_grid_T.nc 
    82 VAR2_U=uocetr_eff  ;   V2U_PREF=sbcmodMV      ;   V2U_SUFF=_1Y_grid_U.nc 
    83 VAR2_V=vocetr_eff  ;   V2V_PREF=sbcmodMV      ;   V2V_SUFF=_1Y_grid_V.nc 
    84 VAR2_ICE=siconc    ;   V2ICE_PREF=sbcmodMV    ;   V2ICE_SUFF=_1M_icemod.nc 
    85 VAR2_Ithick=sithic ;   V2It_PREF=sbcmodMV     ;   V2It_SUFF=_1M_icemod.nc 
    86 VAR2_SNOW=sndept   ;   V2SNOW_PREF=sbcmodMV   ;   V2SNOW_SUFF=_1M_icemod.nc 
    87 VAR2_IvelV=sivelv  ;   V2IvV_PREF=sbcmodMV    ;   V2IvV_SUFF=_1M_icemod.nc 
    88 VAR2_Ivel=sivelo   ;   V2Iv_PREF=sbcmodMV     ;   V2Iv_SUFF=_1M_icemod.nc 
     77VAR2_T=thetao      ;   V2T_PREF=heatbug      ;   V2T_SUFF=_1Y_grid_T.nc 
     78VAR2_S=so          ;   V2S_PREF=heatbug      ;   V2S_SUFF=_1Y_grid_T.nc 
     79VAR2_SSH=zos       ;   V2SSH_PREF=heatbug    ;   V2SSH_SUFF=_1Y_grid_T.nc 
     80VAR2_Q=tohfls      ;   V2Q_PREF=heatbug      ;   V2Q_SUFF=_1Y_grid_T.nc 
     81VAR2_EMP=wfo       ;   V2EMP_PREF=heatbug    ;   V2EMP_SUFF=_1Y_grid_T.nc 
     82VAR2_U=uocetr_eff  ;   V2U_PREF=heatbug      ;   V2U_SUFF=_1Y_grid_U.nc 
     83VAR2_V=vocetr_eff  ;   V2V_PREF=heatbug      ;   V2V_SUFF=_1Y_grid_V.nc 
     84VAR2_ICE=siconc    ;   V2ICE_PREF=heatbug    ;   V2ICE_SUFF=_1M_icemod.nc 
     85VAR2_Ithick=sithic ;   V2It_PREF=heatbug     ;   V2It_SUFF=_1M_icemod.nc 
     86VAR2_SNOW=sndept   ;   V2SNOW_PREF=heatbug   ;   V2SNOW_SUFF=_1M_icemod.nc 
     87VAR2_IvelV=sivelv  ;   V2IvV_PREF=heatbug    ;   V2IvV_SUFF=_1M_icemod.nc 
     88VAR2_Ivel=sivelo   ;   V2Iv_PREF=heatbug     ;   V2Iv_SUFF=_1M_icemod.nc 
    8989# 
    9090######################### Export Variables ############################### 
Note: See TracChangeset for help on using the changeset viewer.