| 39 |
use dimens_m, only: jjm |
use dimens_m, only: jjm |
| 40 |
use dynetat0_m, only: rlatv |
use dynetat0_m, only: rlatv |
| 41 |
use nr_util, only: pi |
use nr_util, only: pi |
| 42 |
use numer_rec_95, only: regr1_step_av, regr3_lint |
use numer_rec_95, only: regr3_lint, regr1_conserv, slopes |
| 43 |
use netcdf95, only: nf95_open, nf95_gw_var, nf95_close, & |
use netcdf95, only: nf95_open, nf95_gw_var, nf95_close, & |
| 44 |
nf95_inq_varid, handle_err, nf95_put_var |
nf95_inq_varid, handle_err, nf95_put_var |
| 45 |
use netcdf, only: nf90_nowrite, nf90_get_var |
use netcdf, only: nf90_nowrite, nf90_get_var |
| 204 |
! Regrid in latitude: |
! Regrid in latitude: |
| 205 |
! We average with respect to sine of latitude, which is |
! We average with respect to sine of latitude, which is |
| 206 |
! equivalent to weighting by cosine of latitude: |
! equivalent to weighting by cosine of latitude: |
| 207 |
v_regr_lat(jjm+1:1:-1, :, 1:12) = regr1_step_av(o3_par_in, & |
call regr1_conserv(o3_par_in, & |
| 208 |
xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/))) |
xs = sin(lat_in_edg), xt = (/- 1., sin(rlatv(jjm:1:-1)), 1./), & |
| 209 |
|
slope = slopes(o3_par_in, sin(lat_in_edg)), & |
| 210 |
|
vt = v_regr_lat(jjm + 1:1:- 1, :, 1:12)) |
| 211 |
! (invert order of indices in "v_regr_lat" because "rlatu" is |
! (invert order of indices in "v_regr_lat" because "rlatu" is |
| 212 |
! decreasing) |
! decreasing) |
| 213 |
|
|
Parent Directory
| 
Revision Log
| 
Patch