00001 subroutine adjust_matter_to_surface(new_rsf)
00002 use phys_constant, only : long, nmpt
00003 use def_matter, only : emd, rs
00004 use def_matter_parameter, only : emdc
00005 use grid_parameter, only : nrf, ntf, npf, &
00006 & ntfeq, ntfxy, npfyzp, npfxzp, npfxzm, &
00007 & ratio, NS_shape, EQ_point, r_surf
00008 implicit none
00009 integer :: it, ip
00010 real(long) :: new_rsf, rat
00011
00012
00013 rat = new_rsf/rs(ntfeq,npfxzp)
00014
00015 rs(0:ntf,0:npf) = rat*rs(0:ntf,0:npf)
00016
00017
00018 end subroutine adjust_matter_to_surface