/[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 107 by guez, Thu Sep 11 15:09:15 2014 UTC trunk/Sources/dyn3d/etat0.f revision 138 by guez, Fri May 22 23:13:19 2015 UTC
# Line 21  contains Line 21  contains
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: rlatu, rlonv, rlonu, rlatv, aire_2d, apoln, apols, &
23           cu_2d, cv_2d, inigeom           cu_2d, cv_2d, inigeom
     use conf_gcm_m, only: dayref, anneeref  
24      use dimens_m, only: iim, jjm, llm, nqmx      use dimens_m, only: iim, jjm, llm, nqmx
25      use dimphy, only: zmasq      use dimphy, only: zmasq
26      use dimsoil, only: nsoilmx      use dimsoil, only: nsoilmx
27      use disvert_m, only: ap, bp, preff, pa, disvert      use disvert_m, only: ap, bp, preff, pa, disvert
28        use dynetat0_m, only: day_ref, annee_ref
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 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
48      use regr_pr_o3_m, only: regr_pr_o3      use regr_pr_o3_m, only: regr_pr_o3
     use serre, only: alphax  
49      use startdyn, only: start_init_dyn      use startdyn, only: start_init_dyn
50      USE start_init_orog_m, only: start_init_orog, mask      USE start_init_orog_m, only: start_init_orog, mask
51      use start_init_phys_m, only: start_init_phys      use start_init_phys_m, only: start_init_phys
52      use start_inter_3d_m, only: start_inter_3d      use start_inter_3d_m, only: start_inter_3d
53      use temps, only: itau_phy, annee_ref, day_ref      use temps, only: itau_phy
54      use test_disvert_m, only: test_disvert      use test_disvert_m, only: test_disvert
55        use unit_nml_m, only: unit_nml
56    
57      REAL, intent(out):: phis(:, :) ! (iim + 1, jjm + 1)      REAL, intent(out):: phis(:, :) ! (iim + 1, jjm + 1)
58      ! surface geopotential, in m2 s-2      ! surface geopotential, in m2 s-2
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 123  contains Line 119  contains
119      ! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)",      ! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)",
120      ! for interface "l")      ! for interface "l")
121    
122        namelist /etat0_nml/ day_ref, annee_ref
123    
124      !---------------------------------      !---------------------------------
125    
126      print *, "Call sequence information: etat0"      print *, "Call sequence information: etat0"
127    
128        print *, "Enter namelist 'etat0_nml'."
129        read(unit=*, nml=etat0_nml)
130        write(unit_nml, nml=etat0_nml)
131    
132      CALL iniconst      CALL iniconst
133    
134      ! Construct a grid:      ! Construct a grid:
# Line 137  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 296  contains Line 298  contains
298      ! Calcul interm\'ediaire :      ! Calcul interm\'ediaire :
299      CALL massdair(p3d, masse)      CALL massdair(p3d, masse)
300    
     print *, 'ALPHAX = ', alphax  
   
301      forall (l = 1:llm)      forall (l = 1:llm)
302         masse(:, 1, l) = SUM(aire_2d(:iim, 1) * masse(:iim, 1, l)) / apoln         masse(:, 1, l) = SUM(aire_2d(:iim, 1) * masse(:iim, 1, l)) / apoln
303         masse(:, jjm + 1, l) = &         masse(:, jjm + 1, l) = &
# Line 307  contains Line 307  contains
307      ! Initialisation pour traceurs:      ! Initialisation pour traceurs:
308      call iniadvtrac      call iniadvtrac
309      itau_phy = 0      itau_phy = 0
     day_ref = dayref  
     annee_ref = anneeref  
310    
311      CALL geopot(teta, pk , pks, phis, phi)      CALL geopot(teta, pk , pks, phis, phi)
312      CALL caldyn0(ucov, vcov, teta, ps, masse, pk, phis, phi, w, pbaru, &      CALL caldyn0(ucov, vcov, teta, ps, masse, pk, phis, phi, w, pbaru, &
313           pbarv)           pbarv)
314      CALL dynredem0("start.nc", dayref, phis)      CALL dynredem0("start.nc", day_ref, phis)
315      CALL dynredem1("start.nc", vcov, ucov, teta, q, masse, ps, itau=0)      CALL dynredem1("start.nc", vcov, ucov, teta, q, masse, ps, itau=0)
316    
317      ! Initialisations :      ! Initialisations :
# Line 350  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.107  
changed lines
  Added in v.138
  ViewVC Help
Powered by ViewVC 1.1.21