Changeset 59 for branches/2016/dev_v3.20_2016_platelet/GRAPHICS/ALL_stocks.pro

Timestamp:

12/12/17 14:33:19 (7 years ago)

Author:

vancop

Message:

graphics for Nice

File:

• branches/2016/dev_v3.20_2016_platelet/GRAPHICS/ALL_stocks.pro (modified) (2 diffs)

branches/2016/dev_v3.20_2016_platelet/GRAPHICS/ALL_stocks.pro

@@ -213 +213 @@
 ; YROSIAE
 ; taken from yrosiae ASPECT log sheets (Fripiat et al, oct 2015)
+   doy_obs    = [ 263 , 299,  318 , 279 , 306 , 327,  292 , 312  ,335 ]
+   Ichla_obs  = [ 1.4683  , 4.6038  ,45.4413, 128.0884,  44.5830 ,219.1416 ,208.6254 ,108.0657 , 87.9038 ]
+   ISi_obs    = [ 22.6631 , 21.1690 , 19.2174,  18.8949 , 18.9943,  18.9517,  20.5963,  16.6855,  13.8374 ]
+   IN_obs     = [ 9.0613  , 5.6130  , 9.4723,  10.0763,  11.4708 , 14.6193 , 14.4107 ,  4.7585 ,  2.1876 ]
+   IP_obs     = [ 0.4976  , 0.4067  , 1.1311,   1.8459,   2.3452 ,  3.7857 ,  3.6948 ,  1.3270 ,  1.2602 ]
+ENDIF
+
+IF ( site EQ 'N_ICE_RL' ) THEN BEGIN
+; N_ICE_RL
+; ALLLLLL WRONG !!!
    doy_obs    = [ 263 , 299,  318 , 279 , 306 , 327,  292 , 312  ,335 ]
    Ichla_obs  = [ 1.4683  , 4.6038  ,45.4413, 128.0884,  44.5830 ,219.1416 ,208.6254 ,108.0657 , 87.9038 ]
@@ -492 +502 @@
    XYOUTS, xmin+ (xmax-xmin)*0.10, ymax*0.80, ztitle
 
-   ; Iron stocks
-   ymin = 0.0000   & ymax = MAX(dFet)*1.2
-
-   PLOT, [xmin, xmax], [ ymin, ymax ], /NODATA, charsize = cs, $
-         XTITLE = 'mon', YTITLE = 'Iron (mumol/m2)', XSTYLE = 1, YSTYLE = 1
-
-   FOR i_run = nruns - 1, 0, -1 DO BEGIN
-      LOADCT, ct
-      OPLOT, mon(i_run,0:ntss(i_run)-1), dFet(0:ntss(i_run)-1), color = 150, thick = 2
-      OPLOT, mon(i_run,0:ntss(i_run)-1), aFet(0:ntss(i_run)-1), color = 0, thick = 2
-      OPLOT, mon(i_run,0:ntss(i_run)-1), eFet(0:ntss(i_run)-1), color = 0, thick = 1
-   ENDFOR
-   LOADCT, 0
-
-   ztitle = 'Max dFe stock: '+STRCOMPRESS(STRING(MAX(dFet), FORMAT='(F6.2)'),/REMOVE_ALL)+' mumol/m2'
-   XYOUTS, xmin+ (xmax-xmin)*0.10, ymax*0.90, ztitle
+;  ; Iron stocks
+;  ymin = 0.0000   & ymax = 500. ;ymax = MAX(dFet)*1.2
+
+;  PLOT, [xmin, xmax], [ ymin, ymax ], /NODATA, charsize = cs, $
+;        XTITLE = 'mon', YTITLE = 'Iron (mumol/m2)', XSTYLE = 1, YSTYLE = 1
+
+;  FOR i_run = nruns - 1, 0, -1 DO BEGIN
+;     LOADCT, ct
+;     OPLOT, mon(i_run,0:ntss(i_run)-1), dFet(0:ntss(i_run)-1), color = 150, thick = 2
+;     OPLOT, mon(i_run,0:ntss(i_run)-1), aFet(0:ntss(i_run)-1), color = 0, thick = 2
+;     OPLOT, mon(i_run,0:ntss(i_run)-1), eFet(0:ntss(i_run)-1), color = 0, thick = 1
+;  ENDFOR
+;  LOADCT, 0
+
+;  ztitle = 'Max dFe stock: '+STRCOMPRESS(STRING(MAX(dFet), FORMAT='(F6.2)'),/REMOVE_ALL)+' mumol/m2'
+;  XYOUTS, xmin+ (xmax-xmin)*0.10, ymax*0.90, ztitle