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 xios, only: xios_initialize, xios_finalize

  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)

  call xios_initialize("LMDZE")
  
  ! 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

  call xios_finalize

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