/[lmdze]/trunk/libf/dyn3d/limit.f90
ViewVC logotype

Diff of /trunk/libf/dyn3d/limit.f90

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

revision 15 by guez, Fri Aug 1 15:24:12 2008 UTC revision 22 by guez, Fri Jul 31 15:18:47 2009 UTC
# Line 27  contains Line 27  contains
27      use numer_rec, only: spline, splint      use numer_rec, only: spline, splint
28      use grid_change, only: dyn_phy      use grid_change, only: dyn_phy
29    
30      use netcdf95, only: handle_err, nf95_get_coord, NF95_CLOSE, NF95_DEF_DIM, &      use netcdf95, only: handle_err, nf95_gw_var, NF95_CLOSE, NF95_DEF_DIM, &
31           nf95_enddef, NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, &           nf95_enddef, NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, &
32           nf95_inq_varid, NF95_OPEN           nf95_inq_varid, NF95_OPEN
33      use netcdf, only: NF90_CLOBBER, nf90_def_var, NF90_FLOAT, NF90_GET_VAR, &      use netcdf, only: NF90_CLOBBER, nf90_def_var, NF90_FLOAT, NF90_GET_VAR, &
# Line 97  contains Line 97  contains
97    
98      ! Read coordinate variables:      ! Read coordinate variables:
99    
100      call nf95_get_coord(ncid, "longitude", dlon_ini)      call nf95_inq_varid(ncid, "longitude", varid)
101        call nf95_gw_var(ncid, varid, dlon_ini)
102      imdep = size(dlon_ini)      imdep = size(dlon_ini)
103    
104      call nf95_get_coord(ncid, "latitude", dlat_ini)      call nf95_inq_varid(ncid, "latitude", varid)
105        call nf95_gw_var(ncid, varid, dlat_ini)
106      jmdep = size(dlat_ini)      jmdep = size(dlat_ini)
107    
108      call nf95_get_coord(ncid, "temps", timeyear)      call nf95_inq_varid(ncid, "temps", varid)
109        call nf95_gw_var(ncid, varid, timeyear)
110      lmdep = size(timeyear)      lmdep = size(timeyear)
111    
112      ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep))      ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep))
# Line 148  contains Line 151  contains
151      PRINT *, 'Processing sea ice...'      PRINT *, 'Processing sea ice...'
152      call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid)      call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid)
153    
154      call nf95_get_coord(ncid, "longitude", dlon_ini)      call nf95_inq_varid(ncid, "longitude", varid)
155        call nf95_gw_var(ncid, varid, dlon_ini)
156      imdep = size(dlon_ini)      imdep = size(dlon_ini)
157    
158      call nf95_get_coord(ncid, "latitude", dlat_ini)      call nf95_inq_varid(ncid, "latitude", varid)
159        call nf95_gw_var(ncid, varid, dlat_ini)
160      jmdep = size(dlat_ini)      jmdep = size(dlat_ini)
161    
162      call nf95_inq_dimid(ncid, "time", dimid)      call nf95_inq_dimid(ncid, "time", dimid)
# Line 245  contains Line 250  contains
250      PRINT *, 'Traitement de la sst'      PRINT *, 'Traitement de la sst'
251      call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid)      call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid)
252    
253      call nf95_get_coord(ncid, "longitude", dlon_ini)      call nf95_inq_varid(ncid, "longitude", varid)
254        call nf95_gw_var(ncid, varid, dlon_ini)
255      imdep = size(dlon_ini)      imdep = size(dlon_ini)
256    
257      call nf95_get_coord(ncid, "latitude", dlat_ini)      call nf95_inq_varid(ncid, "latitude", varid)
258        call nf95_gw_var(ncid, varid, dlat_ini)
259      jmdep = size(dlat_ini)      jmdep = size(dlat_ini)
260    
261      call nf95_inq_dimid(ncid, "time", dimid)      call nf95_inq_dimid(ncid, "time", dimid)
# Line 314  contains Line 321  contains
321      PRINT *, 'Traitement de l albedo'      PRINT *, 'Traitement de l albedo'
322      call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid)      call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid)
323    
324      call nf95_get_coord(ncid, "longitude", dlon_ini)      call nf95_inq_varid(ncid, "longitude", varid)
325        call nf95_gw_var(ncid, varid, dlon_ini)
326      imdep = size(dlon_ini)      imdep = size(dlon_ini)
327    
328      call nf95_get_coord(ncid, "latitude", dlat_ini)      call nf95_inq_varid(ncid, "latitude", varid)
329        call nf95_gw_var(ncid, varid, dlat_ini)
330      jmdep = size(dlat_ini)      jmdep = size(dlat_ini)
331    
332      call nf95_get_coord(ncid, "temps", timeyear)      call nf95_inq_varid(ncid, "temps", varid)
333        call nf95_gw_var(ncid, varid, timeyear)
334      lmdep = size(timeyear)      lmdep = size(timeyear)
335    
336      ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep))      ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep))
Legend:
Removed from v.15  
changed lines
  Added in v.22
  ViewVC Help
Powered by ViewVC 1.1.21