! +-======-+ 
!  Copyright (c) 2003-2007 United States Government as represented by 
!  the Admistrator of the National Aeronautics and Space Administration.  
!  All Rights Reserved.
!  
!  THIS OPEN  SOURCE  AGREEMENT  ("AGREEMENT") DEFINES  THE  RIGHTS  OF USE,
!  REPRODUCTION,  DISTRIBUTION,  MODIFICATION AND REDISTRIBUTION OF CERTAIN 
!  COMPUTER SOFTWARE ORIGINALLY RELEASED BY THE UNITED STATES GOVERNMENT AS 
!  REPRESENTED BY THE GOVERNMENT AGENCY LISTED BELOW ("GOVERNMENT AGENCY").  
!  THE UNITED STATES GOVERNMENT, AS REPRESENTED BY GOVERNMENT AGENCY, IS AN 
!  INTENDED  THIRD-PARTY  BENEFICIARY  OF  ALL  SUBSEQUENT DISTRIBUTIONS OR 
!  REDISTRIBUTIONS  OF THE  SUBJECT  SOFTWARE.  ANYONE WHO USES, REPRODUCES, 
!  DISTRIBUTES, MODIFIES  OR REDISTRIBUTES THE SUBJECT SOFTWARE, AS DEFINED 
!  HEREIN, OR ANY PART THEREOF,  IS,  BY THAT ACTION, ACCEPTING IN FULL THE 
!  RESPONSIBILITIES AND OBLIGATIONS CONTAINED IN THIS AGREEMENT.
!  
!  Government Agency: National Aeronautics and Space Administration
!  Government Agency Original Software Designation: GSC-15354-1
!  Government Agency Original Software Title:  GEOS-5 GCM Modeling Software
!  User Registration Requested.  Please Visit http://opensource.gsfc.nasa.gov
!  Government Agency Point of Contact for Original Software:  
!  			Dale Hithon, SRA Assistant, (301) 286-2691
!  
! +-======-+ 


use MAPL_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