27 |
use interpolation, only: spline, splint |
use interpolation, only: spline, splint |
28 |
use grid_change, only: dyn_phy |
use grid_change, only: dyn_phy |
29 |
|
|
30 |
use netcdf95, only: handle_err, coordin, & |
use netcdf95, only: handle_err, nf95_get_coord, NF95_CLOSE, NF95_DEF_DIM, & |
31 |
NF90_CLOBBER, NF95_CLOSE, NF95_DEF_DIM, nf90_def_var, nf95_enddef, & |
nf95_enddef, NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, & |
32 |
NF90_FLOAT, NF90_GET_VAR, NF90_GLOBAL, NF90_NOWRITE, NF90_PUT_ATT, & |
nf95_inq_varid, NF95_OPEN |
33 |
NF90_PUT_VAR, NF90_UNLIMITED, & |
use netcdf, only: NF90_CLOBBER, nf90_def_var, NF90_FLOAT, NF90_GET_VAR, & |
34 |
NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, nf95_inq_varid, & |
NF90_GLOBAL, NF90_NOWRITE, NF90_PUT_ATT, NF90_PUT_VAR, & |
35 |
nf95_open |
NF90_UNLIMITED |
36 |
|
|
37 |
! Variables local to the procedure: |
! Variables local to the procedure: |
38 |
|
|
95 |
PRINT *, 'Processing rugosity...' |
PRINT *, 'Processing rugosity...' |
96 |
call NF95_OPEN('Rugos.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('Rugos.nc', NF90_NOWRITE, ncid) |
97 |
|
|
98 |
dlon_ini => coordin(ncid, "longitude") |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
99 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
100 |
|
|
101 |
dlat_ini => coordin(ncid, "latitude") |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
102 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
103 |
|
|
104 |
timeyear => coordin(ncid, "temps") |
call nf95_get_coord(ncid, "temps", timeyear) |
105 |
lmdep = size(timeyear) |
lmdep = size(timeyear) |
106 |
|
|
107 |
ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |
ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |
144 |
PRINT *, 'Processing sea ice...' |
PRINT *, 'Processing sea ice...' |
145 |
call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid) |
146 |
|
|
147 |
dlon_ini => coordin(ncid, "longitude") |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
148 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
149 |
|
|
150 |
dlat_ini => coordin(ncid, "latitude") |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
151 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
152 |
|
|
153 |
call nf95_inq_dimid(ncid, "time", dimid) |
call nf95_inq_dimid(ncid, "time", dimid) |
241 |
PRINT *, 'Traitement de la sst' |
PRINT *, 'Traitement de la sst' |
242 |
call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid) |
243 |
|
|
244 |
dlon_ini => coordin(ncid, "longitude") |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
245 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
246 |
|
|
247 |
dlat_ini => coordin(ncid, "latitude") |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
248 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
249 |
|
|
250 |
call nf95_inq_dimid(ncid, "time", dimid) |
call nf95_inq_dimid(ncid, "time", dimid) |
310 |
PRINT *, 'Traitement de l albedo' |
PRINT *, 'Traitement de l albedo' |
311 |
call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid) |
312 |
|
|
313 |
dlon_ini => coordin(ncid, "longitude") |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
314 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
315 |
|
|
316 |
dlat_ini => coordin(ncid, "latitude") |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
317 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
318 |
|
|
319 |
timeyear => coordin(ncid, "temps") |
call nf95_get_coord(ncid, "temps", timeyear) |
320 |
lmdep = size(timeyear) |
lmdep = size(timeyear) |
321 |
|
|
322 |
ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |
ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |