source: trunk/procs/macros/make_wdiv.pro @ 87

;
; make Wind Stress divergence
;
FUNCTION make_wdiv, file_name, ncdf_db, BOXZOOM = boxzoom, TIME_1 = time_1,  TIME_2 =  time_2, ZMTYP = zmtyp, ALL_DATA = all_data

@common
@com_eg
;
;
; Read taux and tauy
;
   ;old_boite = [lon1, lon2, lat1, lat2, prof1, prof2]
   ;domdef
   vargrid = 'U'
   new_file = new_filename(file_name, 'T', vargrid)
   taux = nc_read(new_file,'sozotaux', ncdf_db, BOXZOOM = boxzoom, TIME_1 = time_1,  TIME_2 =  time_2, ALL_DATA = all_data)
   vargrid = 'V'
   new_file = new_filename(file_name, 'T', vargrid)
   tauy = nc_read(new_file,'sometauy', ncdf_db, BOXZOOM = boxzoom, TIME_1 = time_1,  TIME_2 =  time_2, ALL_DATA = all_data)

   wdiv = div(taux.data, tauy.data)
   ;domdef, old_boite

   field = {name: '', data: wdiv, legend: '', units: '', origin: '', dim: 0, direc:''}
   
   field.origin = taux.origin
   field.dim = taux.dim
   field.direc = taux.direc

   return, field
END