34 |
use start_init_phys_m, only: start_init_phys |
use start_init_phys_m, only: start_init_phys |
35 |
use nr_util, only: assert, pi |
use nr_util, only: assert, pi |
36 |
|
|
37 |
REAL, intent(out):: tsol_2d(:, :) |
REAL, intent(in):: tsol_2d(:, :) ! (iim + 1, jjm + 1) |
38 |
REAL, intent(out):: psol(:, :) ! surface pressure, in Pa |
REAL, intent(out):: psol(:, :) ! (iim + 1, jjm + 1) surface pressure, in Pa |
39 |
|
|
40 |
! Local: |
! Local: |
41 |
|
|
101 |
CALL inter_barxy(lon_rad, lat_rad(:jml_dyn -1), var_ana, rlonu(:iim), & |
CALL inter_barxy(lon_rad, lat_rad(:jml_dyn -1), var_ana, rlonu(:iim), & |
102 |
rlatv, tmp_var) |
rlatv, tmp_var) |
103 |
psol = gr_int_dyn(tmp_var) |
psol = gr_int_dyn(tmp_var) |
|
CALL start_init_phys(tsol_2d) |
|
104 |
|
|
105 |
psol(:iim, :) = psol(:iim, :) & |
psol(:iim, :) = psol(:iim, :) & |
106 |
* (1. + (z(:iim, :) - phis(:iim, :)) / 287. / tsol_2d(:iim, :)) |
* (1. + (z(:iim, :) - phis(:iim, :)) / 287. / tsol_2d(:iim, :)) |