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 dimphy, only: klon
  use dynetat0_m, only: read_serre
  use etat0_m, only: etat0
  use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub
  use indicesol, only: nbsrf
  use limit_mod, only: limit
  use unit_nml_m, only: unit_nml, set_unit_nml

  implicit none

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

  REAL pctsrf(klon, nbsrf)
  ! ("pctsrf(i, :)" is the composition of the surface at horizontal
  ! position "i")

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

  call set_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, pctsrf)
  CALL limit(pctsrf)
  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 dimphy, only: klon
9	use dynetat0_m, only: read_serre
10	use etat0_m, only: etat0
11	use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub
12	use indicesol, only: nbsrf
13	use limit_mod, only: limit
14	use unit_nml_m, only: unit_nml, set_unit_nml
15
16	implicit none
17
18	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 set_unit_nml
27	open(unit_nml, file="used_namelists.txt", status="replace", action="write")
28	CALL conf_gcm
29	call read_comdissnew
30	call read_serre
31	CALL etat0(phis, pctsrf)
32	CALL limit(pctsrf)
33	CALL grilles_gcm_netcdf_sub(phis)
34	close(unit_nml)
35	print *, "ce0l: done"
36
37	END PROGRAM ce0l