/[lmdze]/trunk/Sources/phylmd/phytrac.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/phytrac.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 159 by guez, Tue Jul 21 15:29:52 2015 UTC revision 171 by guez, Tue Sep 29 19:48:59 2015 UTC
# Line 38  contains Line 38  contains
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
# Line 194  contains Line 195  contains
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    
# Line 339  contains Line 342  contains
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

Legend:
Removed from v.159  
changed lines
  Added in v.171

  ViewVC Help
Powered by ViewVC 1.1.21