00001 Subroutine mit_ene2prho(ene,pre,rho) 00002 use def_qeos_parameter 00003 implicit none 00004 real(8),intent(inout) :: rho, pre, ene 00005 pre=aq*(ene-enesurf_gcm1) 00006 rho=0.0d0 00007 END Subroutine mit_ene2prho