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