! the latitudes need to be flipped with flip_lat.f
!
!use dissipation_00.0kyrBP_ITfricPD_I0.02_191556.nc
!repeat/i=1:1439 (list/file=ITfricPD_M2_inverse_lat2.txt/append/form=(2881(2x,E14.7))/nohead/order=y diss_m2)
!repeat/i=1:1439 (list/file=ITfricPD_S2_inverse_lat2.txt/append/form=(2881(2x,E14.7))/nohead/order=y diss_s2)
!repeat/i=1:1439 (list/file=ITfricPD_K1_inverse_lat2.txt/append/form=(2881(2x,E14.7))/nohead/order=y diss_k1)
!repeat/i=1:1439 (list/file=ITfricPD_O1_inverse_lat2.txt/append/form=(2881(2x,E14.7))/nohead/order=y diss_o1)

use dissipation_21.0kyrBP_ITfricLGM0_I0.02_583323.nc
let diss = if (mz eq 1) then diss_m2 else 0
repeat/i=1:1439 (list/file=ITfricLGM_M2_inverse_lat2.txt/append/form=(2881(2x,E14.7))/nohead/order=y diss)
let diss = if (mz eq 1) then diss_s2 else 0
repeat/i=1:1439 (list/file=ITfricLGM_S2_inverse_lat2.txt/append/form=(2881(2x,E14.7))/nohead/order=y diss)
let diss = if (mz eq 1) then diss_k1 else 0
repeat/i=1:1439 (list/file=ITfricLGM_K1_inverse_lat2.txt/append/form=(2881(2x,E14.7))/nohead/order=y diss)
let diss = if (mz eq 1) then diss_o1 else 0
repeat/i=1:1439 (list/file=ITfricLGM_O1_inverse_lat2.txt/append/form=(2881(2x,E14.7))/nohead/order=y diss)