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ées verticales hybrides, 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".
  ! Pour Van-Leer plus vapeur d'eau saturée : iadv(1)=4
  ! Pour Van-Leer : iadv=10

  use clesphys2, only: read_clesphys2
  use com_io_dyn, only: histid, histvid, histaveid
  use comconst, only: daysec, cpp, dtvr, g, rad, r, initialize
  use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv
  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
  use dimens_m, only: iim, jjm, llm, nqmx
  use dimphy, only: klon
  use dynetat0_m, only: dynetat0, day_ini
  use dynredem0_m, only: dynredem0
  use grid_change, only: dyn_phy, init_dyn_phy
  use iniadvtrac_m, only: iniadvtrac
  use inidissip_m, only: inidissip
  use inigeom_m, only: inigeom
  use initdynav_m, only: initdynav
  use inithist_m, only: inithist
  USE calendar, only: ioconf_calendar
  use histcom, only: histclo
  use leapfrog_m, only: leapfrog
  use logic, only: iflag_phys
  use paramet_m, only: ip1jm, ip1jmp1
  use suphec_m, only: suphec
  use temps, only: day_ref, annee_ref, day_end, itau_dyn
  use tracstoke, only: istdyn, istphy

  IMPLICIT NONE

  REAL zdtvr ! time step for dynamics, in s

  ! Variables dynamiques :
  REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants
  REAL teta(ip1jmp1, llm) ! température potentielle 
  REAL q(ip1jmp1, llm, nqmx) ! champs advectés
  REAL ps(ip1jmp1) ! pression au sol (Pa)

  REAL masse(ip1jmp1, 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)

  namelist /main_nml/true_calendar

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

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

  ! Initialisations:
  call initialize

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

  ! Lecture des fichiers "gcm.def" ou "run.def" :
  call read_clesphys2
  CALL conf_gcm

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

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

  ! On remet le calendrier à zero si demandé :
  if (annee_ref /= anneeref .or. day_ref /= dayref) then
     print *, 'Attention : les dates initiales lues dans le fichier ' // &
          '"start" ne correspondent pas à celles lues dans "gcm.def".'
     if (raz_date) then
        print *, 'On réinitialise à la date lue dans "gcm.def".'
        annee_ref = anneeref
        day_ref = dayref
        day_ini = dayref
        itau_dyn = 0
        time_0 = 0.
     else
        print *, 'On garde les dates du fichier "start".'
     endif
  ELSE
     raz_date = .false.
  endif

  ! On recalcule éventuellement le pas de temps :
  zdtvr = daysec / REAL(day_step)
  IF (dtvr /= zdtvr) THEN
     print *, 'Warning: the time steps in the ".def" file and in ' // &
          '"start.nc" are different'
     print *, 'dtvr (from "start.nc") = ', dtvr
     print *, 'zdtvr (from ".def") = ', zdtvr
     print *, 'Using the value from the ".def" file.'
     dtvr = zdtvr
  ENDIF

  CALL iniconst 
  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
  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, zdtvr, nqmx, histid, histvid, &
       t_ops = iecri * daysec, t_wrt = iecri * daysec)
  CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, &
       t_ops = iperiod * zdtvr, t_wrt = periodav * daysec)

  ! 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     

  ! Intégration temporelle du modèle :
  CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)

  call histclo
  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	! General circulation model of LMD
6			! Avec coordonnées verticales hybrides, avec nouveaux opérateurs de
7	guez	36	! dissipation "*" (gradiv2, divgrad2, nxgraro2)
8	guez	3	! Possibilité de choisir le schéma pour l'advection de "q", en
9			! modifiant "iadv" dans "traceur.def".
10			! Pour Van-Leer plus vapeur d'eau saturée : iadv(1)=4
11			! Pour Van-Leer : iadv=10
12
13	guez	26	use clesphys2, only: read_clesphys2
14			use com_io_dyn, only: histid, histvid, histaveid
15	guez	23	use comconst, only: daysec, cpp, dtvr, g, rad, r, initialize
16	guez	26	use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv
17	guez	37	use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond
18	guez	3	use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, &
19			nday, raz_date, periodav, conf_gcm
20	guez	26	use dimens_m, only: iim, jjm, llm, nqmx
21			use dimphy, only: klon
22	guez	25	use dynetat0_m, only: dynetat0, day_ini
23	guez	26	use dynredem0_m, only: dynredem0
24	guez	3	use grid_change, only: dyn_phy, init_dyn_phy
25	guez	18	use iniadvtrac_m, only: iniadvtrac
26	guez	26	use inidissip_m, only: inidissip
27			use inigeom_m, only: inigeom
28			use initdynav_m, only: initdynav
29			use inithist_m, only: inithist
30	guez	30	USE calendar, only: ioconf_calendar
31			use histcom, only: histclo
32	guez	3	use leapfrog_m, only: leapfrog
33	guez	26	use logic, only: iflag_phys
34			use paramet_m, only: ip1jm, ip1jmp1
35	guez	37	use suphec_m, only: suphec
36	guez	26	use temps, only: day_ref, annee_ref, day_end, itau_dyn
37			use tracstoke, only: istdyn, istphy
38	guez	3
39			IMPLICIT NONE
40
41			REAL zdtvr ! time step for dynamics, in s
42
43			! Variables dynamiques :
44			REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vents covariants
45			REAL teta(ip1jmp1, llm) ! température potentielle
46			REAL q(ip1jmp1, llm, nqmx) ! champs advectés
47			REAL ps(ip1jmp1) ! pression au sol (Pa)
48
49			REAL masse(ip1jmp1, llm) ! masse d'air
50	guez	7	REAL phis(iim + 1, jjm + 1) ! géopotentiel au sol
51	guez	3
52			! Variables pour le fichier histoire :
53	guez	25	REAL time_0 ! time in day, as a fraction of day, in [0, 1[
54	guez	3
55			! Calendrier :
56			LOGICAL:: true_calendar = .false. ! default value
57
58			logical mask_v(iim + 1, jjm)
59			! (mask for points in the "v" grid, first index is for longitude,
60			! second index is for latitude)
61
62			namelist /main_nml/true_calendar
63
64			!------------------------------------------------------------
65
66			print *, "Enter namelist 'main_nml'."
67			read (unit=*, nml=main_nml)
68			write(unit=*, nml=main_nml)
69
70			! Initialisations:
71			call initialize
72
73			! Choix du calendrier :
74			if (true_calendar) then
75			call ioconf_calendar('gregorian')
76			else
77			call ioconf_calendar('360d')
78			endif
79
80			! Lecture des fichiers "gcm.def" ou "run.def" :
81	guez	13	call read_clesphys2
82			CALL conf_gcm
83	guez	3
84			! Initialisation des traceurs
85			! Choix du schéma pour l'advection dans le fichier "traceur.def" ou via INCA
86	guez	23	call iniadvtrac
87	guez	3
88			! Lecture du fichier "start.nc" :
89			CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)
90
91	guez	36	! On remet le calendrier à zero si demandé :
92	guez	3	if (annee_ref /= anneeref .or. day_ref /= dayref) then
93			print *, 'Attention : les dates initiales lues dans le fichier ' // &
94			'"start" ne correspondent pas à celles lues dans "gcm.def".'
95	guez	15	if (raz_date) then
96	guez	3	print *, 'On réinitialise à la date lue dans "gcm.def".'
97			annee_ref = anneeref
98			day_ref = dayref
99			day_ini = dayref
100			itau_dyn = 0
101			time_0 = 0.
102	guez	15	else
103			print *, 'On garde les dates du fichier "start".'
104	guez	3	endif
105			ELSE
106	guez	15	raz_date = .false.
107	guez	3	endif
108
109	guez	36	! On recalcule éventuellement le pas de temps :
110	guez	3	zdtvr = daysec / REAL(day_step)
111			IF (dtvr /= zdtvr) THEN
112			print *, 'Warning: the time steps in the ".def" file and in ' // &
113			'"start.nc" are different'
114			print *, 'dtvr (from "start.nc") = ', dtvr
115			print *, 'zdtvr (from ".def") = ', zdtvr
116			print *, 'Using the value from the ".def" file.'
117			dtvr = zdtvr
118			ENDIF
119	guez	33
120	guez	36	CALL iniconst
121			CALL inigeom ! initialisation de la géometrie
122			CALL inifilr ! initialisation du filtre
123	guez	27	CALL inidissip
124	guez	3	call init_dyn_phy
125
126			! Initialisation de la physique :
127			IF (iflag_phys == 1) THEN
128	guez	37	rlatd(1)=rlatu(1)
129			rlatd(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.)
130			rlatd(klon)= rlatu(jjm + 1)
131	guez	3
132	guez	37	rlond(1)=0.
133			rlond(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.)
134			rlond(klon)= 0.
135	guez	3
136	guez	37	cuphy = pack(cu_2d, dyn_phy)
137	guez	3
138			! Construct a mask for points in the "v" grid:
139			mask_v = .true.
140			mask_v(2:, 1) = .false.
141			mask_v(iim + 1, 2:) = .false.
142
143	guez	37	cvphy(:klon - 1) = pack(cv_2d, mask_v)
144			cvphy(klon) = cv_2d(1, jjm)
145			! (that value of "cv_2d" is used twice in "cvphy")
146	guez	3
147	guez	37	airephy = pack(aire_2d, dyn_phy)
148			CALL suphec
149	guez	3	ENDIF
150
151			! Initialisation des entrées-sorties :
152			day_end = day_ini + nday
153			print *, "day_ini = ", day_ini
154			print *, "day_end = ", day_end
155
156	guez	5	CALL dynredem0("restart.nc", day_end, phis)
157	guez	3	CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, &
158	guez	27	t_ops = iecri * daysec, t_wrt = iecri * daysec)
159	guez	3	CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, &
160	guez	27	t_ops = iperiod * zdtvr, t_wrt = periodav * daysec)
161	guez	3
162			! Choix des fréquences de stockage pour le hors-ligne :
163	guez	36	istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4
164	guez	3	istphy = istdyn / iphysiq
165
166			! Intégration temporelle du modèle :
167	guez	23	CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)
168	guez	3
169	guez	10	call histclo
170			print *, 'Simulation finished'
171			print *, 'Everything is cool'
172
173	guez	3	END PROGRAM gcm