/[lmdze]/trunk/dyn3d/dynredem0.f
ViewVC logotype

Diff of /trunk/dyn3d/dynredem0.f

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 5 by guez, Mon Mar 3 16:32:04 2008 UTC
# Line 1  Line 1 
1  SUBROUTINE dynredem0(fichnom,iday_end,phis,nq)  SUBROUTINE dynredem0(fichnom, iday_end, phis)
2    
3    ! From dyn3d/dynredem.F, v 1.2 2004/06/22 11:45:30    ! From dyn3d/dynredem.F, v 1.2 2004/06/22 11:45:30
4    
5    ! Ecriture du fichier de redémarrage au format NetCDF (initialisation)    ! Ecriture du fichier de redémarrage au format NetCDF (initialisation)
6    
7    USE IOIPSL, only: ymds2ju, ju2ymds    USE IOIPSL, only: ymds2ju, ju2ymds
8    use dimens_m, only: iim, jjm, llm    use dimens_m, only: iim, jjm, llm, nqmx
9    use paramet_m, only: ip1jmp1, iip1, jjp1, llmp1    use paramet_m, only: ip1jmp1, iip1, jjp1, llmp1
10    use comconst, only: rad, cpp, daysec, dtvr, kappa, g, omeg    use comconst, only: rad, cpp, daysec, dtvr, kappa, g, omeg
11    use comvert, only: pa, bp, ap, nivsigs, preff, presnivs, nivsig    use comvert, only: pa, bp, ap, nivsigs, preff, presnivs, nivsig
# Line 15  SUBROUTINE dynredem0(fichnom,iday_end,ph Line 15  SUBROUTINE dynredem0(fichnom,iday_end,ph
15    use temps, only: annee_ref, day_ref, itaufin, itau_dyn    use temps, only: annee_ref, day_ref, itaufin, itau_dyn
16    use ener    use ener
17    use advtrac_m, only: tname, ttext    use advtrac_m, only: tname, ttext
18      use netcdf95, only: nf95_create, NF90_CLOBBER, nf95_put_att, NF90_GLOBAL, &
19           nf95_def_dim, NF90_UNLIMITED
20    
21    IMPLICIT NONE    IMPLICIT NONE
22    
23    include "netcdf.inc"    CHARACTER(len=*), intent(in):: fichnom
   
   !   Arguments:  
24    INTEGER, intent(in):: iday_end    INTEGER, intent(in):: iday_end
25    REAL, intent(in):: phis(ip1jmp1)    REAL, intent(in):: phis(ip1jmp1)
   CHARACTER(len=*) fichnom  
   INTEGER nq  
26    
27    !   Local:    !   Local:
28    INTEGER iq,l  
29      include "netcdf.inc"
30    
31      INTEGER iq, l
32    INTEGER length    INTEGER length
33    PARAMETER (length = 100)    PARAMETER (length = 100)
34    REAL tab_cntrl(length) ! tableau des parametres du run    REAL tab_cntrl(length) ! tableau des parametres du run
35    INTEGER ierr    INTEGER ierr
36    
37    !   Variables locales pour NetCDF:    !   Variables locales pour NetCDF:
38    !  
39    INTEGER dims2(2), dims3(3), dims4(4)    INTEGER dims2(2), dims3(3), dims4(4)
40    INTEGER idim_index    INTEGER idim_index
41    INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv    INTEGER idim_rlonu, idim_rlonv, idim_rlatu, idim_rlatv
42    INTEGER idim_s, idim_sig    INTEGER idim_s, idim_sig
43    INTEGER idim_tim    INTEGER idim_tim
44    INTEGER nid,nvarid    INTEGER nid, nvarid
45    
46    REAL zjulian,hours    REAL zjulian, hours
47    INTEGER yyears0,jjour0, mmois0    INTEGER yyears0, jjour0, mmois0
48    character(len=30) unites    character(len=30) unites
49    
50    !-----------------------------------------------------------------------    !-----------------------------------------------------------------------
# Line 53  SUBROUTINE dynredem0(fichnom,iday_end,ph Line 54  SUBROUTINE dynredem0(fichnom,iday_end,ph
54    call ymds2ju(annee_ref, 1, iday_end, 0.0, zjulian)    call ymds2ju(annee_ref, 1, iday_end, 0.0, zjulian)
55    call ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)    call ju2ymds(zjulian, yyears0, mmois0, jjour0, hours)
56    
57    DO l=1,length    DO l=1, length
58       tab_cntrl(l) = 0.       tab_cntrl(l) = 0.
59    ENDDO    ENDDO
60    tab_cntrl(1)  = REAL(iim)    tab_cntrl(1)  = REAL(iim)
# Line 75  SUBROUTINE dynredem0(fichnom,iday_end,ph Line 76  SUBROUTINE dynredem0(fichnom,iday_end,ph
76    tab_cntrl(17) = ang0    tab_cntrl(17) = ang0
77    tab_cntrl(18) = pa    tab_cntrl(18) = pa
78    tab_cntrl(19) = preff    tab_cntrl(19) = preff
79    !  
80    !    .....    parametres  pour le zoom      ......      ! Paramètres  pour le zoom :
81    
82    tab_cntrl(20)  = clon    tab_cntrl(20)  = clon
83    tab_cntrl(21)  = clat    tab_cntrl(21)  = clat
84    tab_cntrl(22)  = grossismx    tab_cntrl(22)  = grossismx
85    tab_cntrl(23)  = grossismy    tab_cntrl(23)  = grossismy
86    !  
87    IF ( fxyhypb )   THEN    IF ( fxyhypb )   THEN
88       tab_cntrl(24) = 1.       tab_cntrl(24) = 1.
89       tab_cntrl(25) = dzoomx       tab_cntrl(25) = dzoomx
# Line 102  SUBROUTINE dynredem0(fichnom,iday_end,ph Line 103  SUBROUTINE dynredem0(fichnom,iday_end,ph
103    
104    tab_cntrl(30) = REAL(iday_end)    tab_cntrl(30) = REAL(iday_end)
105    tab_cntrl(31) = REAL(itau_dyn + itaufin)    tab_cntrl(31) = REAL(itau_dyn + itaufin)
106    !  
107    !    .........................................................    call nf95_create(fichnom, NF90_CLOBBER, nid)
108    !    call nf95_put_att(nid, NF90_GLOBAL, "title", &
109    ! Creation du fichier:         "Fichier de démarrage dynamique")
110    !  
   ierr = NF_CREATE(fichnom, NF_CLOBBER, nid)  
   IF (ierr.NE.NF_NOERR) THEN  
      WRITE(6,*)" Pb d ouverture du fichier "//fichnom  
      WRITE(6,*)' ierr = ', ierr  
      stop 1  
   ENDIF  
   !  
   ! Preciser quelques attributs globaux:  
   !  
   ierr = NF_PUT_ATT_TEXT (nid, NF_GLOBAL, "title", 27, &  
        "Fichier demmarage dynamique")  
   !  
111    ! Definir les dimensions du fichiers:    ! Definir les dimensions du fichiers:
112    !  
113    ierr = NF_DEF_DIM (nid, "index", length, idim_index)    call nf95_def_dim(nid, "index", length, idim_index)
114    ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)    call NF95_DEF_DIM(nid, "rlonu", iip1, idim_rlonu)
115    ierr = NF_DEF_DIM (nid, "rlatu", jjp1, idim_rlatu)    call NF95_DEF_DIM(nid, "rlatu", jjp1, idim_rlatu)
116    ierr = NF_DEF_DIM (nid, "rlonv", iip1, idim_rlonv)    call NF95_DEF_DIM(nid, "rlonv", iip1, idim_rlonv)
117    ierr = NF_DEF_DIM (nid, "rlatv", jjm, idim_rlatv)    call NF95_DEF_DIM(nid, "rlatv", jjm, idim_rlatv)
118    ierr = NF_DEF_DIM (nid, "sigs", llm, idim_s)    call NF95_DEF_DIM(nid, "sigs", llm, idim_s)
119    ierr = NF_DEF_DIM (nid, "sig", llmp1, idim_sig)    call NF95_DEF_DIM(nid, "sig", llmp1, idim_sig)
120    ierr = NF_DEF_DIM (nid, "temps", NF_UNLIMITED, idim_tim)    call NF95_DEF_DIM(nid, "temps", NF90_UNLIMITED, idim_tim)
121    !  
122    ierr = NF_ENDDEF(nid) ! sortir du mode de definition    ierr = NF_ENDDEF(nid) ! sortir du mode de definition
123    !  
124    ! Definir et enregistrer certains champs invariants:    ! Definir et enregistrer certains champs invariants:
125    !  
126    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
127    ierr = NF_DEF_VAR (nid,"controle",NF_FLOAT,1,idim_index,nvarid)    ierr = NF_DEF_VAR (nid, "controle", NF_FLOAT, 1, idim_index, nvarid)
128    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22, &    call nf95_put_att(nid, nvarid, "title", "Parametres de controle")
        "Parametres de controle")  
129    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
130    ierr = NF_PUT_VAR_REAL (nid,nvarid,tab_cntrl)    ierr = NF_PUT_VAR_REAL (nid, nvarid, tab_cntrl)
131    !  
132    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
133    ierr = NF_DEF_VAR (nid,"rlonu",NF_FLOAT,1,idim_rlonu,nvarid)    ierr = NF_DEF_VAR (nid, "rlonu", NF_FLOAT, 1, idim_rlonu, nvarid)
134    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23, &    call nf95_put_att(nid, nvarid, "title", "Longitudes des points U")
        "Longitudes des points U")  
135    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
136    ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonu)    ierr = NF_PUT_VAR_REAL (nid, nvarid, rlonu)
137    !  
138    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
139    ierr = NF_DEF_VAR (nid,"rlatu",NF_FLOAT,1,idim_rlatu,nvarid)    ierr = NF_DEF_VAR (nid, "rlatu", NF_FLOAT, 1, idim_rlatu, nvarid)
140    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22, &    call nf95_put_att(nid, nvarid, "title", "Latitudes des points U")
        "Latitudes des points U")  
141    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
142    ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatu)    ierr = NF_PUT_VAR_REAL (nid, nvarid, rlatu)
143    !  
144    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
145    ierr = NF_DEF_VAR (nid,"rlonv",NF_FLOAT,1,idim_rlonv,nvarid)    ierr = NF_DEF_VAR (nid, "rlonv", NF_FLOAT, 1, idim_rlonv, nvarid)
146    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 23, &    call nf95_put_att(nid, nvarid, "title", "Longitudes des points V")
        "Longitudes des points V")  
147    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
148    ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonv)    ierr = NF_PUT_VAR_REAL (nid, nvarid, rlonv)
149    !  
150    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
151    ierr = NF_DEF_VAR (nid,"rlatv",NF_FLOAT,1,idim_rlatv,nvarid)    ierr = NF_DEF_VAR (nid, "rlatv", NF_FLOAT, 1, idim_rlatv, nvarid)
152    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22, &    call nf95_put_att(nid, nvarid, "title", "Latitudes des points V")
        "Latitudes des points V")  
153    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
154    ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatv)    ierr = NF_PUT_VAR_REAL (nid, nvarid, rlatv)
155    !  
156    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
157    ierr = NF_DEF_VAR (nid,"nivsigs",NF_FLOAT,1,idim_s,nvarid)    ierr = NF_DEF_VAR (nid, "nivsigs", NF_FLOAT, 1, idim_s, nvarid)
158    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 28, &    call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches s")
        "Numero naturel des couches s")  
159    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
160    ierr = NF_PUT_VAR_REAL (nid,nvarid,nivsigs)    ierr = NF_PUT_VAR_REAL (nid, nvarid, nivsigs)
161    !  
162    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
163    ierr = NF_DEF_VAR (nid,"nivsig",NF_FLOAT,1,idim_sig,nvarid)    ierr = NF_DEF_VAR (nid, "nivsig", NF_FLOAT, 1, idim_sig, nvarid)
164    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 32, &    call nf95_put_att(nid, nvarid, "title", "Numero naturel des couches sigma")
        "Numero naturel des couches sigma")  
165    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
166    ierr = NF_PUT_VAR_REAL (nid,nvarid,nivsig)    ierr = NF_PUT_VAR_REAL (nid, nvarid, nivsig)
167    !  
168    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
169    ierr = NF_DEF_VAR (nid,"ap",NF_FLOAT,1,idim_sig,nvarid)    ierr = NF_DEF_VAR (nid, "ap", NF_FLOAT, 1, idim_sig, nvarid)
170    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 26, &    call nf95_put_att(nid, nvarid, "title", "Coefficient A pour hybride")
        "Coefficient A pour hybride")  
171    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
172    ierr = NF_PUT_VAR_REAL (nid,nvarid,ap)    ierr = NF_PUT_VAR_REAL (nid, nvarid, ap)
173    !  
174    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
175    ierr = NF_DEF_VAR (nid,"bp",NF_FLOAT,1,idim_sig,nvarid)    ierr = NF_DEF_VAR (nid, "bp", NF_FLOAT, 1, idim_sig, nvarid)
176    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 26, &    call nf95_put_att(nid, nvarid, "title", "Coefficient B pour hybride")
        "Coefficient B pour hybride")  
177    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
178    ierr = NF_PUT_VAR_REAL (nid,nvarid,bp)    ierr = NF_PUT_VAR_REAL (nid, nvarid, bp)
179    !  
180    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
181    ierr = NF_DEF_VAR (nid,"presnivs",NF_FLOAT,1,idim_s,nvarid)    ierr = NF_DEF_VAR (nid, "presnivs", NF_FLOAT, 1, idim_s, nvarid)
182    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
183    ierr = NF_PUT_VAR_REAL (nid,nvarid,presnivs)    ierr = NF_PUT_VAR_REAL (nid, nvarid, presnivs)
184    !  
185    ! Coefficients de passage cov. <-> contra. <--> naturel    ! Coefficients de passage cov. <-> contra. <--> naturel
186    !  
187    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
188    dims2(1) = idim_rlonu    dims2(1) = idim_rlonu
189    dims2(2) = idim_rlatu    dims2(2) = idim_rlatu
190    ierr = NF_DEF_VAR (nid,"cu",NF_FLOAT,2,dims2,nvarid)    ierr = NF_DEF_VAR (nid, "cu", NF_FLOAT, 2, dims2, nvarid)
191    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29, &    call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour U")
        "Coefficient de passage pour U")  
192    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
193    ierr = NF_PUT_VAR_REAL (nid,nvarid,cu)    ierr = NF_PUT_VAR_REAL (nid, nvarid, cu)
194    !  
195    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
196    dims2(1) = idim_rlonv    dims2(1) = idim_rlonv
197    dims2(2) = idim_rlatv    dims2(2) = idim_rlatv
198    ierr = NF_DEF_VAR (nid,"cv",NF_FLOAT,2,dims2,nvarid)    ierr = NF_DEF_VAR (nid, "cv", NF_FLOAT, 2, dims2, nvarid)
199    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 29, &    call nf95_put_att(nid, nvarid, "title", "Coefficient de passage pour V")
        "Coefficient de passage pour V")  
200    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
201    ierr = NF_PUT_VAR_REAL (nid,nvarid,cv)    ierr = NF_PUT_VAR_REAL (nid, nvarid, cv)
202    !  
203    ! Aire de chaque maille:    ! Aire de chaque maille:
204    !  
205    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
206    dims2(1) = idim_rlonv    dims2(1) = idim_rlonv
207    dims2(2) = idim_rlatu    dims2(2) = idim_rlatu
208    ierr = NF_DEF_VAR (nid,"aire",NF_FLOAT,2,dims2,nvarid)    ierr = NF_DEF_VAR (nid, "aire", NF_FLOAT, 2, dims2, nvarid)
209    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 22, &    call nf95_put_att(nid, nvarid, "title", "Aires de chaque maille")
        "Aires de chaque maille")  
210    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
211    ierr = NF_PUT_VAR_REAL (nid,nvarid,aire)    ierr = NF_PUT_VAR_REAL (nid, nvarid, aire)
212    !  
213    ! Geopentiel au sol:    ! Geopentiel au sol:
214    !  
215    ierr = NF_REDEF (nid)    ierr = NF_REDEF (nid)
216    dims2(1) = idim_rlonv    dims2(1) = idim_rlonv
217    dims2(2) = idim_rlatu    dims2(2) = idim_rlatu
218    ierr = NF_DEF_VAR (nid,"phisinit",NF_FLOAT,2,dims2,nvarid)    ierr = NF_DEF_VAR (nid, "phisinit", NF_FLOAT, 2, dims2, nvarid)
219    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19, &    call nf95_put_att(nid, nvarid, "title", "Geopotentiel au sol")
        "Geopotentiel au sol")  
220    ierr = NF_ENDDEF(nid)    ierr = NF_ENDDEF(nid)
221    ierr = NF_PUT_VAR_REAL (nid,nvarid,phis)    ierr = NF_PUT_VAR_REAL (nid, nvarid, phis)
222    !  
223    ! Definir les variables pour pouvoir les enregistrer plus tard:    ! Definir les variables pour pouvoir les enregistrer plus tard:
224    !  
225    ierr = NF_REDEF (nid) ! entrer dans le mode de definition    ierr = NF_REDEF (nid) ! entrer dans le mode de definition
   !  
   ierr = NF_DEF_VAR (nid,"temps",NF_FLOAT,1,idim_tim,nvarid)  
   ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 19, &  
        "Temps de simulation")  
   write(unites,200)yyears0,mmois0,jjour0  
 200 format('days since ',i4,'-',i2.2,'-',i2.2,' 00:00:00')  
   ierr = NF_PUT_ATT_TEXT (nid, nvarid, "units", 30, &  
        unites)  
226    
227    !    ierr = NF_DEF_VAR (nid, "temps", NF_FLOAT, 1, idim_tim, nvarid)
228      call nf95_put_att(nid, nvarid, "title", "Temps de simulation")
229      write(unites, 200)yyears0, mmois0, jjour0
230    200 format('days since ', i4, '-', i2.2, '-', i2.2, ' 00:00:00')
231      call nf95_put_att(nid, nvarid, "units", unites)
232    
233    
234    dims4(1) = idim_rlonu    dims4(1) = idim_rlonu
235    dims4(2) = idim_rlatu    dims4(2) = idim_rlatu
236    dims4(3) = idim_s    dims4(3) = idim_s
237    dims4(4) = idim_tim    dims4(4) = idim_tim
238    ierr = NF_DEF_VAR (nid,"ucov",NF_FLOAT,4,dims4,nvarid)    ierr = NF_DEF_VAR (nid, "ucov", NF_FLOAT, 4, dims4, nvarid)
239    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9, &    call nf95_put_att(nid, nvarid, "title", "Vitesse U")
240         "Vitesse U")  
   !  
241    dims4(1) = idim_rlonv    dims4(1) = idim_rlonv
242    dims4(2) = idim_rlatv    dims4(2) = idim_rlatv
243    dims4(3) = idim_s    dims4(3) = idim_s
244    dims4(4) = idim_tim    dims4(4) = idim_tim
245    ierr = NF_DEF_VAR (nid,"vcov",NF_FLOAT,4,dims4,nvarid)    ierr = NF_DEF_VAR (nid, "vcov", NF_FLOAT, 4, dims4, nvarid)
246    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 9, &    call nf95_put_att(nid, nvarid, "title", "Vitesse V")
247         "Vitesse V")  
   !  
248    dims4(1) = idim_rlonv    dims4(1) = idim_rlonv
249    dims4(2) = idim_rlatu    dims4(2) = idim_rlatu
250    dims4(3) = idim_s    dims4(3) = idim_s
251    dims4(4) = idim_tim    dims4(4) = idim_tim
252    ierr = NF_DEF_VAR (nid,"teta",NF_FLOAT,4,dims4,nvarid)    ierr = NF_DEF_VAR (nid, "teta", NF_FLOAT, 4, dims4, nvarid)
253    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 11, &    call nf95_put_att(nid, nvarid, "title", "Temperature")
254         "Temperature")  
   !  
255    dims4(1) = idim_rlonv    dims4(1) = idim_rlonv
256    dims4(2) = idim_rlatu    dims4(2) = idim_rlatu
257    dims4(3) = idim_s    dims4(3) = idim_s
258    dims4(4) = idim_tim    dims4(4) = idim_tim
259    IF(nq.GE.1) THEN    DO iq=1, nqmx
260       DO iq=1,nq       ierr = NF_DEF_VAR (nid, tname(iq), NF_FLOAT, 4, dims4, nvarid)
261          ierr = NF_DEF_VAR (nid,tname(iq),NF_FLOAT,4,dims4,nvarid)       call nf95_put_att(nid, nvarid, "title", ttext(iq))
262          ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,ttext(iq))    ENDDO
263       ENDDO  
   ENDIF  
   !  
264    dims4(1) = idim_rlonv    dims4(1) = idim_rlonv
265    dims4(2) = idim_rlatu    dims4(2) = idim_rlatu
266    dims4(3) = idim_s    dims4(3) = idim_s
267    dims4(4) = idim_tim    dims4(4) = idim_tim
268    ierr = NF_DEF_VAR (nid,"masse",NF_FLOAT,4,dims4,nvarid)    ierr = NF_DEF_VAR (nid, "masse", NF_FLOAT, 4, dims4, nvarid)
269    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12, &    call nf95_put_att(nid, nvarid, "title", "C est quoi ?")
270         "C est quoi ?")  
   !  
271    dims3(1) = idim_rlonv    dims3(1) = idim_rlonv
272    dims3(2) = idim_rlatu    dims3(2) = idim_rlatu
273    dims3(3) = idim_tim    dims3(3) = idim_tim
274    ierr = NF_DEF_VAR (nid,"ps",NF_FLOAT,3,dims3,nvarid)    ierr = NF_DEF_VAR (nid, "ps", NF_FLOAT, 3, dims3, nvarid)
275    ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 15, &    call nf95_put_att(nid, nvarid, "title", "Pression au sol")
276         "Pression au sol")  
   !  
277    ierr = NF_ENDDEF(nid) ! sortir du mode de definition    ierr = NF_ENDDEF(nid) ! sortir du mode de definition
278    ierr = NF_CLOSE(nid) ! fermer le fichier    ierr = NF_CLOSE(nid) ! fermer le fichier
279    
280    PRINT*,'iim,jjm,llm,iday_end',iim,jjm,llm,iday_end    PRINT*, 'iim, jjm, llm, iday_end', iim, jjm, llm, iday_end
281    PRINT*,'rad,omeg,g,cpp,kappa', &    PRINT*, 'rad, omeg, g, cpp, kappa', rad, omeg, g, cpp, kappa
        rad,omeg,g,cpp,kappa  
282    
283  END SUBROUTINE dynredem0  END SUBROUTINE dynredem0
Legend:
Removed from v.3  
changed lines
  Added in v.5
  ViewVC Help
Powered by ViewVC 1.1.21