Changeset 4004

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

Location:

branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts

Files:

• std_ts_AMOC.pro (modified) (5 diffs)
• std_ts_vardef.sh (modified) (3 diffs)
• std_ts_vardef.sh.l3fwb78 (deleted)
• std_ts_vardef.sh.l3fwb78-500yfwb0 (deleted)

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

branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh

@@ -18 +18 @@
 idl_command=/Applications/itt/idl64/bin/idl
 #
-PS_DIR=$( pwd )/ts_500yfwb0_ps_500y
-PDF_DIR=$( pwd )/ts_500yfwb0_pdf_500y
+PS_DIR=$( pwd )/ts_qsriceok_ps
+PDF_DIR=$( pwd )/ts_qsriceok_pdf
 HTML_DIR=$( pwd )/html
 SAXO_DIR=/Users/sflod/SAXO_DIR
@@ -55 +55 @@
 #DATE1=00010101      ;   DATE2=00101231
 #
-VAR1_T=votemper         ;   V1T_PREF=500yfwb0      ;   V1T_SUFF=_1Y_grid_T.nc
-VAR1_S=vosaline         ;   V1S_PREF=500yfwb0      ;   V1S_SUFF=_1Y_grid_T.nc
-VAR1_SSH=sossheig       ;   V1SSH_PREF=500yfwb0    ;   V1SSH_SUFF=_1Y_grid_T.nc
-VAR1_Q=sohflxdo         ;   V1Q_PREF=500yfwb0      ;   V1Q_SUFF=_1Y_grid_T.nc
-VAR1_EMP=sowaflup       ;   V1EMP_PREF=500yfwb0    ;   V1EMP_SUFF=_1Y_grid_T.nc
-VAR1_U=vozoeftr         ;   V1U_PREF=500yfwb0      ;   V1U_SUFF=_1Y_grid_U.nc
-VAR1_V=vomeeftr         ;   V1V_PREF=500yfwb0      ;   V1V_SUFF=_1Y_grid_V.nc
-VAR1_ICE=iiceconc       ;   V1ICE_PREF=500yfwb0    ;   V1ICE_SUFF=_1M_icemod.nc
-VAR1_Ithick=iicethic    ;   V1It_PREF=500yfwb0     ;   V1It_SUFF=_1M_icemod.nc
-VAR1_SNOW=isnowthi      ;   V1SNOW_PREF=500yfwb0   ;   V1SNOW_SUFF=_1M_icemod.nc
+VAR1_T=votemper         ;   V1T_PREF=qsriceok      ;   V1T_SUFF=_1Y_grid_T.nc
+VAR1_S=vosaline         ;   V1S_PREF=qsriceok      ;   V1S_SUFF=_1Y_grid_T.nc
+VAR1_SSH=sossheig       ;   V1SSH_PREF=qsriceok    ;   V1SSH_SUFF=_1Y_grid_T.nc
+VAR1_Q=sohflxdo         ;   V1Q_PREF=qsriceok      ;   V1Q_SUFF=_1Y_grid_T.nc
+VAR1_EMP=sowaflup       ;   V1EMP_PREF=qsriceok    ;   V1EMP_SUFF=_1Y_grid_T.nc
+VAR1_U=vozoeftr         ;   V1U_PREF=qsriceok      ;   V1U_SUFF=_1Y_grid_U.nc
+VAR1_V=vomeeftr         ;   V1V_PREF=qsriceok      ;   V1V_SUFF=_1Y_grid_V.nc
+VAR1_ICE=iiceconc       ;   V1ICE_PREF=qsriceok    ;   V1ICE_SUFF=_1M_icemod.nc
+VAR1_Ithick=iicethic    ;   V1It_PREF=qsriceok     ;   V1It_SUFF=_1M_icemod.nc
+VAR1_SNOW=isnowthi      ;   V1SNOW_PREF=qsriceok   ;   V1SNOW_SUFF=_1M_icemod.nc
 #
 #===================== EXP2 =====================
@@ -71 +71 @@
 DATE1_2=20010101   ;   DATE2_2=25001231
 #
-VAR2_T=votemper         ;   V2T_PREF=500yfwb0      ;   V2T_SUFF=_1Y_grid_T.nc
-VAR2_S=vosaline         ;   V2S_PREF=500yfwb0      ;   V2S_SUFF=_1Y_grid_T.nc
-VAR2_SSH=sossheig       ;   V2SSH_PREF=500yfwb0    ;   V2SSH_SUFF=_1Y_grid_T.nc
-VAR2_Q=sohflxdo         ;   V2Q_PREF=500yfwb0      ;   V2Q_SUFF=_1Y_grid_T.nc
-VAR2_EMP=sowaflup       ;   V2EMP_PREF=500yfwb0    ;   V2EMP_SUFF=_1Y_grid_T.nc
-VAR2_U=vozoeftr         ;   V2U_PREF=500yfwb0      ;   V2U_SUFF=_1Y_grid_U.nc
-VAR2_V=vomeeftr         ;   V2V_PREF=500yfwb0      ;   V2V_SUFF=_1Y_grid_V.nc
-VAR2_ICE=iiceconc       ;   V2ICE_PREF=500yfwb0    ;   V2ICE_SUFF=_1M_icemod.nc
-VAR2_Ithick=iicethic    ;   V2It_PREF=500yfwb0     ;   V2It_SUFF=_1M_icemod.nc
-VAR2_SNOW=isnowthi      ;   V2SNOW_PREF=500yfwb0   ;   V2SNOW_SUFF=_1M_icemod.nc
+VAR2_T=votemper         ;   V2T_PREF=qsriceok      ;   V2T_SUFF=_1Y_grid_T.nc
+VAR2_S=vosaline         ;   V2S_PREF=qsriceok      ;   V2S_SUFF=_1Y_grid_T.nc
+VAR2_SSH=sossheig       ;   V2SSH_PREF=qsriceok    ;   V2SSH_SUFF=_1Y_grid_T.nc
+VAR2_Q=sohflxdo         ;   V2Q_PREF=qsriceok      ;   V2Q_SUFF=_1Y_grid_T.nc
+VAR2_EMP=sowaflup       ;   V2EMP_PREF=qsriceok    ;   V2EMP_SUFF=_1Y_grid_T.nc
+VAR2_U=vozoeftr         ;   V2U_PREF=qsriceok      ;   V2U_SUFF=_1Y_grid_U.nc
+VAR2_V=vomeeftr         ;   V2V_PREF=qsriceok      ;   V2V_SUFF=_1Y_grid_V.nc
+VAR2_ICE=iiceconc       ;   V2ICE_PREF=qsriceok    ;   V2ICE_SUFF=_1M_icemod.nc
+VAR2_Ithick=iicethic    ;   V2It_PREF=qsriceok     ;   V2It_SUFF=_1M_icemod.nc
+VAR2_SNOW=isnowthi      ;   V2SNOW_PREF=qsriceok   ;   V2SNOW_SUFF=_1M_icemod.nc
 #
 ######################### Export Variables ###############################