! This is lo_lat.jnl
! Find the low latitude ocean inventory in total MP

def symbol inc = $1%1%

list ($itr)

let isfile = `test_opendap("sens_rise/test($itr)/tavg.01951.01.01.nc")`

if `isfile EQ 0` then

  use "sens_rise/test($itr)/tavg.01951.01.01.nc"

  let totmp=o_mp+o_mpa+o_mpp
  let posmp = if totmp GT 0 then totmp else 0
  let lolat = posmp[x=@din,y=40S:40N@din,l=61,k=@din] !particles in lolat ocean

  let prop = lolat/posmp[x=@din,y=@din,z=@din,l=61]

list/format=(F6.3)/nohead/append/file=metric.txt prop

SPAWN grep -Po "pinp=\K[^,]*" sens_rise/test($itr)/control.in >> param1.txt
SPAWN grep -Po "prise=\K[^,]*" sens_rise/test($itr)/control.in >> param2.txt

  canc data 1
endif

def symbol itr = `($itr)+($inc)`