/[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 7 by guez, Mon Mar 31 12:24:17 2008 UTC revision 22 by guez, Fri Jul 31 15:18:47 2009 UTC
# Line 2  program test_disvert Line 2  program test_disvert
2    
3    use dimens_m, only: llm    use dimens_m, only: llm
4    use comvert, only: pa, disvert, ap, bp, preff    use comvert, only: pa, disvert, ap, bp, preff
   use comconst, only: initialize  
5    use new_unit_m, only: new_unit    use new_unit_m, only: new_unit
6    
7    implicit none    implicit none
# Line 14  program test_disvert Line 13  program test_disvert
13    !---------------------    !---------------------
14    
15    print *, "llm = ", llm    print *, "llm = ", llm
   call initialize  
16    pa = 5e4    pa = 5e4
17    call disvert    call disvert
18    sigma = ap / pa + bp    sigma = ap / pa + bp
19    p = ap + bp * preff    p = ap + bp * preff
20    
21    unit = new_unit()    call new_unit(unit)
22    open(unit, file="test_disvert.csv", status="replace", action="write")    open(unit, file="test_disvert.csv", status="replace", action="write")
23    write(unit, fmt=*) '"ap (Pa)" "bp" "sigma" "p (Pa)"' ! title line    write(unit, fmt=*) '"ap (Pa)" "bp" "sigma" "p (Pa)"' ! title line
24    do l = 1, llm    do l = 1, llm
Legend:
Removed from v.7  
changed lines
  Added in v.22
  ViewVC Help
Powered by ViewVC 1.1.21