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