Changes between Version 16 and Version 17 of DevelopmentActivities/ORCHIDEE-MICT-IMBALANCE-P/Evaluation


Ignore:
Timestamp:
2015-09-15T10:21:59+02:00 (9 years ago)
Author:
ajornet
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • DevelopmentActivities/ORCHIDEE-MICT-IMBALANCE-P/Evaluation

    v16 v17  
    208208CHANGE 
    209209{{{ 
    210       !! To get the extra energy used to melt the snowpack 
    211       IF (ok_explicitsnow .AND. temp_sol_new (ji) > tp_00 .AND. SUM(snowdz(ji,:)) .GT. zero .AND. snowcap(ji) .GT. zero) THEN 
    212  
    213          lwup_tmp(ji) = emis(ji) * c_stefan * temp_sol(ji)**4 + & 
    214            &     quatre * emis(ji) * c_stefan * temp_sol(ji)**3 * & 
    215            &     (tp_00 - temp_sol(ji)) 
     210!! To get the extra energy used to melt the snowpack 
     211IF (ok_explicitsnow .AND. temp_sol_new (ji) > tp_00 .AND. SUM(snowdz(ji,:)) .GT. zero .AND. snowcap(ji) .GT. zero) THEN 
     212 
     213lwup_tmp(ji) = emis(ji) * c_stefan * temp_sol(ji)**4 + & 
     214&     quatre * emis(ji) * c_stefan * temp_sol(ji)**3 * & 
     215&     (tp_00 - temp_sol(ji)) 
    216216}}} 
    217217 
    … …  
    219219To  
    220220{{{ 
    221   REAL(r_std), PARAMETER                                      :: EPS =1E-3_r_std         
    222    !! To get the extra energy used to melt the snowpack 
    223       IF (ok_explicitsnow .AND. temp_sol_new (ji) > tp_00 .AND. SUM(snowdz(ji,:)) .GT. zero .AND. snowcap(ji) .GT. EPS) THEN 
    224  
    225          lwup_tmp(ji) = emis(ji) * c_stefan * temp_sol(ji)**4 + & 
    226            &     quatre * emis(ji) * c_stefan * temp_sol(ji)**3 * & 
    227            &     (tp_00 - temp_sol(ji)) 
    228 }}} 
     221REAL(r_std), PARAMETER :: EPS =1E-3_r_std         
     222!! To get the extra energy used to melt the snowpack 
     223IF (ok_explicitsnow .AND. temp_sol_new (ji) > tp_00 .AND. SUM(snowdz(ji,:)) .GT. zero .AND. snowcap(ji) .GT. EPS) THEN 
     224 
     225lwup_tmp(ji) = emis(ji) * c_stefan * temp_sol(ji)**4 + & 
     226&     quatre * emis(ji) * c_stefan * temp_sol(ji)**3 * & 
     227&     (tp_00 - temp_sol(ji)) 
     228}}}