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	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	3	use etat0_mod, only: etat0
9	guez	107	use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub
10	guez	57	use jumble, only: new_unit
11	guez	3	use limit_mod, only: limit
12	guez	139	use read_serre_m, only: read_serre
13	guez	57	use unit_nml_m, only: unit_nml
14	guez	3
15			implicit none
16
17	guez	107	REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2
18
19	guez	3	!-------------------------------------
20
21	guez	57	call new_unit(unit_nml)
22	guez	79	open(unit_nml, file="used_namelists.txt", status="replace", action="write")
23	guez	3	CALL conf_gcm
24	guez	191	call read_comdissnew
25	guez	113	call read_serre
26	guez	107	CALL etat0(phis)
27	guez	3	CALL limit
28	guez	107	CALL grilles_gcm_netcdf_sub(phis)
29	guez	57	close(unit_nml)
30	guez	212	print *, "ce0l: done"
31	guez	3
32	guez	68	END PROGRAM ce0l