Ignore:
Timestamp:
Jan 25, 2002, 5:20:10 PM (22 years ago)
Author:
lmdzadmin
Message:

Adaptation à la version couplée
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/dyn3d/limit_netcdf.F

    r316 r320  
    22C $Header$
    33C
    4       SUBROUTINE limit_netcdf ( interbar, extrap, oldice )
     4      SUBROUTINE limit_netcdf ( interbar, extrap, oldice, masque )
    55c
    66      IMPLICIT none
     
    2828#include "comgeom2.h"
    2929#include "comconst.h"
     30#include "dimphy.h"
    3031c
    3132c-----------------------------------------------------------------------
     
    3334
    3435      INTEGER KIDIA, KFDIA, KLON, KLEV
    35       PARAMETER (KIDIA=1,KFDIA=iim*(jjm-1)+2,
    36      .           KLON=KFDIA-KIDIA+1,KLEV=llm)
    3736c-----------------------------------------------------------------------
    3837      REAL phy_nat(klon,360), phy_nat0(klon)
     
    195194      CALL mask_c_o(imdep, jmdep, dlon_msk, dlat_msk,champ_msk,
    196195     .             iim, jjp1, rlonv, rlatu, champint)
    197       CALL gr_int_dyn(champint, masque, iim, jjp1)
    198       DO i = 1, iim
    199          masque(i,1) = FLOAT(NINT(masque(i,1)))
    200          masque(i,jjp1) = FLOAT(NINT(masque(i,jjp1)))
    201       ENDDO
     196c      CALL gr_int_dyn(champint, masque, iim, jjp1)
     197c      DO i = 1, iim
     198c         masque(i,1) = FLOAT(NINT(masque(i,1)))
     199c         masque(i,jjp1) = FLOAT(NINT(masque(i,jjp1)))
     200c      ENDDO
    202201      DO i = 1, iim
    203202      DO j = 1, jjp1
     
    205204      ENDDO
    206205      ENDDO
    207       CALL gr_dyn_fi(1, iip1, jjp1, klon, masque, phy_nat0)
     206c      CALL gr_dyn_fi(1, iip1, jjp1, klon, masque, phy_nat0)
    208207      ierr = NF_CLOSE(ncid)
    209208c
     
    591590     .                  champan(1,1,k), phy_ice(1,k))
    592591         DO i = 1, klon
    593             phy_nat(i,k) = phy_nat0(i)
     592            phy_nat(i,k) = zmasq(i)
    594593            IF ( (phy_ice(i,k) - 0.5).GE.1.e-5 ) THEN
    595                IF (NINT(phy_nat0(i)).EQ.0) THEN
     594               IF (NINT(zmasq(i)).EQ.0) THEN
    596595                  phy_nat(i,k) = 3.0
    597596               ELSE
Note: See TracChangeset for help on using the changeset viewer.