Changeset 320 for LMDZ.3.3/branches/rel-LF/libf/dyn3d/limit_netcdf.F
- Timestamp:
- Jan 25, 2002, 5:20:10 PM (22 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/dyn3d/limit_netcdf.F
r316 r320 2 2 C $Header$ 3 3 C 4 SUBROUTINE limit_netcdf ( interbar, extrap, oldice )4 SUBROUTINE limit_netcdf ( interbar, extrap, oldice, masque ) 5 5 c 6 6 IMPLICIT none … … 28 28 #include "comgeom2.h" 29 29 #include "comconst.h" 30 #include "dimphy.h" 30 31 c 31 32 c----------------------------------------------------------------------- … … 33 34 34 35 INTEGER KIDIA, KFDIA, KLON, KLEV 35 PARAMETER (KIDIA=1,KFDIA=iim*(jjm-1)+2,36 . KLON=KFDIA-KIDIA+1,KLEV=llm)37 36 c----------------------------------------------------------------------- 38 37 REAL phy_nat(klon,360), phy_nat0(klon) … … 195 194 CALL mask_c_o(imdep, jmdep, dlon_msk, dlat_msk,champ_msk, 196 195 . iim, jjp1, rlonv, rlatu, champint) 197 CALL gr_int_dyn(champint, masque, iim, jjp1)198 DO i = 1, iim199 masque(i,1) = FLOAT(NINT(masque(i,1)))200 masque(i,jjp1) = FLOAT(NINT(masque(i,jjp1)))201 ENDDO196 c CALL gr_int_dyn(champint, masque, iim, jjp1) 197 c DO i = 1, iim 198 c masque(i,1) = FLOAT(NINT(masque(i,1))) 199 c masque(i,jjp1) = FLOAT(NINT(masque(i,jjp1))) 200 c ENDDO 202 201 DO i = 1, iim 203 202 DO j = 1, jjp1 … … 205 204 ENDDO 206 205 ENDDO 207 CALL gr_dyn_fi(1, iip1, jjp1, klon, masque, phy_nat0)206 c CALL gr_dyn_fi(1, iip1, jjp1, klon, masque, phy_nat0) 208 207 ierr = NF_CLOSE(ncid) 209 208 c … … 591 590 . champan(1,1,k), phy_ice(1,k)) 592 591 DO i = 1, klon 593 phy_nat(i,k) = phy_nat0(i)592 phy_nat(i,k) = zmasq(i) 594 593 IF ( (phy_ice(i,k) - 0.5).GE.1.e-5 ) THEN 595 IF (NINT( phy_nat0(i)).EQ.0) THEN594 IF (NINT(zmasq(i)).EQ.0) THEN 596 595 phy_nat(i,k) = 3.0 597 596 ELSE
Note: See TracChangeset
for help on using the changeset viewer.