tstqsat.F90 Source File


This file depends on

sourcefile~~tstqsat.f90~~EfferentGraph sourcefile~tstqsat.f90 tstqsat.F90 sourcefile~base.f90 Base.F90 sourcefile~tstqsat.f90->sourcefile~base.f90

Source Code

program tstqsat

use MAPLBase_Mod


integer i,l
integer i1, i2, rate

real*4, dimension(140:340) :: q1, q2, q3, q4, q5, q6, t, p
real*8, dimension(140:340) :: qd1, qd2, qd3, qd4, qd5, qd6, td, pd


open(20,file='eqsat_verification.dat',form='formatted',status='unknown')

call MAPL_EQsatSet(Formulation=MAPL_UseMurphyKoopQsat)

do i=140,340
   p(i)=100000.
   t(i)=i
   pd(i)=100000.
   td(i)=i
   q3(i)=MAPL_Eqsat(t(i))
   q4(i)=MAPL_Eqsat(t(i),OverIce=.true.)
   q5(i)=MAPL_Eqsat(t(i),P(i))
   q6(i)=MAPL_Eqsat(t(i),P(i),OverIce=.true.)
enddo

call mapl_eqsatset(usetable=.false.)

do i=140,340
   q1(i)=mapl_EQSAT(t(i))
   q2(i)=mapl_EQSAT(t(i),OverIce=.true.)
   write(20,*)  t(i),q1(i),q3(i),q2(i),q4(i),q5(i),q6(i)
enddo


do i=140,340
t(i)=i
enddo

call MAPL_EqsatSet(Formulation=MAPL_UseStarrQsat)
q1=MAPL_Eqsat(T)
q4=MAPL_Eqsat(T,OverIce=.true.)
call MAPL_EqsatSet(Formulation=MAPL_UseCAMQsat)
q2=MAPL_Eqsat(T)
q5=MAPL_Eqsat(T,OverIce=.true.)
call MAPL_EqsatSet(Formulation=MAPL_UseMurphyKoopQsat)
q3=MAPL_Eqsat(T)
q6=MAPL_Eqsat(T,OverIce=.true.)


do i=140,340
   write(20,*)  t(i),q1(i),q2(i),q3(i),q4(i),q5(i),q6(i)
enddo


call MAPL_EqsatSet(Formulation=1)
q1=MAPL_Eqsat(T,p)
q4=MAPL_Eqsat(T,p,OverIce=.true.)
call MAPL_EqsatSet(Formulation=2)
q2=MAPL_Eqsat(T,p)
q5=MAPL_Eqsat(T,p,OverIce=.true.)
call MAPL_EqsatSet(Formulation=3)
q3=MAPL_Eqsat(T,p)
q6=MAPL_Eqsat(T,p,OverIce=.true.)


do i=140,340
   write(20,*)  t(i),q1(i),q2(i),q3(i),q4(i),q5(i),q6(i)
enddo


call mapl_eqsatset(usetable=.true.)


call MAPL_EqsatSet(Formulation=MAPL_UseStarrQsat)
q1=MAPL_Eqsat(T)
q4=MAPL_Eqsat(T,OverIce=.true.)
call MAPL_EqsatSet(Formulation=MAPL_UseCAMQsat)
q2=MAPL_Eqsat(T)
q5=MAPL_Eqsat(T,OverIce=.true.)
call MAPL_EqsatSet(Formulation=MAPL_UseMurphyKoopQsat)
q3=MAPL_Eqsat(T)
q6=MAPL_Eqsat(T,OverIce=.true.)


do i=140,340
   write(20,*)  t(i),q1(i),q2(i),q3(i),q4(i),q5(i),q6(i)
enddo


call MAPL_EqsatSet(Formulation=MAPL_UseStarrQsat)
q1=MAPL_Eqsat(T)
q4=MAPL_Eqsat(T,OverIce=.true.)
qd1=MAPL_Eqsat(Td)
qd4=MAPL_Eqsat(Td,OverIce=.true.)

do i=140,340
   write(20,*)  t(i),q1(i),qd1(i),q4(i),qd4(i)
enddo

call MAPL_EqsatSet(Formulation=MAPL_UseCAMQsat)
q2=MAPL_Eqsat(T)
q5=MAPL_Eqsat(T,OverIce=.true.)
qd2=MAPL_Eqsat(Td)
qd5=MAPL_Eqsat(Td,OverIce=.true.)

do i=140,340
   write(20,*)  t(i),q2(i),qd2(i),q5(i),qd5(i)
enddo

call MAPL_EqsatSet(Formulation=MAPL_UseMurphyKoopQsat)
q3=MAPL_Eqsat(T)
q6=MAPL_Eqsat(T,OverIce=.true.)
qd3=MAPL_Eqsat(Td)
qd6=MAPL_Eqsat(Td,OverIce=.true.)

do i=140,340
   write(20,*)  t(i),q3(i),qd3(i),q6(i),qd6(i)
enddo

q6=MAPL_Eqsat(T,p)

call MAPL_EqsatSet(MixingRatio=.true.)

q2=MAPL_Eqsat(T)
q5=MAPL_Eqsat(T,p)

do i=140,340
   write(20,*)  t(i),q3(i),q2(i),q6(i),q5(i)
enddo


stop

end