Changeset 4004 for branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_AMOC.pro

Timestamp:

2013-08-07T16:32:18+02:00 (11 years ago)

Author:

flavoni

Message:

correct script for ts AMOC, in IDL scripts, see ticket #1115

File:

• branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_AMOC.pro (modified) (5 diffs)

branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_AMOC.pro

@@ -17 +17 @@
   trans = strlowcase(var_name) EQ 'vomeeftr'
 
+  ; ADD definition of mask_atl_nomed to compute AMOC (with no Mediterranean Sea)
+  msk_atl_nomed = read_ncdf( 'atlmsk_nomed', filename = std_file_msksub, /nostruct, _extra = ex )
   FOR i = 0, nfiles-1 DO BEGIN
     var = read_ncdf(var_name, allrecords = allrec, filename = list[i], /nostruct)
@@ -22 +24 @@
 ;
     FOR t = 0, jpt-1 DO BEGIN
+      msfatl = msf(var[*, *, *, t], msk_atl_nomed, TRANSPORT = trans, /nostruct, indexboxzoom = ind)
       msftot = msf(var[*, *, *, t], TRANSPORT = trans, /nostruct, indexboxzoom = ind)
       yaxis = gphit[ind[0], ind[2]:ind[3]]
@@ -28 +31 @@
       indy = where(yaxis gt 40 and yaxis le 50)
       domdef, 0, 3500
-      t45 = [t45, max(msftot[indy, firstzw:lastzw], /NaN)]
+      ;SF commented because we've to compute msftot starting from msf atlantic; i.e. without mediterranean sea.
+      ;SF t45 = [t45, max(msftot[indy, firstzw:lastzw], /NaN)]
+      t45 = [t45, max(msfatl[indy, firstzw:lastzw], /NaN)]
       
 ; computation of max atlantic Antarctic Bottom Water between 80°S and 65°S
@@ -79 +84 @@
 
 ; compute the MSF
-
+;  
   std_ts_read_AMOC, vV1, date1, date2, prefix, suffix, t45, t70, t50
   if prefix NE prefix2 then begin
@@ -115 +120 @@
   if KEYWORD_SET(postscript) then openps, filename+'_3.ps', portrait = 1
 
-  title = prefix+' '+d1_d2+'!C'+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2
+  title = prefix+' '+d1_d2+'!C '+'Max AntArctic Bottom Cell between 65S and 30N '+d1_d2
   pltt, -t50, 't', 5., 30., date1, date2, /REMPLI, /PORTRAIT, window = 3 $
         ,  small = [1, 2, 1],  TITLE = title, YTITLE = varunit, _extra = ex