Changes between Version 43 and Version 44 of DevelopmentActivities/MergeHydro/Martial_notes_on_merge

Timestamp:

2012-01-04T14:29:49+01:00 (12 years ago)

Author:

mmaipsl

Comment:

--

DevelopmentActivities/MergeHydro/Martial_notes_on_merge

@@ -120 +120 @@
  * '''hydrol_canop''' : je me pose la question du '''veget/veget_max''' ici car on définit '''qsintveg et precisol''' dans cette routine. 
    Si l'on passe veget_max, on ne tient plus compte du LAI pour calculer ces deux variables importantes. Je laisse donc veget.
-   C'était bien un bogue de la version LMD.
+   C'était bien un bogue de la version LMD. [[BR]]
+   Voici les corrections que l'on a appliqué pour rendre le modèle conservatif :
+{{{
+    precisol(:,1)=veget_max(:,1)*precip_rain(:)
+    DO jv=2,nvm
+       DO ji = 1, kjpindex
+          zqsintvegnew(ji,jv) = MIN (qsintveg(ji,jv),qsintmax(ji,jv)) 
+          IF ( ok_throughfall_by_pft ) THEN
+             ! correction throughfall Nathalie - Juin 2006
+             precisol(ji,jv) = (veget(ji,jv)*throughfall_by_pft(jv)*precip_rain(ji)) + &
+                  qsintveg(ji,jv) - zqsintvegnew (ji,jv) + &
+                  (veget_max(ji,jv) - veget(ji,jv))*precip_rain(ji)
+          ELSE
+             precisol(ji,jv) = qsintveg(ji,jv) - zqsintvegnew (ji,jv) + &
+                  (veget_max(ji,jv) - veget(ji,jv))*precip_rain(ji)
+          ENDIF
+       ENDDO
+    END DO
+[...]
+    DO jv=2,nvm
+       qsintveg(:,jv) = zqsintvegnew (:,jv)
+    END DO
+}}}
+   On ne fait donc plus aucun calcul de l'interception pour le PFT 1. On prends bien en compte le precip_rain sur les surface fracnobio dans hydrol_snow.
+
 == Formules de convolution/déconvolution ==
 La convolution est définit par les corr_veg_soil ("(:,nvm,nstm) percentage of each veg. type on each soil of each grid point")