# ------------------------------------------------
# T, (u, omg) (赤道断面) の重ね描き
t = ape_new 1
l = ape_new 13

t_con = t.go("ml_t").cut(true,0,true)
t_pa = l.go("ml_t").cut(true,0,true)

u_con = t.go("ml_u").cut(true,0,true)
om_con = t.go("ml_om").cut(true,0,true)

u_pa = l.go("ml_u").cut(true,0,true)
om_pa = l.go("ml_om").cut(true,0,true)

t_anm = t_pa - t_con
t_anm = t_anm.lon_lotate
u = u_pa - u_con
u = u.lon_lotate
om = om_pa -om_con
om = om.lon_lotate

u_data  = NArray.sfloat( u.axis(0).pos.val.size ,
   u.axis(1).pos.val.size ).fill!(0.0)
om_data  = NArray.sfloat( om.axis(0).pos.val.size ,
   om.axis(1).pos.val.size ).fill!(0.0)

num_l = u.axis(0).pos.val.size/40
num_z = u.axis(1).pos.val.size/17

u.axis(0).pos.val.size.times{ |num|
      if (num % num_l) == 0
           u_data[num,true]  =  u.data.val[num,true]
           om_data[num,true]  =  om.data.val[num,true]
      elsif   
           u_data[num,true]  = 0.0
           om_data[num,true]  = 0.0
      end
}  

t.plot(t_anm.rename("ml_tuom_anm"),u_data,-om_data)



# ------------------------------------------------
# T, (u, v) (500 hPa) の重ね描き
t = ape_new 1
l = ape_new 13

t_con = t.go("ml_t")[true,true,5]
t_pa = l.go("ml_t")[true,true,5]

u_con = t.go("ml_u")[true,true,5]
v_con = t.go("ml_v")[true,true,5]

u_pa = l.go("ml_u")[true,true,5]
v_pa = l.go("ml_v")[true,true,5]

t_anm = t_pa - t_con
t_anm = t_anm.lon_lotate
u = u_pa - u_con
u = u.lon_lotate
v = v_pa -v_con
v = v.lon_lotate

u_data  = NArray.sfloat( u.axis(0).pos.val.size ,
   u.axis(1).pos.val.size ).fill!(0.0)
v_data  = NArray.sfloat( v.axis(0).pos.val.size ,
   v.axis(1).pos.val.size ).fill!(0.0)

num_l = u.axis(0).pos.val.size/40
num_z = u.axis(1).pos.val.size/20

u.axis(0).pos.val.size.times{ |num|
      if (num % num_l) == 0
           u_data[num,true]  =  u.data.val[num,true]
           v_data[num,true]  =  v.data.val[num,true]
      elsif   
           u_data[num,true]  = 0.0
           v_data[num,true]  = 0.0
      end
}  

u.axis(1).pos.val.size.times{ |num|
     if ((num+2) % num_z) == 0
           u_data[true,num]  =  u_data[true,num]
           v_data[true,num]  =  v_data[true,num]
     elsif
           u_data[true,num]  = 0.0
           v_data[true,num]  = 0.0
     end
}

t.plot(t_anm.rename("sh_tuv500_anm"),u_data,v_data)




# ------------------------------------------------
# gph, (u, v) (850 hPa) の重ね描き
t = ape_new 1
l = ape_new 13

phi_con = t.go("ml_phi")[true,true,2]
phi_pa = l.go("ml_phi")[true,true,2]

u_con = t.go("ml_u")[true,true,2]
v_con = t.go("ml_v")[true,true,2]

u_pa = l.go("ml_u")[true,true,2]
v_pa = l.go("ml_v")[true,true,2]

phi = phi_pa - phi_con
phi = phi.lon_lotate
u = u_pa - u_con
u = u.lon_lotate
v = v_pa -v_con
v = v.lon_lotate

u_data  = NArray.sfloat( u.axis(0).pos.val.size ,
   u.axis(1).pos.val.size ).fill!(0.0)
v_data  = NArray.sfloat( v.axis(0).pos.val.size ,
   v.axis(1).pos.val.size ).fill!(0.0)

num_l = u.axis(0).pos.val.size/40
num_z = u.axis(1).pos.val.size/20

u.axis(0).pos.val.size.times{ |num|
      if (num % num_l) == 0
           u_data[num,true]  =  u.data.val[num,true]
           v_data[num,true]  =  v.data.val[num,true]
      elsif   
           u_data[num,true]  = 0.0
           v_data[num,true]  = 0.0
      end
}  

u.axis(1).pos.val.size.times{ |num|
     if ((num+2) % num_z) == 0
           u_data[true,num]  =  u_data[true,num]
           v_data[true,num]  =  v_data[true,num]
     elsif
           u_data[true,num]  = 0.0
           v_data[true,num]  = 0.0
     end
}

t.plot(phi.rename("sh_phiuv850_anm"),u_data,v_data)


# ------------------------------------------------
# gph, (u, v) (250 hPa) の重ね描き
t = ape_new 1
l = ape_new 13

phi_con = t.go("ml_phi")[true,true,8]
phi_pa = l.go("ml_phi")[true,true,8]

u_con = t.go("ml_u")[true,true,8]
v_con = t.go("ml_v")[true,true,8]

u_pa = l.go("ml_u")[true,true,8]
v_pa = l.go("ml_v")[true,true,8]

phi = phi_pa - phi_con
phi = phi.lon_lotate
u = u_pa - u_con
u = u.lon_lotate
v = v_pa -v_con
v = v.lon_lotate

u_data  = NArray.sfloat( u.axis(0).pos.val.size ,
   u.axis(1).pos.val.size ).fill!(0.0)
v_data  = NArray.sfloat( v.axis(0).pos.val.size ,
   v.axis(1).pos.val.size ).fill!(0.0)

num_l = u.axis(0).pos.val.size/40
num_z = u.axis(1).pos.val.size/20

u.axis(0).pos.val.size.times{ |num|
      if (num % num_l) == 0
           u_data[num,true]  =  u.data.val[num,true]
           v_data[num,true]  =  v.data.val[num,true]
      elsif   
           u_data[num,true]  = 0.0
           v_data[num,true]  = 0.0
      end
}  

u.axis(1).pos.val.size.times{ |num|
     if ((num+2) % num_z) == 0
           u_data[true,num]  =  u_data[true,num]
           v_data[true,num]  =  v_data[true,num]
     elsif
           u_data[true,num]  = 0.0
           v_data[true,num]  = 0.0
     end
}

t.plot(phi.rename("sh_phiuv250_anm"),u_data,v_data)


# --------------------------------------------------


# ps, (u, v) の重ね描き
t = ape_new 23
l = ape_new 20

ps_con = t.go("sh_ps")
ps_pa = l.go("sh_ps")

t = ape_new 1
l = ape_new 13

u_con = t.go("ml_u")[true,true,0]
v_con = t.go("ml_v")[true,true,0]

u_pa = l.go("ml_u")[true,true,0]
v_pa = l.go("ml_v")[true,true,0]

ps = ps_pa - ps_con
ps = ps.lon_lotate
u = u_pa - u_con
u = u.lon_lotate
v = v_pa -v_con
v = v.lon_lotate

u_data  = NArray.sfloat( u.axis(0).pos.val.size ,
   u.axis(1).pos.val.size ).fill!(0.0)
v_data  = NArray.sfloat( v.axis(0).pos.val.size ,
   v.axis(1).pos.val.size ).fill!(0.0)

num_l = u.axis(0).pos.val.size/40
num_z = u.axis(1).pos.val.size/20

u.axis(0).pos.val.size.times{ |num|
      if (num % num_l) == 0
           u_data[num,true]  =  u.data.val[num,true]
           v_data[num,true]  =  v.data.val[num,true]
      elsif   
           u_data[num,true]  = 0.0
           v_data[num,true]  = 0.0
      end
}  

u.axis(1).pos.val.size.times{ |num|
     if ((num+2) % num_z) == 0
           u_data[true,num]  =  u_data[true,num]
           v_data[true,num]  =  v_data[true,num]
     elsif
           u_data[true,num]  = 0.0
           v_data[true,num]  = 0.0
     end
}

t.plot(ps.rename("sh_psuv_anm"),u_data,v_data)