Sources/dyn3d/gcm.f

PROGRAM gcm

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

  ! 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".

  use comconst, only: daysec, 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, iecri, iphysiq, nday, periodav, &
       conf_gcm, iflag_phys
  use conf_guide_m, only: conf_guide
  use dimens_m, only: iim, jjm, llm, nqmx
  USE disvert_m, ONLY : disvert
  use dynetat0_m, only: dynetat0, day_ini
  use dynredem0_m, only: dynredem0
  use grid_change, only: dyn_phy, init_dyn_phy
  use histclo_m, only: histclo
  use iniadvtrac_m, only: iniadvtrac
  use inidissip_m, only: inidissip
  use inifilr_m, only: inifilr
  use initdynav_m, only: initdynav
  use inithist_m, only: inithist
  use init_dynzon_m, only: init_dynzon
  USE ioconf_calendar_m, only: ioconf_calendar
  use jumble, only: new_unit
  use leapfrog_m, only: leapfrog
  use netcdf95, only: nf95_close
  use suphec_m, only: suphec
  use tracstoke, only: istdyn, istphy
  use unit_nml_m, only: unit_nml
  use yoethf_m, only: yoethf
  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) ! champs advect\'es
  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 new_unit(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 iniadvtrac
  CALL iniconst
  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
     call yoethf
  ENDIF

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

  ! Choix des fr\'equences de stockage pour le hors-ligne :
  istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4
  istphy = istdyn / iphysiq     

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