/[lmdze]/trunk/phylmd/Interface_surf/interfsurf_hq.f
ViewVC logotype

Diff of /trunk/phylmd/Interface_surf/interfsurf_hq.f

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

revision 323 by guez, Mon Dec 3 17:52:21 2018 UTC revision 324 by guez, Wed Feb 6 15:58:03 2019 UTC
# Line 131  contains Line 131  contains
131         cal = RCPD / soilcap         cal = RCPD / soilcap
132    
133         CALL calcul_fluxs(ts, p1lay, cal, beta, cdragh, ps, qsurf, &         CALL calcul_fluxs(ts, p1lay, cal, beta, cdragh, ps, qsurf, &
134              radsol + soilflux, temp_air, q1lay, u1lay, v1lay, tAcoef, &              radsol + soilflux, temp_air, q1lay, u1lay, v1lay, tAcoef, qAcoef, &
135              qAcoef, tBcoef, qBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, &              tBcoef, qBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, &
136              dflux_l, dif_grnd = 0.)              dflux_l, dif_grnd = 0.)
137         CALL fonte_neige(is_ter, rain_fall, snow_fall, snow, qsol, &         CALL fonte_neige(is_ter, rain_fall, snow_fall, snow, qsol, &
138              tsurf_new, evap, fqcalving, ffonte, run_off_lic_0, run_off_lic)              tsurf_new, evap, fqcalving, ffonte, run_off_lic_0, run_off_lic)
# Line 162  contains Line 162  contains
162         DO ii = 1, size(knindex)         DO ii = 1, size(knindex)
163            IF (pctsrf_new_sic(ii) < EPSFRA) then            IF (pctsrf_new_sic(ii) < EPSFRA) then
164               snow(ii) = 0.               snow(ii) = 0.
165               tsurf_new(ii) = RTT - 1.8               tsurf(ii) = RTT - 1.8
166               tsoil(ii, :) = RTT - 1.8               tsoil(ii, :) = RTT - 1.8
167            else            else
168               tsurf_new(ii) = ts(ii)               tsurf(ii) = ts(ii)
169            endif            endif
170         enddo         enddo
171    
172         CALL soil(is_sic, snow, tsurf_new, tsoil, soilcap, soilflux)         CALL soil(is_sic, snow, tsurf, tsoil, soilcap, soilflux)
173         cal = RCPD / soilcap         cal = RCPD / soilcap
        tsurf = tsurf_new  
174         beta = 1.         beta = 1.
175         CALL calcul_fluxs(tsurf, p1lay, cal, beta, cdragh, ps, qsurf, &         CALL calcul_fluxs(tsurf, p1lay, cal, beta, cdragh, ps, qsurf, &
176              radsol + soilflux, temp_air, q1lay, u1lay, v1lay, tAcoef, &              radsol + soilflux, temp_air, q1lay, u1lay, v1lay, tAcoef, qAcoef, &
177              qAcoef, tBcoef, qBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, &              tBcoef, qBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, &
178              dflux_l, dif_grnd = 1. / tau_gl)              dflux_l, dif_grnd = 1. / tau_gl)
179         CALL fonte_neige(is_sic, rain_fall, snow_fall, snow, qsol, &         CALL fonte_neige(is_sic, rain_fall, snow_fall, snow, qsol, &
180              tsurf_new, evap, fqcalving, ffonte, run_off_lic_0, run_off_lic)              tsurf_new, evap, fqcalving, ffonte, run_off_lic_0, run_off_lic)
# Line 189  contains Line 188  contains
188    
189         z0_new = SQRT(0.002**2 + rugoro**2)         z0_new = SQRT(0.002**2 + rugoro**2)
190      case (is_lic)      case (is_lic)
191         ! Surface "glacier continentaux" appel a l'interface avec le sol         ! Surface "glacier continentaux" appel \`a l'interface avec le sol
192    
193         CALL soil(is_lic, snow, ts, tsoil, soilcap, soilflux)         CALL soil(is_lic, snow, ts, tsoil, soilcap, soilflux)
194         cal = RCPD / soilcap         cal = RCPD / soilcap
195         beta = 1.         beta = 1.
196         call calcul_fluxs(ts, p1lay, cal, beta, cdragh, ps, qsurf, &         call calcul_fluxs(ts, p1lay, cal, beta, cdragh, ps, qsurf, &
197              radsol + soilflux, temp_air, q1lay, u1lay, v1lay, tAcoef, &              radsol + soilflux, temp_air, q1lay, u1lay, v1lay, tAcoef, qAcoef, &
198              qAcoef, tBcoef, qBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, &              tBcoef, qBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, &
199              dflux_l, dif_grnd = 0.)              dflux_l, dif_grnd = 0.)
200         call fonte_neige(is_lic, rain_fall, snow_fall, snow, qsol, &         call fonte_neige(is_lic, rain_fall, snow_fall, snow, qsol, &
201              tsurf_new, evap, fqcalving, ffonte, run_off_lic_0, run_off_lic)              tsurf_new, evap, fqcalving, ffonte, run_off_lic_0, run_off_lic)
Legend:
Removed from v.323  
changed lines
  Added in v.324
  ViewVC Help
Powered by ViewVC 1.1.21