27 |
use numer_rec, only: spline, splint |
use numer_rec, only: spline, splint |
28 |
use grid_change, only: dyn_phy |
use grid_change, only: dyn_phy |
29 |
|
|
30 |
use netcdf95, only: handle_err, nf95_get_coord, NF95_CLOSE, NF95_DEF_DIM, & |
use netcdf95, only: handle_err, nf95_gw_var, NF95_CLOSE, NF95_DEF_DIM, & |
31 |
nf95_enddef, NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, & |
nf95_enddef, NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, & |
32 |
nf95_inq_varid, NF95_OPEN |
nf95_inq_varid, NF95_OPEN |
33 |
use netcdf, only: NF90_CLOBBER, nf90_def_var, NF90_FLOAT, NF90_GET_VAR, & |
use netcdf, only: NF90_CLOBBER, nf90_def_var, NF90_FLOAT, NF90_GET_VAR, & |
97 |
|
|
98 |
! Read coordinate variables: |
! Read coordinate variables: |
99 |
|
|
100 |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
call nf95_inq_varid(ncid, "longitude", varid) |
101 |
|
call nf95_gw_var(ncid, varid, dlon_ini) |
102 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
103 |
|
|
104 |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
call nf95_inq_varid(ncid, "latitude", varid) |
105 |
|
call nf95_gw_var(ncid, varid, dlat_ini) |
106 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
107 |
|
|
108 |
call nf95_get_coord(ncid, "temps", timeyear) |
call nf95_inq_varid(ncid, "temps", varid) |
109 |
|
call nf95_gw_var(ncid, varid, timeyear) |
110 |
lmdep = size(timeyear) |
lmdep = size(timeyear) |
111 |
|
|
112 |
ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |
ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |
151 |
PRINT *, 'Processing sea ice...' |
PRINT *, 'Processing sea ice...' |
152 |
call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid) |
153 |
|
|
154 |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
call nf95_inq_varid(ncid, "longitude", varid) |
155 |
|
call nf95_gw_var(ncid, varid, dlon_ini) |
156 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
157 |
|
|
158 |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
call nf95_inq_varid(ncid, "latitude", varid) |
159 |
|
call nf95_gw_var(ncid, varid, dlat_ini) |
160 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
161 |
|
|
162 |
call nf95_inq_dimid(ncid, "time", dimid) |
call nf95_inq_dimid(ncid, "time", dimid) |
250 |
PRINT *, 'Traitement de la sst' |
PRINT *, 'Traitement de la sst' |
251 |
call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid) |
252 |
|
|
253 |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
call nf95_inq_varid(ncid, "longitude", varid) |
254 |
|
call nf95_gw_var(ncid, varid, dlon_ini) |
255 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
256 |
|
|
257 |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
call nf95_inq_varid(ncid, "latitude", varid) |
258 |
|
call nf95_gw_var(ncid, varid, dlat_ini) |
259 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
260 |
|
|
261 |
call nf95_inq_dimid(ncid, "time", dimid) |
call nf95_inq_dimid(ncid, "time", dimid) |
321 |
PRINT *, 'Traitement de l albedo' |
PRINT *, 'Traitement de l albedo' |
322 |
call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid) |
323 |
|
|
324 |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
call nf95_inq_varid(ncid, "longitude", varid) |
325 |
|
call nf95_gw_var(ncid, varid, dlon_ini) |
326 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
327 |
|
|
328 |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
call nf95_inq_varid(ncid, "latitude", varid) |
329 |
|
call nf95_gw_var(ncid, varid, dlat_ini) |
330 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
331 |
|
|
332 |
call nf95_get_coord(ncid, "temps", timeyear) |
call nf95_inq_varid(ncid, "temps", varid) |
333 |
|
call nf95_gw_var(ncid, varid, timeyear) |
334 |
lmdep = size(timeyear) |
lmdep = size(timeyear) |
335 |
|
|
336 |
ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |
ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |