Changeset 777 for trunk/LMDZ.MARS/util/zrecast.F90

Timestamp:

Sep 11, 2012, 3:50:51 PM (12 years ago)

Author:

emillour

Message:

Mars GCM utilities:
Bug correction (for the "above local surface", ie ztype==3, case) in zrecast utility.
EM

File:

• trunk/LMDZ.MARS/util/zrecast.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/util/zrecast.F90

@@ -565 +565 @@
   ! ask for # of points and end values for pressure or above areoid cases
   write(*,*) ""
-  if ((ztype.le.2).or.(ztype.le.4)) then
+  if ((ztype.le.2).or.(ztype.ge.4)) then
     write(*,*) "Enter min and max of vertical coordinate (Pa or m)"
     write(*,*) " (in that order and on the same line)"
@@ -573 +573 @@
       read(*,*) zamin,zamax
     endif
-  endif ! of if ((ztype.le.2).or.(ztype.le.4))
+  endif ! of if ((ztype.le.2).or.(ztype.ge.4))
 
   ! Build corresponding vertical coordinates
@@ -668 +668 @@
       read(*,*) zareoid(i)
     enddo
-  else if (ztype.eq.2) then ! above local surface
+  else if (ztype.eq.3) then ! above local surface
     allocate(zsurface(nblev),stat=ierr)
     if (ierr.ne.0) then