Changeset 22 for branches/2016


Ignore:
Timestamp:
10/05/16 14:26:39 (8 years ago)
Author:
vancop
Message:

Ludivine graphic files

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  
    1717indir  = '' 
    1818outdir = '' 
    19 dt     = 86400 
     19dt     = 3600. 
    2020c_bio_model = 'BFMSI' 
    2121site        = 'TOURNAI' 
    … …  
    8787 
    8888zunit = 12. 
    89 AoCb = AoCb * zunit 
    90 eoCb = eoCb * zunit 
    91 AoCt = AoCt * zunit 
    92 eoCt = eoCt * zunit 
    93 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 
     93syn    = syn       * zunit * dt 
     94lys    = lys       * zunit * dt 
     95rem    = rem       * zunit * dt 
    9696 
    9797ENDIF 
    … …  
    126126 
    127127IF ( site EQ 'YROSIAE' ) THEN BEGIN 
    128    mon = numt / 30.41666666 + 1 
    129 ENDIF 
     128   truc = numt / 24 
     129   mon  = truc / 30.4666666 + 1 
     130ENDIF 
     131 
     132;IF ( site EQ 'YROSIAE' ) THEN BEGIN 
     133;   mon = numt / 30.41666666 + 1 
     134;ENDIF 
    130135 
    131136; 
    … …  
    239244 
    240245; contour parameters 
    241 nlevels = 13 
    242 levels = [ 0,1,2,3,4,5,6,7,8,9,10,11,12] 
     246nlevels = 12 
     247levels = [ 0.,1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,11.,12.] 
    243248colors = [ 241,221,201,181,161,141,121,101,81,61,41,21,1] 
    244249 
    … …  
    247252tmax = MAX(mon) 
    248253 
    249 title = ' Salinity ' 
     254title = ' Salinite totale ' 
    250255 
    251256icecontour, s_i, h_i, h_s, nlevels, colors, levels, ct, z_ip, mon, nlay, nts, $ 
    252257            tmin, tmax, zmin, zmax, title 
    253258 
    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 
     265horpos   = 0.02 + i_plot_h * dh 
     266vertpos  = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
     267format='(i3)' 
     268;'(f5.0)' 
    263269colorbarn, horpos, vertpos, cb_length, 0.02, MIN(levels), MAX(levels), leveling, nlevels+1, levels, colors, format 
    264270 
    … …  
    271277 
    272278; contour parameters 
    273 nlevels = 11 
     279nlevels = 10 
    274280levels = [   0.,   2.,  4.,   5.,   8.,  10.,  15.,20., 25., 30., 50.] 
    275281colors = [   20,  60,  100,  150,  160,  170, 180, 190, 200, 210, 250 ] 
    … …  
    277283 
    278284ct = 33 ; colortable 
    279 title = ' Brine volume ' 
     285title = ' Fraction de saumures ' 
    280286 
    281287icecontour, e_i, h_i, h_s, nlevels, colors, levels, ct, z_ip, mon, nlay, nts, $ 
    … …  
    303309 
    304310ct = 11 ; colortable 
    305 title = ' Rayleigh number ' 
     311title = 'Nombre de Rayleigh ' 
    306312 
    307313; contours 
    … …  
    325331 
    326332; contour parameters 
    327 nlevels = 10 
     333nlevels = 9 
    328334levels = [  0.0,0.015, 0.1,  1.0,  10.,  15.,  20.,  30., 50., 80.,120., 150. ] 
    329335levels = [  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 ] 
     336levels = [  0.0, .1,  2.5,  5., 7.5,  10.,  12.5,  15.,  17.5, 20. ] 
     337colors = [   20, 40,   60,  80, 120,  140,   180,  200,  220,   240 ] 
    332338 
    333339ct = 3 ; colortable 
    … …  
    337343            tmin, tmax, zmin, zmax, title 
    338344 
    339 LOADCT, 1 
    340 oplot, mon_obs, -hi_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR = 150 
    341 oplot, mon_obs, hs_obs, PSYM=1, SYMSIZE=1.0, THICK = 3, COLOR = 150 
    342 LOADCT, ct 
     345;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 
    343349 
    344350; colorbar 
    … …  
    364370 
    365371ct = 3 ; colortable 
    366 title = ' Temperature ' 
     372title = ' Temperature (Celsius)' 
    367373 
    368374icecontour, t_i, h_i, h_s, nlevels, colors, levels, ct, z_ip, mon, nlay, nts, $ 
    … …  
    370376 
    371377; colorbar 
    372 format='(f5.0)' 
     378format='(i3)' 
    373379horpos   = 0.02 + i_plot_h * dh 
    374380vertpos  = 0.015 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
    … …  
    394400 
    395401ct = 1 ; colortable 
    396 title = ' Brine Salinity ' 
     402title = ' Salinite des saumures ' 
    397403 
    398404icecontour, sigma, h_i, h_s, nlevels, colors, levels, ct, z_ip, mon, nlay, nts, $ 
    … …  
    432438ENDIF 
    433439 
    434 ct =  5 ; colortable 
    435 title = ' NO3 (mmol/m3)' 
     440ct =  7 ; colortable 
     441title = ' DIN total (mmolN/m3)' 
    436442 
    437443icecontour, no3b, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 
    … …  
    469475ENDIF 
    470476 
    471 ct = 5 ; colortable 
    472 title = ' PO4 (mmol/m3)' 
     477ct = 7 ; colortable 
     478title = ' DIP total (mmolP/m3)' 
    473479 
    474480icecontour, PO4b, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 
    … …  
    509515ENDIF 
    510516 
    511 ct = 5 ; colortable 
    512 title = ' DSi (mmol/m3)' 
     517ct = 7 ; colortable 
     518title = ' DSi total (mmolSi/m3)' 
    513519 
    514520icecontour, dsib, h_i, h_s, nlevels, colors, levels, ct, z_ib, mon, nlay_bio, nts, $ 
    … …  
    537543 
    538544ct = 3  ; colortable 
    539 title = ' Light limitation ' 
     545title = ' limitation Lumiere ' 
    540546 
    541547; contours 
    … …  
    565571 
    566572ct = 3  ; colortable 
    567 title = ' Temp limitation ' 
     573title = ' limitation Temperature ' 
    568574 
    569575; contours 
    … …  
    592598 
    593599ct = 3  ; colortable 
    594 title = ' Salt limitation ' 
     600title = ' Inhibition Salinite ' 
    595601 
    596602; contours 
    … …  
    619625 
    620626ct = 3 ; colortable 
    621 title = ' NO3 limitation ' 
     627title = ' limitation DIN ' 
    622628 
    623629; contours 
    … …  
    646652 
    647653ct = 3 ; colortable 
    648 title = ' PO4 limitation ' 
     654title = ' limitation DIP ' 
    649655 
    650656; contours 
    … …  
    673679 
    674680ct = 3 ; colortable 
    675 title = ' DSi limitation ' 
     681title = ' limitation DSi ' 
    676682 
    677683; contours 
    … …  
    731737nlevels = 8 
    732738fmin = 0.0 & fmax = MAX(AoCb) 
    733 fmin = 0.0 & fmax = 1000. 
     739;fmin = 0.0 & fmax = 1000. 
    734740cmin = 255. & cmax = 80. 
    735741do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 
    736742 
    737743ct = 8 ; colortable 
    738 title = ' AoCb (mgC/m3)' 
     744title = ' AoC totale (mmolC/m3)' 
    739745 
    740746; contours 
    … …  
    742748            tmin, tmax, zmin, zmax, title 
    743749; colorbar 
    744 format='(f5.0)' 
     750;format='(f5.0)' 
     751format='(i3)' 
    745752horpos   = 0.02 + i_plot_h * dh 
    746753vertpos  = 0.0155 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
    … …  
    759766nlevels = 8 
    760767fmin = 0.0 & fmax = MAX(eOCb) 
    761 fmin = 0.0 & fmax = 1000. 
     768;fmin = 0.0 & fmax = 1000. 
    762769cmin = 255. & cmax = 80. 
    763770do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 
    764771 
    765772ct = 8 ; colortable 
    766 title = ' eoCb (mgC/m3)' 
     773title = ' eoC totale (mmolC/m3)' 
    767774 
    768775; contours 
    … …  
    770777            tmin, tmax, zmin, zmax, title 
    771778; colorbar 
    772 format='(f5.0)' 
     779;format='(f5.0)' 
     780format='(i3)' 
    773781horpos   = 0.02 + i_plot_h * dh 
    774782vertpos  = 0.0155 + ( numplot_y - i_plot_v ) * dv ; position of colorbar. 
    … …  
    787795nlevels = 10 
    788796fmin = 0.0 & fmax = MAX(chla) 
    789 fmin = 0.0 & fmax = 80. 
     797;fmin = 0.0 & fmax = 80. 
    790798cmin = 255. & cmax = 80. 
    791799do_levels, fmin, fmax, cmin, cmax, nlevels, levels, colors 
    … …  
    801809 
    802810ct = 8 ; colortable 
    803 title = ' chla (mg/m3) ' 
     811title = ' chla totale (mg chla/m3) ' 
    804812 
    805813; contours 
    … …  
    807815            tmin, tmax, zmin, zmax, title 
    808816; colorbar 
    809 format='(f5.0)' 
     817;format='(f5.0)' 
     818format='(i3)' 
    810819horpos   = 0.02 + i_plot_h * dh 
    811820vertpos  = 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  
    165165; 
    166166 
     167POC = eoCb + DAFb 
     168 
    167169; observed normalized profiles from Lannuzel et al. (2007) 
    168170 
    … …  
    274276po4_obs(zaddr) = po4_obs(zaddr)  * rhoi / rhow 
    275277 
     278 
     279 
     280; depths for POC 
     281N_POC = REPLICATE(6,9) 
     282depth_POC = FLTARR(N_obs, 6) 
     283depth_POC[0, 0:N_POC(0)-1] = [ 12.5, 37.5, 65.05, 103.5, 128.5, 146. ] / 100. 
     284depth_POC[1, 0:N_POC(1)-1] = [   15, 45  , 75   , 105  , 135  , 155  ] / 100. 
     285depth_POC[2, 0:N_POC(2)-1] = [   16, 47  , 79.5 , 111  , 143  , 164  ] / 100. 
     286depth_POC[3, 0:N_POC(3)-1] = [   16, 48  , 80   , 112  , 144  , 165  ] / 100. 
     287depth_POC[4, 0:N_POC(4)-1] = [   16, 48  , 80.5 , 113.3, 145.8, 167  ] / 100. 
     288depth_POC[5, 0:N_POC(5)-1] = [16.35, 49.05, 71.75, 114.45, 147.15, 168.5 ] / 100. 
     289depth_POC[6, 0:N_POC(6)-1] = [16.05, 48.15, 80.25, 112.35, 144.45, 165.5 ] / 100. 
     290depth_POC[7, 0:N_POC(7)-1] = [   16, 48, 80, 112, 144, 165 ] / 100. 
     291depth_POC[8, 0:N_POC(8)-1] = [ 17, 49.75, 85, 112.75, 144.25, 165 ] / 100. 
     292; POC 
     293POC_obs = FLTARR(N_obs, 6) 
     294POC_obs[0,0:N_POC(0)-1] = [      5.4, miss_val,   24.9,    165.0,    0.1,   99.6 ] 
     295POC_obs[1,0:N_POC(1)-1] = [ miss_val,      9.8,   12.4,     17.8,   15.2,  110.4 ] 
     296POC_obs[2,0:N_POC(2)-1] = [     12.0,    343.0,    9.4, miss_val,    0.1, 3148.7 ] 
     297POC_obs[3,0:N_POC(3)-1] = [ miss_val,     11.3,   19.1,     11.4,    3.8, 2427.6 ] 
     298POC_obs[4,0:N_POC(4)-1] = [      7.7,      1.3,    1.9,     11.8,   31.9, 1726.0 ] 
     299POC_obs[5,0:N_POC(5)-1] = [     19.4,      6.4,    6.0,     12.6,    9.6, 3839.4 ] 
     300POC_obs[6,0:N_POC(6)-1] = [      6.3,      1.5,    2.5,     10.7,   19.9, 1327.0 ] 
     301POC_obs[7,0:N_POC(7)-1] = [      5.1,     10.1,    6.3,      6.5,   23.0, 1445.1 ] 
     302POC_obs[8,0:N_POC(8)-1] = [      3.7,      7.0,    3.5,      3.8,   11.6, 2438.7 ] 
     303 
     304zaddr = WHERE(POC_obs NE miss_val) 
     305POC_obs(zaddr) = POC_obs(zaddr)  * rhoi / rhow 
     306 
     307 
     308 
     309 
     310 
    276311; physical observations 
    277312 
    … …  
    364399nox_br_obs(iaddr) = nox_obs(iaddr) / e_i_nut(iaddr) 
    365400 
     401iaddr = where ( POC_obs NE miss_val ) 
     402POC_br_obs = FLTARR(N_obs, 6) 
     403POC_br_obs(iaddr) = POC_obs(iaddr) / e_i_nut(iaddr) 
     404 
    366405;*** Mean observed Profiles 
    367406;--- interpolate T, S and e on standard depths --- 
    … …  
    429468nox_obs_mean = FLTARR(N_dpt) & nox_obs_std = FLTARR(N_dpt) 
    430469po4_obs_mean = FLTARR(N_dpt) & po4_obs_std = FLTARR(N_dpt) 
     470POC_obs_mean = FLTARR(N_dpt) & POC_obs_std = FLTARR(N_dpt) 
     471 
    431472depth_nut_mean = FLTARR(N_dpt) 
    432473depth_nut_std = FLTARR(N_obs,6) 
     474 
     475depth_POC_mean = FLTARR(N_dpt) 
     476depth_POC_std = FLTARR(N_obs,6) 
     477 
    433478FOR i_obs = 0, N_obs - 1 DO BEGIN 
    434479   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) 
    435481ENDFOR 
    436482 
    … …  
    445491   po4_obs_std(i)  = STDDEV( po4_obs(i_addr,i) ) 
    446492   depth_nut_mean(i) = MEAN( depth_nut_std(i_addr,i) ) 
     493ENDFOR 
     494 
     495FOR 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) ) 
    447500ENDFOR 
    448501 
    … …  
    462515dsi_mod_mean = FLTARR(nlay)  & dsi_mod_mean(*) = 0. 
    463516dsi_mod_std  = FLTARR(nlay)  & dsi_mod_std (*) = 0. 
     517POC_mod_mean = FLTARR(nlay)  & POC_mod_mean(*) = 0. 
     518POC_mod_std = FLTARR(nlay)   & POC_mod_std (*) = 0. 
     519 
     520 
    464521FOR i = 0, nlay - 1 DO BEGIN 
    465522 
    … …  
    473530      nox_mod_mean(i)  = nox_mod_mean(i) + no3b(i,i_mod) / FLOAT(N_obs) 
    474531      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) 
    475533   ENDFOR 
    476534 
    … …  
    482540   zsum_n = 0 
    483541   zsum_p = 0 
     542   zsum_j = 0 
    484543   FOR i_obs = 0, N_obs - 1 DO BEGIN 
    485544      i_mod = WHERE(doy EQ doy_obs(i_obs)-1) 
    … …  
    491550      zsum_n = zsum_n + ( no3b(i,i_mod) - nox_mod_mean(i) ) ^2 / FLOAT(N_obs) 
    492551      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) 
    493553   ENDFOR 
    494554   t_i_mod_std(i)  = SQRT( zsum   ) 
    … …  
    499559   po4_mod_std(i)  = SQRT( zsum_p ) 
    500560   dsi_mod_std(i)  = SQRT( zsum_d ) 
    501  
     561   POC_mod_std(i)  = SQRT( zsum_j )  
    502562ENDFOR 
    503563 
    … …  
    677737   LOADCT, 0 
    678738   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 = 2 
     739   PLOT, [ 0.00, 4.0 ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'mmmol P / m3', XTICKFORMAT='(F1.2)', XMINOR = 2, YMINOR = 2 
    680740   XYOUTS, 0.37, -0.3, ztitle 
    681741 
    … …  
    819879;============================================================================== 
    820880; 
     881numplot_x = 5 
     882numplot_y = 2 
     883 
    821884!P.MULTI=[0,numplot_x, numplot_y] 
    822885 
    … …  
    825888; prepare plot 
    826889LOADCT, 0 
    827 ztitle = "Mean T" 
    828 PLOT, [ -20., 0. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'deg C', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
     890ztitle = "Temperature (C)" 
     891PLOT, [ -20., 0. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
    829892 
    830893; add obs 
    … …  
    844907; prepare plot 
    845908LOADCT, 0 
    846 ztitle = "Mean S" 
    847 PLOT, [0., 12. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = 'g/kg', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
     909ztitle = "Salinite totale" 
     910PLOT, [0., 15. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
    848911 
    849912; add obs 
    … …  
    865928; prepare plot 
    866929LOADCT, 0 
    867 ztitle = "Mean e" 
    868    PLOT, [0., 50. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '%', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
     930ztitle = "Fraction de saumures" 
     931   PLOT, [0., 50. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
    869932 
    870933; add obs 
    … …  
    888951; prepare plot 
    889952LOADCT, 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 = ztitle 
     953ztitle = "chla (mg/m3)" 
     954PLOT, [ 0., 2500. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '',XTICKINTERVAL=1500, XTICKFORMAT='(I4)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
    892955; add obs 
    893956OPLOT, chla_obs_mean, - depth_std, PSYM = 8, THICK = 3, SYMSIZE = 0.8 
    … …  
    897960ENDFOR 
    898961; plot seawater value 
    899 LOADCT, 1 
    900 OPLOT, [ 0.08 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 
     962;LOADCT, 1 
     963;OPLOT, [ 0.08 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 
    901964 
    902965; add model 
    … …  
    911974; prepare plot 
    912975LOADCT, 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 = ztitle 
     976ztitle = "DIN total (mmolN/m3)" 
     977PLOT, [ 0., 12. ], [ -1.0, 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I3)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
    915978; add obs 
    916 OPLOT, nox_obs_mean, - depth_nut_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8 
    917 FOR i = 0, 5 DO BEGIN 
    918    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 ENDFOR 
     979;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 
    921984 
    922985; plot seawater value 
    923 LOADCT, 1 
    924 OPLOT, [ 31.16 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 
     986;LOADCT, 1 
     987;OPLOT, [ 31.16 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 
    925988 
    926989; add model 
    … …  
    935998; prepare plot 
    936999LOADCT, 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 = ztitle 
     1000ztitle = "DIP total (mmolP/m3)" 
     1001PLOT, [ 0., 1. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(F4.1)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
    9391002; add obs 
    940 OPLOT, po4_obs_mean, - depth_nut_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8 
    941 FOR i = 0, N_std - 1 DO BEGIN 
    942    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 ENDFOR 
     1003;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 
    9451008 
    9461009; plot seawater value 
    947 LOADCT, 1 
    948 OPLOT, [ 2.05 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 
     1010;LOADCT, 1 
     1011;OPLOT, [ 2.05 ], [ -1 ], PSYM = 8, THICK = 3, SYMSIZE = 1.6, COLOR = 150 
    9491012 
    9501013; add model 
    … …  
    9591022; prepare plot 
    9601023LOADCT, 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 = ztitle 
     1024ztitle = "DSi total (mmolSi/m3)" 
     1025PLOT, [ 0., 30. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKFORMAT='(I2)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
    9631026; add obs 
    964 OPLOT, dsi_obs_mean, - depth_nut_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8 
    965 FOR i = 0, N_std - 1 DO BEGIN 
    966    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 ENDFOR 
     1027;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 
    9691032 
    9701033; plot seawater value 
    … …  
    9781041OPLOT, dsi_mod_mean+dsi_mod_std, -zz, THICK = 1, COLOR = 150 
    9791042OPLOT, dsi_mod_mean-dsi_mod_std, -zz, THICK = 1, COLOR = 150 
     1043 
     1044 
     1045;--- POC --- 
     1046 
     1047; prepare plot 
     1048LOADCT, 0 
     1049ztitle = "POC total (mmolC/m3)" 
     1050PLOT, [ 0., 3000. ], [ -1., 0. ], CHARSIZE = cs, /NODATA, SUBTITLE = '', XTICKINTERVAL=1500, XTICKFORMAT='(I4)', XMINOR = 2, YMINOR = 2, TITLE = ztitle 
     1051; add obs 
     1052OPLOT, POC_obs_mean, - depth_POC_mean, PSYM = 8, THICK = 3, SYMSIZE = 0.8 
     1053FOR 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) ] 
     1056ENDFOR 
     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 
     1063LOADCT, 9 
     1064OPLOT, POC_mod_mean, -zz, THICK = 3, COLOR = 150 
     1065LOADCT, 0 
     1066OPLOT, POC_mod_mean+POC_mod_std, -zz, THICK = 1, COLOR = 150 
     1067OPLOT, POC_mod_mean-POC_mod_std, -zz, THICK = 1, COLOR = 150 
     1068 
     1069; 
    9801070 
    9811071; 
Note: See TracChangeset for help on using the changeset viewer.