Changeset 326
- Timestamp:
- 2011-07-18T11:42:45+02:00 (13 years ago)
- Location:
- branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/intersurf.f90
r311 r326 2902 2902 control_flags%ok_dgvm = .FALSE. 2903 2903 CALL getin_p('STOMATE_OK_DGVM',control_flags%ok_dgvm) 2904 2905 IF ( control_flags%ok_dgvm ) THEN2906 WRITE(numout,*) 'You try to use LPJ ',control_flags%ok_dgvm, &2907 ' with this version. '2908 WRITE(numout,*) 'It is not possible because it has to be modified ', &2909 ' to give correct values.'2910 CALL ipslerr (2,'intsurf_config', &2911 & 'Use of STOMATE_OK_DGVM is not stable for this version.',&2912 & 'ORCHIDEE should not give correct results with this option activated.', &2913 & 'Please disable DGVM to use this version of ORCHIDEE.')2914 ENDIF2915 !2916 !!$ DS : reading of parameters associated to ok_dgvm2917 IF ( control_flags%ok_dgvm ) THEN2918 CALL getin_dgvm_parameters2919 ENDIF2920 2904 2921 2905 ! … … 4586 4570 & dt, hist_dt_stom, hori_id, hist_PFTaxis_id, & 4587 4571 & hist_pool_10axis_id, hist_pool_100axis_id, & 4588 & hist_pool_11axis_id, hist_pool_101axis_id, control_flags) 4589 ! deforestation axis added as arguments 4572 & hist_pool_11axis_id, hist_pool_101axis_id) 4590 4573 4591 4574 !- end definition … … 4683 4666 & hist_dt, hist_hori_id, hist_PFTaxis_id, & 4684 4667 & hist_pool_10axis_id, hist_pool_100axis_id, & 4685 & hist_pool_11axis_id, hist_pool_101axis_id , control_flags)4668 & hist_pool_11axis_id, hist_pool_101axis_id) 4686 4669 ! deforestation axis added as arguments 4687 4670 … … 4711 4694 INTEGER(i_std),INTENT(in) :: hist_pool_10axis_id,hist_pool_100axis_id 4712 4695 INTEGER(i_std),INTENT(in) :: hist_pool_11axis_id,hist_pool_101axis_id 4713 !- Flags that (de)activate parts of the model4714 TYPE(control_type), INTENT(in) :: control_flags4715 4696 !- 4716 4697 !- 1 local … … 4941 4922 & nvm,1,nvm, hist_PFTaxis_id,32, ave(3), dt, hist_dt) 4942 4923 4943 IF ( control_flags%ok_dgvm ) THEN 4944 ! Adaptation to climate 4945 CALL histdef (hist_id_stom, & 4946 & TRIM("ADAPTATION "), & 4947 & TRIM("Adaptation to climate (DGVM) "), & 4948 & TRIM("- "), iim,jjm, hist_hori_id, & 4949 & nvm,1,nvm, hist_PFTaxis_id,32, ave(3), dt, hist_dt) 4950 4951 ! Probability from regenerative 4952 CALL histdef (hist_id_stom, & 4953 & TRIM("REGENERATION "), & 4954 & TRIM("Probability from regenerative (DGVM) "), & 4955 & TRIM("- "), iim,jjm, hist_hori_id, & 4956 & nvm,1,nvm, hist_PFTaxis_id,32, ave(3), dt, hist_dt) 4957 ENDIF 4924 ! Adaptation to climate 4925 CALL histdef (hist_id_stom, & 4926 & TRIM("ADAPTATION "), & 4927 & TRIM("Adaptation to climate (DGVM) "), & 4928 & TRIM("- "), iim,jjm, hist_hori_id, & 4929 & nvm,1,nvm, hist_PFTaxis_id,32, ave(3), dt, hist_dt) 4930 4931 ! Probability from regenerative 4932 CALL histdef (hist_id_stom, & 4933 & TRIM("REGENERATION "), & 4934 & TRIM("Probability from regenerative (DGVM) "), & 4935 & TRIM("- "), iim,jjm, hist_hori_id, & 4936 & nvm,1,nvm, hist_PFTaxis_id,32, ave(3), dt, hist_dt) 4958 4937 4959 4938 ! crown area of individuals (m**2) -
branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/slowproc.f90
r311 r326 703 703 ! to be in sechiba when teststomate will have disapeared. 704 704 !MM Problem here with dpu which depends on soil type 705 DO jv= 1, nbdl-1705 DO l = 1, nbdl-1 706 706 ! first 2.0 is dpu 707 707 ! second 2.0 is average 708 diaglev( jv) = dpu_cste/(2**(nbdl-1) -1) * ( ( 2**(jv-1) -1) + ( 2**(jv) -1) ) / 2.0708 diaglev(l) = dpu_cste/(2**(nbdl-1) -1) * ( ( 2**(l-1) -1) + ( 2**(l) -1) ) / 2.0 709 709 ENDDO 710 710 diaglev(nbdl) = dpu_cste … … 2579 2579 ! conservation : 2580 2580 rapport = ( sum_veg - sumvAnthro ) / ( sum_veg - sumvAnthro_old ) 2581 veget_next(ib,1) = veget_last(ib,1) * rapport 2582 DO jv = 2, nvm 2583 IF ( .NOT. natural(jv) ) THEN 2581 DO jv = 1, nvm 2582 IF ( natural(jv) ) THEN 2584 2583 veget_next(ib,jv) = veget_last(ib,jv) * rapport 2585 2584 ENDIF 2586 2585 ENDDO 2587 2586 ! test 2588 IF ( ABS( SUM(veget_next(ib,:)) - sum_veg ) > EPSILON(un) ) THEN2587 IF ( ABS( SUM(veget_next(ib,:)) - sum_veg ) > 10*EPSILON(un) ) THEN 2589 2588 WRITE(numout,*) "No conservation of sum of veget for point ",ib,",(",lalo(ib,1),",",lalo(ib,2),")" 2590 2589 WRITE(numout,*) "last sum of veget ",sum_veg," new sum of veget ",SUM(veget_next(ib,:))," error : ",& 2591 2590 & SUM(veget_next(ib,:)) - sum_veg 2592 WRITE(numout,*) "Anthropic modifica ztions : last ",sumvAnthro_old," new ",sumvAnthro2591 WRITE(numout,*) "Anthropic modifications : last ",sumvAnthro_old," new ",sumvAnthro 2593 2592 CALL ipslerr (3,'slowproc_update', & 2594 2593 & 'No conservation of sum of veget_next', &
Note: See TracChangeset
for help on using the changeset viewer.