trunk/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ée verticale
  ! hybride, avec nouveaux opérateurs de dissipation "*" ("gradiv2",
  ! "divgrad2", "nxgraro2"). Possibilité de choisir le schéma pour
  ! l'advection de "q", en modifiant "iadv" dans "traceur.def".

  use comconst, only: daysec, dtvr, iniconst
  use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv, inigeom
  use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond
  use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, &
       nday, raz_date, periodav, conf_gcm, iflag_phys
  use conf_guide_m, only: conf_guide
  use dimens_m, only: iim, jjm, llm, nqmx
  use dimphy, only: klon
  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 temps, only: day_ref, annee_ref, day_end, itau_dyn
  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érature potentielle 
  REAL q(iim + 1, jjm + 1, llm, nqmx) ! champs advectés
  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éopotentiel au sol

  ! Variables pour le fichier histoire :
  REAL time_0 ! time in day, as a fraction of day, in [0, 1[

  ! Calendrier :
  LOGICAL:: true_calendar = .false. ! default value

  logical mask_v(iim + 1, jjm) 
  ! (mask for points in the "v" grid, first index is for longitude,
  ! second index is for latitude)

  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

  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

  ! Initialisation des traceurs
  ! Choix du schéma pour l'advection dans le fichier "traceur.def" ou via INCA
  call iniadvtrac

  CALL iniconst

  ! Lecture du fichier "start.nc" :
  CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)

  ! On remet le calendrier à zéro si demandé :
  if (raz_date) then
     print *, 'On réinitialise à la date lue dans la namelist.'
     annee_ref = anneeref
     day_ref = dayref
     day_ini = dayref
     itau_dyn = 0
  else
     print *, 'On garde les dates du fichier "start".'
  endif

  CALL disvert
  CALL inigeom ! initialisation de la géometrie
  CALL inifilr ! initialisation du filtre
  CALL inidissip
  call init_dyn_phy

  ! Initialisation de la physique :
  IF (iflag_phys == 1) THEN
     rlatd(1)=rlatu(1)
     rlatd(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.)
     rlatd(klon)= rlatu(jjm + 1)

     rlond(1)=0.
     rlond(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.)
     rlond(klon)= 0.

     cuphy = pack(cu_2d, dyn_phy)

     ! Construct a mask for points in the "v" grid:
     mask_v = .true.
     mask_v(2:, 1) = .false.
     mask_v(iim + 1, 2:) = .false.

     cvphy(:klon - 1) = pack(cv_2d, mask_v)
     cvphy(klon) = cv_2d(1, jjm)
     ! (that value of "cv_2d" is used twice in "cvphy")

     airephy = pack(aire_2d, dyn_phy)
     CALL suphec
     call yoethf
  ENDIF

  ! Initialisation des entrées-sorties :
  day_end = day_ini + nday
  print *, "day_ini = ", day_ini
  print *, "day_end = ", day_end

  CALL dynredem0("restart.nc", day_end, phis)
  CALL inithist(day_ref, annee_ref, dtvr, nqmx, t_ops = iecri * daysec, &
       t_wrt = iecri * daysec)
  CALL initdynav(day_ref, annee_ref, dtvr, nqmx, t_ops = iperiod * dtvr, &
       t_wrt = periodav * daysec)
  call init_dynzon(dt_app = dtvr * iperiod)

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

  CALL conf_guide

  ! Intégration temporelle du modèle :
  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	38	! General circulation model of LMD. Avec coordonnée verticale
7			! hybride, avec nouveaux opérateurs de dissipation "*" ("gradiv2",
8			! "divgrad2", "nxgraro2"). Possibilité de choisir le schéma pour
9			! l'advection de "q", en modifiant "iadv" dans "traceur.def".
10
11	guez	79	use comconst, only: daysec, dtvr, iniconst
12	guez	78	use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv, inigeom
13	guez	37	use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond
14	guez	3	use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, &
15	guez	57	nday, raz_date, periodav, conf_gcm, iflag_phys
16	guez	115	use conf_guide_m, only: conf_guide
17	guez	26	use dimens_m, only: iim, jjm, llm, nqmx
18			use dimphy, only: klon
19	guez	79	USE disvert_m, ONLY : disvert
20	guez	25	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 temps, only: day_ref, annee_ref, day_end, itau_dyn
36			use tracstoke, only: istdyn, istphy
37	guez	57	use unit_nml_m, only: unit_nml
38	guez	38	use yoethf_m, only: yoethf
39	guez	110	use createnewfield_m, only: NbField, Ncid
40	guez	3
41			IMPLICIT NONE
42
43			! Variables dynamiques :
44	guez	55	REAL ucov(iim + 1, jjm + 1, llm), vcov(iim + 1, jjm, llm) ! vent covariant
45	guez	43	REAL teta(iim + 1, jjm + 1, llm) ! température potentielle
46	guez	40	REAL q(iim + 1, jjm + 1, llm, nqmx) ! champs advectés
47	guez	39	REAL ps(iim + 1, jjm + 1) ! pression au sol (Pa)
48	guez	55	REAL masse(iim + 1, jjm + 1, llm) ! masse d'air
49	guez	7	REAL phis(iim + 1, jjm + 1) ! géopotentiel au sol
50	guez	3
51			! Variables pour le fichier histoire :
52	guez	25	REAL time_0 ! time in day, as a fraction of day, in [0, 1[
53	guez	3
54			! Calendrier :
55			LOGICAL:: true_calendar = .false. ! default value
56
57			logical mask_v(iim + 1, jjm)
58			! (mask for points in the "v" grid, first index is for longitude,
59			! second index is for latitude)
60
61	guez	108	integer i
62
63	guez	3	namelist /main_nml/true_calendar
64
65			!------------------------------------------------------------
66
67	guez	57	call new_unit(unit_nml)
68	guez	79	open(unit_nml, file="used_namelists.txt", status="replace", action="write")
69	guez	57
70			CALL conf_gcm
71
72	guez	3	print *, "Enter namelist 'main_nml'."
73			read (unit=*, nml=main_nml)
74	guez	57	write(unit_nml, nml=main_nml)
75	guez	3
76			! Choix du calendrier :
77			if (true_calendar) then
78			call ioconf_calendar('gregorian')
79			else
80			call ioconf_calendar('360d')
81			endif
82
83			! Initialisation des traceurs
84			! Choix du schéma pour l'advection dans le fichier "traceur.def" ou via INCA
85	guez	23	call iniadvtrac
86	guez	3
87	guez	79	CALL iniconst
88
89	guez	3	! Lecture du fichier "start.nc" :
90	guez	128	CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)
91	guez	3
92	guez	98	! On remet le calendrier à zéro si demandé :
93			if (raz_date) then
94			print *, 'On réinitialise à la date lue dans la namelist.'
95			annee_ref = anneeref
96			day_ref = dayref
97			day_ini = dayref
98			itau_dyn = 0
99			else
100			print *, 'On garde les dates du fichier "start".'
101	guez	3	endif
102
103	guez	79	CALL disvert
104	guez	36	CALL inigeom ! initialisation de la géometrie
105			CALL inifilr ! initialisation du filtre
106	guez	27	CALL inidissip
107	guez	3	call init_dyn_phy
108
109			! Initialisation de la physique :
110			IF (iflag_phys == 1) THEN
111	guez	37	rlatd(1)=rlatu(1)
112			rlatd(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.)
113			rlatd(klon)= rlatu(jjm + 1)
114	guez	3
115	guez	37	rlond(1)=0.
116			rlond(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.)
117			rlond(klon)= 0.
118	guez	3
119	guez	37	cuphy = pack(cu_2d, dyn_phy)
120	guez	3
121			! Construct a mask for points in the "v" grid:
122			mask_v = .true.
123			mask_v(2:, 1) = .false.
124			mask_v(iim + 1, 2:) = .false.
125
126	guez	37	cvphy(:klon - 1) = pack(cv_2d, mask_v)
127			cvphy(klon) = cv_2d(1, jjm)
128			! (that value of "cv_2d" is used twice in "cvphy")
129	guez	3
130	guez	37	airephy = pack(aire_2d, dyn_phy)
131			CALL suphec
132	guez	38	call yoethf
133	guez	3	ENDIF
134
135			! Initialisation des entrées-sorties :
136			day_end = day_ini + nday
137			print *, "day_ini = ", day_ini
138			print *, "day_end = ", day_end
139
140	guez	5	CALL dynredem0("restart.nc", day_end, phis)
141	guez	79	CALL inithist(day_ref, annee_ref, dtvr, nqmx, t_ops = iecri * daysec, &
142	guez	56	t_wrt = iecri * daysec)
143	guez	79	CALL initdynav(day_ref, annee_ref, dtvr, nqmx, t_ops = iperiod * dtvr, &
144	guez	56	t_wrt = periodav * daysec)
145	guez	57	call init_dynzon(dt_app = dtvr * iperiod)
146	guez	3
147			! Choix des fréquences de stockage pour le hors-ligne :
148	guez	36	istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4
149	guez	3	istphy = istdyn / iphysiq
150
151	guez	115	CALL conf_guide
152
153	guez	3	! Intégration temporelle du modèle :
154	guez	128	CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)
155	guez	3
156	guez	68	close(unit_nml)
157	guez	10	call histclo
158	guez	108
159			do i = 1, nbfield
160	guez	109	call nf95_close(Ncid(i))
161	guez	108	end do
162
163	guez	10	print *, 'Simulation finished'
164			print *, 'Everything is cool'
165
166	guez	3	END PROGRAM gcm