Changeset 705 for LMDZ4/branches/IPSL-CM4_IPCC_patches/libf/phylmd/clmain.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/clmain.F (modified) (8 diffs)

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

@@ -20 +20 @@
 cIM cf JLD    .                  zcoefh,zu1,zv1, t2m, q2m, u10m, v10m)
      .                  zcoefh,zu1,zv1, t2m, q2m, u10m, v10m,
-     .                  fqcalving,ffonte, run_off_lic_0)
+     .                  fqcalving,fqfonte,ffonte, run_off_lic_0)
 cAA .                  itr, tr, flux_surf, d_tr)
 cAA REM:
@@ -73 +73 @@
 c dflux_q derive du flux latent
 c ffonte----Flux thermique utilise pour fondre la neige
+c fqfonte -- quantite d'eau due a la fonte de la calotte
 c fqcalving-Flux d'eau "perdue" par la surface et necessaire pour limiter la
 c           hauteur de neige, en kg/m2/s
@@ -102 +103 @@
       REAL dflux_t(klon), dflux_q(klon)
 cIM cf JLD
-      REAL y_fqcalving(klon), y_ffonte(klon)
-      REAL fqcalving(klon,nbsrf), ffonte(klon,nbsrf)
+      REAL y_fqcalving(klon), y_fqfonte(klon), y_ffonte(klon)
+      REAL fqcalving(klon,nbsrf), fqfonte(klon,nbsrf), ffonte(klon,nbsrf)
       REAL run_off_lic_0(klon), y_run_off_lic_0(klon)
 
@@ -579 +580 @@
 cIM cf JLD    s          y_flux_t, y_flux_q, y_dflux_t, y_dflux_q)
      s          y_flux_t, y_flux_q, y_dflux_t, y_dflux_q,
-     s          y_fqcalving,y_ffonte,y_run_off_lic_0)
+     s          y_fqcalving,y_fqfonte,y_ffonte,y_run_off_lic_0)
 c
 c calculer la longueur de rugosite sur ocean
@@ -643 +644 @@
 cIM cf JLD ??
          fqcalving(i,nsrf) = y_fqcalving(j)        
+         fqfonte(i,nsrf) = y_fqfonte(j)        
          ffonte(i,nsrf) = y_ffonte(j)        
          cdragh(i) = cdragh(i) + ycoefh(j,1)
@@ -824 +826 @@
 cIM cf JLD    s                flux_t, flux_q,dflux_s,dflux_l)
      s                flux_t, flux_q,dflux_s,dflux_l,
-     s                fqcalving,ffonte,run_off_lic_0)
+     s                fqcalving,fqfonte,ffonte,run_off_lic_0)
 
       USE interface_surf
@@ -900 +902 @@
 c hauteur de neige, en kg/m2/s
       REAL fqcalving(klon)
+c Quantite d'eau de fonte de la calotte kg/mé/s
+      REAL fqfonte(klon)
 c======================================================================
       REAL t_grnd  ! temperature de rappel pour glace de mer
@@ -1120 +1124 @@
      s tsol_rad, tsurf_new, alb_new, alblw, emis_new, z0_new, 
 cIM cf JLD    s pctsrf_new, agesno)
-     s pctsrf_new, agesno,fqcalving,ffonte, run_off_lic_0)
+     s pctsrf_new, agesno,fqcalving,fqfonte,ffonte, run_off_lic_0)