Changeset 354
- Timestamp:
- 2011-07-29T16:44:45+02:00 (13 years ago)
- Location:
- branches/ORCHIDEE_EXT/ORCHIDEE
- Files:
-
- 3 edited
-
src_parameters/constantes.f90 (modified) (1 diff)
-
src_sechiba/intersurf.f90 (modified) (10 diffs)
-
src_sechiba/routing.f90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters/constantes.f90
r343 r354 1815 1815 CALL getin_p('MAINT_RESP_COEFF',maint_resp_coeff) 1816 1816 !- 1817 ! soilcarbon parameters 1817 ! soilcarbon parameters 1818 !- 1819 ! 1820 !Config Key = FRAC_CARB_AA 1821 !Config Desc = frac carb coefficients from active pool: depends on clay content 1822 !Config if = OK_STOMATE 1823 !Config Def = 0.0 1824 !Config Help = fraction of the active pool going to the active pool 1825 !Config Units = NONE 1818 1826 CALL getin_p('FRAC_CARB_AA',frac_carb_aa) 1819 CALL getin_p('FRAC_CARB_AP',frac_carb_ap) 1827 ! 1828 !Config Key = FRAC_CARB_AP 1829 !Config Desc = frac carb coefficients from active pool: depends on clay content 1830 !Config if = OK_STOMATE 1831 !Config Def = 0.004 1832 !Config Help = fraction of the active pool going to the passive pool 1833 !Config Units = NONE 1834 CALL getin_p('FRAC_CARB_AP',frac_carb_ap) 1835 ! 1836 !Config Key = FRAC_CARB_SS 1837 !Config Desc = frac_carb_coefficients from slow pool 1838 !Config if = OK_STOMATE 1839 !Config Def = 0.0 1840 !Config Help = fraction of the slow pool going to the slow pool 1841 !Config Units = NONE 1820 1842 CALL getin_p('FRAC_CARB_SS',frac_carb_ss) 1843 ! 1844 !Config Key = FRAC_CARB_SA 1845 !Config Desc = frac_carb_coefficients from slow pool 1846 !Config if = OK_STOMATE 1847 !Config Def = 0.42 1848 !Config Help = fraction of the slow pool going to the active pool 1849 !Config Units = NONE 1821 1850 CALL getin_p('FRAC_CARB_SA',frac_carb_sa) 1851 ! 1852 !Config Key = FRAC_CARB_SP 1853 !Config Desc = frac_carb_coefficients from slow pool 1854 !Config if = OK_STOMATE 1855 !Config Def = 0.03 1856 !Config Help = fraction of the slow pool going to the passive pool 1857 !Config Units = NONE 1822 1858 CALL getin_p('FRAC_CARB_SP',frac_carb_sp) 1859 ! 1860 !Config Key = FRAC_CARB_PP 1861 !Config Desc = frac_carb_coefficients from passive pool 1862 !Config if = OK_STOMATE 1863 !Config Def = 0.0 1864 !Config Help = fraction of the passive pool going to the passive pool 1865 !Config Units = NONE 1823 1866 CALL getin_p('FRAC_CARB_PP',frac_carb_pp) 1867 ! 1868 !Config Key = FRAC_CARB_PA 1869 !Config Desc = frac_carb_coefficients from passive pool 1870 !Config if = OK_STOMATE 1871 !Config Def = 0.45 1872 !Config Help = fraction of the passive pool going to the passive pool 1873 !Config Units = NONE 1824 1874 CALL getin_p('FRAC_CARB_PA',frac_carb_pa) 1875 ! 1876 !Config Key = FRAC_CARB_PS 1877 !Config Desc = frac_carb_coefficients from passive pool 1878 !Config if = OK_STOMATE 1879 !Config Def = 0.0 1880 !Config Help = fraction of the passive pool going to the passive pool 1881 !Config Units = NONE 1825 1882 CALL getin_p('FRAC_CARB_PS',frac_carb_ps) 1826 1883 ! 1884 !Config Key = ACTIVE_TO_PASS_CLAY_FRAC 1885 !Config Desc = 1886 !Config if = OK_STOMATE 1887 !Config Def = .68 1888 !Config Help = 1889 !Config Units = NONE 1827 1890 CALL getin_p('ACTIVE_TO_PASS_CLAY_FRAC',active_to_pass_clay_frac) 1891 ! 1892 !Config Key = CARBON_TAU_IACTIVE 1893 !Config Desc = residence times in carbon pools 1894 !Config if = OK_STOMATE 1895 !Config Def = 0.149 1896 !Config Help = 1897 !Config Units = days (d) 1828 1898 CALL getin_p('CARBON_TAU_IACTIVE',carbon_tau_iactive) 1899 ! 1900 !Config Key = CARBON_TAU_ISLOW 1901 !Config Desc = residence times in carbon pools 1902 !Config if = OK_STOMATE 1903 !Config Def = 5.48 1904 !Config Help = 1905 !Config Units = days (d) 1829 1906 CALL getin_p('CARBON_TAU_ISLOW',carbon_tau_islow) 1907 ! 1908 !Config Key = CARBON_TAU_IPASSIVE 1909 !Config Desc = residence times in carbon pools 1910 !Config if = OK_STOMATE 1911 !Config Def = 241. 1912 !Config Help = 1913 !Config Units = days (d) 1830 1914 CALL getin_p('CARBON_TAU_IPASSIVE',carbon_tau_ipassive) 1915 ! 1916 !Config Key = FLUX_TOT_COEFF 1917 !Config Desc = 1918 !Config if = OK_STOMATE 1919 !Config Def = 1.2, 1.4,.75 1920 !Config Help = 1921 !Config Units = 1831 1922 CALL getin_p('FLUX_TOT_COEFF',flux_tot_coeff) 1832 1923 !- 1833 1924 ! turnover parameters 1925 !- 1834 1926 CALL getin_p('NEW_TURNOVER_TIME_REF',new_turnover_time_ref) 1835 1927 CALL getin_p('DT_TURNOVER_TIME',dt_turnover_time) -
branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/intersurf.f90
r335 r354 501 501 502 502 z0(i,j) = zz0(ik) 503 coastalflow(i,j) = zcoastal(ik)/ 1000.504 riverflow(i,j) = zriver(ik)/ 1000.503 coastalflow(i,j) = zcoastal(ik)/mille 504 riverflow(i,j) = zriver(ik)/mille 505 505 tsol_rad(i,j) = ztsol_rad(ik) 506 506 vevapp(i,j) = zvevapp(ik) … … 520 520 ! Modified fields for variables scattered during the writing 521 521 ! 522 dcoastal(:) = (zcoastal(:))/ 1000.523 driver(:) = (zriver(:))/ 1000.522 dcoastal(:) = (zcoastal(:))/mille 523 driver(:) = (zriver(:))/mille 524 524 ! 525 525 IF ( .NOT. l_first_intersurf) THEN … … 1009 1009 1010 1010 z0(kindex(ik)) = zz0(ik) 1011 coastalflow(kindex(ik)) = zcoastal(ik)/ 1000.1012 riverflow(kindex(ik)) = zriver(ik)/ 1000.1011 coastalflow(kindex(ik)) = zcoastal(ik)/mille 1012 riverflow(kindex(ik)) = zriver(ik)/mille 1013 1013 tsol_rad(kindex(ik)) = ztsol_rad(ik) 1014 1014 vevapp(kindex(ik)) = zvevapp(ik) … … 1028 1028 ! Modified fields for variables scattered during the writing 1029 1029 ! 1030 dcoastal(:) = (zcoastal(:))/ 1000.1031 driver(:) = (zriver(:))/ 1000.1030 dcoastal(:) = (zcoastal(:))/mille 1031 driver(:) = (zriver(:))/mille 1032 1032 ! 1033 1033 IF ( .NOT. l_first_intersurf) THEN … … 1717 1717 1718 1718 z0(ik) = zz0(ik) 1719 coastalflow(ik) = zcoastal(ik)/ 1000.1720 riverflow(ik) = zriver(ik)/ 1000.1719 coastalflow(ik) = zcoastal(ik)/mille 1720 riverflow(ik) = zriver(ik)/mille 1721 1721 tsol_rad(ik) = ztsol_rad(ik) 1722 1722 vevapp(ik) = zvevapp(ik) … … 1747 1747 ! Modified fields for variables scattered during the writing 1748 1748 ! 1749 dcoastal(:) = (zcoastal(:))/ 1000.1750 driver(:) = (zriver(:))/ 1000.1749 dcoastal(:) = (zcoastal(:))/mille 1750 driver(:) = (zriver(:))/mille 1751 1751 ! 1752 1752 IF ( .NOT. l_first_intersurf) THEN … … 2495 2495 2496 2496 z0(ik) = zz0(ik) 2497 coastalflow(ik) = zcoastal(ik)/ 1000.2498 riverflow(ik) = zriver(ik)/ 1000.2497 coastalflow(ik) = zcoastal(ik)/mille 2498 riverflow(ik) = zriver(ik)/mille 2499 2499 tsol_rad(ik) = ztsol_rad(ik) 2500 2500 vevapp(ik) = zvevapp(ik) … … 2527 2527 ! Modified fields for variables scattered during the writing 2528 2528 ! 2529 dcoastal(:) = (zcoastal(:))/ 1000.2530 driver(:) = (zriver(:))/ 1000.2529 dcoastal(:) = (zcoastal(:))/mille 2530 driver(:) = (zriver(:))/mille 2531 2531 ! 2532 2532 IF ( .NOT. l_first_intersurf) THEN … … 2762 2762 ! 2763 2763 !Config Key = NVM 2764 !Config Desc = number of PFTs 2764 !Config Desc = number of PFTs 2765 !Config if = ANYTIME 2766 !Config Def = 13 2767 !Config Help = The number of vegetation types define by the user 2768 !Config Units = NONE 2765 2769 CALL getin_p('NVM',nvm) 2766 2770 WRITE(numout,*)'the number of pfts is : ', nvm … … 2811 2815 ENDIF 2812 2816 2813 !!$ DS : reading of IMPOSE_PARAM 2814 ! Option : do you want to change the values of the parameters 2817 !Config Key = IMPOSE_PARAM 2818 !Config Desc = Do you impose the values of the parameters? 2819 !Config if = ANYTIME 2820 !Config Def = y 2821 !Config Help = This flag can deactivate the reading of some paramters. 2822 ! Useful if you want to use the standard values without commenting the run.def 2823 !Config Units = NONE 2815 2824 CALL getin_p('IMPOSE_PARAM',impose_param) 2825 2826 ! Read and allocate the pft parameters 2816 2827 CALL pft_parameters_main 2828 2817 2829 ! 2818 2830 !Config Key = RIVER_ROUTING -
branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/routing.f90
r311 r354 325 325 ! Convert from kg/dtradia to m^3/dtradia 326 326 ! 327 hydrographs(:) = hydrographs(:)/ 1000.327 hydrographs(:) = hydrographs(:)/mille 328 328 ! 329 329 ENDIF … … 934 934 IF ( route_tobasin(ig,ib) .GT. 0 ) THEN 935 935 ! 936 flow = MIN(fast_reservoir(ig,ib)/((topo_resid(ig,ib)/ 1000.)*fast_tcst*one_day/dt_routing),&936 flow = MIN(fast_reservoir(ig,ib)/((topo_resid(ig,ib)/mille)*fast_tcst*one_day/dt_routing),& 937 937 & fast_reservoir(ig,ib)) 938 938 fast_flow(ig,ib) = flow 939 939 ! 940 flow = MIN(slow_reservoir(ig,ib)/((topo_resid(ig,ib)/ 1000.)*slow_tcst*one_day/dt_routing),&940 flow = MIN(slow_reservoir(ig,ib)/((topo_resid(ig,ib)/mille)*slow_tcst*one_day/dt_routing),& 941 941 & slow_reservoir(ig,ib)) 942 942 slow_flow(ig,ib) = flow 943 943 ! 944 flow = MIN(stream_reservoir(ig,ib)/((topo_resid(ig,ib)/ 1000.)*stream_tcst*one_day/dt_routing),&944 flow = MIN(stream_reservoir(ig,ib)/((topo_resid(ig,ib)/mille)*stream_tcst*one_day/dt_routing),& 945 945 & stream_reservoir(ig,ib)) 946 946 stream_flow(ig,ib) = flow … … 1892 1892 WRITE(numout,*) '-------------> hierarchy ',hierarchy_bx(1,1) 1893 1893 DO jp=1,nbj 1894 WRITE(numout,fmt) INT(hierarchy_bx(1:nbi,jp)/ 1000.)1894 WRITE(numout,fmt) INT(hierarchy_bx(1:nbi,jp)/mille) 1895 1895 ENDDO 1896 1896 WRITE(numout,*) '-------------> topoindex ',topoind_bx(1,1) 1897 1897 DO jp=1,nbj 1898 WRITE(numout,fmt) INT(topoind_bx(1:nbi,jp)/ 1000.)1898 WRITE(numout,fmt) INT(topoind_bx(1:nbi,jp)/mille) 1899 1899 ENDDO 1900 1900 !
Note: See TracChangeset
for help on using the changeset viewer.
