23 |
|
|
24 |
! L. Fairhead, 03/99 |
! L. Fairhead, 03/99 |
25 |
|
|
26 |
USE calendar |
USE comconst |
27 |
|
use conf_gcm_m |
28 |
|
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 |
USE histbeg_totreg_m, ONLY : histbeg_totreg |
32 |
USE histdef_m, ONLY : histdef |
USE histdef_m, ONLY : histdef |
33 |
USE histend_m, ONLY : histend |
USE histend_m, ONLY : histend |
34 |
use histhori_regular_m, only: histhori_regular |
use histhori_regular_m, only: histhori_regular |
35 |
use histsync_m, only: histsync |
use histsync_m, only: histsync |
36 |
USE histvert_m, ONLY : histvert |
USE histvert_m, ONLY : histvert |
|
USE dimens_m |
|
|
USE paramet_m |
|
|
USE comconst |
|
|
USE disvert_m |
|
|
use conf_gcm_m |
|
|
USE comgeom |
|
|
USE serre |
|
|
USE temps, ONLY : annee_ref, day_ref, itau_dyn |
|
|
USE ener |
|
37 |
USE nr_util, ONLY : pi |
USE nr_util, ONLY : pi |
38 |
|
USE paramet_m |
39 |
|
USE temps, ONLY : itau_dyn |
40 |
|
use ymds2ju_m, only: ymds2ju |
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 |
!--------------------------------------------------------- |
!--------------------------------------------------------- |
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', nivsigs, & |
CALL histvert(fileid, 'sig_s', 'Niveaux sigma', 'sigma_level', & |
115 |
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', nivsigs, & |
CALL histvert(filevid, 'sig_s', 'Niveaux sigma', 'sigma_level', & |
118 |
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', nivd, & |
CALL histvert(filedid, 'sig_s', 'Niveaux sigma', 'sigma_level', nivd, & |