trunk/dyn3d/ce0l.f90

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	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	276	use dimphy, only: klon
9			use dynetat0_m, only: read_serre
10			use etat0_m, only: etat0
11	guez	107	use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub
12	guez	276	use indicesol, only: nbsrf
13	guez	3	use limit_mod, only: limit
14	guez	276	use unit_nml_m, only: unit_nml, set_unit_nml
15	guez	3
16			implicit none
17
18	guez	107	REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2
19
20	guez	276	REAL pctsrf(klon, nbsrf)
21			! ("pctsrf(i, :)" is the composition of the surface at horizontal
22			! position "i")
23
24	guez	3	!-------------------------------------
25
26	guez	276	call set_unit_nml
27	guez	79	open(unit_nml, file="used_namelists.txt", status="replace", action="write")
28	guez	3	CALL conf_gcm
29	guez	191	call read_comdissnew
30	guez	113	call read_serre
31	guez	276	CALL etat0(phis, pctsrf)
32			CALL limit(pctsrf)
33	guez	107	CALL grilles_gcm_netcdf_sub(phis)
34	guez	57	close(unit_nml)
35	guez	212	print *, "ce0l: done"
36	guez	3
37	guez	68	END PROGRAM ce0l