Changeset 136

Timestamp:

06/03/09 16:41:57 (15 years ago)

Author:

ericg

Message:

Added update of make_wcurl.pro

File:

• trunk/procs/macros/make_wcurl.pro (modified) (2 diffs)

trunk/procs/macros/make_wcurl.pro

@@ -1 +1 @@
 ;
-; make Wind Stress curl
+; make curl of vector field
 ;
 FUNCTION make_wcurl, file_name, ncdf_db, BOXZOOM = boxzoom, TIME_1 = time_1,  TIME_2 =  time_2, ZMTYP = zmtyp, ALL_DATA = all_data
@@ -7 +7 @@
 @com_eg
 ;
+   IF debug_w THEN print, '   ENTER make_wcurl... '
+
+   var_name1 =  (strsplit(macro_base_fld, ',', ESCAPE = ' ', /EXTRACT))[0]
+   var_name2 =  (strsplit(macro_base_fld, ',', ESCAPE = ' ', /EXTRACT))[1]
+
+;  Build file_name2 if different
+
+   IF strpos (cmd1_back.grid, '#') NE -1 THEN BEGIN 
+      varpos = strpos(file_name,  var_name1)
+      file_name2= strmid(file_name, 0, varpos)+var_name2+'.nc'
+      IF debug_w THEN print, '   file_name2 = ', file_name2
+   ENDIF ELSE file_name2 = file_name
+
 ;
-; 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)
+;  Read the variables in the correspondent netcdf file
+   var_x =  nc_read(file_name, var_name1,  ncdf_db, BOXZOOM = boxzoom, TIME_1 = time_1,  TIME_2 =  time_2)
+   var_y =  nc_read(file_name2, var_name2,  ncdf_db, BOXZOOM = boxzoom, TIME_1 = time_1,  TIME_2 =  time_2)
 
-   wcurl = curl(taux.data, tauy.data)
-   ;domdef, old_boite
+; TO DO: ORCA grid C case
+;   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)
+
+;   for ORCA C grid vcurl = curl(var_x, var_y)
+
+   grille,mask,glam,gphi,gdep,nx,ny,nz,premierx,premiery,premierz,dernierx,derniery,dernierz   
+
+   zu = var_x.data
+   zv = var_y.data
+
+
+   IF jpt EQ 1 THEN BEGIN 
+      wcurl =  (shift(zv, -1, 0)-zv)/e1t[premierx:dernierx, premiery:derniery] + (zu-shift(zu, 0, -1))/e2t[premierx:dernierx, premiery:derniery]
+   ENDIF ELSE BEGIN 
+      e1t3d = reform(reform(e1t[premierx:dernierx, premiery:derniery], nx*ny)#replicate(1, jpt), nx, ny, jpt)
+      e2t3d = reform(reform(e2t[premierx:dernierx, premiery:derniery], nx*ny)#replicate(1, jpt), nx, ny, jpt)
+      wcurl =  (shift(zv, -1, 0, 0)-zv)/e1t3d + (zu-shift(zu, 0, -1, 0))/e2t3d
+   ENDELSE 
 
    field = {name: '', data: wcurl, legend: '', units: '', origin: '', dim: 0, direc:''}
    
-   field.origin = taux.origin
-   field.dim = taux.dim
-   field.direc = taux.direc
+   field.origin = var_x.origin
+   field.dim = var_x.dim
+   field.direc = var_x.direc
 
    return, field