/[lmdze]/trunk/Sources/dyn3d/limit.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/limit.f

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

revision 6 by guez, Wed Feb 27 13:16:39 2008 UTC revision 7 by guez, Mon Mar 31 12:24:17 2008 UTC
# Line 27  contains Line 27  contains
27      use interpolation, only: spline, splint      use interpolation, only: spline, splint
28      use grid_change, only: dyn_phy      use grid_change, only: dyn_phy
29    
30      use netcdf95, only: handle_err, coordin, &      use netcdf95, only: handle_err, nf95_get_coord, NF95_CLOSE, NF95_DEF_DIM, &
31           NF90_CLOBBER, NF95_CLOSE, NF95_DEF_DIM, nf90_def_var, nf95_enddef, &           nf95_enddef, NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, &
32           NF90_FLOAT, NF90_GET_VAR, NF90_GLOBAL, NF90_NOWRITE, NF90_PUT_ATT, &           nf95_inq_varid, NF95_OPEN
33           NF90_PUT_VAR, NF90_UNLIMITED, &      use netcdf, only: NF90_CLOBBER, nf90_def_var, NF90_FLOAT, NF90_GET_VAR, &
34           NF95_CREATE, nf95_inq_dimid, nf95_inquire_dimension, nf95_inq_varid, &           NF90_GLOBAL, NF90_NOWRITE, NF90_PUT_ATT, NF90_PUT_VAR, &
35           nf95_open           NF90_UNLIMITED
36    
37      ! Variables local to the procedure:      ! Variables local to the procedure:
38    
# Line 95  contains Line 95  contains
95      PRINT *, 'Processing rugosity...'      PRINT *, 'Processing rugosity...'
96      call NF95_OPEN('Rugos.nc', NF90_NOWRITE, ncid)      call NF95_OPEN('Rugos.nc', NF90_NOWRITE, ncid)
97    
98      dlon_ini => coordin(ncid, "longitude")      call nf95_get_coord(ncid, "longitude", dlon_ini)
99      imdep = size(dlon_ini)      imdep = size(dlon_ini)
100    
101      dlat_ini => coordin(ncid, "latitude")      call nf95_get_coord(ncid, "latitude", dlat_ini)
102      jmdep = size(dlat_ini)      jmdep = size(dlat_ini)
103    
104      timeyear => coordin(ncid, "temps")      call nf95_get_coord(ncid, "temps", timeyear)
105      lmdep = size(timeyear)      lmdep = size(timeyear)
106    
107      ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep))      ALLOCATE(champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep))
# Line 144  contains Line 144  contains
144      PRINT *, 'Processing sea ice...'      PRINT *, 'Processing sea ice...'
145      call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid)      call NF95_OPEN('amipbc_sic_1x1.nc', NF90_NOWRITE, ncid)
146    
147      dlon_ini => coordin(ncid, "longitude")      call nf95_get_coord(ncid, "longitude", dlon_ini)
148      imdep = size(dlon_ini)      imdep = size(dlon_ini)
149    
150      dlat_ini => coordin(ncid, "latitude")      call nf95_get_coord(ncid, "latitude", dlat_ini)
151      jmdep = size(dlat_ini)      jmdep = size(dlat_ini)
152    
153      call nf95_inq_dimid(ncid, "time", dimid)      call nf95_inq_dimid(ncid, "time", dimid)
# Line 241  contains Line 241  contains
241      PRINT *, 'Traitement de la sst'      PRINT *, 'Traitement de la sst'
242      call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid)      call NF95_OPEN('amipbc_sst_1x1.nc', NF90_NOWRITE, ncid)
243    
244      dlon_ini => coordin(ncid, "longitude")      call nf95_get_coord(ncid, "longitude", dlon_ini)
245      imdep = size(dlon_ini)      imdep = size(dlon_ini)
246    
247      dlat_ini => coordin(ncid, "latitude")      call nf95_get_coord(ncid, "latitude", dlat_ini)
248      jmdep = size(dlat_ini)      jmdep = size(dlat_ini)
249    
250      call nf95_inq_dimid(ncid, "time", dimid)      call nf95_inq_dimid(ncid, "time", dimid)
# Line 310  contains Line 310  contains
310      PRINT *, 'Traitement de l albedo'      PRINT *, 'Traitement de l albedo'
311      call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid)      call NF95_OPEN('Albedo.nc', NF90_NOWRITE, ncid)
312    
313      dlon_ini => coordin(ncid, "longitude")      call nf95_get_coord(ncid, "longitude", dlon_ini)
314      imdep = size(dlon_ini)      imdep = size(dlon_ini)
315    
316      dlat_ini => coordin(ncid, "latitude")      call nf95_get_coord(ncid, "latitude", dlat_ini)
317      jmdep = size(dlat_ini)      jmdep = size(dlat_ini)
318    
319      timeyear => coordin(ncid, "temps")      call nf95_get_coord(ncid, "temps", timeyear)
320      lmdep = size(timeyear)      lmdep = size(timeyear)
321    
322      ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep))      ALLOCATE ( champ(imdep, jmdep), champtime(iim, jjm + 1, lmdep))
Legend:
Removed from v.6  
changed lines
  Added in v.7
  ViewVC Help
Powered by ViewVC 1.1.21