23 |
|
|
24 |
! L. Fairhead, 03/99 |
! L. Fairhead, 03/99 |
25 |
|
|
26 |
USE calendar |
USE comconst |
27 |
USE histcom |
use conf_gcm_m |
28 |
USE dimens_m |
USE dimens_m |
29 |
|
USE disvert_m |
30 |
|
use dynetat0_m, only: day_ref, annee_ref, rlonu, rlatu, rlonv, rlatv |
31 |
|
USE histbeg_totreg_m, ONLY : histbeg_totreg |
32 |
|
USE histdef_m, ONLY : histdef |
33 |
|
USE histend_m, ONLY : histend |
34 |
|
use histhori_regular_m, only: histhori_regular |
35 |
|
use histsync_m, only: histsync |
36 |
|
USE histvert_m, ONLY : histvert |
37 |
|
USE nr_util, ONLY : pi |
38 |
USE paramet_m |
USE paramet_m |
39 |
USE comconst |
USE temps, ONLY : itau_dyn |
40 |
USE comvert |
use ymds2ju_m, only: ymds2ju |
|
USE logic |
|
|
USE comgeom |
|
|
USE serre |
|
|
USE temps, ONLY : annee_ref, day_ref, itau_dyn |
|
|
USE ener |
|
41 |
|
|
42 |
IMPLICIT NONE |
IMPLICIT NONE |
43 |
|
|
51 |
|
|
52 |
! Variables locales |
! Variables locales |
53 |
REAL zjulian |
REAL zjulian |
54 |
CHARACTER*3 str |
CHARACTER(len=3) str |
55 |
CHARACTER*10 ctrac |
CHARACTER(len=10) ctrac |
56 |
INTEGER iq |
INTEGER iq |
57 |
REAL rlong(iip1, jjp1), rlat(iip1, jjp1) |
REAL rlong(iip1, jjp1), rlat(iip1, jjp1) |
58 |
INTEGER uhoriid, vhoriid, thoriid, zvertiid, dhoriid, dvertiid |
INTEGER uhoriid, vhoriid, thoriid, zvertiid, dhoriid, dvertiid |
59 |
INTEGER ii, jj |
INTEGER ii, jj, l |
|
INTEGER zan, idayref |
|
60 |
LOGICAL ok_sync |
LOGICAL ok_sync |
61 |
|
|
62 |
!--------------------------------------------------------- |
!--------------------------------------------------------- |
63 |
|
|
64 |
! Initialisations |
! Initialisations |
|
|
|
|
pi = 4.*atan(1.) |
|
65 |
str = 'q ' |
str = 'q ' |
66 |
ctrac = 'traceur ' |
ctrac = 'traceur ' |
67 |
ok_sync = .TRUE. |
ok_sync = .TRUE. |
68 |
|
|
69 |
! Appel a histbeg: creation du fichier netcdf et initialisations diverses |
! Appel a histbeg: creation du fichier netcdf et initialisations diverses |
70 |
|
|
71 |
zan = annee_ref |
CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) |
|
idayref = day_ref |
|
|
CALL ymds2ju(zan, 1, idayref, 0.0, zjulian) |
|
72 |
|
|
73 |
DO jj = 1, jjp1 |
DO jj = 1, jjp1 |
74 |
DO ii = 1, iip1 |
DO ii = 1, iip1 |
111 |
|
|
112 |
! Appel a histvert pour la grille verticale |
! Appel a histvert pour la grille verticale |
113 |
|
|
114 |
CALL histvert(fileid, 'sig_s', 'Niveaux sigma', 'sigma_level', llm, & |
CALL histvert(fileid, 'sig_s', 'Niveaux sigma', 'sigma_level', & |
115 |
nivsigs, zvertiid) |
(/(real(l), l = 1, llm)/), zvertiid) |
116 |
! Pour le fichier V |
! Pour le fichier V |
117 |
CALL histvert(filevid, 'sig_s', 'Niveaux sigma', 'sigma_level', llm, & |
CALL histvert(filevid, 'sig_s', 'Niveaux sigma', 'sigma_level', & |
118 |
nivsigs, zvertiid) |
(/(real(l), l = 1, llm)/), zvertiid) |
119 |
! pour le fichier def |
! pour le fichier def |
120 |
nivd(1) = 1 |
nivd(1) = 1 |
121 |
CALL histvert(filedid, 'sig_s', 'Niveaux sigma', 'sigma_level', 1, nivd, & |
CALL histvert(filedid, 'sig_s', 'Niveaux sigma', 'sigma_level', nivd, & |
122 |
dvertiid) |
dvertiid) |
123 |
|
|
124 |
! Appels a histdef pour la definition des variables a sauvegarder |
! Appels a histdef pour la definition des variables a sauvegarder |