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

Diff of /trunk/dyn3d/limit.f90

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

revision 323 by guez, Thu Jul 12 14:49:20 2018 UTC revision 324 by guez, Wed Feb 6 15:58:03 2019 UTC
# Line 27  contains Line 27  contains
27           NF90_UNLIMITED           NF90_UNLIMITED
28      use nr_util, only: assert      use nr_util, only: assert
29      use numer_rec_95, only: spline, splint      use numer_rec_95, only: spline, splint
30      use phyetat0_m, only: zmasq      use phyetat0_m, only: masque
31      use start_init_orog_m, only: mask      use start_init_orog_m, only: mask
32      use unit_nml_m, only: unit_nml      use unit_nml_m, only: unit_nml
33    
# Line 249  contains Line 249  contains
249         pctsrf(:, is_sic) = max(phy_ice - pctsrf(:, is_lic), 0.)         pctsrf(:, is_sic) = max(phy_ice - pctsrf(:, is_lic), 0.)
250         ! Il y a des cas o\`u il y a de la glace dans landiceref et         ! Il y a des cas o\`u il y a de la glace dans landiceref et
251         ! pas dans AMIP         ! pas dans AMIP
252         WHERE (1. - zmasq < EPSFRA)         WHERE (1. - masque < EPSFRA)
253            pctsrf(:, is_sic) = 0.            pctsrf(:, is_sic) = 0.
254            pctsrf(:, is_oce) = 0.            pctsrf(:, is_oce) = 0.
255         elsewhere         elsewhere
256            where (pctsrf(:, is_sic) >= 1 - zmasq)            where (pctsrf(:, is_sic) >= 1 - masque)
257               pctsrf(:, is_sic) = 1. - zmasq               pctsrf(:, is_sic) = 1. - masque
258               pctsrf(:, is_oce) = 0.               pctsrf(:, is_oce) = 0.
259            ELSEwhere            ELSEwhere
260               pctsrf(:, is_oce) = 1. - zmasq - pctsrf(:, is_sic)               pctsrf(:, is_oce) = 1. - masque - pctsrf(:, is_sic)
261               where (pctsrf(:, is_oce) < EPSFRA)               where (pctsrf(:, is_oce) < EPSFRA)
262                  pctsrf(:, is_oce) = 0.                  pctsrf(:, is_oce) = 0.
263                  pctsrf(:, is_sic) = 1 - zmasq                  pctsrf(:, is_sic) = 1 - masque
264               end where               end where
265            end where            end where
266         end where         end where
Legend:
Removed from v.323  
changed lines
  Added in v.324
  ViewVC Help
Powered by ViewVC 1.1.21