trunk/dyn3d/gcm.f90

PROGRAM gcm

  ! Authors: P. Le Van, L. Fairhead, F. Hourdin
  ! From "gcm.F", version 1.4, 2006/04/04

  ! General circulation model of LMD. Avec coordonn\'ee verticale
  ! hybride, avec nouveaux op\'erateurs de dissipation "*" ("gradiv2",
  ! "divgrad2", "nxgraro2"). Possibilit\'e de choisir le sch\'ema pour
  ! l'advection de "q", en modifiant "iadv" dans "traceur.def".

  ! Libraries:
  use netcdf95, only: nf95_close

  use comconst, only: dtvr, iniconst
  use comdissnew, only: read_comdissnew
  use comgeom, only:  aire_2d, inigeom
  use comgeomphy, only: airephy
  use conf_gcm_m, only: day_step, iperiod, iphysiq, nday, conf_gcm, iflag_phys
  use conf_guide_m, only: conf_guide
  use dimensions, only: iim, jjm, llm, nqmx
  USE disvert_m, ONLY : disvert
  use dynetat0_m, only: dynetat0, day_ini
  use dynetat0_chosen_m, only: dynetat0_chosen
  use dynredem0_m, only: dynredem0
  use grid_change, only: dyn_phy, init_dyn_phy
  use histclo_m, only: histclo
  use infotrac_init_m, only: infotrac_init
  use inidissip_m, only: inidissip
  use inifilr_m, only: inifilr
  use inithist_m, only: inithist
  use init_dynzon_m, only: init_dynzon
  USE ioconf_calendar_m, only: ioconf_calendar
  use leapfrog_m, only: leapfrog
  use suphec_m, only: suphec
  use unit_nml_m, only: unit_nml, set_unit_nml
  use createnewfield_m, only: NbField, Ncid

  IMPLICIT NONE

  ! Variables dynamiques :
  REAL ucov(iim + 1, jjm + 1, llm), vcov(iim + 1, jjm, llm)  ! vent covariant
  REAL teta(iim + 1, jjm + 1, llm) ! temp\'erature potentielle 
  REAL q(iim + 1, jjm + 1, llm, nqmx) ! mass fraction of advected species
  REAL ps(iim + 1, jjm + 1) ! pression au sol (Pa)
  REAL masse(iim + 1, jjm + 1, llm) ! masse d'air
  REAL phis(iim + 1, jjm + 1) ! g\'eopotentiel au sol

  LOGICAL:: true_calendar = .false. ! default value
  integer i

  namelist /main_nml/true_calendar

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

  call set_unit_nml
  open(unit_nml, file="used_namelists.txt", status="replace", action="write")

  CALL conf_gcm
  call read_comdissnew

  print *, "Enter namelist 'main_nml'."
  read (unit=*, nml=main_nml)
  write(unit_nml, nml=main_nml)

  ! Choix du calendrier :
  if (true_calendar) then
     call ioconf_calendar('gregorian')
  else
     call ioconf_calendar('360d')
  endif

  call infotrac_init
  CALL iniconst
  CALL dynetat0_chosen
  CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)
  CALL disvert
  CALL inigeom ! initialisation de la g\'eometrie
  CALL inifilr ! initialisation du filtre
  CALL inidissip
  call init_dyn_phy

  ! Initialisation de la physique :
  IF (iflag_phys) THEN
     airephy = pack(aire_2d, dyn_phy)
     CALL suphec
  ENDIF

  ! Initialisation des entr\'ees-sorties :
  CALL dynredem0(day_ini + nday, phis)
  CALL inithist(t_ops = dtvr, t_wrt = dtvr)
  call init_dynzon(dt_app = dtvr * iperiod)

  CALL conf_guide
  CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)

  close(unit_nml)
  call histclo

  do i = 1, nbfield
     call nf95_close(Ncid(i))
  end do

  print *, 'Simulation finished'
  print *, 'Everything is cool'

END PROGRAM gcm
1	guez	3	PROGRAM gcm
2
3	guez	36	! Authors: P. Le Van, L. Fairhead, F. Hourdin
4	guez	344	! From "gcm.F", version 1.4, 2006/04/04
5	guez	3
6	guez	139	! General circulation model of LMD. Avec coordonn\'ee verticale
7			! hybride, avec nouveaux op\'erateurs de dissipation "*" ("gradiv2",
8			! "divgrad2", "nxgraro2"). Possibilit\'e de choisir le sch\'ema pour
9	guez	38	! l'advection de "q", en modifiant "iadv" dans "traceur.def".
10
11	guez	344	! Libraries:
12			use netcdf95, only: nf95_close
13
14	guez	263	use comconst, only: dtvr, iniconst
15	guez	191	use comdissnew, only: read_comdissnew
16	guez	168	use comgeom, only: aire_2d, inigeom
17			use comgeomphy, only: airephy
18	guez	266	use conf_gcm_m, only: day_step, iperiod, iphysiq, nday, conf_gcm, iflag_phys
19	guez	115	use conf_guide_m, only: conf_guide
20	guez	265	use dimensions, only: iim, jjm, llm, nqmx
21	guez	79	USE disvert_m, ONLY : disvert
22	guez	168	use dynetat0_m, only: dynetat0, day_ini
23	guez	313	use dynetat0_chosen_m, only: dynetat0_chosen
24	guez	26	use dynredem0_m, only: dynredem0
25	guez	3	use grid_change, only: dyn_phy, init_dyn_phy
26	guez	61	use histclo_m, only: histclo
27	guez	321	use infotrac_init_m, only: infotrac_init
28	guez	26	use inidissip_m, only: inidissip
29	guez	54	use inifilr_m, only: inifilr
30	guez	26	use inithist_m, only: inithist
31	guez	57	use init_dynzon_m, only: init_dynzon
32	guez	92	USE ioconf_calendar_m, only: ioconf_calendar
33	guez	3	use leapfrog_m, only: leapfrog
34	guez	37	use suphec_m, only: suphec
35	guez	276	use unit_nml_m, only: unit_nml, set_unit_nml
36	guez	110	use createnewfield_m, only: NbField, Ncid
37	guez	3
38			IMPLICIT NONE
39
40			! Variables dynamiques :
41	guez	55	REAL ucov(iim + 1, jjm + 1, llm), vcov(iim + 1, jjm, llm) ! vent covariant
42	guez	139	REAL teta(iim + 1, jjm + 1, llm) ! temp\'erature potentielle
43	guez	344	REAL q(iim + 1, jjm + 1, llm, nqmx) ! mass fraction of advected species
44	guez	39	REAL ps(iim + 1, jjm + 1) ! pression au sol (Pa)
45	guez	55	REAL masse(iim + 1, jjm + 1, llm) ! masse d'air
46	guez	139	REAL phis(iim + 1, jjm + 1) ! g\'eopotentiel au sol
47	guez	3
48			LOGICAL:: true_calendar = .false. ! default value
49	guez	108	integer i
50
51	guez	3	namelist /main_nml/true_calendar
52
53			!------------------------------------------------------------
54
55	guez	276	call set_unit_nml
56	guez	79	open(unit_nml, file="used_namelists.txt", status="replace", action="write")
57	guez	57
58			CALL conf_gcm
59	guez	191	call read_comdissnew
60	guez	57
61	guez	3	print *, "Enter namelist 'main_nml'."
62			read (unit=*, nml=main_nml)
63	guez	57	write(unit_nml, nml=main_nml)
64	guez	3
65			! Choix du calendrier :
66			if (true_calendar) then
67			call ioconf_calendar('gregorian')
68			else
69			call ioconf_calendar('360d')
70			endif
71
72	guez	321	call infotrac_init
73	guez	79	CALL iniconst
74	guez	313	CALL dynetat0_chosen
75	guez	128	CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)
76	guez	79	CALL disvert
77	guez	139	CALL inigeom ! initialisation de la g\'eometrie
78	guez	36	CALL inifilr ! initialisation du filtre
79	guez	27	CALL inidissip
80	guez	3	call init_dyn_phy
81
82			! Initialisation de la physique :
83	guez	208	IF (iflag_phys) THEN
84	guez	37	airephy = pack(aire_2d, dyn_phy)
85			CALL suphec
86	guez	3	ENDIF
87
88	guez	139	! Initialisation des entr\'ees-sorties :
89	guez	157	CALL dynredem0(day_ini + nday, phis)
90	guez	266	CALL inithist(t_ops = dtvr, t_wrt = dtvr)
91	guez	57	call init_dynzon(dt_app = dtvr * iperiod)
92	guez	3
93	guez	115	CALL conf_guide
94	guez	128	CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)
95	guez	3
96	guez	68	close(unit_nml)
97	guez	10	call histclo
98	guez	108
99			do i = 1, nbfield
100	guez	109	call nf95_close(Ncid(i))
101	guez	108	end do
102
103	guez	10	print *, 'Simulation finished'
104			print *, 'Everything is cool'
105
106	guez	3	END PROGRAM gcm