! ! $Header$ ! IF (ok_histNMC(3)) THEN c ndex3d = 0 itau_w = itau_phy + itap + start_time * day_step / iphysiq ccc c Champs interpolles sur des niveaux de pression du NMC c c PARAMETER(nout=3) c nout=1 : in=pdtphys, out=mth c nout=2 : in=pdtphys, out=day c nout=3 : in=pdtphys, out=hf ccc CALL histwrite_phy(nid_hfnmc,lNMC,"tnondef",itau_w, $tnondef(:,:,3)) c CALL histwrite_phy(nid_hfnmc,lNMC,"ta",itau_w, $twriteSTD3) c CALL histwrite_phy(nid_hfnmc,lNMC,"zg",itau_w, $phiwriteSTD3) c CALL histwrite_phy(nid_hfnmc,lNMC,"hus",itau_w, $qwriteSTD3) c CALL histwrite_phy(nid_hfnmc,lNMC,"hur",itau_w, $rhwriteSTD3) c CALL histwrite_phy(nid_hfnmc,lNMC,"ua",itau_w, $uwriteSTD3) c CALL histwrite_phy(nid_hfnmc,lNMC,"va",itau_w, $vwriteSTD3) c CALL histwrite_phy(nid_hfnmc,lNMC,"wap",itau_w, $wwriteSTD3) c IF (1.EQ.0) THEN c DO k=1, nlevSTD DO i=1, klon IF(tnondef(i,k,3).NE.missing_val) THEN zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,3))/freq_moyNMC(3) ELSE zx_tmp_fiNC(i,k) = missing_val ENDIF ENDDO ENDDO !k=1, nlevSTD c CALL histwrite_phy(nid_hfnmc,lNMC,"psbg",itau_w, $zx_tmp_fiNC) c CALL histwrite_phy(nid_hfnmc,lNMC,"uv",itau_w, $uvsumSTD(:,:,3)) c CALL histwrite_phy(nid_hfnmc,lNMC,"vq",itau_w, $vqsumSTD(:,:,3)) c CALL histwrite_phy(nid_hfnmc,lNMC,"vT",itau_w, $vTsumSTD(:,:,3)) c CALL histwrite_phy(nid_hfnmc,lNMC,"wq",itau_w, $wqsumSTD(:,:,3)) c CALL histwrite_phy(nid_hfnmc,lNMC,"vphi",itau_w, $vphisumSTD(:,:,3)) c CALL histwrite_phy(nid_hfnmc,lNMC,"wT",itau_w, $wTsumSTD(:,:,3)) c CALL histwrite_phy(nid_hfnmc,lNMC,"uxu",itau_w, $u2sumSTD(:,:,3)) c CALL histwrite_phy(nid_hfnmc,lNMC,"vxv",itau_w, $v2sumSTD(:,:,3)) c CALL histwrite_phy(nid_hfnmc,lNMC,"TxT",itau_w, $T2sumSTD(:,:,3)) c ENDIF !(1.EQ.0) THEN c if (ok_sync) then c$OMP MASTER call histsync(nid_hfnmc) c$OMP END MASTER endif c ENDIF ! (ok_histNMC(3)) THEN