- Timestamp:
- 10/05/16 14:26:39 (8 years ago)
- Location:
- branches/2016/dev_v3.20_2016_gravity_drainage/GRAPHICS
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_v3.20_2016_gravity_drainage/GRAPHICS/ALL_contours.pro
r4 r22 17 17 indir = '' 18 18 outdir = '' 19 dt = 8640019 dt = 3600. 20 20 c_bio_model = 'BFMSI' 21 21 site = 'TOURNAI' … … 87 87 88 88 zunit = 12. 89 AoCb = AoCb * zunit90 eoCb = eoCb * zunit91 AoCt = AoCt * zunit92 eoCt = eoCt * zunit93 syn = syn * zunit * 86400.94 lys = lys * zunit * 86400.95 rem = rem * zunit * 86400.89 ;AoCb = AoCb * zunit 90 ;eoCb = eoCb * zunit 91 ;AoCt = AoCt * zunit 92 ;eoCt = eoCt * zunit 93 syn = syn * zunit * dt 94 lys = lys * zunit * dt 95 rem = rem * zunit * dt 96 96 97 97 ENDIF … … 126 126 127 127 IF ( site EQ 'YROSIAE' ) THEN BEGIN 128 mon = numt / 30.41666666 + 1 129 ENDIF 128 truc = numt / 24 129 mon = truc / 30.4666666 + 1 130 ENDIF 131 132 ;IF ( site EQ 'YROSIAE' ) THEN BEGIN 133 ; mon = numt / 30.41666666 + 1 134 ;ENDIF 130 135 131 136 ; … … 239 244 240 245 ; contour parameters 241 nlevels = 1 3242 levels = [ 0 ,1,2,3,4,5,6,7,8,9,10,11,12]246 nlevels = 12 247 levels = [ 0.,1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,11.,12.] 243 248 colors = [ 241,221,201,181,161,141,121,101,81,61,41,21,1] 244 249 … … 247 252 tmax = MAX(mon) 248 253 249 title = ' Salinit y'254 title = ' Salinite totale ' 250 255 251 256 icecontour, s_i, h_i, h_s, nlevels, colors, levels, ct, z_ip, mon, nlay, nts, $ 252 257 tmin, tmax, zmin, zmax, title 253 258 254 LOADCT, 13 255 oplot, mon_obs, -hi_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR =255 256 oplot, mon_obs, hs_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR = 255 257 LOADCT, ct 258 259 ; colorbar 260 horpos = 0.02 + i_plot_h * dh 261 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 262 format='(f5.0)' 259 ;LOADCT, 13 260 ;oplot, mon_obs, -hi_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR =255 261 ;oplot, mon_obs, hs_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR = 255 262 ;LOADCT, ct 263 264 ; colorbar 265 horpos = 0.02 + i_plot_h * dh 266 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 267 format='(i3)' 268 ;'(f5.0)' 263 269 colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 264 270 … … 271 277 272 278 ; contour parameters 273 nlevels = 1 1279 nlevels = 10 274 280 levels = [ 0., 2., 4., 5., 8., 10., 15.,20., 25., 30., 50.] 275 281 colors = [ 20, 60, 100, 150, 160, 170, 180, 190, 200, 210, 250 ] … … 277 283 278 284 ct = 33 ; colortable 279 title = ' Brine volume'285 title = ' Fraction de saumures ' 280 286 281 287 icecontour, e_i, h_i, h_s, nlevels, colors, levels, ct, z_ip, mon, nlay, nts, $ … … 303 309 304 310 ct = 11 ; colortable 305 title = ' Rayleigh number'311 title = 'Nombre de Rayleigh ' 306 312 307 313 ; contours … … 325 331 326 332 ; contour parameters 327 nlevels = 10333 nlevels = 9 328 334 levels = [ 0.0,0.015, 0.1, 1.0, 10., 15., 20., 30., 50., 80.,120., 150. ] 329 335 levels = [ 0.0,0.015, 0.1, 1.0, 10., 15., 20., 30., 40., 50., 60., 70. ] 330 levels = [ 0.0, 0.05, 0.1, 2.5, 5., 7.5, 10., 12.5, 15., 17.5,20. ]331 colors = [ 20, 40, 60, 80, 140, 160, 180, 190, 200, 210, 220 ]336 levels = [ 0.0, .1, 2.5, 5., 7.5, 10., 12.5, 15., 17.5, 20. ] 337 colors = [ 20, 40, 60, 80, 120, 140, 180, 200, 220, 240 ] 332 338 333 339 ct = 3 ; colortable … … 337 343 tmin, tmax, zmin, zmax, title 338 344 339 LOADCT, 1340 oplot, mon_obs, -hi_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR = 150341 oplot, mon_obs, hs_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR = 150342 LOADCT, ct345 ;LOADCT, 1 346 ;oplot, mon_obs, -hi_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR = 150 347 ;oplot, mon_obs, hs_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR = 150 348 ;LOADCT, ct 343 349 344 350 ; colorbar … … 364 370 365 371 ct = 3 ; colortable 366 title = ' Temperature '372 title = ' Temperature (Celsius)' 367 373 368 374 icecontour, t_i, h_i, h_s, nlevels, colors, levels, ct, z_ip, mon, nlay, nts, $ … … 370 376 371 377 ; colorbar 372 format='( f5.0)'378 format='(i3)' 373 379 horpos = 0.02 + i_plot_h * dh 374 380 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. … … 394 400 395 401 ct = 1 ; colortable 396 title = ' Brine Salinity'402 title = ' Salinite des saumures ' 397 403 398 404 icecontour, sigma, h_i, h_s, nlevels, colors, levels, ct, z_ip, mon, nlay, nts, $ … … 432 438 ENDIF 433 439 434 ct = 5; colortable435 title = ' NO3 (mmol/m3)'440 ct = 7 ; colortable 441 title = ' DIN total (mmolN/m3)' 436 442 437 443 icecontour, no3b, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ … … 469 475 ENDIF 470 476 471 ct = 5; colortable472 title = ' PO4 (mmol/m3)'477 ct = 7 ; colortable 478 title = ' DIP total (mmolP/m3)' 473 479 474 480 icecontour, PO4b, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ … … 509 515 ENDIF 510 516 511 ct = 5; colortable512 title = ' DSi (mmol/m3)'517 ct = 7 ; colortable 518 title = ' DSi total (mmolSi/m3)' 513 519 514 520 icecontour, dsib, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ … … 537 543 538 544 ct = 3 ; colortable 539 title = ' Light limitation'545 title = ' limitation Lumiere ' 540 546 541 547 ; contours … … 565 571 566 572 ct = 3 ; colortable 567 title = ' Temp limitation'573 title = ' limitation Temperature ' 568 574 569 575 ; contours … … 592 598 593 599 ct = 3 ; colortable 594 title = ' Salt limitation'600 title = ' Inhibition Salinite ' 595 601 596 602 ; contours … … 619 625 620 626 ct = 3 ; colortable 621 title = ' NO3 limitation'627 title = ' limitation DIN ' 622 628 623 629 ; contours … … 646 652 647 653 ct = 3 ; colortable 648 title = ' PO4 limitation'654 title = ' limitation DIP ' 649 655 650 656 ; contours … … 673 679 674 680 ct = 3 ; colortable 675 title = ' DSi limitation'681 title = ' limitation DSi ' 676 682 677 683 ; contours … … 731 737 nlevels = 8 732 738 fmin = 0.0 & fmax = MAX(AoCb) 733 fmin = 0.0 & fmax = 1000.739 ;fmin = 0.0 & fmax = 1000. 734 740 cmin = 255. & cmax = 80. 735 741 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 736 742 737 743 ct = 8 ; colortable 738 title = ' AoC b (mgC/m3)'744 title = ' AoC totale (mmolC/m3)' 739 745 740 746 ; contours … … 742 748 tmin, tmax, zmin, zmax, title 743 749 ; colorbar 744 format='(f5.0)' 750 ;format='(f5.0)' 751 format='(i3)' 745 752 horpos = 0.02 + i_plot_h * dh 746 753 vertpos = 0.0155 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. … … 759 766 nlevels = 8 760 767 fmin = 0.0 & fmax = MAX(eOCb) 761 fmin = 0.0 & fmax = 1000.768 ;fmin = 0.0 & fmax = 1000. 762 769 cmin = 255. & cmax = 80. 763 770 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 764 771 765 772 ct = 8 ; colortable 766 title = ' eoC b (mgC/m3)'773 title = ' eoC totale (mmolC/m3)' 767 774 768 775 ; contours … … 770 777 tmin, tmax, zmin, zmax, title 771 778 ; colorbar 772 format='(f5.0)' 779 ;format='(f5.0)' 780 format='(i3)' 773 781 horpos = 0.02 + i_plot_h * dh 774 782 vertpos = 0.0155 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. … … 787 795 nlevels = 10 788 796 fmin = 0.0 & fmax = MAX(chla) 789 fmin = 0.0 & fmax = 80.797 ;fmin = 0.0 & fmax = 80. 790 798 cmin = 255. & cmax = 80. 791 799 do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors … … 801 809 802 810 ct = 8 ; colortable 803 title = ' chla (mg/m3) '811 title = ' chla totale (mg chla/m3) ' 804 812 805 813 ; contours … … 807 815 tmin, tmax, zmin, zmax, title 808 816 ; colorbar 809 format='(f5.0)' 817 ;format='(f5.0)' 818 format='(i3)' 810 819 horpos = 0.02 + i_plot_h * dh 811 820 vertpos = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. -
branches/2016/dev_v3.20_2016_gravity_drainage/GRAPHICS/YROSIAE_profiles.pro
r4 r22 165 165 ; 166 166 167 POC = eoCb + DAFb 168 167 169 ; observed normalized profiles from Lannuzel et al. (2007) 168 170 … … 274 276 po4_obs(zaddr) = po4_obs(zaddr) * rhoi / rhow 275 277 278 279 280 ; depths for POC 281 N_POC = REPLICATE(6,9) 282 depth_POC = FLTARR(N_obs, 6) 283 depth_POC[0, 0:N_POC(0)-1] = [ 12.5, 37.5, 65.05, 103.5, 128.5, 146. ] / 100. 284 depth_POC[1, 0:N_POC(1)-1] = [ 15, 45 , 75 , 105 , 135 , 155 ] / 100. 285 depth_POC[2, 0:N_POC(2)-1] = [ 16, 47 , 79.5 , 111 , 143 , 164 ] / 100. 286 depth_POC[3, 0:N_POC(3)-1] = [ 16, 48 , 80 , 112 , 144 , 165 ] / 100. 287 depth_POC[4, 0:N_POC(4)-1] = [ 16, 48 , 80.5 , 113.3, 145.8, 167 ] / 100. 288 depth_POC[5, 0:N_POC(5)-1] = [16.35, 49.05, 71.75, 114.45, 147.15, 168.5 ] / 100. 289 depth_POC[6, 0:N_POC(6)-1] = [16.05, 48.15, 80.25, 112.35, 144.45, 165.5 ] / 100. 290 depth_POC[7, 0:N_POC(7)-1] = [ 16, 48, 80, 112, 144, 165 ] / 100. 291 depth_POC[8, 0:N_POC(8)-1] = [ 17, 49.75, 85, 112.75, 144.25, 165 ] / 100. 292 ; POC 293 POC_obs = FLTARR(N_obs, 6) 294 POC_obs[0,0:N_POC(0)-1] = [ 5.4, miss_val, 24.9, 165.0, 0.1, 99.6 ] 295 POC_obs[1,0:N_POC(1)-1] = [ miss_val, 9.8, 12.4, 17.8, 15.2, 110.4 ] 296 POC_obs[2,0:N_POC(2)-1] = [ 12.0, 343.0, 9.4, miss_val, 0.1, 3148.7 ] 297 POC_obs[3,0:N_POC(3)-1] = [ miss_val, 11.3, 19.1, 11.4, 3.8, 2427.6 ] 298 POC_obs[4,0:N_POC(4)-1] = [ 7.7, 1.3, 1.9, 11.8, 31.9, 1726.0 ] 299 POC_obs[5,0:N_POC(5)-1] = [ 19.4, 6.4, 6.0, 12.6, 9.6, 3839.4 ] 300 POC_obs[6,0:N_POC(6)-1] = [ 6.3, 1.5, 2.5, 10.7, 19.9, 1327.0 ] 301 POC_obs[7,0:N_POC(7)-1] = [ 5.1, 10.1, 6.3, 6.5, 23.0, 1445.1 ] 302 POC_obs[8,0:N_POC(8)-1] = [ 3.7, 7.0, 3.5, 3.8, 11.6, 2438.7 ] 303 304 zaddr = WHERE(POC_obs NE miss_val) 305 POC_obs(zaddr) = POC_obs(zaddr) * rhoi / rhow 306 307 308 309 310 276 311 ; physical observations 277 312 … … 364 399 nox_br_obs(iaddr) = nox_obs(iaddr) / e_i_nut(iaddr) 365 400 401 iaddr = where ( POC_obs NE miss_val ) 402 POC_br_obs = FLTARR(N_obs, 6) 403 POC_br_obs(iaddr) = POC_obs(iaddr) / e_i_nut(iaddr) 404 366 405 ;*** Mean observed Profiles 367 406 ;--- interpolate T, S and e on standard depths --- … … 429 468 nox_obs_mean = FLTARR(N_dpt) & nox_obs_std = FLTARR(N_dpt) 430 469 po4_obs_mean = FLTARR(N_dpt) & po4_obs_std = FLTARR(N_dpt) 470 POC_obs_mean = FLTARR(N_dpt) & POC_obs_std = FLTARR(N_dpt) 471 431 472 depth_nut_mean = FLTARR(N_dpt) 432 473 depth_nut_std = FLTARR(N_obs,6) 474 475 depth_POC_mean = FLTARR(N_dpt) 476 depth_POC_std = FLTARR(N_obs,6) 477 433 478 FOR i_obs = 0, N_obs - 1 DO BEGIN 434 479 depth_nut_std(i_obs,*) = depth_nut(i_obs,*) / h_nut(i_obs) 480 depth_POC_std(i_obs,*) = depth_POC(i_obs,*) / h_nut(i_obs) 435 481 ENDFOR 436 482 … … 445 491 po4_obs_std(i) = STDDEV( po4_obs(i_addr,i) ) 446 492 depth_nut_mean(i) = MEAN( depth_nut_std(i_addr,i) ) 493 ENDFOR 494 495 FOR i = 0, 5 DO BEGIN 496 i_addr = WHERE( POC_obs(*,i) NE miss_val ) 497 POC_obs_mean(i) = MEAN (POC_obs(i_addr,i) ) 498 POC_obs_std(i) = STDDEV( POC_obs(i_addr,i) ) 499 depth_POC_mean(i) = MEAN( depth_POC_std(i_addr,i) ) 447 500 ENDFOR 448 501 … … 462 515 dsi_mod_mean = FLTARR(nlay) & dsi_mod_mean(*) = 0. 463 516 dsi_mod_std = FLTARR(nlay) & dsi_mod_std (*) = 0. 517 POC_mod_mean = FLTARR(nlay) & POC_mod_mean(*) = 0. 518 POC_mod_std = FLTARR(nlay) & POC_mod_std (*) = 0. 519 520 464 521 FOR i = 0, nlay - 1 DO BEGIN 465 522 … … 473 530 nox_mod_mean(i) = nox_mod_mean(i) + no3b(i,i_mod) / FLOAT(N_obs) 474 531 po4_mod_mean(i) = po4_mod_mean(i) + po4b(i,i_mod) / FLOAT(N_obs) 532 POC_mod_mean(i) = POC_mod_mean(i) + POC(i,i_mod) / FLOAT(N_obs) 475 533 ENDFOR 476 534 … … 482 540 zsum_n = 0 483 541 zsum_p = 0 542 zsum_j = 0 484 543 FOR i_obs = 0, N_obs - 1 DO BEGIN 485 544 i_mod = WHERE(doy EQ doy_obs(i_obs)-1) … … 491 550 zsum_n = zsum_n + ( no3b(i,i_mod) - nox_mod_mean(i) ) ^2 / FLOAT(N_obs) 492 551 zsum_p = zsum_p + ( po4b(i,i_mod) - po4_mod_mean(i) ) ^2 / FLOAT(N_obs) 552 zsum_j = zsum_j + ( POC(i,i_mod) - POC_mod_mean(i) ) ^2 / FLOAT(N_obs) 493 553 ENDFOR 494 554 t_i_mod_std(i) = SQRT( zsum ) … … 499 559 po4_mod_std(i) = SQRT( zsum_p ) 500 560 dsi_mod_std(i) = SQRT( zsum_d ) 501 561 POC_mod_std(i) = SQRT( zsum_j ) 502 562 ENDFOR 503 563 … … 677 737 LOADCT, 0 678 738 ztitle = STRCOMPRESS(STRING(doy_obs(i_obs)), /REMOVE_ALL) 679 PLOT, [ 0.00, 4.0 ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'mmmol P / m3', XTICKFORMAT='( I2)', XMINOR = 2, YMINOR = 2739 PLOT, [ 0.00, 4.0 ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'mmmol P / m3', XTICKFORMAT='(F1.2)', XMINOR = 2, YMINOR = 2 680 740 XYOUTS, 0.37, -0.3, ztitle 681 741 … … 819 879 ;============================================================================== 820 880 ; 881 numplot_x = 5 882 numplot_y = 2 883 821 884 !P.MULTI=[0,numplot_x, numplot_y] 822 885 … … 825 888 ; prepare plot 826 889 LOADCT, 0 827 ztitle = " Mean T"828 PLOT, [ -20., 0. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = ' deg C', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle890 ztitle = "Temperature (C)" 891 PLOT, [ -20., 0. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 829 892 830 893 ; add obs … … 844 907 ; prepare plot 845 908 LOADCT, 0 846 ztitle = " Mean S"847 PLOT, [0., 1 2. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'g/kg', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle909 ztitle = "Salinite totale" 910 PLOT, [0., 15. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 848 911 849 912 ; add obs … … 865 928 ; prepare plot 866 929 LOADCT, 0 867 ztitle = " Mean e"868 PLOT, [0., 50. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = ' %', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle930 ztitle = "Fraction de saumures" 931 PLOT, [0., 50. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 869 932 870 933 ; add obs … … 888 951 ; prepare plot 889 952 LOADCT, 0 890 ztitle = " Mean chla"891 PLOT, [ 0. 01, 40. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'mg chla / m3', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle953 ztitle = "chla (mg/m3)" 954 PLOT, [ 0., 2500. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '',XTICKINTERVAL=1500, XTICKFORMAT='(I4)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 892 955 ; add obs 893 956 OPLOT, chla_obs_mean, - depth_std, PSYM = 8, THICK = 3, SYMSIZE = 0.8 … … 897 960 ENDFOR 898 961 ; plot seawater value 899 LOADCT, 1900 OPLOT, [ 0.08 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150962 ;LOADCT, 1 963 ;OPLOT, [ 0.08 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 901 964 902 965 ; add model … … 911 974 ; prepare plot 912 975 LOADCT, 0 913 ztitle = " Mean NOX"914 PLOT, [ 0. 00, 35. ], [ -1.0, 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'mmmol N / m3', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle976 ztitle = "DIN total (mmolN/m3)" 977 PLOT, [ 0., 12. ], [ -1.0, 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 915 978 ; add obs 916 OPLOT, nox_obs_mean, - depth_nut_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8917 FOR i = 0, 5 DO BEGIN918 OPLOT, [ nox_obs_mean(i) - nox_obs_std(i), nox_obs_mean(i) + nox_obs_std(i) ], $919 [ -depth_nut_mean(i), -depth_nut_mean(i) ]920 ENDFOR979 ;OPLOT, nox_obs_mean, - depth_nut_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8 980 ;FOR i = 0, 5 DO BEGIN 981 ; OPLOT, [ nox_obs_mean(i) - nox_obs_std(i), nox_obs_mean(i) + nox_obs_std(i) ], $ 982 ; [ -depth_nut_mean(i), -depth_nut_mean(i) ] 983 ;ENDFOR 921 984 922 985 ; plot seawater value 923 LOADCT, 1924 OPLOT, [ 31.16 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150986 ;LOADCT, 1 987 ;OPLOT, [ 31.16 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 925 988 926 989 ; add model … … 935 998 ; prepare plot 936 999 LOADCT, 0 937 ztitle = " Mean PO4"938 PLOT, [ 0. 00, 4.0 ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'mmmol P / m3', XTICKFORMAT='(I2)', XMINOR = 2, YMINOR = 2, TITLE = ztitle1000 ztitle = "DIP total (mmolP/m3)" 1001 PLOT, [ 0., 1. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(F4.1)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 939 1002 ; add obs 940 OPLOT, po4_obs_mean, - depth_nut_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8941 FOR i = 0, N_std - 1 DO BEGIN942 OPLOT, [ po4_obs_mean(i) - po4_obs_std(i), po4_obs_mean(i) + po4_obs_std(i) ], $943 [ -depth_nut_mean(i), -depth_nut_mean(i) ]944 ENDFOR1003 ;OPLOT, po4_obs_mean, - depth_nut_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8 1004 ;FOR i = 0, N_std - 1 DO BEGIN 1005 ; OPLOT, [ po4_obs_mean(i) - po4_obs_std(i), po4_obs_mean(i) + po4_obs_std(i) ], $ 1006 ; [ -depth_nut_mean(i), -depth_nut_mean(i) ] 1007 ;ENDFOR 945 1008 946 1009 ; plot seawater value 947 LOADCT, 1948 OPLOT, [ 2.05 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 1501010 ;LOADCT, 1 1011 ;OPLOT, [ 2.05 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 949 1012 950 1013 ; add model … … 959 1022 ; prepare plot 960 1023 LOADCT, 0 961 ztitle = " Mean DSi"962 PLOT, [ 0. 00, 90. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'mmmol N / m3', XTICKFORMAT='(I2)', XMINOR = 2, YMINOR = 2, TITLE = ztitle1024 ztitle = "DSi total (mmolSi/m3)" 1025 PLOT, [ 0., 30. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I2)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 963 1026 ; add obs 964 OPLOT, dsi_obs_mean, - depth_nut_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8965 FOR i = 0, N_std - 1 DO BEGIN966 OPLOT, [ dsi_obs_mean(i) - dsi_obs_std(i), dsi_obs_mean(i) + dsi_obs_std(i) ], $967 [ -depth_nut_mean(i), -depth_nut_mean(i) ]968 ENDFOR1027 ;OPLOT, dsi_obs_mean, - depth_nut_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8 1028 ;FOR i = 0, N_std - 1 DO BEGIN 1029 ; OPLOT, [ dsi_obs_mean(i) - dsi_obs_std(i), dsi_obs_mean(i) + dsi_obs_std(i) ], $ 1030 ; [ -depth_nut_mean(i), -depth_nut_mean(i) ] 1031 ;ENDFOR 969 1032 970 1033 ; plot seawater value … … 978 1041 OPLOT, dsi_mod_mean+dsi_mod_std, -zz, THICK = 1, COLOR = 150 979 1042 OPLOT, dsi_mod_mean-dsi_mod_std, -zz, THICK = 1, COLOR = 150 1043 1044 1045 ;--- POC --- 1046 1047 ; prepare plot 1048 LOADCT, 0 1049 ztitle = "POC total (mmolC/m3)" 1050 PLOT, [ 0., 3000. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKINTERVAL=1500, XTICKFORMAT='(I4)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 1051 ; add obs 1052 OPLOT, POC_obs_mean, - depth_POC_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8 1053 FOR i = 0, N_std - 1 DO BEGIN 1054 OPLOT, [ POC_obs_mean(i) - POC_obs_std(i), POC_obs_mean(i) + POC_obs_std(i) ], $ 1055 [ -depth_POC_mean(i), -depth_POC_mean(i) ] 1056 ENDFOR 1057 1058 ; plot seawater value 1059 ;LOADCT, 1 1060 ;OPLOT, [ 84.35 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 1061 1062 ; add model 1063 LOADCT, 9 1064 OPLOT, POC_mod_mean, -zz, THICK = 3, COLOR = 150 1065 LOADCT, 0 1066 OPLOT, POC_mod_mean+POC_mod_std, -zz, THICK = 1, COLOR = 150 1067 OPLOT, POC_mod_mean-POC_mod_std, -zz, THICK = 1, COLOR = 150 1068 1069 ; 980 1070 981 1071 ;
Note: See TracChangeset
for help on using the changeset viewer.