/[lmdze]/trunk/phylmd/pbl_surface.f
ViewVC logotype

Diff of /trunk/phylmd/pbl_surface.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 282 by guez, Fri Jul 20 16:46:48 2018 UTC revision 283 by guez, Fri Jul 20 17:08:44 2018 UTC
# Line 217  contains Line 217  contains
217      yv = 0.      yv = 0.
218      yt = 0.      yt = 0.
219      yq = 0.      yq = 0.
     y_dflux_t = 0.  
     y_dflux_q = 0.  
220      yrugoro = 0.      yrugoro = 0.
221      d_ts = 0.      d_ts = 0.
222      flux_t = 0.      flux_t = 0.
# Line 353  contains Line 351  contains
351                 ytsoil(:knon, :), yqsol(:knon), mu0, yrugos(:knon), &                 ytsoil(:knon, :), yqsol(:knon), mu0, yrugos(:knon), &
352                 yrugoro(:knon), yu(:knon, 1), yv(:knon, 1), ycoefh(:knon, :), &                 yrugoro(:knon), yu(:knon, 1), yv(:knon, 1), ycoefh(:knon, :), &
353                 ycdragh(:knon), yt(:knon, :), yq(:knon, :), yts(:knon), &                 ycdragh(:knon), yt(:knon, :), yq(:knon, :), yts(:knon), &
354                 ypaprs(:knon, :), ypplay(:knon, :), ydelp, yrads(:knon), &                 ypaprs(:knon, :), ypplay(:knon, :), ydelp(:knon, :), &
355                 yalb(:knon), snow(:knon), yqsurf, yrain_f, ysnow_f, &                 yrads(:knon), yalb(:knon), snow(:knon), yqsurf(:knon), yrain_f, &
356                 yfluxlat(:knon), pctsrf_new_sic, yagesno(:knon), &                 ysnow_f, yfluxlat(:knon), pctsrf_new_sic, yagesno(:knon), &
357                 y_d_t(:knon, :), y_d_q(:knon, :), y_d_ts(:knon), &                 y_d_t(:knon, :), y_d_q(:knon, :), y_d_ts(:knon), &
358                 yz0_new(:knon), y_flux_t(:knon), y_flux_q(:knon), &                 yz0_new(:knon), y_flux_t(:knon), y_flux_q(:knon), &
359                 y_dflux_t(:knon), y_dflux_q(:knon), y_fqcalving(:knon), &                 y_dflux_t(:knon), y_dflux_q(:knon), y_fqcalving(:knon), &
360                 y_ffonte, y_run_off_lic_0)                 y_ffonte, y_run_off_lic_0)
361    
362            ! calculer la longueur de rugosite sur ocean            ! calculer la longueur de rugosite sur ocean
363    
364            yrugm = 0.            yrugm = 0.
365    
366            IF (nsrf == is_oce) THEN            IF (nsrf == is_oce) THEN
367               DO j = 1, knon               DO j = 1, knon
368                  yrugm(j) = 0.018 * ycdragm(j) * (yu(j, 1)**2 + yv(j, 1)**2) &                  yrugm(j) = 0.018 * ycdragm(j) * (yu(j, 1)**2 + yv(j, 1)**2) &
# Line 371  contains Line 371  contains
371                  yrugm(j) = max(1.5E-05, yrugm(j))                  yrugm(j) = max(1.5E-05, yrugm(j))
372               END DO               END DO
373            END IF            END IF
           DO j = 1, knon  
              y_dflux_t(j) = y_dflux_t(j) * ypct(j)  
              y_dflux_q(j) = y_dflux_q(j) * ypct(j)  
           END DO  
374    
375            DO k = 1, klev            DO k = 1, klev
376               DO j = 1, knon               DO j = 1, knon
# Line 414  contains Line 410  contains
410               ffonte(i, nsrf) = y_ffonte(j)               ffonte(i, nsrf) = y_ffonte(j)
411               cdragh(i) = cdragh(i) + ycdragh(j) * ypct(j)               cdragh(i) = cdragh(i) + ycdragh(j) * ypct(j)
412               cdragm(i) = cdragm(i) + ycdragm(j) * ypct(j)               cdragm(i) = cdragm(i) + ycdragm(j) * ypct(j)
413               dflux_t(i) = dflux_t(i) + y_dflux_t(j)               dflux_t(i) = dflux_t(i) + y_dflux_t(j) * ypct(j)
414               dflux_q(i) = dflux_q(i) + y_dflux_q(j)               dflux_q(i) = dflux_q(i) + y_dflux_q(j) * ypct(j)
415            END DO            END DO
416            IF (nsrf == is_ter) THEN            IF (nsrf == is_ter) THEN
417               qsol(ni(:knon)) = yqsol(:knon)               qsol(ni(:knon)) = yqsol(:knon)
Legend:
Removed from v.282  
changed lines
  Added in v.283
  ViewVC Help
Powered by ViewVC 1.1.21