/[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 44 by guez, Wed Apr 13 12:29:18 2011 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 50  C Line 51  C
51    
52        INTEGER histid, nq        INTEGER histid, nq
53        REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm)        REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm)
54        REAL teta(ip1jmp1*llm),phi(ip1jmp1,llm),ppk(ip1jmp1*llm)                          REAL, intent(in):: teta(ip1jmp1*llm)
55          real phi(ip1jmp1,llm),ppk(ip1jmp1*llm)                  
56        REAL ps(ip1jmp1),masse(ip1jmp1,llm)                          REAL ps(ip1jmp1),masse(ip1jmp1,llm)                  
57        REAL phis(ip1jmp1)                          REAL phis(ip1jmp1)                  
58        REAL q(ip1jmp1,llm,nq)        REAL q(ip1jmp1,llm,nq)
# Line 87  C Line 89  C
89  C  Vents U scalaire  C  Vents U scalaire
90  C  C
91        call gr_u_scal(llm, unat, us)        call gr_u_scal(llm, unat, us)
92        call histwrite(histid, 'u', itau_w, us,        call histwrite(histid, 'u', itau_w, us)
      .               iip1*jjp1*llm, ndex3d)  
93  C  C
94  C  Vents V scalaire  C  Vents V scalaire
95  C  C
96        call gr_v_scal(llm, vnat, vs)        call gr_v_scal(llm, vnat, vs)
97        call histwrite(histid, 'v', itau_w, vs,        call histwrite(histid, 'v', itau_w, vs)
      .               iip1*jjp1*llm, ndex3d)  
98  C  C
99  C  Temperature potentielle moyennee  C  Temperature potentielle moyennee
100  C  C
101        call histwrite(histid, 'theta', itau_w, teta,        call histwrite(histid, 'theta', itau_w, teta)
      .                iip1*jjp1*llm, ndex3d)  
102  C  C
103  C  Temperature moyennee  C  Temperature moyennee
104  C  C
105        do ii = 1, ijp1llm        do ii = 1, ijp1llm
106          tm(ii) = teta(ii) * ppk(ii)/cpp          tm(ii) = teta(ii) * ppk(ii)/cpp
107        enddo        enddo
108        call histwrite(histid, 'temp', itau_w, tm,        call histwrite(histid, 'temp', itau_w, tm)
      .                iip1*jjp1*llm, ndex3d)  
109  C  C
110  C  Geopotentiel  C  Geopotentiel
111  C  C
112        call histwrite(histid, 'phi', itau_w, phi,        call histwrite(histid, 'phi', itau_w, phi)
      .                iip1*jjp1*llm, ndex3d)  
113  C  C
114  C  Traceurs  C  Traceurs
115  C  C
116          DO iq=1,nq          DO iq=1,nq
117            call histwrite(histid, ttext(iq), itau_w, q(:,:,iq),            call histwrite(histid, ttext(iq), itau_w, q(:,:,iq))
      .                   iip1*jjp1*llm, ndex3d)  
118          enddo          enddo
119  C  C
120  C  Masse  C  Masse
121  C  C
122         call histwrite(histid, 'masse', itau_w, masse, iip1*jjp1, ndex2d)         call histwrite(histid, 'masse', itau_w, masse)
123  C  C
124  C  Pression au sol  C  Pression au sol
125  C  C
126         call histwrite(histid, 'ps', itau_w, ps, iip1*jjp1, ndex2d)         call histwrite(histid, 'ps', itau_w, ps)
127  C  C
128  C  Geopotentiel au sol  C  Geopotentiel au sol
129  C  C
130         call histwrite(histid, 'phis', itau_w, phis, iip1*jjp1, ndex2d)         call histwrite(histid, 'phis', itau_w, phis)
131  C  C
132  C  Fin  C  Fin
133  C  C
Legend:
Removed from v.3  
changed lines
  Added in v.44
  ViewVC Help
Powered by ViewVC 1.1.21