/[lmdze]/trunk/Sources/dyn3d/Guide/guide.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/Guide/guide.f

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

revision 103 by guez, Fri Aug 29 13:00:05 2014 UTC revision 107 by guez, Thu Sep 11 15:09:15 2014 UTC
# Line 24  CONTAINS Line 24  CONTAINS
24      use dump2d_m, only: dump2d      use dump2d_m, only: dump2d
25      USE exner_hyb_m, ONLY: exner_hyb      USE exner_hyb_m, ONLY: exner_hyb
26      USE inigrads_m, ONLY: inigrads      USE inigrads_m, ONLY: inigrads
27      use netcdf, only: nf90_nowrite, nf90_close, nf90_inq_dimid      use netcdf, only: nf90_nowrite
28      use netcdf95, only: nf95_inquire_dimension, nf95_open      use netcdf95, only: nf95_close, nf95_inq_dimid, nf95_inquire_dimension, &
29             nf95_open
30      use nr_util, only: pi      use nr_util, only: pi
31      USE paramet_m, ONLY: iip1, ip1jmp1, jjp1, llmp1      USE paramet_m, ONLY: iip1, ip1jmp1, jjp1, llmp1
32      USE q_sat_m, ONLY: q_sat      USE q_sat_m, ONLY: q_sat
# Line 75  CONTAINS Line 76  CONTAINS
76      real ztau(iim + 1, jjm + 1)      real ztau(iim + 1, jjm + 1)
77    
78      INTEGER ij, l      INTEGER ij, l
79      INTEGER ncidpl, status      INTEGER ncidpl
80      INTEGER rcod, rid      INTEGER rid
81      REAL tau      REAL tau
82      INTEGER, SAVE:: nlev      INTEGER, SAVE:: nlev
83    
# Line 122  CONTAINS Line 123  CONTAINS
123            CALL dump2d(iip1, jjp1, alpha_u, 'COEFF U ')            CALL dump2d(iip1, jjp1, alpha_u, 'COEFF U ')
124            CALL dump2d(iip1, jjp1, alpha_t, 'COEFF T ')            CALL dump2d(iip1, jjp1, alpha_t, 'COEFF T ')
125         ELSE         ELSE
126            ! Cas ou on force exactement par les variables analysees            ! Cas où on force exactement par les variables analysées
127            alpha_t = 0.            alpha_t = 0.
128            alpha_u = 0.            alpha_u = 0.
129            alpha_v = 0.            alpha_v = 0.
# Line 140  CONTAINS Line 141  CONTAINS
141         if (guide_Q) call nf95_open('hur.nc',nf90_nowrite, ncidpl)         if (guide_Q) call nf95_open('hur.nc',nf90_nowrite, ncidpl)
142    
143         IF (ncep) THEN         IF (ncep) THEN
144            status = nf90_inq_dimid(ncidpl, 'LEVEL', rid)            call nf95_inq_dimid(ncidpl, 'LEVEL', rid)
145         ELSE         ELSE
146            status = nf90_inq_dimid(ncidpl, 'PRESSURE', rid)            call nf95_inq_dimid(ncidpl, 'PRESSURE', rid)
147         END IF         END IF
148         call nf95_inquire_dimension(ncidpl, rid, nclen=nlev)         call nf95_inquire_dimension(ncidpl, rid, nclen=nlev)
149         PRINT *, 'nlev', nlev         PRINT *, 'nlev', nlev
150         rcod = nf90_close(ncidpl)         call nf95_close(ncidpl)
151         ! Lecture du premier etat des reanalyses.         ! Lecture du premier etat des reanalyses.
152         CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, &         CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, &
153              masserea2, nlev)              masserea2, nlev)
Legend:
Removed from v.103  
changed lines
  Added in v.107
  ViewVC Help
Powered by ViewVC 1.1.21