/[lmdze]/trunk/Sources/dyn3d/startdyn.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/startdyn.f

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

revision 36 by guez, Thu Dec 2 17:11:04 2010 UTC revision 43 by guez, Fri Apr 8 12:43:31 2011 UTC
# Line 27  CONTAINS Line 27  CONTAINS
27      use comgeom, only: aire_2d, apoln, apols      use comgeom, only: aire_2d, apoln, apols
28      use conf_dat2d_m, only: conf_dat2d      use conf_dat2d_m, only: conf_dat2d
29      use inter_barxy_m, only: inter_barxy      use inter_barxy_m, only: inter_barxy
     use comconst, only: pi  
30      use comgeom, only: rlonu, rlatv      use comgeom, only: rlonu, rlatv
31      use dimens_m, only: iim, jjm      use dimens_m, only: iim, jjm
32      use gr_int_dyn_m, only: gr_int_dyn      use gr_int_dyn_m, only: gr_int_dyn
33      use start_init_orog_m, only: phis      use start_init_orog_m, only: phis
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      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    
42      REAL date, dt      REAL date, dt
43      INTEGER itau(1)      INTEGER itau(1)
     INTEGER i, j  
44      REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:)      REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:)
45    
46      REAL, ALLOCATABLE:: lon_dyn(:, :), lat_dyn(:, :)      REAL, ALLOCATABLE:: lon_dyn(:, :), lat_dyn(:, :)
# Line 56  CONTAINS Line 54  CONTAINS
54    
55      print *, "Call sequence information: start_init_dyn"      print *, "Call sequence information: start_init_dyn"
56      call assert((/size(tsol_2d, 1), size(psol, 1)/) == iim + 1, &      call assert((/size(tsol_2d, 1), size(psol, 1)/) == iim + 1, &
57           "start_init_phys size 1")           "start_init_dyn size 1")
58      call assert((/size(tsol_2d, 2), size(psol, 2)/) == jjm + 1, &      call assert((/size(tsol_2d, 2), size(psol, 2)/) == jjm + 1, &
59           "start_init_phys size 2")           "start_init_dyn size 2")
60      CALL flininfo('ECDYN.nc', iml_dyn, jml_dyn, llm_dyn, ttm_dyn, fid_dyn)      CALL flininfo('ECDYN.nc', iml_dyn, jml_dyn, llm_dyn, ttm_dyn, fid_dyn)
61      print *, "iml_dyn = ", iml_dyn, ", jml_dyn = ", jml_dyn, &      print *, "iml_dyn = ", iml_dyn, ", jml_dyn = ", jml_dyn, &
62           ", llm_dyn = ", llm_dyn, ", ttm_dyn = ", ttm_dyn           ", llm_dyn = ", llm_dyn, ", ttm_dyn = ", ttm_dyn
# Line 103  CONTAINS Line 101  CONTAINS
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, :))
Legend:
Removed from v.36  
changed lines
  Added in v.43
  ViewVC Help
Powered by ViewVC 1.1.21