Changeset 710 for LMDZ4/branches/IPSL-CM4_IPCC_branch/libf/phylmd/physiq.F

Timestamp:

Sep 20, 2006, 12:54:27 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Integration des modifications d'OM et JLD pour une meilleure prise en compte
de la conservation du flux d'eau
LF

File:

• LMDZ4/branches/IPSL-CM4_IPCC_branch/libf/phylmd/physiq.F (modified) (8 diffs)

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

@@ -3 +3 @@
 !
 C
-c $Header$
-c
       SUBROUTINE physiq (nlon,nlev,nqmax,
      .            debut,lafin,rjourvrai,gmtime,pdtphys,
@@ -592 +590 @@
 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 +971 @@
 c Declaration des constantes et des fonctions thermodynamiques
 c
+         SAVE fqfonte
 #include "YOMCST.h"
 #include "YOETHF.h"
@@ -1026 +1026 @@
        IF (debut) THEN
 C
+         fqfonte(:,:)=0.
          IF (if_ebil.ge.1) d_h_vcol_phy=0.
 c
@@ -1510 +1511 @@
      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 +1578 @@
          zxffonte(i) = 0.0
          zxfqcalving(i) = 0.0
+         zxfqfonte(i) = 0.0
 c
          IF ( abs( pctsrf(i, is_ter) + pctsrf(i, is_lic) + 
@@ -1603 +1605 @@
             zxfqcalving(i) = zxfqcalving(i) + 
      .                      fqcalving(i,nsrf)*pctsrf(i,nsrf)
+            zxfqfonte(i) = zxfqfonte(i) + 
+     .                      fqfonte(i,nsrf)*pctsrf(i,nsrf)
 c        ENDIF 
         ENDDO
@@ -1622 +1626 @@
           IF (pctsrf(i,nsrf) .LT. epsfra) 
      .    fqcalving(i,nsrf) = zxfqcalving(i)
+          IF (pctsrf(i,nsrf) .LT. epsfra) 
+     .    fqfonte(i,nsrf) = zxfqfonte(i)
         ENDDO
       ENDDO