18 |
! From "etat0_netcdf.F", version 1.3 2005/05/25 13:10:09 |
! From "etat0_netcdf.F", version 1.3 2005/05/25 13:10:09 |
19 |
|
|
20 |
use caldyn0_m, only: caldyn0 |
use caldyn0_m, only: caldyn0 |
21 |
use comconst, only: dtvr, daysec, cpp, kappa |
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 |
24 |
use conf_gcm_m, only: day_step, iphysiq, dayref, anneeref |
use conf_gcm_m, only: dayref, anneeref |
25 |
use dimens_m, only: iim, jjm, llm, nqmx |
use dimens_m, only: iim, jjm, llm, nqmx |
26 |
use dimphy, only: zmasq |
use dimphy, only: zmasq |
27 |
use dimsoil, only: nsoilmx |
use dimsoil, only: nsoilmx |
28 |
use disvert_m, only: ap, bp, preff, pa |
use disvert_m, only: ap, bp, preff, pa, disvert |
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 |
109 |
REAL phi(iim + 1, jjm + 1, llm) |
REAL phi(iim + 1, jjm + 1, llm) |
110 |
REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm) |
REAL pbaru(ip1jmp1, llm), pbarv(ip1jm, llm) |
111 |
REAL w(ip1jmp1, llm) |
REAL w(ip1jmp1, llm) |
|
REAL phystep |
|
112 |
|
|
113 |
real sig1(klon, llm) ! section adiabatic updraft |
real sig1(klon, llm) ! section adiabatic updraft |
114 |
real w01(klon, llm) ! vertical velocity within adiabatic updraft |
real w01(klon, llm) ! vertical velocity within adiabatic updraft |
117 |
|
|
118 |
print *, "Call sequence information: etat0" |
print *, "Call sequence information: etat0" |
119 |
|
|
120 |
dtvr = daysec / real(day_step) |
CALL iniconst |
|
print *, 'dtvr = ', dtvr |
|
121 |
|
|
122 |
! Construct a grid: |
! Construct a grid: |
123 |
|
|
124 |
pa = 5e4 |
pa = 5e4 |
125 |
CALL iniconst |
CALL disvert |
126 |
CALL inigeom |
CALL inigeom |
127 |
CALL inifilr |
CALL inifilr |
128 |
|
|
311 |
CALL dynredem0("start.nc", dayref, phis) |
CALL dynredem0("start.nc", dayref, phis) |
312 |
CALL dynredem1("start.nc", vcov, ucov, teta, q, masse, ps, itau=0) |
CALL dynredem1("start.nc", vcov, ucov, teta, q, masse, ps, itau=0) |
313 |
|
|
|
! Ecriture état initial physique: |
|
|
print *, "iphysiq = ", iphysiq |
|
|
phystep = dtvr * REAL(iphysiq) |
|
|
print *, 'phystep = ', phystep |
|
|
|
|
314 |
! Initialisations : |
! Initialisations : |
315 |
tsolsrf(:, is_ter) = tsol |
tsolsrf(:, is_ter) = tsol |
316 |
tsolsrf(:, is_lic) = tsol |
tsolsrf(:, is_lic) = tsol |