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

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

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

trunk/phylmd/Interface_surf/interfsur_lim.f revision 104 by guez, Thu Sep 4 10:05:52 2014 UTC trunk/Sources/phylmd/Interface_surf/interfsur_lim.f revision 155 by guez, Wed Jul 8 17:03:45 2015 UTC
# Line 4  module interfsur_lim_m Line 4  module interfsur_lim_m
4    
5  contains  contains
6    
7    SUBROUTINE interfsur_lim(itime, dtime, jour, nisurf, knon, knindex, debut, &    SUBROUTINE interfsur_lim(itime, dtime, jour, knindex, debut, alblw, z0_new)
        alb_new, z0_new)  
8    
9      ! Cette routine sert d'interface entre le modèle atmosphérique et      ! Cette routine sert d'interface entre le mod\`ele atmosph\'erique et
10      ! un fichier de conditions aux limites.      ! un fichier de conditions aux limites.
11    
12      ! Laurent FAIRHEAD, February 2000      ! Laurent FAIRHEAD, February 2000
# Line 21  contains Line 20  contains
20      real, intent(IN):: dtime ! pas de temps de la physique (en s)      real, intent(IN):: dtime ! pas de temps de la physique (en s)
21      integer, intent(IN):: jour ! jour a lire dans l'annee      integer, intent(IN):: jour ! jour a lire dans l'annee
22    
23      integer, intent(IN):: nisurf      integer, intent(in):: knindex(:) ! (knon)
24      ! index de la surface à traiter (1 = sol continental)      ! index des points de la surface \`a traiter
25    
26      integer, intent(IN):: knon ! nombre de points dans le domaine a traiter      logical, intent(IN):: debut ! premier appel \`a la physique (initialisation)
27        real, intent(out):: alblw(:) ! (knon) albedo lu
28      integer, intent(in):: knindex(:) ! (klon)      real, intent(out):: z0_new(:) ! (klon) longueur de rugosit\'e lue
     ! index des points de la surface à traiter  
   
     logical, intent(IN):: debut ! premier appel à la physique (initialisation)  
     real, intent(out):: alb_new(:) ! (klon) albedo lu  
     real, intent(out):: z0_new(:) ! (klon) longueur de rugosité lue  
29    
30      ! Local:      ! Local:
31    
32        integer knon ! nombre de points dans le domaine a traiter
33    
34      integer, save:: lmt_pas ! frequence de lecture des conditions limites      integer, save:: lmt_pas ! frequence de lecture des conditions limites
35      ! (en pas de physique)      ! (en pas de physique)
36    
37      logical, save:: deja_lu_sur      logical, save:: deja_lu_sur
38      ! jour à lire déjà lu pour une surface précédente      ! jour \`a lire d\'ej\`a lu pour une surface pr\'ec\'edente
39    
40      integer, save:: jour_lu_sur      integer, save:: jour_lu_sur
41    
# Line 50  contains Line 46  contains
46    
47      !------------------------------------------------------------      !------------------------------------------------------------
48    
49        knon = size(knindex)
50    
51      if (debut) then      if (debut) then
52         lmt_pas = nint(86400./dtime * 1.0) ! pour une lecture une fois par jour         lmt_pas = nint(86400./dtime * 1.0) ! pour une lecture une fois par jour
53         jour_lu_sur = jour - 1         jour_lu_sur = jour - 1
# Line 67  contains Line 65  contains
65         call NF95_INQ_VARID(ncid, 'ALB', varid)         call NF95_INQ_VARID(ncid, 'ALB', varid)
66         call NF95_GET_VAR(ncid, varid, alb_lu, start=(/1, jour/))         call NF95_GET_VAR(ncid, varid, alb_lu, start=(/1, jour/))
67    
68         ! Lecture rugosité         ! Lecture rugosit\'e
69         call NF95_INQ_VARID(ncid, 'RUG', varid)         call NF95_INQ_VARID(ncid, 'RUG', varid)
70         call NF95_GET_VAR(ncid, varid, rug_lu, start=(/1, jour/))         call NF95_GET_VAR(ncid, varid, rug_lu, start=(/1, jour/))
71    
# Line 77  contains Line 75  contains
75      endif      endif
76    
77      ! Recopie des variables dans les champs de sortie      ! Recopie des variables dans les champs de sortie
78      alb_new(:knon) = alb_lu(knindex(:knon))      alblw = alb_lu(knindex)
79      z0_new(:knon) = rug_lu(knindex(:knon))      z0_new(:knon) = rug_lu(knindex)
     alb_new(knon + 1:) = 999999.  
80      z0_new(knon + 1:) = 999999.      z0_new(knon + 1:) = 999999.
81    
82    END SUBROUTINE interfsur_lim    END SUBROUTINE interfsur_lim
Legend:
Removed from v.104  
changed lines
  Added in v.155
  ViewVC Help
Powered by ViewVC 1.1.21