Changeset 522 for LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F

Timestamp:

May 6, 2004, 10:57:36 AM (20 years ago)

Author:

lmdzadmin

Message:

Modifications pour la fermeture en eau (fonte des glaciers) JLD, OM
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F (modified) (8 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F

@@ -18 +18 @@
 cIM cf JLD    .                  zcoefh,zu1,zv1, t2m, q2m, u10m, v10m)
      .                  zcoefh,zu1,zv1, t2m, q2m, u10m, v10m,
-     .                  fqcalving,ffonte)
+     .                  fqcalving,ffonte, run_off_lic_0)
 cAA .                  itr, tr, flux_surf, d_tr)
 cAA REM:
@@ -101 +101 @@
       REAL y_fqcalving(klon), y_ffonte(klon)
       REAL fqcalving(klon,nbsrf), ffonte(klon,nbsrf)
+      REAL run_off_lic_0(klon), y_run_off_lic_0(klon)
 
       REAL flux_u(klon,klev, nbsrf), flux_v(klon,klev, nbsrf)
@@ -442 +443 @@
         yrads(j) =  ysolsw(j)+ ysollw(j)
         ypaprs(j,klev+1) = paprs(i,klev+1)
+        y_run_off_lic_0(j) = run_off_lic_0(i)
       END DO
 C
@@ -557 +559 @@
 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)
+     s          y_fqcalving,y_ffonte,y_run_off_lic_0)
 c
 c calculer la longueur de rugosite sur ocean
@@ -635 +637 @@
       END DO
       END IF 
+      IF ( nsrf .eq. is_lic ) THEN 
+        DO j = 1, knon
+          i = ni(j)
+          run_off_lic_0(i) = y_run_off_lic_0(j)
+        END DO
+      END IF 
 c$$$ PB ajout pour soil
       ftsoil(:,:,nsrf) = 0.
@@ -793 +801 @@
 cIM cf JLD    s                flux_t, flux_q,dflux_s,dflux_l)
      s                flux_t, flux_q,dflux_s,dflux_l,
-     s                fqcalving,ffonte)
+     s                fqcalving,ffonte,run_off_lic_0)
 
       USE interface_surf
@@ -835 +843 @@
       REAL agesno(klon)
       REAL rugoro(klon)
+      REAL run_off_lic_0(klon)! runof glacier au pas de temps precedent
       integer jour            ! jour de l'annee en cours
       real rmu0(klon)         ! cosinus de l'angle solaire zenithal
@@ -1078 +1087 @@
      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)
+     s pctsrf_new, agesno,fqcalving,ffonte, run_off_lic_0)