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

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

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

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 138 by guez, Fri May 22 23:13:19 2015 UTC
# Line 32  contains Line 32  contains
32      use geopot_m, only: geopot      use geopot_m, only: geopot
33      use grid_atob, only: grille_m      use grid_atob, only: grille_m
34      use grid_change, only: init_dyn_phy, dyn_phy      use grid_change, only: init_dyn_phy, dyn_phy
     use histclo_m, only: histclo  
35      use indicesol, only: is_oce, is_sic, is_ter, is_lic, epsfra      use indicesol, only: is_oce, is_sic, is_ter, is_lic, epsfra
36      use iniadvtrac_m, only: iniadvtrac      use iniadvtrac_m, only: iniadvtrac
37      use inifilr_m, only: inifilr      use inifilr_m, only: inifilr
# Line 42  contains Line 41  contains
41           nf95_inq_varid, nf95_open           nf95_inq_varid, nf95_open
42      use nr_util, only: pi, assert      use nr_util, only: pi, assert
43      use paramet_m, only: ip1jm, ip1jmp1      use paramet_m, only: ip1jm, ip1jmp1
44        use phyetat0_m, only: rlat, rlon
45      use phyredem_m, only: phyredem      use phyredem_m, only: phyredem
46      use q_sat_m, only: q_sat      use q_sat_m, only: q_sat
47      use regr_lat_time_coefoz_m, only: regr_lat_time_coefoz      use regr_lat_time_coefoz_m, only: regr_lat_time_coefoz
# Line 59  contains Line 59  contains
59    
60      ! Local:      ! Local:
61    
     REAL latfi(klon), lonfi(klon)  
     ! (latitude and longitude of a point of the scalar grid identified  
     ! by a simple index, in degrees)  
   
62      REAL, dimension(iim + 1, jjm + 1, llm):: ucov, t3d, teta      REAL, dimension(iim + 1, jjm + 1, llm):: ucov, t3d, teta
63      REAL vcov(iim + 1, jjm, llm)      REAL vcov(iim + 1, jjm, llm)
64    
# Line 143  contains Line 139  contains
139      CALL inigeom      CALL inigeom
140      CALL inifilr      CALL inifilr
141    
142      latfi(1) = 90.      rlat(1) = 90.
143      latfi(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.) * 180. / pi      rlat(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.) * 180. / pi
144      ! (with conversion to degrees)      ! (with conversion to degrees)
145      latfi(klon) = - 90.      rlat(klon) = - 90.
146    
147      lonfi(1) = 0.      rlon(1) = 0.
148      lonfi(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.) * 180. / pi      rlon(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.) * 180. / pi
149      ! (with conversion to degrees)      ! (with conversion to degrees)
150      lonfi(klon) = 0.      rlon(klon) = 0.
151    
152      call start_init_orog(phis, zmea_2d, zstd_2d, zsig_2d, zgam_2d, zthe_2d, &      call start_init_orog(phis, zmea_2d, zstd_2d, zsig_2d, zgam_2d, zthe_2d, &
153           zpic_2d, zval_2d) ! also compute "mask"           zpic_2d, zval_2d) ! also compute "mask"
# Line 352  contains Line 348  contains
348      sig1 = 0.      sig1 = 0.
349      w01 = 0.      w01 = 0.
350    
351      call phyredem("startphy.nc", latfi, lonfi, pctsrf, tsoil(:, 1, :), tsoil, &      call phyredem("startphy.nc", pctsrf, tsoil(:, 1, :), tsoil, &
352           tsoil(:, 1, is_oce), seaice, qsolsrf, pack(qsol_2d, dyn_phy), snsrf, &           tsoil(:, 1, is_oce), seaice, qsolsrf, pack(qsol_2d, dyn_phy), snsrf, &
353           albe, alblw, evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, &           albe, alblw, evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, &
354           frugs, agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, &           frugs, agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, &
355           q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01)           q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01)
     CALL histclo  
356    
357    END SUBROUTINE etat0    END SUBROUTINE etat0
358    
Legend:
Removed from v.134  
changed lines
  Added in v.138
  ViewVC Help
Powered by ViewVC 1.1.21