/[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

trunk/dyn3d/etat0.f revision 129 by guez, Fri Feb 13 18:22:38 2015 UTC trunk/Sources/dyn3d/etat0.f revision 139 by guez, Tue May 26 17:46:03 2015 UTC
# Line 19  contains Line 19  contains
19    
20      use caldyn0_m, only: caldyn0      use caldyn0_m, only: caldyn0
21      use comconst, only: cpp, kappa, iniconst      use comconst, only: cpp, kappa, iniconst
22      use comgeom, only: rlatu, rlonv, rlonu, rlatv, aire_2d, apoln, apols, &      use comgeom, only: aire_2d, apoln, apols, cu_2d, cv_2d, inigeom
          cu_2d, cv_2d, inigeom  
23      use dimens_m, only: iim, jjm, llm, nqmx      use dimens_m, only: iim, jjm, llm, nqmx
24      use dimphy, only: zmasq      use dimphy, only: zmasq
25      use dimsoil, only: nsoilmx      use dimsoil, only: nsoilmx
26      use disvert_m, only: ap, bp, preff, pa, disvert      use disvert_m, only: ap, bp, preff, pa, disvert
27      use dynetat0_m, only: day_ref, annee_ref      use dynetat0_m, only: day_ref, annee_ref, xprimp025, xprimm025, rlatu1, &
28             rlatu2, rlatu, rlatv, yprimu1, yprimu2, rlonu, rlonv, xprimu, xprimv
29      use dynredem0_m, only: dynredem0      use dynredem0_m, only: dynredem0
30      use dynredem1_m, only: dynredem1      use dynredem1_m, only: dynredem1
31      use exner_hyb_m, only: exner_hyb      use exner_hyb_m, only: exner_hyb
32        use fxhyp_m, only: fxhyp
33        use fyhyp_m, only: fyhyp
34      use geopot_m, only: geopot      use geopot_m, only: geopot
35      use grid_atob, only: grille_m      use grid_atob, only: grille_m
36      use grid_change, only: init_dyn_phy, dyn_phy      use grid_change, only: init_dyn_phy, dyn_phy
     use histclo_m, only: histclo  
37      use indicesol, only: is_oce, is_sic, is_ter, is_lic, epsfra      use indicesol, only: is_oce, is_sic, is_ter, is_lic, epsfra
38      use iniadvtrac_m, only: iniadvtrac      use iniadvtrac_m, only: iniadvtrac
39      use inifilr_m, only: inifilr      use inifilr_m, only: inifilr
# Line 42  contains Line 43  contains
43           nf95_inq_varid, nf95_open           nf95_inq_varid, nf95_open
44      use nr_util, only: pi, assert      use nr_util, only: pi, assert
45      use paramet_m, only: ip1jm, ip1jmp1      use paramet_m, only: ip1jm, ip1jmp1
46        use phyetat0_m, only: rlat, rlon
47      use phyredem_m, only: phyredem      use phyredem_m, only: phyredem
48      use q_sat_m, only: q_sat      use q_sat_m, only: q_sat
49      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 61  contains
61    
62      ! Local:      ! Local:
63    
     REAL latfi(klon), lonfi(klon)  
     ! (latitude and longitude of a point of the scalar grid identified  
     ! by a simple index, in degrees)  
   
64      REAL, dimension(iim + 1, jjm + 1, llm):: ucov, t3d, teta      REAL, dimension(iim + 1, jjm + 1, llm):: ucov, t3d, teta
65      REAL vcov(iim + 1, jjm, llm)      REAL vcov(iim + 1, jjm, llm)
66    
# Line 140  contains Line 138  contains
138      pa = 5e4      pa = 5e4
139      CALL disvert      CALL disvert
140      call test_disvert      call test_disvert
141    
142        CALL fyhyp(rlatu, rlatv, rlatu2, yprimu2, rlatu1, yprimu1)
143        CALL fxhyp(xprimm025, rlonv, xprimv, rlonu, xprimu, xprimp025)
144    
145        rlatu(1) = pi / 2.
146        rlatu(jjm + 1) = -rlatu(1)
147    
148      CALL inigeom      CALL inigeom
149      CALL inifilr      CALL inifilr
150    
151      latfi(1) = 90.      rlat(1) = 90.
152      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
153      ! (with conversion to degrees)      ! (with conversion to degrees)
154      latfi(klon) = - 90.      rlat(klon) = - 90.
155    
156      lonfi(1) = 0.      rlon(1) = 0.
157      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
158      ! (with conversion to degrees)      ! (with conversion to degrees)
159      lonfi(klon) = 0.      rlon(klon) = 0.
160    
161      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, &
162           zpic_2d, zval_2d) ! also compute "mask"           zpic_2d, zval_2d) ! also compute "mask"
# Line 352  contains Line 357  contains
357      sig1 = 0.      sig1 = 0.
358      w01 = 0.      w01 = 0.
359    
360      call phyredem("startphy.nc", latfi, lonfi, pctsrf, tsoil(:, 1, :), tsoil, &      call phyredem("startphy.nc", pctsrf, tsoil(:, 1, :), tsoil, &
361           tsoil(:, 1, is_oce), seaice, qsolsrf, pack(qsol_2d, dyn_phy), snsrf, &           tsoil(:, 1, is_oce), seaice, qsolsrf, pack(qsol_2d, dyn_phy), snsrf, &
362           albe, alblw, evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, &           albe, alblw, evap, rain_fall, snow_fall, solsw, sollw, fder, radsol, &
363           frugs, agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, &           frugs, agesno, zmea, zstd, zsig, zgam, zthe, zpic, zval, t_ancien, &
364           q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01)           q_ancien, rnebcon, ratqs, clwcon, run_off_lic_0, sig1, w01)
     CALL histclo  
365    
366    END SUBROUTINE etat0    END SUBROUTINE etat0
367    
Legend:
Removed from v.129  
changed lines
  Added in v.139
  ViewVC Help
Powered by ViewVC 1.1.21