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

Timestamp:

Jul 26, 2000, 2:58:36 PM (24 years ago)

Author:

lmdzadmin

Message:

Deplacement des dernieres lignes concernant le sol hors de la physique
LF

File:

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

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

@@ -1 +1 @@
       SUBROUTINE clmain(dtime,pctsrf,t,q,u,v,
+     .                  jour, rmu0,
      .                  ok_veget,ts,
      .                  paprs,pplay,radsol,snow,qsol,evap,albe,
@@ -85 +86 @@
       REAL rugmer(klon)
       REAL cdragh(klon), cdragm(klon)
+      integer jour            ! jour de l'annee en cours
+      real rmu0(klon)         ! cosinus de l'angle solaire zenithal
       LOGICAL debut, lafin, ok_veget
 cAA      INTEGER itr
@@ -112 +115 @@
       EXTERNAL clqh, clvent, coefkz, calbeta, cltrac
 c======================================================================
-      REAL yts(klon), yrugos(klon), ypct(klon)
+      REAL yts(klon), yrugos(klon), ypct(klon), yz0_new(klon)
       REAL ycal(klon), ybeta(klon), ydif(klon), yalb(klon),yevap(klon)
       REAL yu1(klon), yv1(klon)
@@ -329 +332 @@
 c calculer la diffusion de "q" et de "h"
       CALL clqh(knon, dtime, nsrf, ni, pctsrf, rlon, rlat,
-     e          yu1, yv1,
-     e          ycoefh,yt,yq,yts,ypaprs,ypplay,ydelp,yrads,
-     e          yevap,yalb, ysnow, yqsol, yrain_f, ysnow_f,
-     e          yfder, ytaux, ytauy, ysollw, ysolsw,
+     e          jour, rmu0,
+     e          yu1, yv1, ycoefh,
+     e          yt,yq,yts,ypaprs,ypplay,
+     e          ydelp,yrads, yevap,yalb, ysnow, yqsol, 
+     e          yrain_f, ysnow_f, yfder, ytaux, ytauy, 
+     e          ysollw, ysolsw,
      s          pctsrf_new,
-     s          y_d_t, y_d_q, y_d_ts, 
+     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)
 c
@@ -391 +396 @@
          snow(i,nsrf) = ysnow(j)
          qsol(i,nsrf) = yqsol(j)
+         rugos(i,nsrf) = yz0_new(j)
          rugmer(i) = yrugm(j)
          cdragh(i) = cdragh(i) + ycoefh(j,1)
@@ -445 +451 @@
       END
       SUBROUTINE clqh(knon,dtime,nisurf,knindex,pctsrf, rlon, rlat,
+     e                jour, rmu0
      e                u1lay,v1lay,coef,
      e                t,q,ts,paprs,pplay,
@@ -451 +458 @@
      e                lwdown, swdown, 
      s                pctsrf_new,
-     s                d_t, d_q, d_ts, flux_t, flux_q,dflux_s,dflux_l)
+     s                d_t, d_q, d_ts, z0_new, 
+     s                flux_t, flux_q,dflux_s,dflux_l)
 
       USE interface_surf
@@ -491 +499 @@
       REAL qsol(klon)         ! humidite de la surface
       real precip_rain(klon), precip_snow(klon)
+      integer jour            ! jour de l'annee en cours
+      real rmu0(klon)         ! cosinus de l'angle solaire zenithal
       integer knindex(klon)
       real pctsrf(klon,nbsrf)
@@ -554 +564 @@
 c Rajout pour l'interface
       integer itime
-      integer jour
       integer nisurf
       logical debut, lafin, ok_veget
@@ -699 +708 @@
       ccanopy = 0.
 
-      CALL interfsurf(itime, dtime, jour,
-     . klon, iim, jjm, nisurf, knon, knindex, pctsrf, rlon, rlat,
-     . debut, lafin, ok_veget,
-     . zlev1,  u1lay, v1lay, temp_air, spechum, hum_air, ccanopy, 
-     . tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef,
-     . precip_rain, precip_snow, lwdown, swnet, swdown,
-     . fder, taux, tauy,
-     . albedo, snow, qsol,
-     . ts, p1lay, psref, radsol,
-     . ocean,zmasq
-     . evap, fluxsens, fluxlat, dflux_l, dflux_s,              
-     . tsol_rad, tsurf_new, alb_new, emis_new, z0_new, pctsrf_new)
+      CALL interfsurf(itime, dtime, jour, rmu0,
+     e klon, iim, jjm, nisurf, knon, knindex, pctsrf, rlon, rlat,
+     e debut, lafin, ok_veget,
+     e zlev1,  u1lay, v1lay, temp_air, spechum, hum_air, ccanopy, 
+     e tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef,
+     e precip_rain, precip_snow, lwdown, swnet, swdown,
+     e fder, taux, tauy,
+     e albedo, snow, qsol,
+     e ts, p1lay, psref, radsol,
+     e ocean,zmasq,
+     s evap, fluxsens, fluxlat, dflux_l, dflux_s,              
+     s tsol_rad, tsurf_new, alb_new, emis_new, z0_new, pctsrf_new)
 
       flux_t(:,1) = fluxsens