/[lmdze]/trunk/libf/test_disvert.f90
ViewVC logotype

Diff of /trunk/libf/test_disvert.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 37 by guez, Tue Dec 21 15:45:48 2010 UTC revision 38 by guez, Thu Jan 6 17:52:19 2011 UTC
# Line 29  program test_disvert Line 29  program test_disvert
29    call new_unit(unit)    call new_unit(unit)
30    open(unit, file="test_disvert.csv", status="replace", action="write")    open(unit, file="test_disvert.csv", status="replace", action="write")
31    ! Title line:    ! Title line:
32    write(unit, fmt=*) '"ap (hPa)" "bp" "s" "p (hPa)", "z (km)", "delta z (km)"'    write(unit, fmt=*) '"ap (hPa)" "bp" "s" "p (hPa)" "z (km)"'
33    do l = 1, llm - 1    do l = 1, llm
34       write(unit, fmt=*) ap(l) / 100., bp(l), s(l), p(l) / 100., z(l), &       write(unit, fmt=*) ap(l) / 100., bp(l), s(l), p(l) / 100., z(l)
           z(l+1) - z(l)  
35    end do    end do
   write(unit, fmt=*) ap(llm) / 100., bp(llm), s(llm), p(llm) / 100., z(llm)  
36    close(unit)    close(unit)
37    print *, 'The file "test_disvert.csv" has been created.'    print *, 'The file "test_disvert.csv" has been created.'
38    
39      open(unit, file="delta_z.csv", status="replace", action="write")
40      ! Title line:
41      write(unit, fmt=*) '"z middle (km)" "delta z (km)"'
42      do l = 1, llm - 1
43         write(unit, fmt=*) (z(l) + z(l+1)) / 2, z(l+1) - z(l)
44      end do
45      close(unit)
46      print *, 'The file "delta_z.csv" has been created.'
47    
48  end program test_disvert  end program test_disvert
Legend:
Removed from v.37  
changed lines
  Added in v.38
  ViewVC Help
Powered by ViewVC 1.1.21