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

Diff of /trunk/dyn3d/etat0.f

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

trunk/Sources/dyn3d/etat0.f revision 212 by guez, Thu Jan 12 12:31:31 2017 UTC trunk/dyn3d/etat0.f revision 265 by guez, Tue Mar 20 09:35:59 2018 UTC
# Line 21  contains Line 21  contains
21      use comconst, only: cpp, kappa, iniconst      use comconst, only: cpp, kappa, iniconst
22      use comgeom, only: aire_2d, apoln, apols, cu_2d, cv_2d, inigeom      use comgeom, only: aire_2d, apoln, apols, cu_2d, cv_2d, inigeom
23      use conf_gcm_m, only: nday      use conf_gcm_m, only: nday
24      use dimens_m, only: iim, jjm, llm, nqmx      use dimensions, only: iim, jjm, llm, nqmx
25      use dimphy, only: zmasq      use dimphy, only: zmasq
26      use dimsoil, only: nsoilmx      use dimsoil, only: nsoilmx
27      use disvert_m, only: ap, bp, preff, pa, disvert      use disvert_m, only: ap, bp, preff, pa, disvert
# Line 91  contains Line 91  contains
91      ! D\'eclarations pour lecture glace de mer :      ! D\'eclarations pour lecture glace de mer :
92      INTEGER iml_lic, jml_lic      INTEGER iml_lic, jml_lic
93      INTEGER ncid, varid      INTEGER ncid, varid
94      REAL, pointer:: dlon_lic(:), dlat_lic(:)      REAL, ALLOCATABLE:: dlon_lic(:), dlat_lic(:)
95      REAL, ALLOCATABLE:: fraclic(:, :) ! fraction land ice      REAL, ALLOCATABLE:: fraclic(:, :) ! fraction land ice
96      REAL flic_tmp(iim + 1, jjm + 1) ! fraction land ice temporary      REAL flic_tmp(iim + 1, jjm + 1) ! fraction land ice temporary
97    
# Line 261  contains Line 261  contains
261           rlatu)           rlatu)
262      flic_tmp(iim + 1, :) = flic_tmp(1, :)      flic_tmp(iim + 1, :) = flic_tmp(1, :)
263    
     deallocate(dlon_lic, dlat_lic) ! pointers  
   
264      ! Passage sur la grille physique      ! Passage sur la grille physique
265      pctsrf = 0.      pctsrf = 0.
266      pctsrf(:, is_lic) = pack(flic_tmp, dyn_phy)      pctsrf(:, is_lic) = pack(flic_tmp, dyn_phy)
267      ! Ad\'equation avec le maque terre/mer      ! Ad\'equation avec le maque terre/mer
268      WHERE (pctsrf(:, is_lic) < EPSFRA) pctsrf(:, is_lic) = 0.      WHERE (pctsrf(:, is_lic) < EPSFRA) pctsrf(:, is_lic) = 0.
269      WHERE (zmasq < EPSFRA) pctsrf(:, is_lic) = 0.      WHERE (zmasq < EPSFRA) pctsrf(:, is_lic) = 0.
270      pctsrf(:, is_ter) = zmasq      where (zmasq <= EPSFRA) pctsrf(:, is_ter) = zmasq
271      where (zmasq > EPSFRA)      where (zmasq > EPSFRA)
272         where (pctsrf(:, is_lic) >= zmasq)         where (pctsrf(:, is_lic) >= zmasq)
273            pctsrf(:, is_lic) = zmasq            pctsrf(:, is_lic) = zmasq

Legend:
Removed from v.212  
changed lines
  Added in v.265

  ViewVC Help
Powered by ViewVC 1.1.21