1 |
guez |
68 |
PROGRAM ce0l |
2 |
guez |
3 |
|
3 |
|
|
! This program sets the initial and boundary values. |
4 |
|
|
|
5 |
guez |
191 |
use comdissnew, only: read_comdissnew |
6 |
guez |
3 |
use conf_gcm_m, only: conf_gcm |
7 |
guez |
265 |
use dimensions, only: iim, jjm |
8 |
guez |
276 |
use dimphy, only: klon |
9 |
|
|
use dynetat0_m, only: read_serre |
10 |
|
|
use etat0_m, only: etat0 |
11 |
guez |
107 |
use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub |
12 |
guez |
276 |
use indicesol, only: nbsrf |
13 |
guez |
3 |
use limit_mod, only: limit |
14 |
guez |
276 |
use unit_nml_m, only: unit_nml, set_unit_nml |
15 |
guez |
3 |
|
16 |
|
|
implicit none |
17 |
|
|
|
18 |
guez |
107 |
REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2 |
19 |
|
|
|
20 |
guez |
276 |
REAL pctsrf(klon, nbsrf) |
21 |
|
|
! ("pctsrf(i, :)" is the composition of the surface at horizontal |
22 |
|
|
! position "i") |
23 |
|
|
|
24 |
guez |
3 |
!------------------------------------- |
25 |
|
|
|
26 |
guez |
276 |
call set_unit_nml |
27 |
guez |
79 |
open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
28 |
guez |
3 |
CALL conf_gcm |
29 |
guez |
191 |
call read_comdissnew |
30 |
guez |
113 |
call read_serre |
31 |
guez |
276 |
CALL etat0(phis, pctsrf) |
32 |
|
|
CALL limit(pctsrf) |
33 |
guez |
107 |
CALL grilles_gcm_netcdf_sub(phis) |
34 |
guez |
57 |
close(unit_nml) |
35 |
guez |
212 |
print *, "ce0l: done" |
36 |
guez |
3 |
|
37 |
guez |
68 |
END PROGRAM ce0l |