--- trunk/Sources/phylmd/phytrac.f	2015/07/21 15:29:52	159
+++ trunk/Sources/phylmd/phytrac.f	2015/09/29 19:48:59	171
@@ -38,6 +38,7 @@
     use indicesol, only: nbsrf
     use initrrnpb_m, only: initrrnpb
     use minmaxqfi_m, only: minmaxqfi
+    use netcdf, only: NF90_FILL_float
     use netcdf95, only: nf95_inq_varid, nf95_get_var, nf95_put_var
     use nflxtr_m, only: nflxtr
     use nr_util, only: assert
@@ -194,6 +195,8 @@
 
        call nf95_inq_varid(ncid_startphy, "trs", varid)
        call nf95_get_var(ncid_startphy, varid, trs(:, 1))
+       if (any(trs(:, 1) == NF90_FILL_float)) call abort_gcm("phytrac", &
+            "some missing values in trs(:, 1)")
 
        ! Initialisation de la fraction d'aerosols lessivee
 
@@ -339,7 +342,7 @@
        ! Ozone as a tracer:
        if (mod(itap - 1, lmt_pas) == 0) then
           ! Once per day, update the coefficients for ozone chemistry:
-          call regr_pr_comb_coefoz(julien)
+          call regr_pr_comb_coefoz(julien, paprs, pplay)
        end if
        call o3_chem(julien, gmtime, t_seri, zmasse, pdtphys, tr_seri(:, :, 3))
     end if