Changes between Version 36 and Version 37 of DevelopmentActivities/MergeHydro/Martial_notes_on_merge

Timestamp:

2011-12-02T12:04:14+01:00 (13 years ago)

Author:

mmaipsl

Comment:

--

DevelopmentActivities/MergeHydro/Martial_notes_on_merge

@@ -122 +122 @@
    C'était bien un bogue de la version LMD.
 == Formules de convolution/déconvolution ==
-La convolution est définit par les corr_veg_soil
+La convolution est définit par les corr_veg_soil ("(:,nvm,nstm) percentage of each veg. type on each soil of each grid point")
+et les cvs_over_veg ("(:,nvm,nstm) old value of corr_veg_soil/veget_max kept from diag to next split" : commentaire sûrement faux d'après la définition).
 {{{
        ! somme(corr_veg_soil / vegtot / veget_max ) = 1 - vegtot = frac_nobio
@@ -145 +146 @@
        ENDDO
 }}}
+et 
+{{{
+    cvs_over_veg(:,:,:) = zero
+    DO jv=1,nvm
+       DO ji=1,kjpindex
+          IF(veget_max(ji,jv).GT.min_sechiba) THEN
+             DO jst=1,nstm
+                cvs_over_veg(ji,jv,jst) = corr_veg_soil(ji,jv,jst)/vegtot(ji) / veget_max(ji,jv)
+             ENDDO
+          ENDIF
+       END DO
+    END DO
+}}}
+
 Donc pour chaque point de terre :
 {{{
 #!formula
 #density=100
-\displaystyle \sum_{jv=1,nvm} \sum_{jst=1,nstm}\left( \frac {corr\_veg\_soil(jv,jst)}{vegtot \times veget\_max(jv)} \right) = 1 - vegtot= frac\_nobio 
+\[
+\begin{array}{lcl}
+\displaystyle \sum_{jv=1,nvm} \sum_{jst=1,nstm}\frac {corr\_veg\_soil(jv,jst)}{vegtot \times veget\_max(jv)} &=& \\
+\displaystyle \sum_{jv=1,nvm} \sum_{jst=1,nstm} cvs\_over\_veg(jv,jst)  &=& 1 - vegtot \\
+&=& frac\_nobio  
+\]
 }}}
 soit 
@@ -155 +175 @@
 #!formula
 #density=100
-\displaystyle \sum_{jv=1,nvm} \sum_{jst=1,nstm}\left( \frac {corr\_veg\_soil(jv,jst)}{veget\_max(jv)} \right) = vegtot\times (1 - vegtot) = vegtot \times frac\_nobio 
-}}}
+\displaystyle \sum_{jv=1,nvm} \sum_{jst=1,nstm}\frac {corr\_veg\_soil(jv,jst)}{veget\_max(jv)} = vegtot\times (1 - vegtot) = vegtot \times frac\_nobio 
+}}}
+Ce qui semble assez étrange, 
+
 La formule de convolution est alors :
 {{{
+    !
+    !
+    ! split 2d variables into 3d variables, per soil type
+    !
+    precisol_ns(:,:)=zero
     DO jv=1,nvm
        DO jst=1,nstm
@@ -169 +196 @@
        END DO
     END DO
-
+}}}
+et la fomule de déconvolution est donnée par le test :
+{{{
     !
     ! Now we check if the deconvolution is correct and conserves the fluxes:
@@ -220 +249 @@
 
     ENDIF
-
-}}}
+}}}
+
 == traitement des corrections de Nathalie ==
 === Gestion du throughfall_by_pft ===