Changeset 59
- Timestamp:
- 12/12/17 14:33:19 (7 years ago)
- Location:
- branches/2016/dev_v3.20_2016_platelet/GRAPHICS
- Files:
-
- 1 added
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_v3.20_2016_platelet/GRAPHICS/ALL_contours.pro
r6 r59 129 129 ENDIF 130 130 131 IF ( site EQ 'N_ICE_RL' ) THEN BEGIN 132 mon = numt * dt / 86400. / 30.41666666 + 1 133 ENDIF 134 131 135 ; 132 136 ;------------------------------------------------------------------------------ … … 199 203 ENDIF 200 204 205 IF ( site EQ 'N_ICE_RL' ) THEN BEGIN ; all wrong here, just a silly copy paste 206 day_obs = [ 120, 120, 120, 121, 121, 121, 124, 124, 124, 124, 125, 125, 125, 125, 125, 125, 125, 125, 125, 126, 126, 126, 127, 127, 127, 127, 130, 132, 134, 136, 138, 140, 142, 144, 146, 149, 152, 154, 127, 130, 132, 134, 136, 138, 140, 142, 144, 146, 149, 152, 154, 127, 130, 132, 134, 136, 138, 140, 142, 144, 146, 149, 152, 154, 127, 130, 132, 134, 136, 138, 140, 142, 144, 146, 149, 152, 154, 127, 130, 132, 134, 136, 138, 140, 142, 144, 146, 149, 152, 154 ] 207 208 hi_obs = [ 16, 17, 13, 17, 17, 17, 16, 16, 16, 17, 16, 17.5, 20, 14, 17, 14, 18, 30, 19, 15, 15, 15, 19, 18, 19, 17, 19.35, 21.25, 22, 22.1, 22.2, 22, 21.5, 22, 22, 23, 21, 19, 24, 22.25, 23.25, 25.25, 26, 26, 26.5, 24.5, 24, 24, 23, 23, 20.5, 21.75, 23.25, 24.75, 25, 26.5, 27, 26, 25, 25, 25, 24, 23, 21, 21.5, 22.75, 24.5, 25.4, 25.75, 25.7, 25.5, 24, 24.5, 25, 24.5, 23, 21.5, 20.25, 22.25, 23.88, 24.75, 25, 24.7, 25, 24, 24.5, 24.5, 24, 23, 20 ] / 100. 209 hi_obs_std = REPLICATE( 0., N_ELEMENTS(day_obs) ) 210 hs_obs = REPLICATE( 0., N_ELEMENTS(day_obs) ) 211 hs_obs_std = REPLICATE( 0., N_ELEMENTS(day_obs) ) 212 mon_obs = day_obs / 30.467 + 1 213 ENDIF 214 201 215 ; 202 216 ;------------------------------------------------------------------------------ … … 231 245 ENDIF 232 246 247 IF ( site EQ 'N_ICE_RL' ) THEN BEGIN 248 zmin = -0.5 249 zmax = 0.0 250 ENDIF 251 233 252 ;------------------- 234 253 ; Salinity contours … … 416 435 ; contour parameters 417 436 nlevels = 10 437 418 438 colors = [ 20, 40, 60, 80, 100, 120, 140, 160, 180, 200, 220 ] 419 439 IF ( site EQ 'BARROW' ) THEN BEGIN 420 440 levels = [ 0., 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0] 421 441 ENDIF 442 422 443 IF ( site EQ 'ISPOL' ) THEN BEGIN 423 444 levels = [ 0., 2., 4., 6., 8., 10., 12.,14., 16., 18., 20.] … … 426 447 colors = [ 20, 60, 100, 140, 180, 220, 240 ] 427 448 ENDIF 428 IF ( site EQ 'YROSIAE' ) THEN BEGIN 449 450 IF ( site EQ 'YROSIAE' OR site EQ 'N_ICE_RL' ) THEN BEGIN 429 451 nlevels = 7. 430 452 levels = [ 0., 1.6, 4., 6., 8., 10., 12., 20.] … … 463 485 colors = [ 20, 60, 100, 140, 180, 220, 240 ] 464 486 ENDIF 465 IF ( site EQ 'YROSIAE' ) THEN BEGIN487 IF ( site EQ 'YROSIAE' OR site EQ 'N_ICE_RL' ) THEN BEGIN 466 488 nlevels = 7 467 489 levels = [ 0., 0.24, 0.4, 0.6, 0.8, 1., 2. ] … … 503 525 colors = [ 20, 60, 100, 140, 180, 220, 240 ] 504 526 ENDIF 505 IF ( site EQ 'YROSIAE' ) THEN BEGIN527 IF ( site EQ 'YROSIAE' OR site EQ 'N_ICE_RL' ) THEN BEGIN 506 528 nlevels = 7. 507 529 levels = [ 0., 3.9, 8., 12., 16., 20., 24., 40.] … … 1225 1247 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 1226 1248 1227 ;1228 ;--------------1229 ; dFe1230 ;--------------1231 ;1232 !P.MULTI=[0,numplot_x, numplot_y]1233 i_plot = 01234 i_plot_h = i_plot MOD numplot_x1235 i_plot_v = i_plot / numplot_x + 11236 1237 ; contour parameters1238 nlevels = 101239 nlevels = 101240 fmin = 0.0 & fmax = MAX(dFeb) ;1241 cmin = 000 & cmax = 2551242 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors1243 1244 ct = 3 ; colortable1245 title = ' dFe (mumol/m3)'1246 1247 icecontour, dFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $1248 tmin, tmax, zmin, zmax, title1249 1250 ; colorbar1251 format='(f5.1)'1252 horpos = 0.02 + i_plot_h * dh1253 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar.1254 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format1255 1256 ;1257 ;--------------1258 ; aFe1259 ;--------------1260 ;1261 i_plot = i_plot + 11262 i_plot_h = i_plot MOD numplot_x1263 i_plot_v = i_plot / numplot_x + 11264 1265 ; contour parameters1266 nlevels = 101267 nlevels = 101268 fmin = 0.0 & fmax = MAX(aFeb) ;1269 cmin = 000 & cmax = 2551270 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors1271 1272 ct = 3 ; colortable1273 title = ' aFe (mumol/m3)'1274 1275 icecontour, aFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $1276 tmin, tmax, zmin, zmax, title1277 1278 ; colorbar1279 format='(f5.1)'1280 horpos = 0.02 + i_plot_h * dh1281 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar.1282 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format1283 1284 ;1285 ;--------------1286 ; eFe1287 ;--------------1288 ;1289 i_plot = i_plot + 11290 i_plot_h = i_plot MOD numplot_x1291 i_plot_v = i_plot / numplot_x + 11292 1293 ; contour parameters1294 nlevels = 101295 nlevels = 101296 fmin = 0.0 & fmax = MAX(eFeb) ;1297 cmin = 000 & cmax = 2551298 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors1299 1300 ct = 3 ; colortable1301 title = ' eFe (mumol/m3)'1302 1303 icecontour, eFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $1304 tmin, tmax, zmin, zmax, title1305 1306 ; colorbar1307 format='(f5.1)'1308 horpos = 0.02 + i_plot_h * dh1309 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar.1310 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format1311 1312 ;1313 ;------------------------------------------------------------------------------1314 ; End of the script1315 ;------------------------------------------------------------------------------1316 ;1317 1249 IF ( device EQ 'PS' ) THEN BEGIN 1318 1250 DEVICE, /CLOSE … … 1320 1252 !P.MULTI=[0,2,2] 1321 1253 ENDIF 1254 1255 STOP 1256 1257 ; 1258 ;-------------- 1259 ; dFe 1260 ;-------------- 1261 ; 1262 !P.MULTI=[0,numplot_x, numplot_y] 1263 i_plot = 0 1264 i_plot_h = i_plot MOD numplot_x 1265 i_plot_v = i_plot / numplot_x + 1 1266 1267 ; contour parameters 1268 nlevels = 10 1269 nlevels = 10 1270 fmin = 0.0 & fmax = 60. 1271 cmin = 000 & cmax = 255 1272 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 1273 1274 ct = 3 ; colortable 1275 title = ' dFe (mumol/m3)' 1276 1277 icecontour, dFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 1278 tmin, tmax, zmin, zmax, title 1279 1280 ; colorbar 1281 format='(f5.1)' 1282 horpos = 0.02 + i_plot_h * dh 1283 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 1284 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 1285 1286 ; 1287 ;-------------- 1288 ; aFe 1289 ;-------------- 1290 ; 1291 i_plot = i_plot + 1 1292 i_plot_h = i_plot MOD numplot_x 1293 i_plot_v = i_plot / numplot_x + 1 1294 1295 ; contour parameters 1296 nlevels = 10 1297 nlevels = 10 1298 fmin = 0.0 & fmax = MAX(aFeb) ; 1299 cmin = 000 & cmax = 255 1300 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 1301 1302 ct = 3 ; colortable 1303 title = ' aFe (mumol/m3)' 1304 1305 icecontour, aFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 1306 tmin, tmax, zmin, zmax, title 1307 1308 ; colorbar 1309 format='(f5.1)' 1310 horpos = 0.02 + i_plot_h * dh 1311 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 1312 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 1313 1314 ; 1315 ;-------------- 1316 ; eFe 1317 ;-------------- 1318 ; 1319 i_plot = i_plot + 1 1320 i_plot_h = i_plot MOD numplot_x 1321 i_plot_v = i_plot / numplot_x + 1 1322 1323 ; contour parameters 1324 nlevels = 10 1325 nlevels = 10 1326 fmin = 0.0 & fmax = MAX(eFeb) ; 1327 cmin = 000 & cmax = 255 1328 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 1329 1330 ct = 3 ; colortable 1331 title = ' eFe (mumol/m3)' 1332 1333 icecontour, eFeb, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 1334 tmin, tmax, zmin, zmax, title 1335 1336 ; colorbar 1337 format='(f5.1)' 1338 horpos = 0.02 + i_plot_h * dh 1339 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 1340 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 1341 1342 ; 1343 ;------------------------------------------------------------------------------ 1344 ; End of the script 1345 ;------------------------------------------------------------------------------ 1346 ; 1347 IF ( device EQ 'PS' ) THEN BEGIN 1348 DEVICE, /CLOSE 1349 SET_PLOT, "X" 1350 !P.MULTI=[0,2,2] 1351 ENDIF 1322 1352 END 1323 1353 -
branches/2016/dev_v3.20_2016_platelet/GRAPHICS/ALL_stocks.pro
r6 r59 213 213 ; YROSIAE 214 214 ; taken from yrosiae ASPECT log sheets (Fripiat et al, oct 2015) 215 doy_obs = [ 263 , 299, 318 , 279 , 306 , 327, 292 , 312 ,335 ] 216 Ichla_obs = [ 1.4683 , 4.6038 ,45.4413, 128.0884, 44.5830 ,219.1416 ,208.6254 ,108.0657 , 87.9038 ] 217 ISi_obs = [ 22.6631 , 21.1690 , 19.2174, 18.8949 , 18.9943, 18.9517, 20.5963, 16.6855, 13.8374 ] 218 IN_obs = [ 9.0613 , 5.6130 , 9.4723, 10.0763, 11.4708 , 14.6193 , 14.4107 , 4.7585 , 2.1876 ] 219 IP_obs = [ 0.4976 , 0.4067 , 1.1311, 1.8459, 2.3452 , 3.7857 , 3.6948 , 1.3270 , 1.2602 ] 220 ENDIF 221 222 IF ( site EQ 'N_ICE_RL' ) THEN BEGIN 223 ; N_ICE_RL 224 ; ALLLLLL WRONG !!! 215 225 doy_obs = [ 263 , 299, 318 , 279 , 306 , 327, 292 , 312 ,335 ] 216 226 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 493 503 494 495 ymin = 0.0000 &ymax = MAX(dFet)*1.2496 497 498 499 500 501 502 503 504 505 506 507 508 509 504 ; ; Iron stocks 505 ; ymin = 0.0000 & ymax = 500. ;ymax = MAX(dFet)*1.2 506 507 ; PLOT, [xmin, xmax], [ ymin, ymax ], /NODATA, charsize = cs, $ 508 ; XTITLE = 'mon', YTITLE = 'Iron (mumol/m2)', XSTYLE = 1, YSTYLE = 1 509 510 ; FOR i_run = nruns - 1, 0, -1 DO BEGIN 511 ; LOADCT, ct 512 ; OPLOT, mon(i_run,0:ntss(i_run)-1), dFet(0:ntss(i_run)-1), color = 150, thick = 2 513 ; OPLOT, mon(i_run,0:ntss(i_run)-1), aFet(0:ntss(i_run)-1), color = 0, thick = 2 514 ; OPLOT, mon(i_run,0:ntss(i_run)-1), eFet(0:ntss(i_run)-1), color = 0, thick = 1 515 ; ENDFOR 516 ; LOADCT, 0 517 518 ; ztitle = 'Max dFe stock: '+STRCOMPRESS(STRING(MAX(dFet), FORMAT='(F6.2)'),/REMOVE_ALL)+' mumol/m2' 519 ; XYOUTS, xmin+ (xmax-xmin)*0.10, ymax*0.90, ztitle 510 520 511 521
Note: See TracChangeset
for help on using the changeset viewer.