3 |
! This program sets the initial and boundary values. |
! This program sets the initial and boundary values. |
4 |
|
|
5 |
use conf_gcm_m, only: conf_gcm |
use conf_gcm_m, only: conf_gcm |
6 |
|
use dimens_m, only: iim, jjm |
7 |
use etat0_mod, only: etat0 |
use etat0_mod, only: etat0 |
8 |
|
use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub |
9 |
use jumble, only: new_unit |
use jumble, only: new_unit |
10 |
use limit_mod, only: limit |
use limit_mod, only: limit |
11 |
|
use serre, only: read_serre |
12 |
use unit_nml_m, only: unit_nml |
use unit_nml_m, only: unit_nml |
13 |
|
|
14 |
implicit none |
implicit none |
15 |
|
|
16 |
|
REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2 |
17 |
|
|
18 |
!------------------------------------- |
!------------------------------------- |
19 |
|
|
20 |
call new_unit(unit_nml) |
call new_unit(unit_nml) |
21 |
open(unit_nml, file="used_namelists", status="replace", action="write") |
open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
22 |
CALL conf_gcm |
CALL conf_gcm |
23 |
CALL etat0 |
call read_serre |
24 |
|
CALL etat0(phis) |
25 |
CALL limit |
CALL limit |
26 |
|
CALL grilles_gcm_netcdf_sub(phis) |
27 |
close(unit_nml) |
close(unit_nml) |
28 |
|
|
29 |
END PROGRAM ce0l |
END PROGRAM ce0l |