Changeset 156 for LMDZ.3.3/branches/rel-LF/libf

Timestamp:

Nov 10, 2000, 11:45:53 AM (24 years ago)

Author:

lmdzadmin

Message:

Menage et initialisation a zero des variables locales
LF

File:

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

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

@@ -43 +43 @@
 c pplay----input-R- pression au milieu de couche (Pa)
 c radsol---input-R- flux radiatif net (positif vers le sol) en W/m**2
-c capsol---input-R- inversion de l'effective capacite du sol (J/m2/K)
-c beta-----input-R- coefficient de l'evaporation reelle (0 a 1)
-c dif_grnd-input-R- coeff. de diffusion (chaleur) vers le sol profond
 c rlat-----input-R- latitude en degree
 c rugos----input-R- longeur de rugosite (en m)
@@ -74 +71 @@
 #include "indicesol.h"
 c
-      LOGICAL soil_model
 c
       REAL dtime
@@ -121 +117 @@
 c======================================================================
       REAL yts(klon), yrugos(klon), ypct(klon), yz0_new(klon)
-      REAL ycal(klon), ybeta(klon), ydif(klon), yalb(klon),yevap(klon)
+      REAL yalb(klon),yevap(klon)
       REAL yu1(klon), yv1(klon)
       real ysnow(klon), yqsol(klon)
       real yrain_f(klon), ysnow_f(klon)
-      real ysollw(klon), ysolsw(klon), ysolswnet(klon)
+      real ysollw(klon), ysolsw(klon)
       real yfder(klon), ytaux(klon), ytauy(klon)
       REAL yrugm(klon), yrads(klon)
@@ -135 +131 @@
       REAL y_dflux_t(klon), y_dflux_q(klon)
       REAL ycoefh(klon,klev), ycoefm(klon,klev)
-      REAL ygamt(klon,2:klev) ! contre-gradient pour temperature
-      REAL ygamq(klon,2:klev) ! contre-gradient pour humidite
       REAL yu(klon,klev), yv(klon,klev)
       REAL yt(klon,klev), yq(klon,klev)
@@ -149 +143 @@
       LOGICAL ok_nonloc
       PARAMETER (ok_nonloc=.FALSE.)
-      REAL y_cd_h(klon), y_cd_m(klon)
       REAL ycoefm0(klon,klev), ycoefh0(klon,klev)
 c
@@ -155 +148 @@
       REAL u1lay(klon), v1lay(klon)
       REAL delp(klon,klev)
-      REAL capsol(klon), beta(klon), dif_grnd(klon)
-      REAL cal(klon)
-      REAL soilcap(klon,nbsrf), soilflux(klon,nbsrf)
       REAL totalflu(klon)
       INTEGER i, k, nsrf 
@@ -194 +184 @@
          zv1(i) = 0.0
       ENDDO
+      ypct = 0.0
+      yts = 0.0
+      ysnow = 0.0
+      yevap = 0.0
+      yqsol = 0.0
+      yalb = 0.0
+      yrain_f = 0.0
+      ysnow_f = 0.0
+      yfder = 0.0
+      ytaux = 0.0
+      ytauy = 0.0
+      ysolsw = 0.0
+      ysollw = 0.0
+      yrugos = 0.0
+      yu1 = 0.0
+      yv1 = 0.0
+      yrads = 0.0
+      ypaprs = 0.0
+      ypaprs = 0.0
+      ypplay = 0.0
+      ydelp = 0.0
+      yu = 0.0
+      yv = 0.0
+      yt = 0.0
+      yq = 0.0
+      pctsrf_new = 0.0
+      y_flux_u = 0.0
+      y_flux_v = 0.0
+
       DO nsrf = 1, nbsrf
       DO i = 1, klon
@@ -238 +257 @@
 
       DO 99999 nsrf = 1, nbsrf
-c
-c prescrire les proprietes du sol:
-c      CALL calbeta(dtime,nsrf,snow,qsol, beta,capsol,dif_grnd)
-c      IF (.NOT.soil_model) THEN
-c         DO i = 1, klon
-c            cal(i) = RCPD * capsol(i)
-c            totalflu(i) = radsol(i)
-c         ENDDO
-c      ELSE
-c         DO i = 1, klon
-c            totalflu(i) = soilflux(i,nsrf) + radsol(i)
-c            IF (nsrf.EQ.is_oce) THEN
-c               cal(i) = 0.0
-c            ELSE
-c               cal(i) = RCPD / soilcap(i,nsrf)
-c            ENDIF
-c         ENDDO
-c      ENDIF
-c
       totalflu = radsol
 
@@ -339 +339 @@
      e          knon, nsrf, ni, pctsrf,
      e          ok_veget, ocean, npas, nexca,
-     e          rmu0,
+     e          rmu0, yrugos,
      e          yu1, yv1, ycoefh,
      e          yt,yq,yts,ypaprs,ypplay,
@@ -437 +437 @@
      e                knon, nisurf, knindex, pctsrf,
      e                ok_veget, ocean, npas, nexca,
-     e                rmu0,
+     e                rmu0, rugos,
      e                u1lay,v1lay,coef,
      e                t,q,ts,paprs,pplay,
@@ -469 +469 @@
 c                               vent (dV/dz); la premiere valeur
 c                               indique la valeur de Cdrag (sans unite)
-      REAL cal(klon)          ! Cp/cal indique la capacite calorifique
-c                               surfacique du sol
-      REAL beta(klon)         ! evap. reelle / evapotranspiration
-      REAL dif_grnd(klon)     ! coeff. diffusion vers le sol profond
       REAL t(klon,klev)       ! temperature (K)
       REAL q(klon,klev)       ! humidite specifique (kg/kg)
@@ -488 +484 @@
       integer jour            ! jour de l'annee en cours
       real rmu0(klon)         ! cosinus de l'angle solaire zenithal
+      real rugos(klon)        ! rugosite
       integer knindex(klon)
       real pctsrf(klon,nbsrf)
@@ -505 +502 @@
       REAL dflux_s(klon) ! derivee du flux sensible dF/dTs
       REAL dflux_l(klon) ! derivee du flux latent dF/dTs
-      REAL dflux_g(klon) ! derivee du flux du sol profond dF/dTs
 c======================================================================
       REAL t_grnd  ! temperature de rappel pour glace de mer
@@ -513 +509 @@
 c======================================================================
       INTEGER i, k
-      REAL zx_a
-      REAL zx_b
-      REAL zx_qs
-      REAL zx_dq_s_dh
-      REAL zx_h_grnd
-      REAL zx_cq0(klon)
-      REAL zx_dq0(klon)
       REAL zx_cq(klon,klev)
       REAL zx_dq(klon,klev)
@@ -527 +516 @@
       REAL zx_buf2(klon)
       REAL zx_coef(klon,klev)
-      REAL zx_q_0(klon)
-      REAL zx_h_ts(klon)
-      REAL zx_sl(klon)
       REAL local_h(klon,klev) ! enthalpie potentielle
       REAL local_q(klon,klev)
@@ -543 +529 @@
       REAL zdelz
 c======================================================================
-C Variables intermediaires pour le calcul des fluxs a la surface
-      real zx_mh(klon), zx_nh(klon), zx_oh(klon)
-      real zx_mq(klon), zx_nq(klon), zx_oq(klon)
-      real zx_k1(klon), zx_dq_s_dt(klon)
-      real zx_qsat(klon)
-c======================================================================
-      REAL zcor, zdelta, zcvm5
       logical contreg
       parameter (contreg=.true.)
@@ -563 +542 @@
       real tq_cdrag(klon), petAcoef(klon), peqAcoef(klon)
       real petBcoef(klon), peqBcoef(klon)
-      real lwdown(klon), swnet(klon), swdown(klon), ps(klon)
+      real lwdown(klon), swnet(klon), swdown(klon)
       real p1lay(klon)
-      real coef1lay(klon)
 
 ! Parametres de sortie
@@ -703 +681 @@
      e tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef,
      e precip_rain, precip_snow, lwdown, swnet, swdown,
-     e fder, taux, tauy,
+     e fder, taux, tauy, rugos,
      e albedo, snow, qsol,
      e ts, p1lay, psref, radsol,
@@ -891 +869 @@
 c ts-------input-R- temperature du sol (en Kelvin)
 c rugos----input-R- longeur de rugosite (en m)
-c xlat-----input-R- latitude en degree
 c u--------input-R- vitesse u
 c v--------input-R- vitesse v
@@ -913 +890 @@
       REAL u(klon,klev), v(klon,klev), t(klon,klev), q(klon,klev)
       REAL rugos(klon)
-      REAL xlat(klon)
 c
       REAL pcfm(klon,klev), pcfh(klon,klev)
@@ -1293 +1269 @@
 c Calculer quelques parametres pour appliquer la couche limite
 c ------------------------------------------------------------
-!#include "dimensions.h"
-!#include "dimphy.h"
+#include "dimensions.h"
+#include "dimphy.h"
 #include "YOMCST.h"
 #include "indicesol.h"
@@ -1311 +1287 @@
 c
       REAL dtime
-      REAL snow(knon,nbsrf), qsol(knon,nbsrf)
+      REAL snow(klon), qsol(klon)
       INTEGER indice, knon
 C
-      REAL vbeta(knon)
-      REAL vcal(knon)
-      REAL vdif(knon)
-C
+      REAL vbeta(klon)
+      REAL vcal(klon)
+      REAL vdif(klon)
+C
+
       IF (indice.EQ.is_oce) THEN
       DO i = 1, knon
@@ -1329 +1306 @@
       DO i = 1, knon
           vcal(i) = calice
-          IF (snow(i,is_sic) .GT. 0.0) vcal(i) = calsno
+          IF (snow(i) .GT. 0.0) vcal(i) = calsno
           vbeta(i)  = 1.0
           vdif(i) = 1.0/tau_gl
@@ -1339 +1316 @@
       DO i = 1, knon
           vcal(i) = calsol
-          IF (snow(i,is_ter) .GT. 0.0) vcal(i) = calsno
-          vbeta(i)  = MIN(2.0*qsol(i,is_ter)/mx_eau_sol, 1.0)
+          IF (snow(i) .GT. 0.0) vcal(i) = calsno
+          vbeta(i)  = MIN(2.0*qsol(i)/mx_eau_sol, 1.0)
           vdif(i) = 0.0
       ENDDO
@@ -1348 +1325 @@
       DO i = 1, knon
           vcal(i) = calice
-          IF (snow(i,is_lic) .GT. 0.0) vcal(i) = calsno
+          IF (snow(i) .GT. 0.0) vcal(i) = calsno
           vbeta(i)  = 1.0
           vdif(i) = 0.0