[524] | 1 | ! |
---|
[1403] | 2 | ! $Id: write_histmthNMC.h 1665 2012-10-09 13:35:26Z emillour $ |
---|
[524] | 3 | ! |
---|
[1374] | 4 | IF (ok_histNMC(1)) THEN |
---|
[524] | 5 | c |
---|
[644] | 6 | ndex3d = 0 |
---|
[1665] | 7 | itau_w = itau_phy + itap + start_time * day_step / iphysiq |
---|
[644] | 8 | ccc |
---|
| 9 | c Champs interpolles sur des niveaux de pression du NMC |
---|
[524] | 10 | c |
---|
[1352] | 11 | c PARAMETER(nout=3) |
---|
| 12 | c nout=1 : in=pdtphys, out=mth |
---|
| 13 | c nout=2 : in=pdtphys, out=day |
---|
| 14 | c nout=3 : in=pdtphys, out=hf |
---|
[644] | 15 | ccc |
---|
[1539] | 16 | CALL histwrite_phy(nid_mthnmc,lNMC,"tnondef",itau_w, |
---|
| 17 | $tnondef(:,:,1)) |
---|
[524] | 18 | c |
---|
[1539] | 19 | CALL histwrite_phy(nid_mthnmc,lNMC,"ta",itau_w, |
---|
| 20 | $twriteSTD(:,:,1)) |
---|
[524] | 21 | c |
---|
[1539] | 22 | CALL histwrite_phy(nid_mthnmc,lNMC,"zg",itau_w, |
---|
| 23 | $phiwriteSTD(:,:,1)) |
---|
[524] | 24 | c |
---|
[1539] | 25 | CALL histwrite_phy(nid_mthnmc,lNMC,"hus",itau_w, |
---|
| 26 | $qwriteSTD(:,:,1)) |
---|
[524] | 27 | c |
---|
[1539] | 28 | CALL histwrite_phy(nid_mthnmc,lNMC,"hur",itau_w, |
---|
| 29 | $rhwriteSTD(:,:,1)) |
---|
[644] | 30 | c |
---|
[1539] | 31 | CALL histwrite_phy(nid_mthnmc,lNMC,"ua",itau_w, |
---|
| 32 | $uwriteSTD(:,:,1)) |
---|
[644] | 33 | c |
---|
[1539] | 34 | CALL histwrite_phy(nid_mthnmc,lNMC,"va",itau_w, |
---|
| 35 | $vwriteSTD(:,:,1)) |
---|
[644] | 36 | c |
---|
[1539] | 37 | CALL histwrite_phy(nid_mthnmc,lNMC,"wap",itau_w, |
---|
| 38 | $wwriteSTD(:,:,1)) |
---|
[644] | 39 | c |
---|
| 40 | DO k=1, nlevSTD |
---|
| 41 | DO i=1, klon |
---|
[1352] | 42 | IF(tnondef(i,k,1).NE.missing_val) THEN |
---|
[1398] | 43 | zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,1))/freq_moyNMC(1) |
---|
[644] | 44 | ELSE |
---|
[1352] | 45 | zx_tmp_fiNC(i,k) = missing_val |
---|
[524] | 46 | ENDIF |
---|
[644] | 47 | ENDDO |
---|
| 48 | ENDDO !k=1, nlevSTD |
---|
[524] | 49 | c |
---|
[1539] | 50 | CALL histwrite_phy(nid_mthnmc,lNMC,"psbg",itau_w, |
---|
| 51 | $zx_tmp_fiNC) |
---|
[524] | 52 | c |
---|
[1539] | 53 | CALL histwrite_phy(nid_mthnmc,lNMC,"uv",itau_w, |
---|
| 54 | $uvsumSTD(:,:,1)) |
---|
[524] | 55 | c |
---|
[1539] | 56 | CALL histwrite_phy(nid_mthnmc,lNMC,"vq",itau_w, |
---|
| 57 | $vqsumSTD(:,:,1)) |
---|
[524] | 58 | c |
---|
[1539] | 59 | CALL histwrite_phy(nid_mthnmc,lNMC,"vT",itau_w, |
---|
| 60 | $vTsumSTD(:,:,1)) |
---|
[524] | 61 | c |
---|
[1539] | 62 | CALL histwrite_phy(nid_mthnmc,lNMC,"wq",itau_w, |
---|
| 63 | $wqsumSTD(:,:,1)) |
---|
[524] | 64 | c |
---|
[1539] | 65 | CALL histwrite_phy(nid_mthnmc,lNMC,"vphi",itau_w, |
---|
| 66 | $vphisumSTD(:,:,1)) |
---|
[524] | 67 | c |
---|
[1539] | 68 | CALL histwrite_phy(nid_mthnmc,lNMC,"wT",itau_w, |
---|
| 69 | $wTsumSTD(:,:,1)) |
---|
[524] | 70 | c |
---|
[1539] | 71 | CALL histwrite_phy(nid_mthnmc,lNMC,"uxu",itau_w, |
---|
| 72 | $u2sumSTD(:,:,1)) |
---|
[644] | 73 | c |
---|
[1539] | 74 | CALL histwrite_phy(nid_mthnmc,lNMC,"vxv",itau_w, |
---|
| 75 | $v2sumSTD(:,:,1)) |
---|
[644] | 76 | c |
---|
[1539] | 77 | CALL histwrite_phy(nid_mthnmc,lNMC,"TxT",itau_w, |
---|
| 78 | $T2sumSTD(:,:,1)) |
---|
[644] | 79 | c |
---|
[1398] | 80 | DO k=1, nlevSTD |
---|
| 81 | DO i=1, klon |
---|
| 82 | IF(O3sumSTD(i,k,1).NE.missing_val) THEN |
---|
| 83 | zx_tmp_fiNC(i,k) = O3sumSTD(i,k,1) * 1.e+9 |
---|
| 84 | ELSE |
---|
| 85 | zx_tmp_fiNC(i,k) = missing_val |
---|
| 86 | ENDIF |
---|
| 87 | ENDDO |
---|
| 88 | ENDDO !k=1, nlevSTD |
---|
[1539] | 89 | CALL histwrite_phy(nid_mthnmc,lNMC,"tro3",itau_w, |
---|
[1398] | 90 | $ zx_tmp_fiNC) |
---|
| 91 | c |
---|
| 92 | if (read_climoz == 2) THEN |
---|
| 93 | DO k=1, nlevSTD |
---|
| 94 | DO i=1, klon |
---|
| 95 | IF(O3daysumSTD(i,k,1).NE.missing_val) THEN |
---|
| 96 | zx_tmp_fiNC(i,k) = O3daysumSTD(i,k,1) * 1.e+9 |
---|
| 97 | ELSE |
---|
| 98 | zx_tmp_fiNC(i,k) = missing_val |
---|
| 99 | ENDIF |
---|
| 100 | ENDDO |
---|
| 101 | ENDDO !k=1, nlevSTD |
---|
| 102 | c |
---|
[1539] | 103 | CALL histwrite_phy(nid_mthnmc,lNMC,"tro3_daylight", |
---|
| 104 | $itau_w, zx_tmp_fiNC) |
---|
[1398] | 105 | endif |
---|
| 106 | c |
---|
[524] | 107 | if (ok_sync) then |
---|
[766] | 108 | c$OMP MASTER |
---|
[1352] | 109 | call histsync(nid_mthnmc) |
---|
[766] | 110 | c$OMP END MASTER |
---|
[524] | 111 | endif |
---|
[1374] | 112 | c |
---|
| 113 | ENDIF !(ok_histNMC(1)) THEN |
---|