Changeset 326 for branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/slowproc.f90
- Timestamp:
- 2011-07-18T11:42:45+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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.