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

Timestamp:

May 14, 2014, 11:00:12 AM (11 years ago)

Author:

emillour

Message:

zrecast utility: minor correction for the very specific "mcd" case,
no need to read in zamin and zamax.
EM

File:

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

trunk/LMDZ.MARS/util/zrecast.F90

@@ -141 +141 @@
                             ! has to be computed; false if these are given
                             ! by the user
-logical :: auto_mcd_levels ! Flag: specific case for MCD automatic above
-                           ! local surface levels
+logical :: auto_mcd_levels=.false. ! Flag: specific case for MCD automatic 
+                           ! above local surface levels
 integer,dimension(4) :: edges,corner ! needed to write variables for big files
 
@@ -610 +610 @@
   elseif(trim(text).eq."no") then
     auto_vert_levels=.false.
+    auto_mcd_levels=.false.
   elseif(trim(text).eq."mcd") then
     auto_vert_levels=.true.
@@ -618 +619 @@
 if (auto_vert_levels) then
   ! ask for # of points and end values for pressure or above areoid cases
-  write(*,*) ""
-  write(*,*) "Enter min and max of vertical coordinate (Pa or m)"
-  write(*,*) " (in that order and on the same line)"
-  if (ztype.eq.1) then ! pressure coordinate
-    read(*,*) pmin,pmax
-  else ! altitude coordinate
-     read(*,*) zamin,zamax
+  ! except in MCD case
+  if (.not.auto_mcd_levels) then
+    write(*,*) ""
+    write(*,*) "Enter min and max of vertical coordinate (Pa or m)"
+    write(*,*) " (in that order and on the same line)"
+    if (ztype.eq.1) then ! pressure coordinate
+      read(*,*) pmin,pmax
+    else ! altitude coordinate, except in MCD case
+      read(*,*) zamin,zamax
+    endif
   endif