--- trunk/libf/bibio/writedynav.f 2008/02/27 13:16:39 3 +++ trunk/libf/bibio/writedynav.f 2011/04/13 12:29:18 44 @@ -4,7 +4,6 @@ subroutine writedynav( histid, nq, time, vcov, , ucov,teta,ppk,phi,q,masse,ps,phis) - USE ioipsl C C Ecriture du fichier histoire au format IOIPSL C @@ -31,6 +30,8 @@ C ===================================================================== C C Declarations + USE histwrite_m + use histcom use dimens_m use paramet_m use comconst @@ -40,7 +41,7 @@ use serre use temps use ener - use advtrac_m + use iniadvtrac_m implicit none @@ -50,7 +51,8 @@ INTEGER histid, nq REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm) - REAL teta(ip1jmp1*llm),phi(ip1jmp1,llm),ppk(ip1jmp1*llm) + REAL, intent(in):: teta(ip1jmp1*llm) + real phi(ip1jmp1,llm),ppk(ip1jmp1*llm) REAL ps(ip1jmp1),masse(ip1jmp1,llm) REAL phis(ip1jmp1) REAL q(ip1jmp1,llm,nq) @@ -87,51 +89,45 @@ C Vents U scalaire C call gr_u_scal(llm, unat, us) - call histwrite(histid, 'u', itau_w, us, - . iip1*jjp1*llm, ndex3d) + call histwrite(histid, 'u', itau_w, us) C C Vents V scalaire C call gr_v_scal(llm, vnat, vs) - call histwrite(histid, 'v', itau_w, vs, - . iip1*jjp1*llm, ndex3d) + call histwrite(histid, 'v', itau_w, vs) C C Temperature potentielle moyennee C - call histwrite(histid, 'theta', itau_w, teta, - . iip1*jjp1*llm, ndex3d) + call histwrite(histid, 'theta', itau_w, teta) C C Temperature moyennee C do ii = 1, ijp1llm tm(ii) = teta(ii) * ppk(ii)/cpp enddo - call histwrite(histid, 'temp', itau_w, tm, - . iip1*jjp1*llm, ndex3d) + call histwrite(histid, 'temp', itau_w, tm) C C Geopotentiel C - call histwrite(histid, 'phi', itau_w, phi, - . iip1*jjp1*llm, ndex3d) + call histwrite(histid, 'phi', itau_w, phi) C C Traceurs C DO iq=1,nq - call histwrite(histid, ttext(iq), itau_w, q(:,:,iq), - . iip1*jjp1*llm, ndex3d) + call histwrite(histid, ttext(iq), itau_w, q(:,:,iq)) enddo C C Masse C - call histwrite(histid, 'masse', itau_w, masse, iip1*jjp1, ndex2d) + call histwrite(histid, 'masse', itau_w, masse) C C Pression au sol C - call histwrite(histid, 'ps', itau_w, ps, iip1*jjp1, ndex2d) + call histwrite(histid, 'ps', itau_w, ps) C C Geopotentiel au sol C - call histwrite(histid, 'phis', itau_w, phis, iip1*jjp1, ndex2d) + call histwrite(histid, 'phis', itau_w, phis) C C Fin C