Changeset 3893 for trunk/LMDZ.GENERIC/libf/phystd/radioactive_tracers.F90
- Timestamp:
- Aug 13, 2025, 7:53:49 PM (7 days ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/phystd/radioactive_tracers.F90
r3866 r3893 1 2 1 subroutine radioactive_tracers(ngrid,nlayer,nq,ptimestep,pq,zdqradio) 3 2 … … 44 43 if (half_life(iq)/10 .le. ptimestep) then ! sanity test 45 44 write(*,*)"timestep is more than half-life / 10 of tracer ", noms(iq) 46 call abort 45 call abort_physic("radioactive_tracers","timestep is more than half-life / 10 of tracer",1) 47 46 end if ! dt > 0.1 * half_life 48 47 zdqradio(1:ngrid,1:nlayer,iq) = zdqradio(1:ngrid,1:nlayer,iq) - LOG(2.)/half_life(iq)*pq(1:ngrid,1:nlayer,iq) ! <== decay
Note: See TracChangeset
for help on using the changeset viewer.