Changeset 2751 for trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_msf.pro

Timestamp:

2011-04-26T10:45:52+02:00 (13 years ago)

Author:

smasson

Message:

large update of IDL_scripts for Is-ENES

File:

• trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_msf.pro (moved) (moved from trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/plt_msf.pro) (1 diff)

trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_msf.pro

@@ -1 @@
-pro plt_msf, V1, V2, sEXP1, sEXP2, start_end, POSTSCRIPT = postscript, SUBBASIN = subbasin, MASK_FILENAME = mask_filename,_extra = ex
+pro std_plot_msf, V1, V2, SUBBASIN = subbasin, POSTSCRIPT = postscript, _extra = ex
 
   compile_opt idl2, strictarrsubs
 
-  @common
+@common  
+@std_common
 
-  if KEYWORD_SET(SUBBASIN) then subname = subbasin 
-
-
-  if KEYWORD_SET(POSTSCRIPT) then begin
-     openps, '024_msf_'+subname+'_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.ps',  portrait = 1
-  endif
-
+  cdti3 = string(cnt, format = '(i3.3)')
+  print, cdti3 + ') ' + blabla
+  if KEYWORD_SET(SUBBASIN) then subname = subbasin else subname = "Glo"
+  filename = cdti3 + '_MSF_'+subname+'_'+std_file1_V
+  if std_file1_V NE std_file2_V then filename = filename + '_'+std_file2_V
+  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
+;
   if KEYWORD_SET(SUBBASIN) then begin
-     
-     CASE 1 of 
-        subbasin eq 'Atl' : var = 'atlmsk_nomed'
-        subbasin eq 'Ind' : var = 'indmsk'
-        subbasin eq 'Pac' : var = 'indpacmsk'
-        subbasin eq 'Glo' : var = 'glomsk_nomed'
-     ENDCASE
-                                ;
-     msk = read_ncdf( var, filename = mask_filename, /nostruct, _extra = ex )
-
-  endif else msk = tmask[*,*,0]
+    CASE subname of 
+      'Atl' : var = 'atlmsk_nomed'
+      'Ind' : var = 'indmsk'
+      'IndoPac' : var = 'indpacmsk'
+      'GloNoMed' : var = 'glomsk_nomed'
+    ENDCASE
+    msk = read_ncdf( var, filename = std_file_msksub, /nostruct, _extra = ex )
+  endif else msk = tmask[*, *, 0]
 ;
-
-  
   CASE subname OF
-     'Glo':lat_ext=[-80,90]
-     'Atl':lat_ext=[-30,90]
-     'Ind':lat_ext=[-30,30]
-     'Pac':lat_ext=[-30,70]
-     
+    'GloNoMed':lat_ext = [-80, 90]
+    'Glo':lat_ext = [-80, 90]
+    'Atl':lat_ext = [-30, 90]
+    'Ind':lat_ext = [-30, 30]
+    'IndoPac':lat_ext = [-30, 70]
   ENDCASE
-
- if subbasin eq 'Pac' then subname = "IndoPac"
   
   domdef, 0, 6000
 
-  mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma)
-  mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma)
-  title=sEXP1+' Meridional Stream Function '+subname+'/ year '+start_end
-  pltz, mm1, -20., 20., int = 1., boxzoom = [ind[0:1],lat_ext, 0, 5500], /xindex  , FORMAT = '(I3)' $
-        , small=[1,2,1], COAST_THICK=2 , zoom = 5500, maskdta = ma, /no_partial, TITLE = title, style = 'so0so'
+  mm1 = msf(V1.arr, indexboxzoom = ind, maskout = ma)
+  title = 'Meridional Stream Function, '+subname+'!C'+std_file1_V
+  pltz, mm1, 'yz', -20., 20., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)', /portrait $
+        , small = [1, 2, 1], COAST_THICK = 2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, style = 'so0so'
 
-  if sEXP1 ne sEXP2 then begin
-     
-     title=sEXP1+' - '+sEXP2+subname+'/ year '+start_end
-     pltz, mm1.arr-mm2.arr, -10., 10., int = 1., boxzoom = [ind[0:1],lat_ext, 0, 5500], /xindex,  FORMAT = '(I3)'  $
-           , small=[1,2,2], COAST_THICK=2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, /NOERASE, style = 'so0so'
-     
+  if std_file1_V NE std_file2_V then begin
+    mm2 = msf(V2.arr, indexboxzoom = ind, maskout = ma)
+    title = title+' - '+std_file2_V
+    pltz, mm1.arr-mm2.arr, 'yz', -10., 10., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex,  FORMAT = '(I3)'  $
+          , small = [1, 2, 2], COAST_THICK = 2, zoom = 5500, maskdta = ma, /no_partial, TITLE = title, /NOERASE, style = 'so0so'
   endif
   
-  title = + subname+ 'Meridional Stream Function'
-  text = [   title, '<img width="80%" src=/Users/sflod/TOOLS/my_IDL/PS/'+sEXP1+'_'+sEXP2+'/ps/024_msf_'+subname+'_'+sEXP1+'_'+sEXP2+'_y'+start_end+'.png  />  ' ]
-  putfile, '/tmp/prova.txt', text
-
-  if KEYWORD_SET(POSTSCRIPT) then begin
-     closeps
-  endif
-
+  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
+  if KEYWORD_SET(postscript) then closeps
 
   return