Changeset 13843

Timestamp:

2020-11-20T16:00:36+01:00 (3 years ago)

Author:

laurent

Message:

Improving plots over ice in STATION_ASF...

File:

• NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_ICE.py (modified) (6 diffs)

NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/tests/STATION_ASF/EXPREF/plot_station_asf_ICE.py

@@ -28 +28 @@
 rDPI=100.
 
-l_color = [ '0.85'    ,  '#ffed00' , '#008ab8' , '0.4'  ] ; # colors to differentiate algos on the plot
-l_width = [  4        ,     3      ,    2      ,  1     ] ; # line-width to differentiate algos on the plot
-l_style = [ '-'       ,    '-'     ,   '-'     , '--'   ] ; # line-style
+#ffed00: yellow ON
+#E8A727: ornage
+
+l_color = [ '0.3' , '#E8A727', '0.1'  , '#008ab8'  ] ; # colors to differentiate algos on the plot
+l_width = [   2   ,    2     ,  1.5   ,    2       ] ; # line-width to differentiate algos on the plot
+l_style = [  '-'  ,    '-'   , '--'   ,   '-'      ] ; # line-style
 
 
@@ -43 +46 @@
 L_PREC = [     3     ,     3     ,      0     ,      0        ,      0      ,      0        ,     0       ,     2       ]
 L_ANOM = [   False   ,   False   ,   True      ,    True       ,    True      ,    True     ,    True     ,   True      ]
+L_MAXT = [  10000.   ,     10000.,  10000.   ,  10000.      ,  10000.    ,  10000.      , 10000.     ,      1.5    ]
+L_MINT = [    0.001  ,     0.001 , -10000.   , -10000.      , -10000.    , -10000.      ,-10000.     ,   -10000.   ]
+
 
 # About STATION_ASF output files to read:
@@ -92 +98 @@
 cf_in    = []
 list_exp = []
+list_ext = [] ; # for title
 list_frc = []
 for fn in ldir:
@@ -100 +107 @@
             clab = fn[np:-nc-ns]
             [ cexp, cfrc ] = str.split(clab, '_', 1)
-            print('  ===> Experiment = '+cexp+', Forcing = '+cfrc)
+            cexp = cexp.split('-')[1] ; # removing air-sea algo name...
+            cext = cexp
+            if   cexp == "LG15":
+                cext ="Lüpkes & Gryanik, 2015"
+            elif cexp == "LU12":
+                cext ="Lüpkes et al., 2012"
+            elif cexp == "AN05":
+                cext ="Andreas et al., 2005"
+            elif cexp == "CSTC":
+                cext ="Constant coefficients"
+            print('  ===> Experiment = '+cexp+' ('+cext+'), Forcing = '+cfrc)
             list_exp.append(cexp)
+            list_ext.append(cext)
             list_frc.append(cfrc)
             cf_in.append(fpn)
 nbf = len( set(list_frc) )
+
+
 
 if nbf == 0:
@@ -175 +195 @@
         #
         id_in = Dataset(cf_in[ja])
-        xF[:,ja] = id_in.variables[L_VNEM[jv]][:,1,1] # only the center point of the 3x3 spatial domain!
+        xF[:,ja] = id_in.variables[L_VNEM[jv]][:,1,1]  ; # only the center point of the 3x3 spatial domain!        
         if ja == 0: cvar_lnm = id_in.variables[L_VNEM[jv]].long_name
         id_in.close()
+        #
+        id_toolarge, = nmp.where( xF[:,ja] > L_MAXT[jv] ) # 
+        xF[id_toolarge,ja] = L_MAXT[jv]
+        id_toosmall, = nmp.where( xF[:,ja] < L_MINT[jv] ) ; print("id_toosmall =", id_toosmall)
+        xF[id_toosmall,ja] = L_MINT[jv]
 
     idx_okay = nmp.where( nmp.abs(xF) < 1.e+10 )
@@ -191 +216 @@
         fplot = nmp.ma.masked_where( xF[:,ja]==0., xF[:,ja] )
         plt.plot(vtime, fplot, '-', color=l_color[ja], \
-                 linestyle=l_style[ja], linewidth=l_width[ja], label=list_exp[ja], zorder=10+ja)
+                 linestyle=l_style[ja], linewidth=l_width[ja], label=list_ext[ja], alpha=0.6 ) #zorder=10+ja)
 
     fmin, fmax = round_bounds( nmp.min(xF[idx_okay]) , nmp.max(xF[idx_okay]), base=L_BASE[jv], prec=L_PREC[jv])