trunk/dyn3d/ce0l.f

PROGRAM ce0l

  ! This program sets the initial and boundary values.

  use comdissnew, only: read_comdissnew
  use conf_gcm_m, only: conf_gcm
  use dimensions, only: iim, jjm
  use etat0_mod, only: etat0
  use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub
  use jumble, only: new_unit
  use limit_mod, only: limit
  use read_serre_m, only: read_serre
  use unit_nml_m, only: unit_nml

  implicit none

  REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2

  !-------------------------------------

  call new_unit(unit_nml)
  open(unit_nml, file="used_namelists.txt", status="replace", action="write")
  CALL conf_gcm
  call read_comdissnew
  call read_serre
  CALL etat0(phis)
  CALL limit
  CALL grilles_gcm_netcdf_sub(phis)
  close(unit_nml)
  print *, "ce0l: done"

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