Changeset 2361 for LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd

Timestamp:

Sep 7, 2015, 5:50:29 PM (10 years ago)

Author:

dcugnet

Message:

In etat0dyn: removed few useless lines: "masque" is always known because etat0dyn is called after etat0phys.
In grid_atob: shape error in routine fine2coarse now fixed: "msk" argument and local variable mask must have the dimensions of the output grid. Working unit of dist_sphe is no longer degree, but radian.

File:

• LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd/etat0dyn_netcdf.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/dynlonlat_phylonlat/phylmd/etat0dyn_netcdf.F90

@@ -36 +36 @@
   USE ioipsl,         ONLY: flininfo, flinopen, flinget, flinclo, histclo
   USE assert_eq_m,    ONLY: assert_eq
-  USE indice_sol_mod, ONLY: epsfra
   IMPLICIT NONE
 
@@ -87 +86 @@
   INTEGER            :: i, j, l, ji, itau, iday
   REAL               :: xpn, xps, time, phystep
-  REAL, DIMENSION(iip1,jjp1)       :: psol, masque_tmp
+  REAL, DIMENSION(iip1,jjp1)       :: psol
   REAL, DIMENSION(iip1,jjp1,llm+1) :: p3d
   REAL, DIMENSION(iip1,jjp1,llm)   :: uvent, t3d, tpot, qsat, qd
@@ -101 +100 @@
 
   deg2rad = pi/180.0
-
-! Compute ground geopotential and possibly the mask.
-!*******************************************************************************
-  masque_tmp(:,:)=masque(:,:)
-  WRITE(fmt,"(i4,'i1)')")iip1 ; fmt='('//ADJUSTL(fmt)
-  IF(ALL(masque==-99999.)) THEN                         !--- KEEP NEW MASK
-    masque=masque_tmp
-    IF(prt_level>=1) THEN
-      WRITE(lunout,*)'BUILT MASK :'
-      WRITE(lunout,fmt) NINT(masque)
-    END IF
-    WHERE(   masque(:,:)<EPSFRA) masque(:,:)=0.
-    WHERE(1.-masque(:,:)<EPSFRA) masque(:,:)=1.
-  END IF
 
 ! Compute psol AND tsol, knowing phis.