Changeset 280 for LMDZ.3.3/branches

Timestamp:

Oct 19, 2001, 12:29:20 PM (23 years ago)

Author:

lmdzadmin

Message:

Modif dans coefkz pour remplacer l'humidite du niveau 1 par qsat(tsol)
Passage des deux albedos de surface (vis et nir)
LF

File:

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

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

@@ -8 +8 @@
      .                  ok_veget, ocean, npas, nexca, ts,
      .                  soil_model,ftsoil,
-     .                  paprs,pplay,radsol,snow,qsol,evap,albe,fluxlat,
+     .                  paprs,pplay,radsol,snow,qsol,evap,albe,alblw,
+     .                  fluxlat,
      .                  rain_f, snow_f, solsw, sollw, sollwdown, fder,
      .                  rlon, rlat, cufi, cvfi, rugos,
@@ -112 +113 @@
       REAL evap(klon,nbsrf)
       REAL albe(klon,nbsrf)
+      REAL alblw(klon,nbsrf)
 c$$$ PB
       REAL fluxlat(klon,nbsrf)
@@ -135 +137 @@
       REAL yts(klon), yrugos(klon), ypct(klon), yz0_new(klon)
       REAL yalb(klon)
+      REAL yalblw(klon)
       REAL yu1(klon), yv1(klon)
       real ysnow(klon), yqsol(klon), yagesno(klon)
@@ -258 +261 @@
       yqsol = 0.0
       yalb = 0.0
+      yalblw = 0.0
       yrain_f = 0.0
       ysnow_f = 0.0
@@ -370 +374 @@
         yqsol(j) = qsol(i,nsrf)
         yalb(j) = albe(i,nsrf)
+        yalblw(j) = alblw(i,nsrf)
         yrain_f(j) = rain_f(i)
         ysnow_f(j) = snow_f(i)
@@ -452 +457 @@
      e          yu1, yv1, ycoefh,
      e          yt,yq,yts,ypaprs,ypplay,
-     e          ydelp,yrads,yalb, ysnow, yqsol, 
+     e          ydelp,yrads,yalb, yalblw, ysnow, yqsol, 
      e          yrain_f, ysnow_f, yfder, ytaux, ytauy, 
 c$$$     e          ysollw, ysolsw,
@@ -499 +504 @@
 c
       albe(:, nsrf) = 0.
+      alblw(:, nsrf) = 0.
       snow(:, nsrf) = 0.
       qsol(:, nsrf) = 0.
@@ -507 +513 @@
          d_ts(i,nsrf) = y_d_ts(j)
          albe(i,nsrf) = yalb(j)
+         alblw(i,nsrf) = yalblw(j)
          snow(i,nsrf) = ysnow(j)
          qsol(i,nsrf) = yqsol(j)
@@ -572 +579 @@
      e                u1lay,v1lay,coef,
      e                t,q,ts,paprs,pplay,
-     e                delp,radsol,albedo,snow,qsol, 
+     e                delp,radsol,albedo,alblw,snow,qsol, 
      e                precip_rain, precip_snow, fder, taux, tauy,
      $                sollw, sollwdown, swdown,fluxlat, 
@@ -612 +619 @@
       REAL radsol(klon)       ! ray. net au sol (Solaire+IR) W/m2
       REAL albedo(klon)       ! albedo de la surface
+      REAL alblw(klon)
       REAL snow(klon)         ! hauteur de neige
       REAL qsol(klon)         ! humidite de la surface
@@ -825 +833 @@
      e ocean, npas, nexca, zmasq,
      s evap, fluxsens, fluxlat, dflux_l, dflux_s,              
-     s tsol_rad, tsurf_new, alb_new, emis_new, z0_new, 
+     s tsol_rad, tsurf_new, alb_new, alblw, emis_new, z0_new, 
      s pctsrf_new, agesno)
 
@@ -1080 +1088 @@
 c contre-gradient pour la chaleur sensible: Kelvin/metre
       REAL gamt(2:klev)
+c essai qsurf
+      real qsurf(klon)
+      real friv, frih
 c
       LOGICAL appel1er
@@ -1116 +1127 @@
          itop(i) = 0
       ENDDO
+
+      do i = 1, knon
+        qsurf(i) = qsatl(ts(i))/paprs(i,1)
+      enddo
+
 c
 c Prescrire la valeur de contre-gradient
@@ -1149 +1165 @@
          zdu2=max(cepdu2,u(i,1)**2+v(i,1)**2)
          zdphi=zgeop(i,1)
-         ztsolv = ts(i) * (1.0+RETV*q(i,1)) ! qsol approx = q(i,1)
+c         ztsolv = ts(i) * (1.0+RETV*q(i,1)) ! qsol approx = q(i,1)
+         ztsolv = ts(i) * (1.0+RETV*qsurf(i))
          ztvd=(t(i,1)+zdphi/RCPD/(1.+RVTMP2*q(i,1)))
      .       *(1.+RETV*q(i,1))
@@ -1157 +1174 @@
            IF (.NOT.zxli) THEN
            zscf=SQRT(1.+cd*ABS(zri(i)))
-           zcfm1(i) = zcdn/(1.+2.0*cb*zri(i)/zscf)
-           zcfh1(i) = zcdn/(1.+3.0*cb*zri(i)*zscf)
+           FRIV = AMAX1(1. / (1.+2.*CB*zri(i)/ZSCF), 0.1)
+!           zcfm1(i) = zcdn/(1.+2.0*cb*zri(i)/ zscf)
+           zcfm1(i) = zcdn * FRIV
+           FRIH = AMAX1(1./ (1.+3.*CB*zri(i)*ZSCF), 0.1 )
+!           zcfh1(i) = zcdn/(1.+3.0*cb*zri(i)*zscf)
+           zcfh1(i) = zcdn * FRIH
            pcfm(i,1) = zcfm1(i)
            pcfh(i,1) = zcfh1(i)
@@ -1169 +1190 @@
            zucf=1./(1.+3.0*cb*cc*zcdn*SQRT(ABS(zri(i))
      .              *(1.0+zgeop(i,1)/(RG*rugos(i)))))
-           zcfm2(i) = zcdn*(1.-2.0*cb*zri(i)*zucf)
-           zcfh2(i) = zcdn*(1.-3.0*cb*zri(i)*zucf)
+           zcfm2(i) = zcdn*amax1((1.-2.0*cb*zri(i)*zucf),0.1)
+           zcfh2(i) = zcdn*amax1((1.-3.0*cb*zri(i)*zucf),0.1)
            pcfm(i,1) = zcfm2(i)
            pcfh(i,1) = zcfh2(i)