/[lmdze]/trunk/libf/bibio/writedynav.f90
ViewVC logotype

Diff of /trunk/libf/bibio/writedynav.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 30 by guez, Thu Apr 1 09:07:28 2010 UTC
# Line 4  Line 4 
4        subroutine writedynav( histid, nq, time, vcov,        subroutine writedynav( histid, nq, time, vcov,
5       ,                          ucov,teta,ppk,phi,q,masse,ps,phis)       ,                          ucov,teta,ppk,phi,q,masse,ps,phis)
6    
       USE ioipsl  
7  C  C
8  C   Ecriture du fichier histoire au format IOIPSL  C   Ecriture du fichier histoire au format IOIPSL
9  C  C
# Line 31  C Line 30  C
30  C =====================================================================  C =====================================================================
31  C  C
32  C   Declarations  C   Declarations
33          USE histwrite_m
34          use histcom
35        use dimens_m        use dimens_m
36        use paramet_m        use paramet_m
37        use comconst        use comconst
# Line 40  C   Declarations Line 41  C   Declarations
41        use serre        use serre
42        use temps        use temps
43        use ener        use ener
44        use advtrac_m        use iniadvtrac_m
45        implicit none        implicit none
46    
47    
# Line 87  C Line 88  C
88  C  Vents U scalaire  C  Vents U scalaire
89  C  C
90        call gr_u_scal(llm, unat, us)        call gr_u_scal(llm, unat, us)
91        call histwrite(histid, 'u', itau_w, us,        call histwrite(histid, 'u', itau_w, us)
      .               iip1*jjp1*llm, ndex3d)  
92  C  C
93  C  Vents V scalaire  C  Vents V scalaire
94  C  C
95        call gr_v_scal(llm, vnat, vs)        call gr_v_scal(llm, vnat, vs)
96        call histwrite(histid, 'v', itau_w, vs,        call histwrite(histid, 'v', itau_w, vs)
      .               iip1*jjp1*llm, ndex3d)  
97  C  C
98  C  Temperature potentielle moyennee  C  Temperature potentielle moyennee
99  C  C
100        call histwrite(histid, 'theta', itau_w, teta,        call histwrite(histid, 'theta', itau_w, teta)
      .                iip1*jjp1*llm, ndex3d)  
101  C  C
102  C  Temperature moyennee  C  Temperature moyennee
103  C  C
104        do ii = 1, ijp1llm        do ii = 1, ijp1llm
105          tm(ii) = teta(ii) * ppk(ii)/cpp          tm(ii) = teta(ii) * ppk(ii)/cpp
106        enddo        enddo
107        call histwrite(histid, 'temp', itau_w, tm,        call histwrite(histid, 'temp', itau_w, tm)
      .                iip1*jjp1*llm, ndex3d)  
108  C  C
109  C  Geopotentiel  C  Geopotentiel
110  C  C
111        call histwrite(histid, 'phi', itau_w, phi,        call histwrite(histid, 'phi', itau_w, phi)
      .                iip1*jjp1*llm, ndex3d)  
112  C  C
113  C  Traceurs  C  Traceurs
114  C  C
115          DO iq=1,nq          DO iq=1,nq
116            call histwrite(histid, ttext(iq), itau_w, q(:,:,iq),            call histwrite(histid, ttext(iq), itau_w, q(:,:,iq))
      .                   iip1*jjp1*llm, ndex3d)  
117          enddo          enddo
118  C  C
119  C  Masse  C  Masse
120  C  C
121         call histwrite(histid, 'masse', itau_w, masse, iip1*jjp1, ndex2d)         call histwrite(histid, 'masse', itau_w, masse)
122  C  C
123  C  Pression au sol  C  Pression au sol
124  C  C
125         call histwrite(histid, 'ps', itau_w, ps, iip1*jjp1, ndex2d)         call histwrite(histid, 'ps', itau_w, ps)
126  C  C
127  C  Geopotentiel au sol  C  Geopotentiel au sol
128  C  C
129         call histwrite(histid, 'phis', itau_w, phis, iip1*jjp1, ndex2d)         call histwrite(histid, 'phis', itau_w, phis)
130  C  C
131  C  Fin  C  Fin
132  C  C
Legend:
Removed from v.3  
changed lines
  Added in v.30
  ViewVC Help
Powered by ViewVC 1.1.21