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_ts_AMOC.pro

    r2751 r4320  
    1414  ts_Time = 0. 
    1515 
    16   trans = strlowcase(var_name) EQ 'vocetr_eff' 
     16  ;SF trans = strlowcase(var_name) EQ 'vocetr_eff' 
     17  trans = strlowcase(var_name) EQ 'vomeeftr' 
    1718 
     19  ; ADD definition of mask_atl_nomed to compute AMOC (with no Mediterranean Sea) 
     20  msk_atl_nomed = read_ncdf( 'atlmsk_nomed', filename = std_file_msksub, /nostruct, _extra = ex ) 
    1821  FOR i = 0, nfiles-1 DO BEGIN 
    1922    var = read_ncdf(var_name, allrecords = allrec, filename = list[i], /nostruct) 
     
    2124; 
    2225    FOR t = 0, jpt-1 DO BEGIN 
     26      msfatl = msf(var[*, *, *, t], msk_atl_nomed, TRANSPORT = trans, /nostruct, indexboxzoom = ind) 
    2327      msftot = msf(var[*, *, *, t], TRANSPORT = trans, /nostruct, indexboxzoom = ind) 
    2428      yaxis = gphit[ind[0], ind[2]:ind[3]] 
     
    2731      indy = where(yaxis gt 40 and yaxis le 50) 
    2832      domdef, 0, 3500 
    29       t45 = [t45, max(msftot[indy, firstzw:lastzw], /NaN)] 
     33      ;SF commented because we've to compute msftot starting from msf atlantic; i.e. without mediterranean sea. 
     34      ;SF t45 = [t45, max(msftot[indy, firstzw:lastzw], /NaN)] 
     35      t45 = [t45, max(msfatl[indy, firstzw:lastzw], /NaN)] 
    3036       
    3137; computation of max atlantic Antarctic Bottom Water between 80°S and 65°S 
     
    7884 
    7985; compute the MSF 
    80  
     86 
    8187  std_ts_read_AMOC, vV1, date1, date2, prefix, suffix, t45, t70, t50 
    8288  if prefix NE prefix2 then begin 
     
    114120  if KEYWORD_SET(postscript) then openps, filename+'_3.ps', portrait = 1 
    115121 
    116   title = prefix+' '+d1_d2+'!C'+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2 
     122  title = prefix+' '+d1_d2+'!C '+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2 
    117123  pltt, -t50, 't', 5., 30., date1, date2, /REMPLI, /PORTRAIT, window = 3 $ 
    118124        ,  small = [1, 2, 1],  TITLE = title, YTITLE = varunit, _extra = ex 
Note: See TracChangeset for help on using the changeset viewer.