/[lmdze]/trunk/phylmd/Mobidic/regr_pr_o3.f
ViewVC logotype

Diff of /trunk/phylmd/Mobidic/regr_pr_o3.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 96 by guez, Wed Mar 5 14:57:53 2014 UTC revision 97 by guez, Fri Apr 25 14:58:31 2014 UTC
# Line 4  module regr_pr_o3_m Line 4  module regr_pr_o3_m
4    
5  contains  contains
6    
7    subroutine regr_pr_o3(o3_mob_regr)    subroutine regr_pr_o3(p3d, o3_mob_regr)
8    
9      ! "regr_pr_o3" stands for "regrid pressure ozone".      ! "regr_pr_o3" stands for "regrid pressure ozone".
10      ! This procedure reads Mobidic ozone mole fraction from      ! This procedure reads Mobidic ozone mole fraction from
# Line 30  contains Line 30  contains
30      use nr_util, only: assert      use nr_util, only: assert
31      use grid_change, only: dyn_phy      use grid_change, only: dyn_phy
32      use numer_rec_95, only: regr1_step_av      use numer_rec_95, only: regr1_step_av
33      use pressure_var, only: p3d  
34        REAL, intent(in):: p3d(:, :, :) ! (iim + 1, jjm + 1, llm+1)
35        ! pressure at layer interfaces, in Pa
36        ! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)",
37        ! for interface "l")
38    
39      real, intent(out):: o3_mob_regr(:, :, :) ! (iim + 1, jjm + 1, llm)      real, intent(out):: o3_mob_regr(:, :, :) ! (iim + 1, jjm + 1, llm)
40      ! (ozone mole fraction from Mobidic adapted to the LMDZ grid)      ! (ozone mole fraction from Mobidic adapted to the LMDZ grid)
# Line 57  contains Line 61  contains
61      !------------------------------------------------------------      !------------------------------------------------------------
62    
63      print *, "Call sequence information: regr_pr_o3"      print *, "Call sequence information: regr_pr_o3"
64      call assert(shape(o3_mob_regr) == (/iim + 1, jjm + 1, llm/), "regr_pr_o3")  
65        call assert(shape(o3_mob_regr) == (/iim + 1, jjm + 1, llm/), &
66             "regr_pr_o3 o3_mob_regr")
67        call assert(shape(p3d) == (/iim + 1, jjm + 1, llm + 1/), &
68             "regr_pr_o3 p3d")
69    
70      call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid)      call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid)
71    
Legend:
Removed from v.96  
changed lines
  Added in v.97
  ViewVC Help
Powered by ViewVC 1.1.21