[3331] | 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) |
---|