Changeset 705 for LMDZ4/branches/IPSL-CM4_IPCC_patches/libf/phylmd/physiq.F

Timestamp:

Sep 19, 2006, 4:23:21 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Pour la conservation du flux d'eau OM, JLD
LF

File:

• LMDZ4/branches/IPSL-CM4_IPCC_patches/libf/phylmd/physiq.F (modified) (7 diffs)

LMDZ4/branches/IPSL-CM4_IPCC_patches/libf/phylmd/physiq.F

@@ -592 +592 @@
 c                               !et necessaire pour limiter la
 c                               !hauteur de neige, en kg/m2/s
-      REAL zxffonte(klon), zxfqcalving(klon)
+      REAL fqfonte(klon,nbsrf)  !Quantite d'eau de fonte des glaciers
+      REAL zxffonte(klon), zxfqcalving(klon),zxfqfonte(klon)
 
       LOGICAL offline           ! Controle du stockage ds "physique"
@@ -972 +973 @@
 c Declaration des constantes et des fonctions thermodynamiques
 c
+         SAVE fqfonte
 #include "YOMCST.h"
 #include "YOETHF.h"
@@ -1026 +1028 @@
        IF (debut) THEN
 C
+         fqfonte(:,:)=0.
          IF (if_ebil.ge.1) d_h_vcol_phy=0.
 c
@@ -1510 +1513 @@
      s            dsens, devap,
      s            ycoefh,yu1,yv1, t2m, q2m, u10m, v10m,
-     s            fqcalving, ffonte, run_off_lic_0) 
+     s            fqcalving, fqfonte,ffonte, run_off_lic_0) 
 c
 CXXX PB
@@ -1577 +1580 @@
          zxffonte(i) = 0.0
          zxfqcalving(i) = 0.0
+         zxfqfonte(i) = 0.0
 c
          IF ( abs( pctsrf(i, is_ter) + pctsrf(i, is_lic) + 
@@ -1603 +1607 @@
             zxfqcalving(i) = zxfqcalving(i) + 
      .                      fqcalving(i,nsrf)*pctsrf(i,nsrf)
+            zxfqfonte(i) = zxfqfonte(i) + 
+     .                      fqfonte(i,nsrf)*pctsrf(i,nsrf)
 c        ENDIF 
         ENDDO
@@ -1622 +1628 @@
           IF (pctsrf(i,nsrf) .LT. epsfra) 
      .    fqcalving(i,nsrf) = zxfqcalving(i)
+          IF (pctsrf(i,nsrf) .LT. epsfra) 
+     .    fqfonte(i,nsrf) = zxfqfonte(i)
         ENDDO
       ENDDO