Changeset 257 for branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/intersurf.f90
- Timestamp:
- 2011-06-17T14:02:17+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/intersurf.f90
r116 r257 7 7 !! 8 8 !! @call sechiba_main 9 !! @Version : $Revision: 1.85 $, $Date: 2010/07/29 15:58:19$9 !! @Version : $Revision: 221 $, $Date: 2011-05-16 17:26:17 +0200 (Mon, 16 May 2011) $ 10 10 !! 11 11 !! @author Marie-Alice Foujols and Jan Polcher 12 12 !! 13 !! $Header: /home/ssipsl/CVSREP/ORCHIDEE/src_sechiba/intersurf.f90,v 1.85 2010/07/29 15:58:19 ssipsl Exp $ 13 !< $HeadURL: http://forge.ipsl.jussieu.fr/orchidee/svn/trunk/ORCHIDEE/src_sechiba/intersurf.f90 $ 14 !< $Date: 2011-05-16 17:26:17 +0200 (Mon, 16 May 2011) $ 15 !< $Author: martial.mancip $ 16 !< $Revision: 221 $ 14 17 !! IPSL (2006) 15 18 !! This software is governed by the CeCILL licence see ORCHIDEE/ORCHIDEE_CeCILL.LIC … … 181 184 ! 182 185 CALL ipslnlf(new_number=numout,old_number=old_fileout) 183 184 186 ! 185 187 IF (l_first_intersurf) THEN … … 225 227 IF ( ok_watchout ) THEN 226 228 IF (is_root_prc) THEN 227 zlev_mean = 0.229 zlev_mean = zero 228 230 DO ik=1, nbp_glo 229 231 j = ((index_g(ik)-1)/iim_g) + 1 … … 391 393 !!$ dt_split_watch,dt_watch,one_day 392 394 !!$ CALL solarang (julian_watch, julian0, iim, jjm, lon, lat, sinang) 393 !!$ WHERE ( sinang(:,:) .LT. EPSILON( 1.) )395 !!$ WHERE ( sinang(:,:) .LT. EPSILON(un) ) 394 396 !!$ isinang(:,:) = isinang(:,:) - 1 395 397 !!$ ENDWHERE … … 529 531 CALL histwrite (hist_id, 'riverflow',itau_sechiba, driver, kjpindex, kindex) 530 532 ! 531 CALL histwrite (hist_id, 'temp_sol', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)532 CALL histwrite (hist_id, 'tsol_max', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)533 CALL histwrite (hist_id, 'tsol_min', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)534 CALL histwrite (hist_id, 'fluxsens', itau_sechiba, fluxsens, iim*jjm, kindex)535 CALL histwrite (hist_id, 'fluxlat', itau_sechiba, fluxlat, iim*jjm, kindex)536 CALL histwrite (hist_id, 'swnet', itau_sechiba, dswnet, iim*jjm, kindex)537 CALL histwrite (hist_id, 'swdown', itau_sechiba, dswdown, iim*jjm, kindex)538 CALL histwrite (hist_id, 'alb_vis', itau_sechiba, albedo(:,:,1), iim*jjm, kindex)539 CALL histwrite (hist_id, 'alb_nir', itau_sechiba, albedo(:,:,2), iim*jjm, kindex)540 CALL histwrite (hist_id, 'tair', itau_sechiba, temp_air, iim*jjm, kindex)541 CALL histwrite (hist_id, 'qair', itau_sechiba, qair, iim*jjm, kindex)533 CALL histwrite (hist_id, 'temp_sol', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 534 CALL histwrite (hist_id, 'tsol_max', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 535 CALL histwrite (hist_id, 'tsol_min', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 536 CALL histwrite (hist_id, 'fluxsens', itau_sechiba, fluxsens, kjpindex, kindex) 537 CALL histwrite (hist_id, 'fluxlat', itau_sechiba, fluxlat, kjpindex, kindex) 538 CALL histwrite (hist_id, 'swnet', itau_sechiba, dswnet, kjpindex, kindex) 539 CALL histwrite (hist_id, 'swdown', itau_sechiba, dswdown, kjpindex, kindex) 540 CALL histwrite (hist_id, 'alb_vis', itau_sechiba, albedo(:,:,1), kjpindex, kindex) 541 CALL histwrite (hist_id, 'alb_nir', itau_sechiba, albedo(:,:,2), kjpindex, kindex) 542 CALL histwrite (hist_id, 'tair', itau_sechiba, temp_air, kjpindex, kindex) 543 CALL histwrite (hist_id, 'qair', itau_sechiba, qair, kjpindex, kindex) 542 544 ! Ajout Nathalie - Juin 2006 - on conserve q2m/t2m 543 CALL histwrite (hist_id, 'q2m', itau_sechiba, qair, iim*jjm, kindex)544 CALL histwrite (hist_id, 't2m', itau_sechiba, temp_air, iim*jjm, kindex)545 CALL histwrite (hist_id, 'q2m', itau_sechiba, qair, kjpindex, kindex) 546 CALL histwrite (hist_id, 't2m', itau_sechiba, temp_air, kjpindex, kindex) 545 547 IF ( hist2_id > 0 ) THEN 546 548 CALL histwrite (hist2_id, 'evap', itau_sechiba, zvevapp, kjpindex, kindex) … … 548 550 CALL histwrite (hist2_id, 'riverflow',itau_sechiba, driver, kjpindex, kindex) 549 551 ! 550 CALL histwrite (hist2_id, 'temp_sol', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)551 CALL histwrite (hist2_id, 'tsol_max', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)552 CALL histwrite (hist2_id, 'tsol_min', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)553 CALL histwrite (hist2_id, 'fluxsens', itau_sechiba, fluxsens, iim*jjm, kindex)554 CALL histwrite (hist2_id, 'fluxlat', itau_sechiba, fluxlat, iim*jjm, kindex)555 CALL histwrite (hist2_id, 'swnet', itau_sechiba, dswnet, iim*jjm, kindex)556 CALL histwrite (hist2_id, 'swdown', itau_sechiba, dswdown, iim*jjm, kindex)557 CALL histwrite (hist2_id, 'alb_vis', itau_sechiba, albedo(:,:,1), iim*jjm, kindex)558 CALL histwrite (hist2_id, 'alb_nir', itau_sechiba, albedo(:,:,2), iim*jjm, kindex)559 CALL histwrite (hist2_id, 'tair', itau_sechiba, temp_air, iim*jjm, kindex)560 CALL histwrite (hist2_id, 'qair', itau_sechiba, qair, iim*jjm, kindex)561 CALL histwrite (hist2_id, 'q2m', itau_sechiba, qair, iim*jjm, kindex)562 CALL histwrite (hist2_id, 't2m', itau_sechiba, temp_air, iim*jjm, kindex)552 CALL histwrite (hist2_id, 'temp_sol', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 553 CALL histwrite (hist2_id, 'tsol_max', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 554 CALL histwrite (hist2_id, 'tsol_min', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 555 CALL histwrite (hist2_id, 'fluxsens', itau_sechiba, fluxsens, kjpindex, kindex) 556 CALL histwrite (hist2_id, 'fluxlat', itau_sechiba, fluxlat, kjpindex, kindex) 557 CALL histwrite (hist2_id, 'swnet', itau_sechiba, dswnet, kjpindex, kindex) 558 CALL histwrite (hist2_id, 'swdown', itau_sechiba, dswdown, kjpindex, kindex) 559 CALL histwrite (hist2_id, 'alb_vis', itau_sechiba, albedo(:,:,1), kjpindex, kindex) 560 CALL histwrite (hist2_id, 'alb_nir', itau_sechiba, albedo(:,:,2), kjpindex, kindex) 561 CALL histwrite (hist2_id, 'tair', itau_sechiba, temp_air, kjpindex, kindex) 562 CALL histwrite (hist2_id, 'qair', itau_sechiba, qair, kjpindex, kindex) 563 CALL histwrite (hist2_id, 'q2m', itau_sechiba, qair, kjpindex, kindex) 564 CALL histwrite (hist2_id, 't2m', itau_sechiba, temp_air, kjpindex, kindex) 563 565 ENDIF 564 566 ELSE 565 567 CALL histwrite (hist_id, 'Evap', itau_sechiba, zvevapp, kjpindex, kindex) 566 CALL histwrite (hist_id, 'SWnet', itau_sechiba, dswnet, iim*jjm, kindex)567 CALL histwrite (hist_id, 'Qh', itau_sechiba, fluxsens, iim*jjm, kindex)568 CALL histwrite (hist_id, 'Qle', itau_sechiba, fluxlat, iim*jjm, kindex)569 CALL histwrite (hist_id, 'AvgSurfT', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)570 CALL histwrite (hist_id, 'RadT', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)568 CALL histwrite (hist_id, 'SWnet', itau_sechiba, dswnet, kjpindex, kindex) 569 CALL histwrite (hist_id, 'Qh', itau_sechiba, fluxsens, kjpindex, kindex) 570 CALL histwrite (hist_id, 'Qle', itau_sechiba, fluxlat, kjpindex, kindex) 571 CALL histwrite (hist_id, 'AvgSurfT', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 572 CALL histwrite (hist_id, 'RadT', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 571 573 IF ( hist2_id > 0 ) THEN 572 574 CALL histwrite (hist2_id, 'Evap', itau_sechiba, zvevapp, kjpindex, kindex) 573 CALL histwrite (hist2_id, 'SWnet', itau_sechiba, dswnet, iim*jjm, kindex)574 CALL histwrite (hist2_id, 'Qh', itau_sechiba, fluxsens, iim*jjm, kindex)575 CALL histwrite (hist2_id, 'Qle', itau_sechiba, fluxlat, iim*jjm, kindex)576 CALL histwrite (hist2_id, 'AvgSurfT', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)577 CALL histwrite (hist2_id, 'RadT', itau_sechiba, temp_sol_NEW, iim*jjm, kindex)575 CALL histwrite (hist2_id, 'SWnet', itau_sechiba, dswnet, kjpindex, kindex) 576 CALL histwrite (hist2_id, 'Qh', itau_sechiba, fluxsens, kjpindex, kindex) 577 CALL histwrite (hist2_id, 'Qle', itau_sechiba, fluxlat, kjpindex, kindex) 578 CALL histwrite (hist2_id, 'AvgSurfT', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 579 CALL histwrite (hist2_id, 'RadT', itau_sechiba, temp_sol_NEW, kjpindex, kindex) 578 580 ENDIF 579 581 ENDIF … … 780 782 ! 781 783 IF ( ok_watchout ) THEN 782 zlev_mean = 0.784 zlev_mean = zero 783 785 DO ik=1, kjpindex 784 786 … … 905 907 !!$ julian_watch = date0_shifted+((itau_sechiba-0.5)/dt_split_watch)*dt_watch/one_day 906 908 !!$ CALL solarang (julian_watch, julian0, iim, jjm, lon, lat, sinang) 907 !!$ WHERE ( sinang(:,:) .LT. EPSILON( 1.) )909 !!$ WHERE ( sinang(:,:) .LT. EPSILON(un) ) 908 910 !!$ isinang(:,:) = isinang(:,:) - 1 909 911 !!$ ENDWHERE … … 1448 1450 IF ( ok_watchout ) THEN 1449 1451 IF (is_root_prc) THEN 1450 zlev_mean = 0.1452 zlev_mean = zero 1451 1453 DO ik=1, nbp_glo 1452 1454 j = ((index_g(ik)-1)/iim_g) + 1 … … 1602 1604 !!$ julian_watch = date0_shifted+((itau_sechiba-0.5)/dt_split_watch)*dt_watch/one_day 1603 1605 !!$ CALL solarang (julian_watch, julian0, iim, jjm, tmp_lon, tmp_lat, sinang) 1604 !!$ WHERE ( sinang(:,:) .LT. EPSILON( 1.) )1606 !!$ WHERE ( sinang(:,:) .LT. EPSILON(un) ) 1605 1607 !!$ isinang(:,:) = isinang(:,:) - 1 1606 1608 !!$ ENDWHERE … … 2178 2180 IF ( ok_watchout ) THEN 2179 2181 IF (is_root_prc) THEN 2180 zlev_mean = 0.2182 zlev_mean = zero 2181 2183 DO ik=1, nbp_glo 2182 2184 j = ((index_g(ik)-1)/iim_g) + 1 … … 2332 2334 !!$ julian_watch = date0_shifted+((itau_sechiba-0.5)/dt_split_watch)*dt_watch/one_day 2333 2335 !!$ CALL solarang (julian_watch, julian0, iim, jjm, tmp_lon, tmp_lat, sinang) 2334 !!$ WHERE ( sinang(:,:) .LT. EPSILON( 1.) )2336 !!$ WHERE ( sinang(:,:) .LT. EPSILON(un) ) 2335 2337 !!$ isinang(:,:) = isinang(:,:) - 1 2336 2338 !!$ ENDWHERE … … 2586 2588 CALL tlen2itau('1Y',dt,date0,year_length) 2587 2589 IF ( TRIM(calendar_str) .EQ. 'gregorian' ) THEN 2588 year_spread= 1.02590 year_spread=un 2589 2591 ELSE 2590 2592 year_spread = one_year/365.2425 … … 2610 2612 ! Real date 2611 2613 CALL ju2ymds (in_julian, year, month, day, sec) 2612 !!$ jur= 0.2614 !!$ jur=zero 2613 2615 !!$ julian_diff = in_julian 2614 2616 !!$ month_len = ioget_mon_len (year,month) … … 2630 2632 ENDIF 2631 2633 ELSE 2632 !!$ in_julian = itau2date(istp-1, 0., dt)2634 !!$ in_julian = itau2date(istp-1, zero, dt) 2633 2635 !!$ CALL ju2ymds (in_julian, year, month, day, sec) 2634 !!$ jur= 0.2636 !!$ jur=zero 2635 2637 !!$ julian_diff = in_julian 2636 2638 !!$ month_len = ioget_mon_len (year,month) … … 2693 2695 CALL getin_p('NVM',nvm) 2694 2696 WRITE(numout,*)'the number of pfts is : ', nvm 2695 !!$DS Debug 28/01/20112696 2697 ! 2697 2698 !Config Key = LONGPRINT … … 2723 2724 ! 2724 2725 dt_watch = dt 2725 CALL getin ('DT_WATCHOUT',dt_watch)2726 CALL getin_p('DT_WATCHOUT',dt_watch) 2726 2727 dt_split_watch = dt_watch / dt 2727 2728 ! … … 2740 2741 ENDIF 2741 2742 2742 2743 2743 !!$ DS : reading of IMPOSE_PARAM 2744 2744 ! Option : do you want to change the values of the parameters 2745 2745 CALL getin_p('IMPOSE_PARAM',impose_param) 2746 ! Calling pft_parameters2747 2746 CALL pft_parameters_main 2748 2747 ! … … 2784 2783 IF ( control_flags%hydrol_cwrr ) THEN 2785 2784 CALL getin_hydrol_cwrr_parameters 2785 ELSE 2786 CALL getin_hydrolc_parameters 2787 ! we read the parameters for the choisnel hydrology 2786 2788 ENDIF 2787 2789 … … 2800 2802 CALL getin_co2_parameters 2801 2803 ENDIF 2802 2803 2804 2805 !!$ DS : reading of IMPOSE_PARAM2806 !!$ ! Option : do you want to change the values of the parameters2807 !!$ CALL getin_p('IMPOS_PARAM',impos_param)2808 !!$ ! Calling pft_parameters2809 !!$ CALL pft_main2810 2804 2811 2805 ! … … 2844 2838 WRITE(numout,*) 'It is not possible because it has to be modified ', & 2845 2839 ' to give correct values.' 2846 CALL ipslerr ( 3,'intsurf_config', &2847 & 'Use of STOMATE_OK_DGVM not allowed withthis version.',&2848 & 'ORCHIDEE will stop.', &2840 CALL ipslerr (2,'intsurf_config', & 2841 & 'Use of STOMATE_OK_DGVM is not stable for this version.',& 2842 & 'ORCHIDEE should not give correct results with this option activated.', & 2849 2843 & 'Please disable DGVM to use this version of ORCHIDEE.') 2850 2844 ENDIF … … 2965 2959 CALL getin_p('SECHIBA_reset_time', overwrite_time) 2966 2960 ! 2967 lev(:) = 0.2961 lev(:) = zero 2968 2962 itau_dep = istp 2969 2963 in_julian = itau2date(istp, date0, dt) … … 3186 3180 !Config Key = WRITE_STEP 3187 3181 !Config Desc = Frequency in seconds at which to WRITE output 3188 !Config Def = 86400.03182 !Config Def = one_day 3189 3183 !Config Help = This variables gives the frequency the output of 3190 3184 !Config the model should be written into the netCDF file. … … 3198 3192 ! 3199 3193 veg(1:nvm) = (/ (REAL(i,r_std),i=1,nvm) /) 3200 !$$ DS DEBUG3201 WRITE(numout,*)'nvm : = ', nvm3202 WRITE(numout,*)'veg : =', veg3203 !$$ nvm =13 (put the calling to getin before)3204 3194 sol(1:ngrnd) = (/ (REAL(i,r_std),i=1,ngrnd) /) 3205 3195 soltyp(1:nstm) = (/ (REAL(i,r_std),i=1,nstm) /) … … 3216 3206 WRITE(flux_sc,'("ave(X*",F8.1,")")') one_day/dt 3217 3207 !WRITE(flux_sc,'("(ave(X)*",F8.1,")")') one_day/dt 3218 WRITE(flux_insec,'("ave(X*",F8.6,")")') 1.0/dt3219 WRITE(flux_scinsec,'("ave(scatter(X*",F8.6,"))")') 1.0/dt3208 WRITE(flux_insec,'("ave(X*",F8.6,")")') un/dt 3209 WRITE(flux_scinsec,'("ave(scatter(X*",F8.6,"))")') un/dt 3220 3210 WRITE(numout,*) flux_op, one_day/dt, dt, dw 3221 3211 !- … … 3371 3361 & iim,jjm, hori_id, 1,1,1, -99, 32, once(1), dt,dw) 3372 3362 ENDIF 3373 IF ( control_flags%ok_stomate .OR. control_flags%stomate_watchout ) THEN3374 CALL histdef (hist_id,'CO2FLUX','Total output CO2 flux', 'gC/day/(m^2 tot)', &3375 & iim,jjm, hori_id, nvm, 1, nvm, vegax_id, 32, avescatter(1), dt, dw)3376 ENDIF3377 3363 !- 3378 3364 !- SECHIBA_HISTLEVEL = 2 … … 3692 3678 CALL histdef(hist_id, 'nobiofrac', 'Fraction of other surface types', '1', & 3693 3679 & iim,jjm, hori_id, nnobio, 1, nnobio, nobioax_id, 32, avescatter(3), dt,dw) 3694 IF ( control_flags%ok_stomate .OR. control_flags%stomate_watchout ) THEN3695 ! Total output CO2 flux3696 CALL histdef (hist_id,'CO2FLUX','Total output CO2 flux', 'gC/day/(m^2 tot)', &3697 & iim,jjm, hori_id, nvm, 1, nvm, vegax_id, 32, avescatter(1), dt, dw)3698 ENDIF3699 3680 !- 3700 3681 !- General energy balance … … 4033 4014 CALL histdef(hist2_id, 'emis', 'Surface emissivity', '?', & 4034 4015 & iim,jjm, hori_id2, 1,1,1, -99, 32, fluxop2(2), dt, dw2) 4035 IF ( control_flags%ok_stomate .OR. control_flags%stomate_watchout ) THEN4036 CALL histdef (hist2_id,'CO2FLUX','Total output CO2 flux', 'gC/day/(m^2 tot)', &4037 & iim,jjm, hori_id2, nvm, 1, nvm, vegax_id2, 32, avescatter2(2), dt, dw2)4038 ENDIF4039 4016 !- 4040 4017 !- SECHIBA_HISTLEVEL2 = 3 … … 4298 4275 CALL histdef(hist2_id, 'nobiofrac', 'Fraction of other surface types', '1', & 4299 4276 & iim,jjm, hori_id2, nnobio, 1, nnobio, nobioax_id2, 32, avescatter2(3), dt, dw2) 4300 IF ( control_flags%ok_stomate .OR. control_flags%stomate_watchout ) THEN4301 CALL histdef (hist2_id,'CO2FLUX','Total output CO2 flux', 'gC/day/(m^2 tot)', &4302 & iim,jjm, hori_id2, nvm, 1, nvm, vegax_id2, 32, avescatter2(1), dt, dw2)4303 ENDIF4304 4277 !- 4305 4278 !- General energy balance … … 4465 4438 hist_days_stom = 10. 4466 4439 CALL getin_p('STOMATE_HIST_DT', hist_days_stom) 4467 IF ( hist_days_stom == -1.) THEN4468 hist_dt_stom = -1.4440 IF ( hist_days_stom == moins_un ) THEN 4441 hist_dt_stom = moins_un 4469 4442 WRITE(numout,*) 'output frequency for STOMATE history file (d): one month.' 4470 4443 ELSE … … 4477 4450 dt_slow_ = one_day 4478 4451 CALL getin_p('DT_SLOW', dt_slow_) 4479 IF ( hist_days_stom /= -1.) THEN4452 IF ( hist_days_stom /= moins_un ) THEN 4480 4453 IF (dt_slow_ > hist_dt_stom) THEN 4481 4454 WRITE(numout,*) "DT_SLOW = ",dt_slow_," , STOMATE_HIST_DT = ",hist_dt_stom … … 4567 4540 !Config Help = Time step of the STOMATE IPCC history file 4568 4541 !- 4569 hist_days_stom_ipcc = 0.4542 hist_days_stom_ipcc = zero 4570 4543 CALL getin_p('STOMATE_IPCC_HIST_DT', hist_days_stom_ipcc) 4571 IF ( hist_days_stom_ipcc == -1.) THEN4572 hist_dt_stom_ipcc = -1.4544 IF ( hist_days_stom_ipcc == moins_un ) THEN 4545 hist_dt_stom_ipcc = moins_un 4573 4546 WRITE(numout,*) 'output frequency for STOMATE IPCC history file (d): one month.' 4574 4547 ELSE … … 4581 4554 dt_slow_ = one_day 4582 4555 CALL getin_p('DT_SLOW', dt_slow_) 4583 IF ( hist_days_stom_ipcc > 0.) THEN4556 IF ( hist_days_stom_ipcc > zero ) THEN 4584 4557 IF (dt_slow_ > hist_dt_stom_ipcc) THEN 4585 4558 WRITE(numout,*) "DT_SLOW = ",dt_slow_," , STOMATE_IPCC_HIST_DT = ",hist_dt_stom_ipcc … … 4822 4795 & 1,1,1, -99,32, ave(5), dt, hist_dt) 4823 4796 4824 ! MonthlyCO2 flux4825 CALL histdef (hist_id_stom, & 4826 & TRIM("CO2FLUX _MONTHLY"), &4827 & TRIM(" Monthly CO2 flux"), &4797 ! CO2 flux 4798 CALL histdef (hist_id_stom, & 4799 & TRIM("CO2FLUX "), & 4800 & TRIM("CO2 flux "), & 4828 4801 & TRIM("gC/m^2/pft/mth "), iim,jjm, hist_hori_id, & 4829 4802 & nvm,1,nvm, hist_PFTaxis_id,32, ave(1), dt, hist_dt) 4830 4803 4831 CALL histdef(hist_id_stom, &4832 & TRIM("CO2FLUX_MONTHLY_SUM "), &4833 & TRIM("Monthly CO2 flux"), &4834 & TRIM("PgC/m^2/mth "), 1,1, hist_hori_id, &4835 & 1,1,1, -99, 32, ave(1), dt, hist_dt)4804 !!$ CALL histdef(hist_id_stom, & 4805 !!$ & TRIM("CO2FLUX_MONTHLY_SUM "), & 4806 !!$ & TRIM("Monthly CO2 flux Sum "), & 4807 !!$ & TRIM("PgC/m^2/mth "), iim,jjm, hist_hori_id, & 4808 !!$ & 1,1,1, -99, 32, 'inst(scatter(X))', dt, hist_dt) 4836 4809 4837 4810 ! Output CO2 flux from fire … … 5121 5094 & TRIM("1/day "), iim,jjm, hist_hori_id, & 5122 5095 & nvm,1,nvm, hist_PFTaxis_id,32, ave(6), dt, hist_dt) 5096 5097 ! Establish tree 5098 CALL histdef (hist_id_stom, & 5099 & TRIM("ESTABTREE "), & 5100 & TRIM("Rate of tree establishement "), & 5101 & TRIM("1/day "), iim,jjm, hist_hori_id, & 5102 & 1,1,1, -99,32, ave(6), dt, hist_dt) 5103 5104 ! Establish grass 5105 CALL histdef (hist_id_stom, & 5106 & TRIM("ESTABGRASS "), & 5107 & TRIM("Rate of grass establishement "), & 5108 & TRIM("1/day "), iim,jjm, hist_hori_id, & 5109 & 1,1,1, -99,32, ave(6), dt, hist_dt) 5123 5110 5124 5111 ! Fraction of plants that dies (light competition)
Note: See TracChangeset
for help on using the changeset viewer.