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

Timestamp:

2011-12-02T11:53:05+01:00 (13 years ago)

Author:

mmaipsl

Comment:

--

DevelopmentActivities/MergeHydro/Martial_notes_on_merge

@@ -121 +121 @@
    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.
+== Formules de convolution/déconvolution ==
+La convolution est définit par les corr_veg_soil
+{{{
+       ! somme(corr_veg_soil / vegtot / veget_max ) = 1 - vegtot = frac_nobio
+       DO ji=1,kjpindex
+          tot_corr_veg_soil(ji)=zero
+          DO jst = 1, nstm
+             DO jv = 1,nvm
+                tot_corr_veg_soil(ji)=tot_corr_veg_soil(ji)+corr_veg_soil(ji,jv,jst) / vegtot(ji) / veget_max(ji,jv)
+             ENDDO
+          ENDDO
+       ENDDO
+
+       DO ji=1,kjpindex
+          IF ( ABS( tot_corr_veg_soil(ji) - (1 - vegtot(ji)) ) .GT. EPS1 ) THEN
+             WRITE(numout,*) 'corr_veg_soil SPLIT FALSE:ji=',ji,tot_corr_veg_soil(ji),(1 - vegtot(ji))
+             WRITE(numout,*) 'err',ABS( tot_corr_veg_soil(ji) - (1 - vegtot(ji)) )
+             WRITE(numout,*) 'vegtot',vegtot(ji)
+             DO jv=1,nvm
+                WRITE(numout,*) 'jv,veget_max,corr_veg_soil',jv,veget_max(ji,jv),corr_veg_soil(ji,jv,:)
+             END DO
+          ENDIF
+       ENDDO
+}}}
+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 
+}}}
+soit 
+{{{
+#!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 
+}}}
+La formule de convolution est alors :
+{{{
+    DO jv=1,nvm
+       DO jst=1,nstm
+          DO ji=1,kjpindex
+             IF(veget_max(ji,jv).GT.min_sechiba) THEN
+                precisol_ns(ji,jst)=precisol_ns(ji,jst)+precisol(ji,jv)* &
+                     & corr_veg_soil(ji,jv,jst) / vegtot(ji) / veget_max(ji,jv))
+             ENDIF
+          END DO
+       END DO
+    END DO
+
+    !
+    ! Now we check if the deconvolution is correct and conserves the fluxes:
+
+    IF (check_cwrr) THEN
+
+
+       tmp_check1(:)=zero
+       tmp_check2(:)=zero  
+
+       ! First we check the precisol and evapnu
+
+       DO jst=1,nstm
+          DO ji=1,kjpindex
+             tmp_check1(ji)=tmp_check1(ji) + &
+                  & precisol_ns(ji,jst)*soiltile(ji,jst)*vegtot(ji)
+          END DO
+       END DO
+
+       DO jv=1,nvm
+          DO ji=1,kjpindex
+             tmp_check2(ji)=tmp_check2(ji) + precisol(ji,jv)
+          END DO
+       END DO
+
+       DO ji=1,kjpindex   
+
+          IF(ABS(tmp_check1(ji)- tmp_check2(ji)).GT.allowed_err) THEN
+             WRITE(numout,*) 'PRECISOL SPLIT FALSE:ji=',ji,tmp_check1(ji),tmp_check2(ji)
+             WRITE(numout,*) 'err',ABS(tmp_check1(ji)- tmp_check2(ji))
+             WRITE(numout,*) 'vegtot',vegtot(ji)
+
+             DO jv=1,nvm
+                WRITE(numout,*) 'jv,veget_max, precisol',jv,veget_max(ji,jv),precisol(ji,jv)
+                DO jst=1,nstm
+                   WRITE(numout,*) 'corr_veg_soil:jst',jst,corr_veg_soil(ji,jv,jst)
+                END DO
+             END DO
+
+             DO jst=1,nstm
+                WRITE(numout,*) 'jst,precisol_ns',jst,precisol_ns(ji,jst)
+                WRITE(numout,*) 'soiltile', soiltile(ji,jst)
+             END DO
+             waterbal_error=.TRUE.
+             CALL ipslerr(2, 'hydrol_split_soil', 'We will STOP after hydrol_split_soil.',&
+                  & 'check_CWRR','PRECISOL SPLIT FALSE')
+          ENDIF
+
+       END DO
+
+    ENDIF
+
+}}}
 == traitement des corrections de Nathalie ==
 === Gestion du throughfall_by_pft ===