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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 98 - (hide annotations)
Tue May 13 17:23:16 2014 UTC (10 years, 1 month ago) by guez
File size: 2642 byte(s) 
Split inter_barxy.f : one procedure per module, one module per
file. Grouped the files into a directory.

Split orbite.f.

Value of raz_date read from the namelist is taken into account
(resetting the step counter) even if annee_ref == anneeref and day_ref
== dayref. raz_date is no longer modified by gcm main unit. (Following
LMDZ.)

Removed argument klon of interfsur_lim. Renamed arguments lmt_alb,
lmt_rug to alb_new, z0_new (same name as corresponding actual
arguments in interfsurf_hq).

Removed argument klon of interfsurf_hq.

Removed arguments qs and d_qs of diagetpq. Were always
zero. Downgraded arguments d_qw, d_ql of diagetpq to local variables,
they were not used in physiq. Removed all computations for solid water
in diagetpq, was just zero.


Downgraded arguments fs_bound, fq_bound of diagphy to local variables,
they were not used in physiq. Encapsulated in a test on iprt all
computations in diagphy.

Removed parameter nbtr of module dimphy. Replaced it everywhere in the
program by nqmx - 2.

Removed parameter rnpb of procedure physiq. Kept the true case in
physiq and phytrac. Could not work with false case anyway.

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