Changeset 370 for LMDZ.3.3/branches/rel-LF/libf/dyn3d/grid_noro.F

Timestamp:

Jul 10, 2002, 5:07:19 PM (23 years ago)

Author:

lmdzadmin

Message:

Prise en compte du fichier de trait de cote (si il existe) et du masque
correspondant dans les parametrisation de FLott
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/dyn3d/grid_noro.F (modified) (5 diffs)

LMDZ.3.3/branches/rel-LF/libf/dyn3d/grid_noro.F

@@ -1 @@
-
+c
+c $Header$
+c
       SUBROUTINE grid_noro(imdep, jmdep, xdata, ydata, zdata,
      .             imar, jmar, x, y,
@@ -82 +84 @@
 c
       REAL a(2200),b(2200),c(1100),d(1100)
+      logical masque_lu
 c
       print *,' parametres de l orographie a l echelle sous maille' 
@@ -87 +90 @@
       rad    = 6 371 229.
       zdeltay=2.*xpi/float(jusn)*rad
+c
+c utilise-t'on un masque lu?
+c
+      masque_lu = .true.
+      if (maxval(mask) == -99999 .and. minval(mask) == -99999) then
+        masque_lu= .false.
+        masque = 0.0
+      endif
+      write(*,*)'Masque lu', masque_lu
 c
 c  quelques tests de dimensions:
@@ -263 +275 @@
 c$$$             mask(ii,jj)=0
 c$$$           ENDIF
-             mask(ii,jj) = num_lan(ii,jj)/num_tot(ii,jj)
+             if (.not. masque_lu) then
+               mask(ii,jj) = num_lan(ii,jj)/num_tot(ii,jj)
+             endif
 c  Mean Orography:
            zmea (ii,jj)=zmea (ii,jj)/weight(ii,jj)
@@ -357 +371 @@
        ENDDO
        ENDDO
-
       print *,'  MEAN ORO:',zllmmea
       print *,'  ST. DEV.:',zllmstd