- Timestamp:
- May 6, 2015, 12:14:12 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_lat_time_coefoz_m.F90
r3819 r3825 49 49 !use comconst_phy_mod, only : pi 50 50 use nrtype, only: pi 51 use comgeomphy, only: rlatv51 use regular_lonlat_mod, only: boundslat_reg, south 52 52 ! Variables local to the procedure: 53 53 … … 211 211 ! equivalent to weighting by cosine of latitude: 212 212 v_regr_lat(nbp_lat:1:-1, :, 1:12) = regr1_step_av(o3_par_in, & 213 xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(nbp_lat-1:1:-1), pi / 2/)))213 xs=sin(lat_in_edg), xt=sin((/- pi / 2, boundslat_reg(nbp_lat-1:1:-1,south), pi / 2/))) 214 214 ! (invert order of indices in "v_regr_lat" because "rlatu" is 215 215 ! in descending order) … … 251 251 !use comconst_phy_mod, only : pi 252 252 use nrtype, only: pi 253 use comgeomphy, only: rlatu 254 253 use regular_lonlat_mod, only : lat_reg 255 254 integer, intent(in):: ncid_in, varid_in(:), n_plev 256 255 character(len=*), intent(in):: name_out(:) ! of NetCDF variables … … 326 325 327 326 ! Write one of the coordinate variables: 328 call nf95_put_var(ncid_out, varid_rlatu, rlatu(nbp_lat:1:-1) / pi * 180.)327 call nf95_put_var(ncid_out, varid_rlatu, lat_reg(nbp_lat:1:-1) / pi * 180.) 329 328 ! (convert from rad to degrees and sort in ascending order) 330 329
Note: See TracChangeset
for help on using the changeset viewer.