38 |
use indicesol, only: nbsrf |
use indicesol, only: nbsrf |
39 |
use initrrnpb_m, only: initrrnpb |
use initrrnpb_m, only: initrrnpb |
40 |
use minmaxqfi_m, only: minmaxqfi |
use minmaxqfi_m, only: minmaxqfi |
41 |
|
use netcdf, only: NF90_FILL_float |
42 |
use netcdf95, only: nf95_inq_varid, nf95_get_var, nf95_put_var |
use netcdf95, only: nf95_inq_varid, nf95_get_var, nf95_put_var |
43 |
use nflxtr_m, only: nflxtr |
use nflxtr_m, only: nflxtr |
44 |
use nr_util, only: assert |
use nr_util, only: assert |
195 |
|
|
196 |
call nf95_inq_varid(ncid_startphy, "trs", varid) |
call nf95_inq_varid(ncid_startphy, "trs", varid) |
197 |
call nf95_get_var(ncid_startphy, varid, trs(:, 1)) |
call nf95_get_var(ncid_startphy, varid, trs(:, 1)) |
198 |
|
if (any(trs(:, 1) == NF90_FILL_float)) call abort_gcm("phytrac", & |
199 |
|
"some missing values in trs(:, 1)") |
200 |
|
|
201 |
! Initialisation de la fraction d'aerosols lessivee |
! Initialisation de la fraction d'aerosols lessivee |
202 |
|
|
342 |
! Ozone as a tracer: |
! Ozone as a tracer: |
343 |
if (mod(itap - 1, lmt_pas) == 0) then |
if (mod(itap - 1, lmt_pas) == 0) then |
344 |
! Once per day, update the coefficients for ozone chemistry: |
! Once per day, update the coefficients for ozone chemistry: |
345 |
call regr_pr_comb_coefoz(julien) |
call regr_pr_comb_coefoz(julien, paprs, pplay) |
346 |
end if |
end if |
347 |
call o3_chem(julien, gmtime, t_seri, zmasse, pdtphys, tr_seri(:, :, 3)) |
call o3_chem(julien, gmtime, t_seri, zmasse, pdtphys, tr_seri(:, :, 3)) |
348 |
end if |
end if |