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

Annotation of /trunk/dyn3d/gcm.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 129 - (hide annotations)
Fri Feb 13 18:22:38 2015 UTC (9 years, 3 months ago) by guez
File size: 4163 byte(s) 
Removed arguments day0, anne0 of procedures initdynav and
inithist. Use directly day_ref, annee_ref instead.

Moved variables annee_ref, day_ref of module temps to module
dynetat0_m. Deleted variables dayref and anneeref of module conf_gcm_m
and removed them from namelist conf_gcm_nml. These variables were
troubling intermediary on the way to annee_ref and day_ref. Gave as
default values to annee_ref and day_ref the default values of dayref
and anneeref. Moved the test on raz_date from main unit gcm to
procedure dynetat0. Created namelist dynetat0_nml. Read annee_ref and
day_ref from standard input in dynetat0 and redefine them from
"start.nc" if not raz_date. Rationale: 1 - Choose the best programming
from the point of view of program gcm only, forgetting program ce0l. 2
- The normal case is to define annee_ref and day_ref from "start.nc"
so put them in module dynetat0_m rather than in conf_gcm_m. 3 - Try to
always read the same namelists in the same order regardless of choices
in previous namelists. Downsides: 1 -We now need the file "dynetat0.f"
for the program ce0l, because dynetat0_m is used by dynredem0. 2 - We
need to define annee_ref and day_ref from procedure etat0.

Removed useless variable day_end of module temps.
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 129 use conf_gcm_m, only: day_step, iperiod, iecri, iphysiq, nday, periodav, &
15     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 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 43 REAL teta(iim + 1, jjm + 1, 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 55 REAL masse(iim + 1, jjm + 1, llm) ! masse d'air
48 guez 7 REAL phis(iim + 1, jjm + 1) ! géopotentiel au sol
49 guez 3
50     ! Calendrier :
51     LOGICAL:: true_calendar = .false. ! default value
52    
53     logical mask_v(iim + 1, jjm)
54     ! (mask for points in the "v" grid, first index is for longitude,
55     ! second index is for latitude)
56    
57 guez 108 integer i
58    
59 guez 3 namelist /main_nml/true_calendar
60    
61     !------------------------------------------------------------
62    
63 guez 57 call new_unit(unit_nml)
64 guez 79 open(unit_nml, file="used_namelists.txt", status="replace", action="write")
65 guez 57
66     CALL conf_gcm
67    
68 guez 3 print *, "Enter namelist 'main_nml'."
69     read (unit=*, nml=main_nml)
70 guez 57 write(unit_nml, nml=main_nml)
71 guez 3
72     ! Choix du calendrier :
73     if (true_calendar) then
74     call ioconf_calendar('gregorian')
75     else
76     call ioconf_calendar('360d')
77     endif
78    
79 guez 23 call iniadvtrac
80 guez 79 CALL iniconst
81 guez 128 CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)
82 guez 79 CALL disvert
83 guez 36 CALL inigeom ! initialisation de la géometrie
84     CALL inifilr ! initialisation du filtre
85 guez 27 CALL inidissip
86 guez 3 call init_dyn_phy
87    
88     ! Initialisation de la physique :
89     IF (iflag_phys == 1) THEN
90 guez 37 rlatd(1)=rlatu(1)
91     rlatd(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.)
92     rlatd(klon)= rlatu(jjm + 1)
93 guez 3
94 guez 37 rlond(1)=0.
95     rlond(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.)
96     rlond(klon)= 0.
97 guez 3
98 guez 37 cuphy = pack(cu_2d, dyn_phy)
99 guez 3
100     ! Construct a mask for points in the "v" grid:
101     mask_v = .true.
102     mask_v(2:, 1) = .false.
103     mask_v(iim + 1, 2:) = .false.
104    
105 guez 37 cvphy(:klon - 1) = pack(cv_2d, mask_v)
106     cvphy(klon) = cv_2d(1, jjm)
107     ! (that value of "cv_2d" is used twice in "cvphy")
108 guez 3
109 guez 37 airephy = pack(aire_2d, dyn_phy)
110     CALL suphec
111 guez 38 call yoethf
112 guez 3 ENDIF
113    
114     ! Initialisation des entrées-sorties :
115 guez 129 CALL dynredem0("restart.nc", day_ini + nday, phis)
116     CALL inithist(dtvr, nqmx, t_ops = iecri * daysec, t_wrt = iecri * daysec)
117     CALL initdynav(dtvr, nqmx, t_ops = iperiod * dtvr, t_wrt = periodav * daysec)
118 guez 57 call init_dynzon(dt_app = dtvr * iperiod)
119 guez 3
120     ! Choix des fréquences de stockage pour le hors-ligne :
121 guez 36 istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4
122 guez 3 istphy = istdyn / iphysiq
123    
124 guez 115 CALL conf_guide
125 guez 128 CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)
126 guez 3
127 guez 68 close(unit_nml)
128 guez 10 call histclo
129 guez 108
130     do i = 1, nbfield
131 guez 109 call nf95_close(Ncid(i))
132 guez 108 end do
133    
134 guez 10 print *, 'Simulation finished'
135     print *, 'Everything is cool'
136    
137 guez 3 END PROGRAM gcm
  ViewVC Help
Powered by ViewVC 1.1.21