15 |
|
|
16 |
! From LMDZ4/libf/dyn3d/bilan_dyn.F, version 1.5 2005/03/16 10:12:17 |
! From LMDZ4/libf/dyn3d/bilan_dyn.F, version 1.5 2005/03/16 10:12:17 |
17 |
|
|
18 |
|
USE comgeom, ONLY: rlatv |
19 |
USE conf_gcm_m, ONLY: day_step, iperiod, periodav |
USE conf_gcm_m, ONLY: day_step, iperiod, periodav |
20 |
|
USE dimens_m, ONLY: jjm, llm |
21 |
|
USE disvert_m, ONLY: presnivs |
22 |
|
use dynetat0_m, only: day_ref, annee_ref |
23 |
USE histbeg_totreg_m, ONLY: histbeg_totreg |
USE histbeg_totreg_m, ONLY: histbeg_totreg |
24 |
USE histdef_m, ONLY: histdef |
USE histdef_m, ONLY: histdef |
25 |
USE histend_m, ONLY: histend |
USE histend_m, ONLY: histend |
26 |
USE histvert_m, ONLY: histvert |
USE histvert_m, ONLY: histvert |
|
USE calendar, ONLY: ymds2ju |
|
|
USE dimens_m, ONLY: jjm, llm |
|
|
USE disvert_m, ONLY: presnivs |
|
|
USE comgeom, ONLY: rlatv |
|
|
USE temps, ONLY: annee_ref, day_ref, itau_dyn |
|
27 |
USE nr_util, ONLY: pi |
USE nr_util, ONLY: pi |
28 |
|
USE temps, ONLY: itau_dyn |
29 |
|
USE ymds2ju_m, ONLY: ymds2ju |
30 |
|
|
31 |
real, intent(in):: dt_app |
real, intent(in):: dt_app |
32 |
|
|
48 |
|
|
49 |
integer horiid, vertiid |
integer horiid, vertiid |
50 |
real julian |
real julian |
|
integer an, dayref |
|
51 |
real rlong(jjm), rlatg(jjm) |
real rlong(jjm), rlatg(jjm) |
52 |
|
|
53 |
!----------------------------------------------------------------- |
!----------------------------------------------------------------- |
61 |
|
|
62 |
! Initialisation du fichier contenant les moyennes zonales |
! Initialisation du fichier contenant les moyennes zonales |
63 |
|
|
64 |
an = annee_ref |
CALL ymds2ju(annee_ref, 1, day_ref, 0.0, julian) |
|
dayref = day_ref |
|
|
CALL ymds2ju(an, 1, dayref, 0.0, julian) |
|
65 |
|
|
66 |
rlong = 0. |
rlong = 0. |
67 |
rlatg = rlatv * 180. / pi |
rlatg = rlatv * 180. / pi |
71 |
|
|
72 |
! Appel à histvert pour la grille verticale |
! Appel à histvert pour la grille verticale |
73 |
|
|
74 |
call histvert(fileid, 'presnivs', 'Niveaux sigma', 'mb', llm, presnivs, & |
call histvert(fileid, 'presnivs', 'Niveaux sigma', 'mb', presnivs, vertiid) |
|
vertiid) |
|
75 |
|
|
76 |
! Appels à histdef pour la définition des variables à sauvegarder |
! Appels à histdef pour la définition des variables à sauvegarder |
77 |
do iQ = 1, nQ |
do iQ = 1, nQ |