/[lmdze]/trunk/dyn3d/gcm.f
ViewVC logotype

Annotation of /trunk/dyn3d/gcm.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 40 - (hide annotations)
Tue Feb 22 13:49:36 2011 UTC (13 years, 2 months ago) by guez
Original Path: trunk/libf/dyn3d/gcm.f90
File size: 5397 byte(s)
"alpha" useless, always 0, in "exner_hyb".

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

  ViewVC Help
Powered by ViewVC 1.1.21