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

Annotation of /trunk/Sources/phylmd/Interface_surf/interfsur_lim.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 178 - (hide annotations)
Fri Mar 11 18:47:26 2016 UTC (8 years, 3 months ago) by guez
File size: 2479 byte(s) 
Moved variables date0, deltat, datasz_max, ncvar_ids, point, buff_pos,
buffer, regular from module histcom_var to modules where they are
defined.

Removed procedure ioipslmpp, useless for a sequential program.

Added argument datasz_max to histwrite_real (to avoid circular
dependency with histwrite).

Removed useless variables and computations everywhere.

Changed real litteral constants from default kind to double precision
in lwb, lwu, lwvn, sw1s, swtt, swtt1, swu.

Removed unused arguments: paer of sw, sw1s, sw2s, swclr; pcldsw of
sw1s, sw2s; pdsig, prayl of swr; co2_ppm of clmain, clqh; tsol of
transp_lay; nsrf of screenp; kcrit and kknu of gwstress; pstd of
orosetup.

Added output of relative humidity.
1 guez 54 module interfsur_lim_m
2    
3     implicit none
4    
5     contains
6    
7 guez 156 SUBROUTINE interfsur_lim(itime, dtime, jour, knindex, debut, albedo, z0_new)
8 guez 54
9 guez 150 ! Cette routine sert d'interface entre le mod\`ele atmosph\'erique et
10 guez 54 ! un fichier de conditions aux limites.
11    
12 guez 104 ! Laurent FAIRHEAD, February 2000
13 guez 54
14 guez 98 USE dimphy, ONLY: klon
15     use netcdf, only: NF90_NOWRITE
16     use netcdf95, only: NF95_close, NF95_GET_VAR, NF95_INQ_VARID, NF95_OPEN
17 guez 54
18 guez 98 integer, intent(IN):: itime ! numero du pas de temps courant
19     real, intent(IN):: dtime ! pas de temps de la physique (en s)
20     integer, intent(IN):: jour ! jour a lire dans l'annee
21 guez 54
22 guez 106 integer, intent(in):: knindex(:) ! (knon)
23 guez 150 ! index des points de la surface \`a traiter
24 guez 98
25 guez 150 logical, intent(IN):: debut ! premier appel \`a la physique (initialisation)
26 guez 156 real, intent(out):: albedo(:) ! (knon) albedo lu
27 guez 150 real, intent(out):: z0_new(:) ! (klon) longueur de rugosit\'e lue
28 guez 98
29     ! Local:
30    
31 guez 106 integer knon ! nombre de points dans le domaine a traiter
32    
33 guez 98 integer, save:: lmt_pas ! frequence de lecture des conditions limites
34 guez 54 ! (en pas de physique)
35    
36 guez 98 logical, save:: deja_lu_sur
37 guez 150 ! jour \`a lire d\'ej\`a lu pour une surface pr\'ec\'edente
38 guez 54
39 guez 98 integer, save:: jour_lu_sur
40 guez 54
41 guez 98 ! Champs lus dans le fichier de conditions aux limites :
42     real, allocatable, save:: alb_lu(:), rug_lu(:)
43    
44     integer ncid, varid
45    
46 guez 54 !------------------------------------------------------------
47    
48 guez 106 knon = size(knindex)
49    
50 guez 54 if (debut) then
51     lmt_pas = nint(86400./dtime * 1.0) ! pour une lecture une fois par jour
52     jour_lu_sur = jour - 1
53     allocate(alb_lu(klon))
54     allocate(rug_lu(klon))
55     endif
56    
57 guez 98 if (jour - jour_lu_sur /= 0) deja_lu_sur = .false.
58 guez 54
59     ! Tester d'abord si c'est le moment de lire le fichier
60 guez 98 if (mod(itime - 1, lmt_pas) == 0 .and. .not. deja_lu_sur) then
61     call NF95_OPEN('limit.nc', NF90_NOWRITE, ncid)
62 guez 54
63     ! Lecture Albedo
64 guez 98 call NF95_INQ_VARID(ncid, 'ALB', varid)
65     call NF95_GET_VAR(ncid, varid, alb_lu, start=(/1, jour/))
66 guez 54
67 guez 150 ! Lecture rugosit\'e
68 guez 98 call NF95_INQ_VARID(ncid, 'RUG', varid)
69     call NF95_GET_VAR(ncid, varid, rug_lu, start=(/1, jour/))
70 guez 54
71 guez 98 call NF95_CLOSE(ncid)
72 guez 54 deja_lu_sur = .true.
73     jour_lu_sur = jour
74     endif
75    
76     ! Recopie des variables dans les champs de sortie
77 guez 156 albedo = alb_lu(knindex)
78 guez 155 z0_new(:knon) = rug_lu(knindex)
79 guez 98 z0_new(knon + 1:) = 999999.
80 guez 54
81     END SUBROUTINE interfsur_lim
82    
83     end module interfsur_lim_m
  ViewVC Help
Powered by ViewVC 1.1.21