/[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 223 by guez, Fri Apr 28 13:22:36 2017 UTC revision 225 by guez, Mon Oct 16 12:35:41 2017 UTC
# Line 45  contains Line 45  contains
45    
46      REAL, intent(inout):: tsoil(:, :) ! (knon, nsoilmx)      REAL, intent(inout):: tsoil(:, :) ! (knon, nsoilmx)
47    
48      REAL, intent(INOUT):: qsol(klon)      REAL, intent(INOUT):: qsol(:) ! (knon)
49      ! column-density of water in soil, in kg m-2      ! column-density of water in soil, in kg m-2
50    
51      real, dimension(klon), intent(IN):: u1_lay, v1_lay      real, intent(IN):: u1_lay(:), v1_lay(:) ! (knon) vitesse 1ere couche
52      ! u1_lay vitesse u 1ere couche  
     ! v1_lay vitesse v 1ere couche  
53      real, dimension(klon), intent(IN):: temp_air, spechum      real, dimension(klon), intent(IN):: temp_air, spechum
54      ! temp_air temperature de l'air 1ere couche      ! temp_air temperature de l'air 1ere couche
55      ! spechum humidite specifique 1ere couche      ! spechum humidite specifique 1ere couche
# Line 166  contains Line 165  contains
165         call interfsur_lim(dtime, julien, knindex, debut, albedo, z0_new)         call interfsur_lim(dtime, julien, knindex, debut, albedo, z0_new)
166    
167         ! Calcul de snow et qsurf, hydrologie adapt\'ee         ! Calcul de snow et qsurf, hydrologie adapt\'ee
168         CALL calbeta(is_ter, snow, qsol(:knon), beta, capsol(:knon), &         CALL calbeta(is_ter, snow, qsol, beta, capsol(:knon), dif_grnd(:knon))
             dif_grnd(:knon))  
169    
170         IF (soil_model) THEN         IF (soil_model) THEN
171            CALL soil(dtime, is_ter, snow, ts, tsoil, soilcap, soilflux)            CALL soil(dtime, is_ter, snow, ts, tsoil, soilcap, soilflux)
# Line 179  contains Line 177  contains
177    
178         CALL calcul_fluxs(dtime, ts, p1lay(:knon), cal, beta, tq_cdrag(:knon), &         CALL calcul_fluxs(dtime, ts, p1lay(:knon), cal, beta, tq_cdrag(:knon), &
179              ps(:knon), qsurf(:knon), radsol, dif_grnd(:knon), &              ps(:knon), qsurf(:knon), radsol, dif_grnd(:knon), &
180              temp_air(:knon), spechum(:knon), u1_lay(:knon), v1_lay(:knon), &              temp_air(:knon), spechum(:knon), u1_lay, v1_lay, &
181              petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &              petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &
182              peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)              peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)
183         CALL fonte_neige(is_ter, dtime, precip_rain(:knon), &         CALL fonte_neige(is_ter, dtime, precip_rain(:knon), &
184              precip_snow(:knon), snow, qsol(:knon), tsurf_new, evap, &              precip_snow(:knon), snow, qsol, tsurf_new, evap, &
185              fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))              fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))
186    
187         call albsno(dtime, agesno, alb_neig, precip_snow(:knon))         call albsno(dtime, agesno, alb_neig, precip_snow(:knon))
# Line 200  contains Line 198  contains
198         dif_grnd = 0.         dif_grnd = 0.
199         call calcul_fluxs(dtime, tsurf, p1lay(:knon), cal, beta, &         call calcul_fluxs(dtime, tsurf, p1lay(:knon), cal, beta, &
200              tq_cdrag(:knon), ps(:knon), qsurf(:knon), radsol, &              tq_cdrag(:knon), ps(:knon), qsurf(:knon), radsol, &
201              dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay(:knon), &              dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay, &
202              v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &              v1_lay, petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &
203              peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)              peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)
204         agesno = 0.         agesno = 0.
205         albedo = alboc_cd(rmu0(knindex)) * fmagic         albedo = alboc_cd(rmu0(knindex)) * fmagic
# Line 219  contains Line 217  contains
217            endif            endif
218         enddo         enddo
219    
220         CALL calbeta(is_sic, snow, qsol(:knon), beta, capsol(:knon), &         CALL calbeta(is_sic, snow, qsol, beta, capsol(:knon), dif_grnd(:knon))
             dif_grnd(:knon))  
221    
222         IF (soil_model) THEN         IF (soil_model) THEN
223            CALL soil(dtime, is_sic, snow, tsurf_new, tsoil, soilcap, &            CALL soil(dtime, is_sic, snow, tsurf_new, tsoil, soilcap, &
# Line 238  contains Line 235  contains
235    
236         CALL calcul_fluxs(dtime, tsurf, p1lay(:knon), cal, beta, &         CALL calcul_fluxs(dtime, tsurf, p1lay(:knon), cal, beta, &
237              tq_cdrag(:knon), ps(:knon), qsurf(:knon), radsol, &              tq_cdrag(:knon), ps(:knon), qsurf(:knon), radsol, &
238              dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay(:knon), &              dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay, &
239              v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &              v1_lay, petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &
240              peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)              peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)
241         CALL fonte_neige(is_sic, dtime, precip_rain(:knon), &         CALL fonte_neige(is_sic, dtime, precip_rain(:knon), &
242              precip_snow(:knon), snow, qsol(:knon), tsurf_new, evap, &              precip_snow(:knon), snow, qsol, tsurf_new, evap, &
243              fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))              fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))
244    
245         ! Compute the albedo:         ! Compute the albedo:
# Line 269  contains Line 266  contains
266    
267         call calcul_fluxs(dtime, ts, p1lay(:knon), cal, beta, tq_cdrag(:knon), &         call calcul_fluxs(dtime, ts, p1lay(:knon), cal, beta, tq_cdrag(:knon), &
268              ps(:knon), qsurf(:knon), radsol, dif_grnd(:knon), &              ps(:knon), qsurf(:knon), radsol, dif_grnd(:knon), &
269              temp_air(:knon), spechum(:knon), u1_lay(:knon), v1_lay(:knon), &              temp_air(:knon), spechum(:knon), u1_lay, v1_lay, &
270              petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &              petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &
271              peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)              peqBcoef(:knon), tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)
272         call fonte_neige(is_lic, dtime, precip_rain(:knon), &         call fonte_neige(is_lic, dtime, precip_rain(:knon), &
273              precip_snow(:knon), snow, qsol(:knon), tsurf_new, evap, &              precip_snow(:knon), snow, qsol, tsurf_new, evap, &
274              fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))              fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))
275    
276         ! calcul albedo         ! calcul albedo
Legend:
Removed from v.223  
changed lines
  Added in v.225
  ViewVC Help
Powered by ViewVC 1.1.21