- Timestamp:
- Jan 31, 2014, 1:42:35 PM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/physiq.F90
r1955 r1956 1577 1577 ! 1578 1578 if (read_climoz >= 1) then 1579 ! 1580 ! 1581 ro3i = int((days_elapsed + jh_cur - jh_1jan) &1582 / ioget_year_len(year_cur)* 360.) + 11579 ! Ozone from a file 1580 ! Update required ozone index: 1581 ro3i = int((days_elapsed + jh_cur - jh_1jan) / ioget_year_len(year_cur) & 1582 * 360.) + 1 1583 1583 if (ro3i == 361) ro3i = 360 1584 ! 1585 ! 1584 ! (This should never occur, except perhaps because of roundup 1585 ! error. See documentation.) 1586 1586 if (ro3i /= co3i) then 1587 ! 1587 ! Update ozone field: 1588 1588 if (read_climoz == 1) then 1589 1589 call regr_pr_av(ncid_climoz, (/"tro3"/), julien=ro3i, & 1590 1590 press_in_edg=press_climoz, paprs=paprs, v3=wo) 1591 1591 else 1592 ! read_climoz == 2 1593 call regr_pr_av(ncid_climoz, & 1594 (/"tro3 ", "tro3_daylight"/), & 1595 julien=ro3i, press_in_edg=press_climoz, paprs=paprs, & 1596 v3=wo) 1592 ! read_climoz == 2 1593 call regr_pr_av(ncid_climoz, (/"tro3 ", "tro3_daylight"/), & 1594 julien=ro3i, press_in_edg=press_climoz, paprs=paprs, v3=wo) 1597 1595 end if 1598 ! 1599 ! 1600 forall (l = 1: read_climoz) wo(:, :, l) = wo(:, :, l) &1601 * rmo3 / rmd *zmasse / dobson_u / 1e31602 ! 1603 ! 1604 ! 1605 ! 1596 ! Convert from mole fraction of ozone to column density of ozone in a 1597 ! cell, in kDU: 1598 forall (l = 1: read_climoz) wo(:, :, l) = wo(:, :, l) * rmo3 / rmd & 1599 * zmasse / dobson_u / 1e3 1600 ! (By regridding ozone values for LMDZ only once every 360th of 1601 ! year, we have already neglected the variation of pressure in one 1602 ! 360th of year. So do not recompute "wo" at each time step even if 1603 ! "zmasse" changes a little.) 1606 1604 co3i = ro3i 1607 1605 end if 1608 1606 ELSEIF (MOD(itap-1,lmt_pas) == 0) THEN 1609 ! 1607 ! Once per day, update ozone from Royer: 1610 1608 1611 1609 IF (solarlong0<-999.) then
Note: See TracChangeset
for help on using the changeset viewer.