Changeset 1638 for trunk/LMDZ.COMMON/libf/dynphy_lonlat

Timestamp:

Dec 6, 2016, 11:45:04 AM (8 years ago)

Author:

slebonnois

Message:

SL: corrections for Venus newstart

Location:

trunk/LMDZ.COMMON/libf/dynphy_lonlat/phyvenus

Files:

• grid_noro.F (moved) (moved from trunk/LMDZ.VENUS/libf/phyvenus/grid_noro.F) (4 diffs)
• startvar.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/dynphy_lonlat/phyvenus/grid_noro.F

@@ -52 +52 @@
 C=======================================================================
 
-      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
       IMPLICIT none 
       
@@ -79 +78 @@
 C INTERMEDIATE FIELDS  (CORRELATIONS OF OROGRAPHY GRADIENT)
 
-      REAL ztz(nbp_lon+1,nbp_lat),zxtzx(nbp_lon+1,nbp_lat)
-      REAL zytzy(nbp_lon+1,nbp_lat),zxtzy(nbp_lon+1,nbp_lat)
-      REAL weight(nbp_lon+1,nbp_lat)
+      REAL ztz(imar+1,jmar),zxtzx(imar+1,jmar)
+      REAL zytzy(imar+1,jmar),zxtzy(imar+1,jmar)
+      REAL weight(imar+1,jmar)
 
 C CORRELATIONS OF USN OROGRAPHY GRADIENTS
@@ -107 +106 @@
 c    
 c
-      if(nbp_lon.ne.imar) STOP 'Problem dim. x'
-      if(nbp_lat-1.ne.jmar-1) STOP 'Problem dim. y'
       IF (imar.GT.2200 .OR. jmar.GT.1100) THEN
          PRINT*, 'imar or jmar too big', imar, jmar
          CALL ABORT
       ENDIF
-
-      IF(imar+1.ne.nbp_lon+1.or.jmar.ne.nbp_lat)THEN
-        print *,' imar or jmar bad dimensions:',imar,jmar
-        call abort
-      ENDIF
-
 
 c      print *,'xdata:',xdata
@@ -297 +288 @@
 C  FIRST FILTER, MOVING AVERAGE OVER 9 POINTS.
 
-       CALL MVA9(zmea,nbp_lon+1,nbp_lat)
-       CALL MVA9(zstd,nbp_lon+1,nbp_lat)
-       CALL MVA9(zpic,nbp_lon+1,nbp_lat)
-       CALL MVA9(zval,nbp_lon+1,nbp_lat)
-       CALL MVA9(zxtzx,nbp_lon+1,nbp_lat)
-       CALL MVA9(zxtzy,nbp_lon+1,nbp_lat) 
-       CALL MVA9(zytzy,nbp_lon+1,nbp_lat)
+       CALL MVA9(zmea,imar+1,jmar)
+       CALL MVA9(zstd,imar+1,jmar)
+       CALL MVA9(zpic,imar+1,jmar)
+       CALL MVA9(zval,imar+1,jmar)
+       CALL MVA9(zxtzx,imar+1,jmar)
+       CALL MVA9(zxtzy,imar+1,jmar) 
+       CALL MVA9(zytzy,imar+1,jmar)
 
        DO ii = 1, imar

trunk/LMDZ.COMMON/libf/dynphy_lonlat/phyvenus/startvar.F90

@@ -428 +428 @@
   IF(check) WRITE(lunout,*)'Compute surface roughness induced by the orography'
   ALLOCATE(rugo   (iml  ,jml))
-  ALLOCATE(tmp_var(iml-1,jml))
   CALL rugsoro(lon_rad, lat_rad, relief_hi,      &
-       lon_in, lat_in, tmp_var)
-  rugo(1:iml-1,:)=tmp_var; rugo(iml,:)=tmp_var(1,:)
-  DEALLOCATE(relief_hi,tmp_var,lon_rad,lat_rad)
+       lon_in, lat_in, rugo)
+  DEALLOCATE(relief_hi,lon_rad,lat_rad)
   RETURN