Changeset 312 for trunk/LMDZ.MARS/util/extract.F90

Timestamp:

Oct 2, 2011, 5:10:54 PM (13 years ago)

Author:

emillour

Message:

Mars GCM utilities: Minor fix on extract (entering a blank line in interactive mode to signal end of list now works as expected).
EM

File:

• trunk/LMDZ.MARS/util/extract.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/util/extract.F90

@@ -37 +37 @@
 character(len=256) :: outfile
 
-character (len=64) :: text ! to store some text
+character (len=256) :: text ! to store some text
 character (len=64) :: varname ! to store the name of the variable to retreive
 
@@ -339 +339 @@
 do 
   ! 2.1 read coordinates and do some sanity checks
+  text="" !initialize text to empty string
+  read(*,'(a)') text ! store input as text (to spot empty input line)
+  if (len_trim(adjustl(text)).eq.0) exit
+  
   if (extract_mode==1) then
-    read(*,*,iostat=status) lon,lat,alt,Ls,LT
+    read(text,*,iostat=status) lon,lat,alt,Ls,LT
     ! Ls in degrees, LT (local true solar time) in hours, i.e.  in [0:24]
   else ! extract_mode==2 , read alt only
-    read(*,*,iostat=status) alt
+    read(text,*,iostat=status) alt
   endif
   if (status.ne.0) exit
@@ -612 +616 @@
 
 end subroutine extraction
+
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 subroutine ls2sol(ls,sol)
@@ -653 +658 @@
 
 end subroutine ls2sol
+