Changeset 326


Ignore:
Timestamp:
2011-07-18T11:42:45+02:00 (13 years ago)
Author:
didier.solyga
Message:

Synchronize intersurf.f90 and slowproc.f90 with the revisions 314 and 317 of the trunk

Location:
branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/intersurf.f90

    r311 r326  
    29022902    control_flags%ok_dgvm = .FALSE. 
    29032903    CALL getin_p('STOMATE_OK_DGVM',control_flags%ok_dgvm) 
    2904  
    2905     IF ( control_flags%ok_dgvm ) THEN 
    2906        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     ENDIF 
    2915     ! 
    2916 !!$    DS : reading of parameters associated to ok_dgvm 
    2917     IF ( control_flags%ok_dgvm ) THEN 
    2918        CALL getin_dgvm_parameters 
    2919     ENDIF     
    29202904 
    29212905    ! 
     
    45864570            & dt, hist_dt_stom, hori_id, hist_PFTaxis_id, & 
    45874571            & 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) 
    45904573 
    45914574       !- end definition 
     
    46834666       &  hist_dt, hist_hori_id, hist_PFTaxis_id, & 
    46844667       & 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) 
    46864669    ! deforestation axis added as arguments 
    46874670 
     
    47114694    INTEGER(i_std),INTENT(in) :: hist_pool_10axis_id,hist_pool_100axis_id 
    47124695    INTEGER(i_std),INTENT(in) :: hist_pool_11axis_id,hist_pool_101axis_id 
    4713     !- Flags that (de)activate parts of the model 
    4714     TYPE(control_type), INTENT(in)             :: control_flags 
    47154696    !- 
    47164697    !- 1 local 
     
    49414922         &               nvm,1,nvm, hist_PFTaxis_id,32, ave(3), dt, hist_dt) 
    49424923 
    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) 
    49584937 
    49594938    ! crown area of individuals (m**2) 
  • branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/slowproc.f90

    r311 r326  
    703703       ! to be in sechiba when teststomate will have disapeared. 
    704704!MM Problem here with dpu which depends on soil type 
    705     DO jv = 1, nbdl-1 
     705    DO l = 1, nbdl-1 
    706706       ! first 2.0 is dpu  
    707707       ! second 2.0 is average 
    708        diaglev(jv) = dpu_cste/(2**(nbdl-1) -1) * ( ( 2**(jv-1) -1) + ( 2**(jv) -1) ) / 2.0 
     708       diaglev(l) = dpu_cste/(2**(nbdl-1) -1) * ( ( 2**(l-1) -1) + ( 2**(l) -1) ) / 2.0 
    709709    ENDDO 
    710710    diaglev(nbdl) = dpu_cste 
     
    25792579             ! conservation : 
    25802580             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 
    25842583                   veget_next(ib,jv) = veget_last(ib,jv) * rapport 
    25852584                ENDIF 
    25862585             ENDDO 
    25872586             ! test 
    2588              IF ( ABS( SUM(veget_next(ib,:)) - sum_veg ) > EPSILON(un) ) THEN 
     2587             IF ( ABS( SUM(veget_next(ib,:)) - sum_veg ) > 10*EPSILON(un) ) THEN 
    25892588                WRITE(numout,*) "No conservation of sum of veget for point ",ib,",(",lalo(ib,1),",",lalo(ib,2),")"  
    25902589                WRITE(numout,*) "last sum of veget ",sum_veg," new sum of veget ",SUM(veget_next(ib,:))," error : ",& 
    25912590                     &                         SUM(veget_next(ib,:)) - sum_veg 
    2592                 WRITE(numout,*) "Anthropic modificaztions : last ",sumvAnthro_old," new ",sumvAnthro                 
     2591                WRITE(numout,*) "Anthropic modifications : last ",sumvAnthro_old," new ",sumvAnthro   
    25932592                CALL ipslerr (3,'slowproc_update', & 
    25942593                     &          'No conservation of sum of veget_next', & 
Note: See TracChangeset for help on using the changeset viewer.