Changeset 454

Timestamp:

Mar 7, 2003, 5:35:57 PM (22 years ago)

Author:

lmdzadmin

Message:

Ajout diagnostics de fonte glace et calving IM/JLD

File:

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

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

@@ -16 +16 @@
      .                  flux_t,flux_q,flux_u,flux_v,cdragh,cdragm,
      .                  dflux_t,dflux_q,
-     .                  zcoefh,zu1,zv1, t2m, q2m, u10m, v10m)
+cIM cf JLD    .                  zcoefh,zu1,zv1, t2m, q2m, u10m, v10m)
+     .                  zcoefh,zu1,zv1, t2m, q2m, u10m, v10m,
+     .                  fqcalving,ffonte)
 cAA .                  itr, tr, flux_surf, d_tr)
 cAA REM:
@@ -68 +70 @@
 c dflux_t derive du flux sensible
 c dflux_q derive du flux latent
+c ffonte----Flux thermique utilise pour fondre la neige
+c fqcalving-Flux d'eau "perdue" par la surface et necessaire pour limiter la
+c           hauteur de neige, en kg/m2/s
 cAA on rajoute en output yu1 et yv1 qui sont les vents dans 
 cAA la premiere couche
@@ -93 +98 @@
       REAL flux_t(klon,klev, nbsrf), flux_q(klon,klev, nbsrf)
       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 flux_u(klon,klev, nbsrf), flux_v(klon,klev, nbsrf)
       REAL rugmer(klon), agesno(klon,nbsrf),rugoro(klon)
@@ -175 +184 @@
       REAL u1lay(klon), v1lay(klon)
       REAL delp(klon,klev)
-      REAL totalflu(klon)
       INTEGER i, k, nsrf 
 cAA   INTEGER it
@@ -368 +376 @@
 
       DO 99999 nsrf = 1, nbsrf
-       totalflu = radsol
 
 c chercher les indices:
@@ -419 +426 @@
         ytaux(j) = flux_u(i,1,nsrf)
         ytauy(j) = flux_v(i,1,nsrf)
-c$$$    ysolsw(j) = solsw(i)
-cIM cf. JLD
         ysolsw(j) = solsw(i,nsrf)
-c       ysolsw(j) = (1 - albe(i,nsrf)) 
-c    $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter) 
-c    $      - pctsrf(i, is_lic) *albe(i,is_lic) 
-c    $      - pctsrf(i, is_oce) *albe(i,is_oce)
-c    $      - pctsrf(i, is_sic) *albe(i,is_sic) 
-c    $      ) * solsw(i)
-cIM cf. JLD     ysollw(j) = sollw(i)
         ysollw(j) = sollw(i,nsrf)
         ysollwdown(j) = sollwdown(i)
@@ -435 +433 @@
         yu1(j) = u1lay(i)
         yv1(j) = v1lay(i)
-c$$$        yrads(j) = totalflu(i)
-cIM cf. JLD
         yrads(j) =  ysolsw(j)+ ysollw(j)
-c       yrads(j) = (1 - albe(i,nsrf)) 
-c    $      /(1 - pctsrf(i,is_ter) * albe(i,is_ter) 
-c    $      - pctsrf(i, is_lic) *albe(i,is_lic) 
-c    $      - pctsrf(i, is_oce) *albe(i,is_oce)
-c    $      - pctsrf(i, is_sic) *albe(i,is_sic) 
-cIM cf. JLD  $      ) * solsw(i) + sollw(i)
-c    $      ) * solsw(i) + ysollw(j)
         ypaprs(j,klev+1) = paprs(i,klev+1)
       END DO
-      IF ( nsrf .eq. is_ter ) THEN 
+C
+C     IF bucket model for continent, copy soil water content
+      IF ( nsrf .eq. is_ter .and. .not. ok_veget ) THEN 
           DO j = 1, knon
             i = ni(j)
@@ -489 +480 @@
       ENDDO
       ENDDO
-
 c
 c
@@ -522 +512 @@
      s          pctsrf_new, yagesno,
      s          y_d_t, y_d_q, y_d_ts, yz0_new,
-     s          y_flux_t, y_flux_q, y_dflux_t, y_dflux_q)
+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)
 c
 c calculer la longueur de rugosite sur ocean
@@ -583 +575 @@
            rugos(i,nsrf) = yrugm(j)
          endif  
+cIM cf JLD ??
+         fqcalving(i,nsrf) = y_fqcalving(j)        
+         ffonte(i,nsrf) = y_ffonte(j)        
          cdragh(i) = cdragh(i) + ycoefh(j,1)
          cdragm(i) = cdragm(i) + ycoefm(j,1)
@@ -751 +746 @@
      s                pctsrf_new, agesno,
      s                d_t, d_q, d_ts, z0_new, 
-     s                flux_t, flux_q,dflux_s,dflux_l)
+cIM cf JLD    s                flux_t, flux_q,dflux_s,dflux_l)
+     s                flux_t, flux_q,dflux_s,dflux_l,
+     s                fqcalving,ffonte)
 
       USE interface_surf
@@ -812 +809 @@
       REAL dflux_s(klon) ! derivee du flux sensible dF/dTs
       REAL dflux_l(klon) ! derivee du flux latent dF/dTs
+cIM cf JLD
+c Flux thermique utiliser pour fondre la neige
+      REAL ffonte(klon)
+c Flux d'eau "perdue" par la surface et nécessaire pour que limiter la
+c hauteur de neige, en kg/m2/s
+      REAL fqcalving(klon)
 c======================================================================
       REAL t_grnd  ! temperature de rappel pour glace de mer
@@ -1028 +1031 @@
      s evap, fluxsens, fluxlat, dflux_l, dflux_s,              
      s tsol_rad, tsurf_new, alb_new, alblw, emis_new, z0_new, 
-     s pctsrf_new, agesno)
+cIM cf JLD    s pctsrf_new, agesno)
+     s pctsrf_new, agesno,fqcalving,ffonte)