1 |
module limit_mod |
module limit_mod |
2 |
|
|
|
! This module is clean: no C preprocessor directive, no include line. |
|
|
|
|
3 |
IMPLICIT none |
IMPLICIT none |
4 |
|
|
5 |
contains |
contains |
12 |
! It uses files with climatological data. |
! It uses files with climatological data. |
13 |
! Both grids must be regular. |
! Both grids must be regular. |
14 |
|
|
15 |
|
use comgeom, only: rlonu, rlatv |
16 |
|
use conf_dat2d_m, only: conf_dat2d |
17 |
use dimens_m, only: iim, jjm |
use dimens_m, only: iim, jjm |
|
use comconst, only: daysec, dtvr |
|
|
use indicesol, only: epsfra, nbsrf, is_ter, is_oce, is_lic, is_sic |
|
18 |
use dimphy, only: klon, zmasq |
use dimphy, only: klon, zmasq |
|
use conf_gcm_m, only: day_step |
|
|
use comgeom, only: rlonu, rlatv |
|
19 |
use etat0_mod, only: pctsrf |
use etat0_mod, only: pctsrf |
|
use start_init_orog_m, only: masque |
|
|
use conf_dat2d_m, only: conf_dat2d |
|
|
use inter_barxy_m, only: inter_barxy |
|
|
use interpolation, only: spline, splint |
|
20 |
use grid_change, only: dyn_phy |
use grid_change, only: dyn_phy |
21 |
|
use indicesol, only: epsfra, nbsrf, is_ter, is_oce, is_lic, is_sic |
22 |
use netcdf95, only: handle_err, nf95_get_coord, NF95_CLOSE, NF95_DEF_DIM, & |
use inter_barxy_m, only: inter_barxy |
23 |
|
use netcdf95, only: handle_err, nf95_gw_var, NF95_CLOSE, NF95_DEF_DIM, & |
24 |
nf95_enddef, NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, & |
nf95_enddef, NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, & |
25 |
nf95_inq_varid, NF95_OPEN |
nf95_inq_varid, NF95_OPEN |
26 |
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, & |
27 |
NF90_GLOBAL, NF90_NOWRITE, NF90_PUT_ATT, NF90_PUT_VAR, & |
NF90_GLOBAL, NF90_NOWRITE, NF90_PUT_ATT, NF90_PUT_VAR, & |
28 |
NF90_UNLIMITED |
NF90_UNLIMITED |
29 |
|
use numer_rec, only: spline, splint |
30 |
|
use start_init_orog_m, only: mask |
31 |
|
use unit_nml_m, only: unit_nml |
32 |
|
|
33 |
! Variables local to the procedure: |
! Variables local to the procedure: |
34 |
|
|
80 |
|
|
81 |
print *, "Enter namelist 'limit_nml'." |
print *, "Enter namelist 'limit_nml'." |
82 |
read (unit=*, nml=limit_nml) |
read (unit=*, nml=limit_nml) |
83 |
write(unit=*, nml=limit_nml) |
write(unit_nml, nml=limit_nml) |
|
|
|
|
! Initializations: |
|
|
dtvr = daysec / real(day_step) |
|
|
CALL inigeom |
|
84 |
|
|
85 |
! Process rugosity: |
! Process rugosity: |
86 |
|
|
87 |
PRINT *, 'Processing rugosity...' |
PRINT *, 'Processing rugosity...' |
88 |
call NF95_OPEN('Rugos.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('Rugos.nc', NF90_NOWRITE, ncid) |
89 |
|
|
90 |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
! Read coordinate variables: |
91 |
|
|
92 |
|
call nf95_inq_varid(ncid, "longitude", varid) |
93 |
|
call nf95_gw_var(ncid, varid, dlon_ini) |
94 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
95 |
|
|
96 |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
call nf95_inq_varid(ncid, "latitude", varid) |
97 |
|
call nf95_gw_var(ncid, varid, dlat_ini) |
98 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
99 |
|
|
100 |
call nf95_get_coord(ncid, "temps", timeyear) |
call nf95_inq_varid(ncid, "temps", varid) |
101 |
|
call nf95_gw_var(ncid, varid, timeyear) |
102 |
lmdep = size(timeyear) |
lmdep = size(timeyear) |
103 |
|
|
104 |
ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |
ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |
105 |
allocate(dlon(imdep), dlat(jmdep)) |
allocate(dlon(imdep), dlat(jmdep)) |
106 |
call NF95_INQ_VARID(ncid, 'RUGOS', varid) |
call NF95_INQ_VARID(ncid, 'RUGOS', varid) |
107 |
|
|
108 |
! Compute "champtime": |
! Read the primary variable day by day and regrid horizontally, |
109 |
|
! result in "champtime": |
110 |
DO l = 1, lmdep |
DO l = 1, lmdep |
111 |
ierr = NF90_GET_VAR(ncid, varid, champ, start=(/1, 1, l/)) |
ierr = NF90_GET_VAR(ncid, varid, champ, start=(/1, 1, l/)) |
112 |
call handle_err("NF90_GET_VAR", ierr) |
call handle_err("NF90_GET_VAR", ierr) |
115 |
CALL inter_barxy(dlon, dlat(:jmdep -1), LOG(champ), rlonu(:iim), & |
CALL inter_barxy(dlon, dlat(:jmdep -1), LOG(champ), rlonu(:iim), & |
116 |
rlatv, champtime(:, :, l)) |
rlatv, champtime(:, :, l)) |
117 |
champtime(:, :, l) = EXP(champtime(:, :, l)) |
champtime(:, :, l) = EXP(champtime(:, :, l)) |
118 |
where (nint(masque(:iim, :)) /= 1) champtime(:, :, l) = 0.001 |
where (nint(mask(:iim, :)) /= 1) champtime(:, :, l) = 0.001 |
119 |
end do |
end do |
120 |
|
|
121 |
call NF95_CLOSE(ncid) |
call NF95_CLOSE(ncid) |
123 |
DEALLOCATE(dlon, dlat, champ, dlon_ini, dlat_ini) |
DEALLOCATE(dlon, dlat, champ, dlon_ini, dlat_ini) |
124 |
allocate(yder(lmdep)) |
allocate(yder(lmdep)) |
125 |
|
|
126 |
|
! Interpolate monthly values to daily values, at each horizontal position: |
127 |
DO j = 1, jjm + 1 |
DO j = 1, jjm + 1 |
128 |
DO i = 1, iim |
DO i = 1, iim |
129 |
yder(:) = SPLINE(timeyear, champtime(i, j, :)) |
yder(:) = SPLINE(timeyear, champtime(i, j, :)) |
143 |
PRINT *, 'Processing sea ice...' |
PRINT *, 'Processing sea ice...' |
144 |
call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid) |
145 |
|
|
146 |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
call nf95_inq_varid(ncid, "longitude", varid) |
147 |
|
call nf95_gw_var(ncid, varid, dlon_ini) |
148 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
149 |
|
|
150 |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
call nf95_inq_varid(ncid, "latitude", varid) |
151 |
|
call nf95_gw_var(ncid, varid, dlat_ini) |
152 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
153 |
|
|
154 |
call nf95_inq_dimid(ncid, "time", dimid) |
call nf95_inq_dimid(ncid, "time", dimid) |
155 |
call NF95_INQuire_DIMension(ncid, dimid, len=lmdep) |
call NF95_INQuire_DIMension(ncid, dimid, nclen=lmdep) |
156 |
print *, 'lmdep = ', lmdep |
print *, 'lmdep = ', lmdep |
157 |
! PM 28/02/2002 : nouvelle coordonnée temporelle, fichiers AMIP |
! PM 28/02/2002 : nouvelle coordonnée temporelle, fichiers AMIP |
158 |
! pas en jours |
! pas en jours |
242 |
PRINT *, 'Traitement de la sst' |
PRINT *, 'Traitement de la sst' |
243 |
call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid) |
244 |
|
|
245 |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
call nf95_inq_varid(ncid, "longitude", varid) |
246 |
|
call nf95_gw_var(ncid, varid, dlon_ini) |
247 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
248 |
|
|
249 |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
call nf95_inq_varid(ncid, "latitude", varid) |
250 |
|
call nf95_gw_var(ncid, varid, dlat_ini) |
251 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
252 |
|
|
253 |
call nf95_inq_dimid(ncid, "time", dimid) |
call nf95_inq_dimid(ncid, "time", dimid) |
254 |
call NF95_INQuire_DIMension(ncid, dimid, len=lmdep) |
call NF95_INQuire_DIMension(ncid, dimid, nclen=lmdep) |
255 |
print *, 'lmdep = ', lmdep |
print *, 'lmdep = ', lmdep |
256 |
!PM28/02/2002 : nouvelle coord temporelle fichiers AMIP pas en jours |
!PM28/02/2002 : nouvelle coord temporelle fichiers AMIP pas en jours |
257 |
! Ici on suppose qu'on a 12 mois (de 30 jours). |
! Ici on suppose qu'on a 12 mois (de 30 jours). |
313 |
PRINT *, 'Traitement de l albedo' |
PRINT *, 'Traitement de l albedo' |
314 |
call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid) |
call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid) |
315 |
|
|
316 |
call nf95_get_coord(ncid, "longitude", dlon_ini) |
call nf95_inq_varid(ncid, "longitude", varid) |
317 |
|
call nf95_gw_var(ncid, varid, dlon_ini) |
318 |
imdep = size(dlon_ini) |
imdep = size(dlon_ini) |
319 |
|
|
320 |
call nf95_get_coord(ncid, "latitude", dlat_ini) |
call nf95_inq_varid(ncid, "latitude", varid) |
321 |
|
call nf95_gw_var(ncid, varid, dlat_ini) |
322 |
jmdep = size(dlat_ini) |
jmdep = size(dlat_ini) |
323 |
|
|
324 |
call nf95_get_coord(ncid, "temps", timeyear) |
call nf95_inq_varid(ncid, "temps", varid) |
325 |
|
call nf95_gw_var(ncid, varid, timeyear) |
326 |
lmdep = size(timeyear) |
lmdep = size(timeyear) |
327 |
|
|
328 |
ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |
ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep)) |