rest new_s1_grav.sav ini xvel 0 yvel 0 ini xdis 0 ydis 0 ini udcoe c_udcoe ini fmod c_fmod f2mod c_f2mod ; interior nwwell C_nwdisc i=i_in j=j_in ;;;;;;;;;;;;;;;;;;; set solve config def c_nw while_stepping ;c_se = (sat(355,35)-c_rsat)/(1.0-c_rsat) ;c_s = sat(355,35) ;c_se = (sat(355,35)-c_rsat)/(1.0-c_rsat) ;c_nwpp = pp(355,35)+(c_vgp0_FENSSCROSS * (c_seˆ(-1.0/c_vga)-1.0)ˆ(1.0-c_vga)) c_nwpp1 = nwgpp(355,35) c_nwpp2 = nwgpp(355,36) c_nwpp3 = nwgpp(355,37) c_nwpp4 = nwgpp(355,38) c_nwpp5 = nwgpp(355,39) c_nwpp6 = nwgpp(355,40) end c_nw set flow on set sratio= 1.0E-8 hist 1 gwtime hist 5 gpp i=355 j=35 hist 6 gpp i=355 j=37 hist 7 gpp i=355 j=42 hist 8 gpp i=355 j=67 hist 9 gpp i=3429 j=8 hist 10 gpp i=3456 j=7 hist 11 gpp i=3456 j=6 hist 12 gpp i=3456 j=3 hist 13 ydisp i=3429 j=8 hist 14 c_nwpp1 hist 15 c_nwpp2 hist 16 c_nwpp3 hist 17 c_nwpp4 hist 18 c_nwpp5 hist 19 c_nwpp6 solve age age_sol_25year