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

Timestamp:

Sep 13, 2011, 9:19:35 AM (13 years ago)

Author:

emillour

Message:

Update of Mars GCM tool 'extract': extracted values are now sent to an ASCII output file ('infile_var.dat', where 'inifil' is the input file name and 'var' is the extracted variable).
EM

File:

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

trunk/LMDZ.MARS/util/extract.F90

@@ -18 +18 @@
 !         coordinate in the input file), one per line, at which values are
 !         sought
-! - output values are given on standard output; with extraction mode 1 as
-!   lines of "lon lat alt Ls LT value" an for a profile (extraction mode 2)
+! - output values are sent to (ASCII) output file 'infile_var_.dat', where
+!   'infile' is the input file name (without trailing '.nc') and
+!   'var' is the sought variable, for extraction mode 1 as
+!   lines of "lon lat alt Ls LT value" and for a profile (extraction mode 2)
 !   as lines of "alt value"
 !
@@ -33 +35 @@
 ! Input file:
 character(len=256) :: infile
+character(len=256) :: outfile
 
 character (len=64) :: text ! to store some text
@@ -310 +313 @@
 ! 2. Process and interpolate
 !===============================================================================
+
+!===============================================================================
+! 2.1 Create output file
+!===============================================================================
+
+outfile=trim(infile(1:index(infile,".nc",back=.true.)-1))//"_"//&
+        trim(varname)//".dat"
+open(42,file=outfile,form="formatted")
+write(*,*) "Output file is: ",trim(outfile)
+
+!===============================================================================
+! 2.2 Extract values
+!===============================================================================
+
 if (extract_mode==1) then ! pointwise extraction
  write(*,*) " Enter values of: lon lat alt Ls LT (all on the same line,"
@@ -370 +387 @@
   ! 2.4 Write value to output
   if (extract_mode==1) then ! pointwise extraction
-    write(*,'(6(1x,1pe12.5))')lon,lat,alt,Ls,LT,value
+    write(42,'(6(1x,1pe12.5))')lon,lat,alt,Ls,LT,value
   else ! profile 
-    write(*,'(6(1x,1pe12.5))')alt,value
+    write(42,'(6(1x,1pe12.5))')alt,value
   endif
 enddo ! of do while
 
+! close output file
+close(42)
 
 end program extract