Changeset 4320 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_AMOC.pro
- Timestamp:
- 2013-11-28T12:41:27+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_AMOC.pro
r2751 r4320 14 14 ts_Time = 0. 15 15 16 trans = strlowcase(var_name) EQ 'vocetr_eff' 16 ;SF trans = strlowcase(var_name) EQ 'vocetr_eff' 17 trans = strlowcase(var_name) EQ 'vomeeftr' 17 18 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 ) 18 21 FOR i = 0, nfiles-1 DO BEGIN 19 22 var = read_ncdf(var_name, allrecords = allrec, filename = list[i], /nostruct) … … 21 24 ; 22 25 FOR t = 0, jpt-1 DO BEGIN 26 msfatl = msf(var[*, *, *, t], msk_atl_nomed, TRANSPORT = trans, /nostruct, indexboxzoom = ind) 23 27 msftot = msf(var[*, *, *, t], TRANSPORT = trans, /nostruct, indexboxzoom = ind) 24 28 yaxis = gphit[ind[0], ind[2]:ind[3]] … … 27 31 indy = where(yaxis gt 40 and yaxis le 50) 28 32 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)] 30 36 31 37 ; computation of max atlantic Antarctic Bottom Water between 80°S and 65°S … … 78 84 79 85 ; compute the MSF 80 86 ; 81 87 std_ts_read_AMOC, vV1, date1, date2, prefix, suffix, t45, t70, t50 82 88 if prefix NE prefix2 then begin … … 114 120 if KEYWORD_SET(postscript) then openps, filename+'_3.ps', portrait = 1 115 121 116 title = prefix+' '+d1_d2+'!C '+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2122 title = prefix+' '+d1_d2+'!C '+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2 117 123 pltt, -t50, 't', 5., 30., date1, date2, /REMPLI, /PORTRAIT, window = 3 $ 118 124 , small = [1, 2, 1], TITLE = title, YTITLE = varunit, _extra = ex
Note: See TracChangeset
for help on using the changeset viewer.