New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2751 for trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_T.pro – NEMO

Ignore:
Timestamp:
2011-04-26T10:45:52+02:00 (13 years ago)
Author:
smasson
Message:

large update of IDL_scripts for Is-ENES

File:
1 moved

Legend:

Unmodified
Added
Removed

  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_T.pro

    r2740 r2751  
    1 pro ts_mean_T, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex 
     1pro std_ts_T, masknp, s_iodir_data, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
    44   
    5   @common 
     5@common 
     6@std_common 
     7 
     8; get exp1 info 
     9  vtemp = getenv('VAR1_T')     &   prefix = getenv('V1T_PREF')         &   suffix = getenv('V1T_SUFF') 
     10  vssh = getenv('VAR1_SSH')    &   sshprefix = getenv('V1SSH_PREF')    &   sshsuffix = getenv('V1SSH_SUFF') 
     11; get exp2 info 
     12  vtemp2 = getenv('VAR2_T')    &   prefix2 = getenv('V2T_PREF')        &   suffix2 = getenv('V2T_SUFF') 
     13  vssh2 = getenv('VAR2_SSH')   &   sshprefix2 = getenv('V2SSH_PREF')   &   sshsuffix2 = getenv('V2SSH_SUFF') 
     14; 
     15  cdti3 = string(cnt, format = '(i3.3)') 
     16  print, cdti3 + ') ' + blabla 
     17  filename = cdti3 + '_ts_T_'+prefix 
     18  if prefix NE prefix2 then filename = filename + '_'+prefix2 
     19  if KEYWORD_SET(postscript) then openps, filename+'_1.ps', portrait = 1 
     20 
     21  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 
     22 
     23; read levitus data 
     24  std_file_Levitus_T = isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato) 
     25  Lev = read_ncdf(getenv('VAR_TEMP_3D'), filename = std_file_Levitus_T, /nostruct )  
     26  Lev_xyz = moyenne(Lev, 'xyz', mask2d = masknp) 
     27  levz = moyenne(temporary(Lev), 'xy', mask2d = masknp, /KEEPBOTTOM) 
     28; 
     29  iodir = std_iodir_data 
     30; 
     31 
     32; read exp1 data 
     33  std_ts_read, vtemp, date1, date2, prefix, suffix, ts_Temp, ts_z $ 
     34               , WITHSSH = vssh, SSHPREFIX = sshprefix, SSHSUFFIX = sshsuffix, LEVZ = levz 
     35 
     36  title = prefix+' '+d1_d2+'!C'+blabla 
     37  pltt, ts_Temp, 't', 1., 4.5, date1, date2, /REMPLI $ 
     38        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex 
     39 
     40  title = prefix+' '+d1_d2+' - Levitus!C'+blabla 
     41  pltt, ts_Temp - Lev_xyz, 't', -1., 1., date1, date2, /REMPLI $ 
     42        , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /NOERASE, _extra = ex 
     43 
     44  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_1.png  />  ' ] 
     45  if KEYWORD_SET(postscript) then closeps 
     46  if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
     47 
     48  title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-1500m) '+blabla 
     49  pltt, ts_z, 'zt', -2., 2., inter = .2, date1, date2, /REMPLI, style = 'so0so' $  
     50        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex 
    651   
    7   @initorca2 
     52  title = prefix+' '+d1_d2+' - Levitus!C ZT-plot (0-6000m) '+blabla 
     53  pltt, ts_z, 'zt', -2., 2., inter = .2, date1, date2, /REMPLI, style = 'so0so' $  
     54        , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 
    855   
    9   std_iodir_climato = '/Users/sflod/idl_PLOTS/CLIMATOLOGIES/' 
    10 ;  std_iodir_interan = '/Users/sflod/idl_PLOTS/INTERAN/' 
    11   std_iodir_data = '/Users/sflod/idl_PLOTS/INTERAN/mercure/' 
     56  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ] 
     57  if KEYWORD_SET(postscript) then closeps 
    1258 
     59  IF prefix NE prefix2 THEN BEGIN 
    1360 
    14   if KEYWORD_SET(POSTSCRIPT) then begin 
    15      openps, sEXP1+'_Temperature_mean.ps', portrait = 1 
    16   endif 
    17    
    18   masknp = read_ncdf('tmaskutil', file = 'mesh_mask_partial.nc', /nostruct) 
     61    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 
     62; read exp2 data 
     63    tsave = time 
     64    std_ts_read, vtemp2, date1_2, date2_2, prefix2, suffix2, ts_Temp2, ts_z2 $ 
     65                 , WITHSSH = vssh2, SSHPREFIX = sshprefix2, SSHSUFFIX = sshsuffix2, LEVZ = levz 
     66    time = tsave   &   IF n_elements(time) NE jpt THEN stop 
    1967 
    20   filename = std_iodir_climato+'data_1y_potential_temperature_nomask.nc' 
    21   Lev = read_ncdf('votemper', filename = filename, /nostruct ) 
    22   Lev_xyz = moyenne(Lev, 'xyz', mask2d = masknp) 
    23   Lev_z = moyenne(Lev, 'xy', mask2d = masknp,/KEEPBOTTOM) 
     68    if KEYWORD_SET(postscript) then openps, filename+'_3.ps', portrait = 1 
    2469 
    25   ts_Temp = 0. 
    26   ts_Temp_2 = 0. 
    27   ts_Temp_Lev = 0. 
    28   ts_Temp_Lev_2 = 0. 
    29   ts_Time = 0. 
    30   ts_z = fltarr(jpk) 
    31   ts_z_2 = fltarr(jpk) 
    32    
    33   for year=year1,year2 do begin 
    34       
    35      date1 = year * 10000L + 101L 
    36      date2 = year * 10000L + 1231L 
    37   
    38 print,s_iodir_data 
     70    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 
     71    pltt, ts_Temp - ts_Temp2, 't', -1., 1., date1, date2, /REMPLI $ 
     72          , COLOR = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, window = 2, _extra = ex 
    3973 
    40      filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1) 
    41      filename2 = getname(s_iodir_data,sEXP2,'1Y','grid_T',date1) 
     74    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C ZT-plot (0-6000m) '+blabla 
     75    pltt, ts_z - ts_z2, 'zt', -1., 1., inter = .1, date1, date2, /REMPLI, style = 'so0so' $  
     76          , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 
    4277 
    43      if filename eq '' then stop  
    44      Temp   = read_ncdf('thetao',date1 ,date2 , filename = filename, /nostruct) 
    45      Temp_2 = read_ncdf('thetao',date1 ,date2 , filename = filename2, /nostruct) 
    46      SSH    =  read_ncdf('zos',date1 ,date2 , filename = filename, /nostruct) 
    47      SSH_2  =  read_ncdf('zos',date1 ,date2 , filename = filename2, /nostruct) 
     78    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_3.png  />  ' ] 
     79    if KEYWORD_SET(postscript) then closeps 
     80     
     81  ENDIF 
    4882 
    49      ts_Temp = [ ts_Temp, moyenne(Temp, 'xyz', mask2d = masknp,ssh = SSH) ] 
    50      ts_Temp_2 = [ ts_Temp_2, moyenne(Temp_2, 'xyz', mask2d = masknp,ssh = SSH_2) ] 
    51      ts_Temp_Lev = [ ts_Temp_Lev, moyenne(Temp, 'xyz', mask2d = masknp,ssh = SSH) - Lev_xyz ] 
    52      ts_Temp_Lev_2 = [ ts_Temp_Lev_2, moyenne(Temp_2, 'xyz', mask2d = masknp,ssh = SSH_2) - Lev_xyz ] 
    53      ts_Time = [ ts_Time, Time] ; Time = date 
    54  
    55      ts_z = [ [ts_z], [moyenne(Temp, 'xy', mask2d = masknp,/KEEPBOTTOM, ssh = SSH)- Lev_z] ] 
    56      ts_z_2 = [ [ts_z_2], [moyenne(Temp_2, 'xy', mask2d = masknp,/KEEPBOTTOM, ssh = SSH_2)- Lev_z] ] 
    57  
    58   endfor 
    59  
    60   ts_Temp = ts_Temp[1:*]        ; remove first record of 0  
    61   ts_Temp_2 = ts_Temp_2[1:*]        ; remove first record of 0 
    62   ts_Temp_Lev = ts_Temp_Lev[1:*] ; remove first record of 0 
    63   ts_Temp_Lev_2 = ts_Temp_Lev_2[1:*] ; remove first record of 0 
    64   ts_z = ts_z[*,1:*]             ; remove first record of 0 
    65   ts_z_2 = ts_z_2[*,1:*]             ; remove first record of 0 
    66    
    67   Time = ts_Time[1:*]           ; remove first record of 0 
    68   jpt = n_elements(Time) 
    69    
    70   title = sEXP1+' Global Annual Mean Temperature --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    71   pltt, ts_Temp, 't', 1., 4.5, INTER = .5, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $ 
    72        , small = [1,2,1], YTITLE='Deg C', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex 
    73  
    74   title = sEXP1+' Global Annual Mean Temperature - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    75   pltt,ts_Temp_Lev, 't', -1., 1., INTER = .25,  year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $ 
    76        , small = [1,2,2], YTITLE='Deg C', TITLE = title, FORMAT = '(I3)', /NOERASE, _extra = ex 
    77  
    78   ;;;;;;;;;;;;;;;;; new page  ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 
    79  
    80 ERASE 
    81  
    82   title = sEXP1+' ZT-plot Temperature Mean (0-1500m) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    83   pltt,ts_z ,'zt', -2., 2., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $  
    84        , small = [1,2,1], YTITLE='Deg C', TITLE = title, boxzoom = 1500., /NOERASE, _extra = ex 
    85    
    86   title = sEXP1+' ZT-plot Temperature Mean (0-6000m) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    87   pltt,ts_z ,'zt', -2., 2., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $  
    88        , small = [1,2,2], YTITLE='Deg C', TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 
    89    
    90 if sEXP2 ne sEXP1 then begin 
    91    
    92 ERASE 
    93   
    94    title = sEXP1+' - '+sEXP2+' Global Annual Mean Temperature --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    95    pltt, ts_Temp - ts_Temp_2, 't', -1., 1., year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $ 
    96          , COLOR = 250, small = [1,2,1], YTITLE='Deg C', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex 
    97     
    98 endif 
    99  
    100    
    101   if KEYWORD_SET(POSTSCRIPT) then begin 
    102      closeps 
    103   endif 
    104    
    105    
    10683  return 
    10784end 
Note: See TracChangeset for help on using the changeset viewer.