/[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

trunk/phylmd/Interface_surf/interfsurf_hq.f revision 116 by guez, Thu Dec 4 16:35:03 2014 UTC trunk/Sources/phylmd/Interface_surf/interfsurf_hq.f revision 150 by guez, Thu Jun 18 13:49:26 2015 UTC
# Line 12  contains Line 12  contains
12         tsurf_new, alb_new, alblw, z0_new, pctsrf_new, agesno, fqcalving, &         tsurf_new, alb_new, alblw, z0_new, pctsrf_new, agesno, fqcalving, &
13         ffonte, run_off_lic_0, flux_o, flux_g)         ffonte, run_off_lic_0, flux_o, flux_g)
14    
15      ! Cette routine sert d'aiguillage entre l'atmosphère et la surface      ! Cette routine sert d'aiguillage entre l'atmosph\`ere et la surface
16      ! en général (sols continentaux, océans, glaces) pour les flux de      ! en g\'en\'eral (sols continentaux, oc\'eans, glaces) pour les flux de
17      ! chaleur et d'humidité.      ! chaleur et d'humidit\'e.
18    
19      ! Laurent Fairhead, February 2000      ! Laurent Fairhead, February 2000
20    
21      USE abort_gcm_m, ONLY: abort_gcm      USE abort_gcm_m, ONLY: abort_gcm
22        use alboc_m, only: alboc
23      USE albsno_m, ONLY: albsno      USE albsno_m, ONLY: albsno
24      use calbeta_m, only: calbeta      use calbeta_m, only: calbeta
25      USE calcul_fluxs_m, ONLY: calcul_fluxs      USE calcul_fluxs_m, ONLY: calcul_fluxs
# Line 83  contains Line 84  contains
84      ! rugos rugosite      ! rugos rugosite
85      ! rugoro rugosite orographique      ! rugoro rugosite orographique
86      real, intent(INOUT):: snow(klon), qsurf(klon)      real, intent(INOUT):: snow(klon), qsurf(klon)
87      real, intent(IN):: tsurf(:) ! (knon) température de surface      real, intent(IN):: tsurf(:) ! (knon) temp\'erature de surface
88      real, dimension(klon), intent(IN):: p1lay      real, dimension(klon), intent(IN):: p1lay
89      ! p1lay pression 1er niveau (milieu de couche)      ! p1lay pression 1er niveau (milieu de couche)
90      real, dimension(klon), intent(IN):: ps      real, dimension(klon), intent(IN):: ps
# Line 95  contains Line 96  contains
96      ! fluxsens flux de chaleur sensible      ! fluxsens flux de chaleur sensible
97      ! fluxlat flux de chaleur latente      ! fluxlat flux de chaleur latente
98      real, dimension(klon), intent(OUT):: dflux_l, dflux_s      real, dimension(klon), intent(OUT):: dflux_l, dflux_s
99      real, intent(OUT):: tsurf_new(knon) ! température au sol      real, intent(OUT):: tsurf_new(knon) ! temp\'erature au sol
100      real, intent(OUT):: alb_new(klon) ! albedo      real, intent(OUT):: alb_new(klon) ! albedo
101      real, dimension(klon), intent(OUT):: alblw      real, dimension(klon), intent(OUT):: alblw
102      real, dimension(klon), intent(OUT):: z0_new      real, dimension(klon), intent(OUT):: z0_new
# Line 104  contains Line 105  contains
105      ! pctsrf_new nouvelle repartition des surfaces      ! pctsrf_new nouvelle repartition des surfaces
106      real, dimension(klon), intent(INOUT):: agesno      real, dimension(klon), intent(INOUT):: agesno
107    
108      ! Flux d'eau "perdue" par la surface et nécessaire pour que limiter la      ! Flux d'eau "perdue" par la surface et n\'ecessaire pour que limiter la
109      ! hauteur de neige, en kg/m2/s      ! hauteur de neige, en kg/m2/s
110      !jld a rajouter real, dimension(klon), intent(INOUT):: fqcalving      !jld a rajouter real, dimension(klon), intent(INOUT):: fqcalving
111      real, dimension(klon), intent(INOUT):: fqcalving      real, dimension(klon), intent(INOUT):: fqcalving
# Line 206  contains Line 207  contains
207    
208         ! calcul albedo: lecture albedo fichier boundary conditions         ! calcul albedo: lecture albedo fichier boundary conditions
209         ! puis ajout albedo neige         ! puis ajout albedo neige
210         call interfsur_lim(itime, dtime, jour, nisurf, knindex, debut, &         call interfsur_lim(itime, dtime, jour, knindex, debut, alb_new, z0_new)
             alb_new, z0_new)  
211    
212         ! calcul snow et qsurf, hydrol adapté         ! calcul snow et qsurf, hydrol adapt\'e
213         CALL calbeta(nisurf, snow(:knon), qsol(:knon), beta(:knon), &         CALL calbeta(nisurf, snow(:knon), qsol(:knon), beta(:knon), &
214              capsol(:knon), dif_grnd(:knon))              capsol(:knon), dif_grnd(:knon))
215    
# Line 245  contains Line 245  contains
245         ! Remplissage des pourcentages de surface         ! Remplissage des pourcentages de surface
246         pctsrf_new(:, nisurf) = pctsrf(:, nisurf)         pctsrf_new(:, nisurf) = pctsrf(:, nisurf)
247      case (is_oce)      case (is_oce)
248         ! Surface "ocean" appel à l'interface avec l'océan         ! Surface "ocean" appel \`a l'interface avec l'oc\'ean
249         ! lecture conditions limites         ! lecture conditions limites
250         call interfoce_lim(itime, dtime, jour, knindex, debut, tsurf_temp, &         call interfoce_lim(itime, dtime, jour, knindex, debut, tsurf_temp, &
251              pctsrf_new)              pctsrf_new)
# Line 275  contains Line 275  contains
275    
276         ! calcul albedo         ! calcul albedo
277         if (minval(rmu0) == maxval(rmu0) .and. minval(rmu0) == -999.999) then         if (minval(rmu0) == maxval(rmu0) .and. minval(rmu0) == -999.999) then
278            CALL alboc(FLOAT(jour), rlat, alb_eau)            CALL alboc(jour, rlat, alb_eau)
279         else ! cycle diurne         else ! cycle diurne
280            CALL alboc_cd(rmu0, alb_eau)            CALL alboc_cd(rmu0, alb_eau)
281         endif         endif
Legend:
Removed from v.116  
changed lines
  Added in v.150
  ViewVC Help
Powered by ViewVC 1.1.21