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 |
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 |
|
|
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. |
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) |