1 | !IM for NMC files |
---|
2 | ! real twriteSTD(klon,nlevSTD,nfiles) |
---|
3 | ! real qwriteSTD(klon,nlevSTD,nfiles) |
---|
4 | ! real rhwriteSTD(klon,nlevSTD,nfiles) |
---|
5 | ! real phiwriteSTD(klon,nlevSTD,nfiles) |
---|
6 | ! real uwriteSTD(klon,nlevSTD,nfiles) |
---|
7 | ! real vwriteSTD(klon,nlevSTD,nfiles) |
---|
8 | ! real wwriteSTD(klon,nlevSTD,nfiles) |
---|
9 | |
---|
10 | real twriteSTD3(klon,nlevSTD3) |
---|
11 | real qwriteSTD3(klon,nlevSTD3) |
---|
12 | real rhwriteSTD3(klon,nlevSTD3) |
---|
13 | real phiwriteSTD3(klon,nlevSTD3) |
---|
14 | real uwriteSTD3(klon,nlevSTD3) |
---|
15 | real vwriteSTD3(klon,nlevSTD3) |
---|
16 | real wwriteSTD3(klon,nlevSTD3) |
---|
17 | |
---|
18 | real tnondefSTD8(klon,nlevSTD8) |
---|
19 | real twriteSTD8(klon,nlevSTD8) |
---|
20 | real qwriteSTD8(klon,nlevSTD8) |
---|
21 | real rhwriteSTD8(klon,nlevSTD8) |
---|
22 | real phiwriteSTD8(klon,nlevSTD8) |
---|
23 | real uwriteSTD8(klon,nlevSTD8) |
---|
24 | real vwriteSTD8(klon,nlevSTD8) |
---|
25 | real wwriteSTD8(klon,nlevSTD8) |
---|
26 | |
---|
27 | real, save :: rlevSTD(nlevSTD) |
---|
28 | DATA rlevSTD/100000., 92500., 85000., 70000., & |
---|
29 | 60000., 50000., 40000., 30000., 25000., 20000., & |
---|
30 | 15000., 10000., 7000., 5000., 3000., 2000., 1000./ |
---|
31 | !$OMP THREADPRIVATE(rlevstd) |
---|
32 | |
---|
33 | CHARACTER*4, SAVE :: clevSTD(nlevSTD) |
---|
34 | DATA clevSTD/'1000','925 ','850 ','700 ','600 ', & |
---|
35 | '500 ','400 ','300 ','250 ','200 ','150 ','100 ', & |
---|
36 | '70 ','50 ','30 ','20 ','10 '/ |
---|
37 | !$OMP THREADPRIVATE(clevSTD) |
---|
38 | |
---|
39 | real, save :: rlevSTD3(nlevSTD3) |
---|
40 | DATA rlevSTD3/85000., 50000., 25000./ |
---|
41 | !$OMP THREADPRIVATE(rlevSTD3) |
---|
42 | |
---|
43 | real, save :: rlevSTD8(nlevSTD8) |
---|
44 | DATA rlevSTD8/100000., 85000., 70000., 50000., 25000., 10000., & |
---|
45 | 5000., 1000./ |
---|
46 | !$OMP THREADPRIVATE(rlevSTD8) |
---|
47 | ! |
---|
48 | REAL geo500(klon) |
---|
49 | |
---|
50 | ! nout : niveau de output des variables a une pression donnee |
---|
51 | logical oknondef(klon,nlevSTD,nout) |
---|
52 | ! |
---|
53 | ! les produits uvSTD, vqSTD, .., T2STD sont calcules |
---|
54 | ! a partir des valeurs instantannees toutes les 6 h |
---|
55 | ! qui sont moyennees sur le mois |
---|
56 | |
---|
57 | REAL zx_tmp_fiNC(klon,nlevSTD) |
---|
58 | |
---|
59 | ! REAL missing_val |
---|
60 | REAL, SAVE :: freq_moyNMC(nout) |
---|
61 | !$OMP THREADPRIVATE(freq_moyNMC) |
---|