1 | c IF (ok_mensuelNMC) THEN |
---|
2 | IF (ok_mensuel) THEN |
---|
3 | c |
---|
4 | IF(region.EQ.'3d') THEN |
---|
5 | zsto = dtime |
---|
6 | cIM 150304 BEG |
---|
7 | zout = dtime * ecrit_mth |
---|
8 | c zout = dtime * REAL(NINT(86400./dtime*ecrit_hf)) |
---|
9 | cIM 150304 END |
---|
10 | cIM 120304 BEG |
---|
11 | c zout = dtime * REAL(NINT(86400./dtime*ecrit_hf)) |
---|
12 | c zsto = zout |
---|
13 | c zsto = dtime |
---|
14 | cIM 120304 END |
---|
15 | ELSEIF(region.EQ.'1d') THEN |
---|
16 | zsto = dtime |
---|
17 | c zout = dtime |
---|
18 | zout = dtime *4. |
---|
19 | ENDIF |
---|
20 | c |
---|
21 | idayref = day_ref |
---|
22 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
---|
23 | c |
---|
24 | CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) |
---|
25 | DO i = 1, iim |
---|
26 | zx_lon(i,1) = rlon(i+1) |
---|
27 | zx_lon(i,jjmp1) = rlon(i+1) |
---|
28 | ENDDO |
---|
29 | DO ll=1,klev |
---|
30 | znivsig(ll)=float(ll) |
---|
31 | ENDDO |
---|
32 | CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) |
---|
33 | CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), |
---|
34 | c071103 |
---|
35 | c . 1,iim,1,jjmp1, |
---|
36 | . imin_debut,nbpti, |
---|
37 | . jmin_debut,nbptj, |
---|
38 | c071103 |
---|
39 | . itau_phy, zjulian, dtime, |
---|
40 | . nhori, nid_nmc) |
---|
41 | write(*,*)'Mensuel NMC ', itau_phy, zjulian |
---|
42 | CALL histvert(nid_nmc, "presnivs", "Vertical levels", "mb", |
---|
43 | . klev, presnivs, nvert) |
---|
44 | c call histvert(nid_nmc, 'sig_s', 'Niveaux sigma','-', |
---|
45 | c . klev, znivsig, nvert) |
---|
46 | c |
---|
47 | c Champs 2D: |
---|
48 | |
---|
49 | c Champs interpolles sur des niveaux de pression du NMC |
---|
50 | c IMIMIM 110304 BEG |
---|
51 | c |
---|
52 | DO k=1, nlevSTD |
---|
53 | c |
---|
54 | bb=clevSTD(k) |
---|
55 | c |
---|
56 | IF(k.GE.2) THEN |
---|
57 | aa=clevSTD(k) |
---|
58 | bb=aa(1:lnblnk1(aa)) |
---|
59 | ENDIF |
---|
60 | c |
---|
61 | CALL histdef(nid_nmc, "t"//bb, |
---|
62 | . "Temperature"//bb//"mb","K", |
---|
63 | . iim,jjmp1,nhori, 1,1,1, -99, 32, |
---|
64 | . "ave(X)", zsto,zout) |
---|
65 | c |
---|
66 | CALL histdef(nid_nmc, "phi"//bb, |
---|
67 | . "Geopotential"//bb//"mb", "m2/s2", |
---|
68 | . iim,jjmp1,nhori, 1,1,1, -99, 32, |
---|
69 | . "ave(X)", zsto,zout) |
---|
70 | c |
---|
71 | CALL histdef(nid_nmc, "q"//bb, |
---|
72 | . "Specific humidity"//bb//"mb","kg/kg", |
---|
73 | . iim,jjmp1,nhori, 1,1,1, -99, 32, |
---|
74 | . "ave(X)", zsto,zout) |
---|
75 | c |
---|
76 | CALL histdef(nid_nmc, "rh"//bb, |
---|
77 | . "Relative humidity"//bb//"mb", "%", |
---|
78 | . iim,jjmp1,nhori, 1,1,1, -99, 32, |
---|
79 | . "ave(X)", zsto,zout) |
---|
80 | c |
---|
81 | CALL histdef(nid_nmc, "u"//bb, |
---|
82 | . "Zonal wind"//bb//"mb","K", |
---|
83 | . iim,jjmp1,nhori, 1,1,1, -99, 32, |
---|
84 | . "ave(X)", zsto,zout) |
---|
85 | c |
---|
86 | CALL histdef(nid_nmc, "v"//bb, |
---|
87 | . "Meridional wind"//bb//"mb","K", |
---|
88 | . iim,jjmp1,nhori, 1,1,1, -99, 32, |
---|
89 | . "ave(X)", zsto,zout) |
---|
90 | c |
---|
91 | ENDDO !nlevSTD |
---|
92 | c IMIMIM 110304 END |
---|
93 | |
---|
94 | CALL histend(nid_nmc) |
---|
95 | c |
---|
96 | c ndex2d = 0 |
---|
97 | c |
---|
98 | ENDIF ! fin de test sur ok_mensuel |
---|
99 | |
---|