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

Timestamp:

2015-09-15T10:20:24+02:00 (9 years ago)

Author:

ajornet

Comment:

--

DevelopmentActivities/ORCHIDEE-MICT-IMBALANCE-P/Evaluation

@@ -200 +200 @@
  - GRAZING_MAP should not be read if grassland management is not activated.
  - If OK_EXPLICITSNOW is FALSE then snowdz is not initialized but used anyway (this is probably the case of other variables, to be checked).
+
+
+== Tao ==
+
+There is a small bug in our MICT version that might induce the failure of some point simulations. In subroutine enerbil_flux of enerbil.f90:
+
+CHANGE
+{{{
+      !! To get the extra energy used to melt the snowpack
+      IF (ok_explicitsnow .AND. temp_sol_new (ji) > tp_00 .AND. SUM(snowdz(ji,:)) .GT. zero .AND. snowcap(ji) .GT. zero) THEN
+
+         lwup_tmp(ji) = emis(ji) * c_stefan * temp_sol(ji)**4 + &
+           &     quatre * emis(ji) * c_stefan * temp_sol(ji)**3 * &
+           &     (tp_00 - temp_sol(ji))
+}}}
+
+
+To 
+{{{
+  REAL(r_std), PARAMETER                                       :: EPS =1E-3_r_std        
+   !! To get the extra energy used to melt the snowpack
+      IF (ok_explicitsnow .AND. temp_sol_new (ji) > tp_00 .AND. SUM(snowdz(ji,:)) .GT. zero .AND. snowcap(ji) .GT. EPS) THEN
+
+         lwup_tmp(ji) = emis(ji) * c_stefan * temp_sol(ji)**4 + &
+           &     quatre * emis(ji) * c_stefan * temp_sol(ji)**3 * &
+           &     (tp_00 - temp_sol(ji))
+}}}