equal
deleted
inserted
replaced
72 ;-------------------------------- |
72 ;-------------------------------- |
73 ; read data: time series, model |
73 ; read data: time series, model |
74 |
74 |
75 fm = addfile (dirm+film7,"r") |
75 fm = addfile (dirm+film7,"r") |
76 |
76 |
|
77 ; FMH: This was Jeff's, but I got an error |
77 data_mod = fm->COL_FIRE_CLOSS(18:25,:,:,:) |
78 data_mod = fm->COL_FIRE_CLOSS(18:25,:,:,:) |
|
79 ; FMH: So I'll try this |
|
80 ;data_mod = fm->COL_FIRE_CLOSS(17:24,:,:,:) |
78 |
81 |
79 delete (fm) |
82 delete (fm) |
80 |
83 |
81 ; Units for these variables are: |
84 ; Units for these variables are: |
82 ; g C/m^2/s |
85 ; g C/m^2/s |
190 delete (good) |
193 delete (good) |
191 |
194 |
192 ;---------------------------------------------------- |
195 ;---------------------------------------------------- |
193 ; compute M_global |
196 ; compute M_global |
194 |
197 |
195 score_max = 1. |
198 score_max = 2.5 |
196 |
199 |
197 Mscore1 = cc_space * cc_space * score_max |
200 Mscore1 = cc_space * cc_space * score_max |
198 |
201 |
199 M_global = sprintf("%.2f", Mscore1) |
202 M_global = sprintf("%.2f", Mscore1) |
200 |
203 |
540 ; data for table : per biome |
543 ; data for table : per biome |
541 |
544 |
542 ; unit change from TgC/month to PgC/month |
545 ; unit change from TgC/month to PgC/month |
543 unit_factor = 1.e-3 |
546 unit_factor = 1.e-3 |
544 |
547 |
545 score_max = 1. |
548 score_max = 2.5 |
546 |
549 |
547 tmp_ob = new((/ntime/),float) |
550 tmp_ob = new((/ntime/),float) |
548 tmp_mod = new((/ntime/),float) |
551 tmp_mod = new((/ntime/),float) |
549 |
552 |
550 total_ob = new((/n_biome/),float) |
553 total_ob = new((/n_biome/),float) |
605 "rm "+plot_name+"."+plot_type) |
608 "rm "+plot_name+"."+plot_type) |
606 |
609 |
607 ;------------------------------------------ |
610 ;------------------------------------------ |
608 ; data for table : global |
611 ; data for table : global |
609 |
612 |
610 score_max = 1. |
613 score_max = 2.5 |
611 |
614 |
612 tmp_ob = ndtooned(yvalues_g(:,0,:)) |
615 tmp_ob = ndtooned(yvalues_g(:,0,:)) |
613 tmp_mod = ndtooned(yvalues_g(:,1,:)) |
616 tmp_mod = ndtooned(yvalues_g(:,1,:)) |
614 |
617 |
615 cc_time = esccr(tmp_mod,tmp_ob,0) |
618 cc_time = esccr(tmp_mod,tmp_ob,0) |