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

Annotation of /trunk/dyn3d/gcm.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 54 - (hide annotations)
Tue Dec 6 15:07:04 2011 UTC (12 years, 5 months ago) by guez
Original Path: trunk/libf/dyn3d/gcm.f90
File size: 5437 byte(s) 
Removed Numerical Recipes procedure "ran1". Replaced calls to "ran1"
in "inidissip" by calls to intrinsic procedures.

Split file "interface_surf.f90" into a file with a module containing
only variables, "interface_surf", and single-procedure files. Gathered
files into directory "Interface_surf".

Added argument "cdivu" to "gradiv" and "gradiv2", "cdivh" to
"divgrad2" and "divgrad", and "crot" to "nxgraro2" and
"nxgrarot". "dissip" now uses variables "cdivu", "cdivh" and "crot"
from module "inidissip_m", so it can pass them to "gradiv2",
etc. Thanks to this modification, we avoid a circular dependency
betwwen "inidissip.f90" and "gradiv2.f90", etc. The value -1. used by
"gradiv2", for instance, during computation of eigenvalues is not the
value "cdivu" computed by "inidissip".

Extracted procedure "start_inter_3d" from module "startdyn", to its
own module.

In "inidissip", unrolled loop on "ii". I find it clearer now.

Moved variables "matriceun", "matriceus", "matricevn", "matricevs",
"matrinvn" and "matrinvs" from module "parafilt" to module
"inifilr_m". Moved variables "jfiltnu", "jfiltnv", "jfiltsu",
"jfiltsv" from module "coefils" to module "inifilr_m".
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 guez 54 use inifilr_m, only: inifilr
26 guez 26 use inigeom_m, only: inigeom
27     use initdynav_m, only: initdynav
28     use inithist_m, only: inithist
29 guez 30 USE calendar, only: ioconf_calendar
30     use histcom, only: histclo
31 guez 3 use leapfrog_m, only: leapfrog
32 guez 26 use logic, only: iflag_phys
33     use paramet_m, only: ip1jm, ip1jmp1
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 38 use yoethf_m, only: yoethf
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 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 3
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     ! Choix du calendrier :
71     if (true_calendar) then
72     call ioconf_calendar('gregorian')
73     else
74     call ioconf_calendar('360d')
75     endif
76    
77     ! Lecture des fichiers "gcm.def" ou "run.def" :
78 guez 13 call read_clesphys2
79     CALL conf_gcm
80 guez 3
81     ! Initialisation des traceurs
82     ! Choix du schéma pour l'advection dans le fichier "traceur.def" ou via INCA
83 guez 23 call iniadvtrac
84 guez 3
85     ! Lecture du fichier "start.nc" :
86     CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)
87    
88 guez 36 ! On remet le calendrier à zero si demandé :
89 guez 3 if (annee_ref /= anneeref .or. day_ref /= dayref) then
90     print *, 'Attention : les dates initiales lues dans le fichier ' // &
91     '"start" ne correspondent pas à celles lues dans "gcm.def".'
92 guez 15 if (raz_date) then
93 guez 3 print *, 'On réinitialise à la date lue dans "gcm.def".'
94     annee_ref = anneeref
95     day_ref = dayref
96     day_ini = dayref
97     itau_dyn = 0
98     time_0 = 0.
99 guez 15 else
100     print *, 'On garde les dates du fichier "start".'
101 guez 3 endif
102     ELSE
103 guez 15 raz_date = .false.
104 guez 3 endif
105    
106 guez 36 ! On recalcule éventuellement le pas de temps :
107 guez 3 zdtvr = daysec / REAL(day_step)
108     IF (dtvr /= zdtvr) THEN
109     print *, 'Warning: the time steps in the ".def" file and in ' // &
110     '"start.nc" are different'
111     print *, 'dtvr (from "start.nc") = ', dtvr
112     print *, 'zdtvr (from ".def") = ', zdtvr
113     print *, 'Using the value from the ".def" file.'
114     dtvr = zdtvr
115     ENDIF
116 guez 33
117 guez 36 CALL iniconst
118     CALL inigeom ! initialisation de la géometrie
119     CALL inifilr ! initialisation du filtre
120 guez 27 CALL inidissip
121 guez 3 call init_dyn_phy
122    
123     ! Initialisation de la physique :
124     IF (iflag_phys == 1) THEN
125 guez 37 rlatd(1)=rlatu(1)
126     rlatd(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.)
127     rlatd(klon)= rlatu(jjm + 1)
128 guez 3
129 guez 37 rlond(1)=0.
130     rlond(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.)
131     rlond(klon)= 0.
132 guez 3
133 guez 37 cuphy = pack(cu_2d, dyn_phy)
134 guez 3
135     ! Construct a mask for points in the "v" grid:
136     mask_v = .true.
137     mask_v(2:, 1) = .false.
138     mask_v(iim + 1, 2:) = .false.
139    
140 guez 37 cvphy(:klon - 1) = pack(cv_2d, mask_v)
141     cvphy(klon) = cv_2d(1, jjm)
142     ! (that value of "cv_2d" is used twice in "cvphy")
143 guez 3
144 guez 37 airephy = pack(aire_2d, dyn_phy)
145     CALL suphec
146 guez 38 call yoethf
147 guez 3 ENDIF
148    
149     ! Initialisation des entrées-sorties :
150     day_end = day_ini + nday
151     print *, "day_ini = ", day_ini
152     print *, "day_end = ", day_end
153    
154 guez 5 CALL dynredem0("restart.nc", day_end, phis)
155 guez 3 CALL inithist(day_ref, annee_ref, zdtvr, nqmx, histid, histvid, &
156 guez 27 t_ops = iecri * daysec, t_wrt = iecri * daysec)
157 guez 3 CALL initdynav(day_ref, annee_ref, zdtvr, nqmx, histaveid, &
158 guez 27 t_ops = iperiod * zdtvr, t_wrt = periodav * daysec)
159 guez 3
160     ! Choix des fréquences de stockage pour le hors-ligne :
161 guez 36 istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4
162 guez 3 istphy = istdyn / iphysiq
163    
164     ! Intégration temporelle du modèle :
165 guez 23 CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)
166 guez 3
167 guez 10 call histclo
168     print *, 'Simulation finished'
169     print *, 'Everything is cool'
170    
171 guez 3 END PROGRAM gcm
  ViewVC Help
Powered by ViewVC 1.1.21