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_S.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_S.pro

    r2740 r2751  
    1 pro ts_mean_S, sEXP1, sEXP2, year1, year2, s_iodir_data, POSTSCRIPT = postscript, _extra = ex 
     1pro std_ts_S, 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  vsal = getenv('VAR1_S')      &   prefix = getenv('V1S_PREF')         &   suffix = getenv('V1S_SUFF') 
     10  vssh = getenv('VAR1_SSH')    &   sshprefix = getenv('V1SSH_PREF')    &   sshsuffix = getenv('V1SSH_SUFF') 
     11; get exp2 info 
     12  vsal2 = getenv('VAR2_S')     &   prefix2 = getenv('V2S_PREF')        &   suffix2 = getenv('V2S_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_S_'+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_S =  isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato) 
     25  Lev = read_ncdf(getenv('VAR_SAL_3D'), filename = std_file_Levitus_S, /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, vsal, date1, date2, prefix, suffix, ts_Sal, ts_z $ 
     34               , WITHSSH = vssh, SSHPREFIX = sshprefix, SSHSUFFIX = sshsuffix, LEVZ = levz 
     35 
     36  title = prefix+'!C'+blabla+' '+d1_d2 
     37  pltt, ts_Sal, 't', 34.54, 34.76, date1, date2, /REMPLI $ 
     38        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, _extra = ex 
     39 
     40  title = prefix+' - Levitus!C'+blabla+' '+d1_d2 
     41  pltt, ts_Sal - Lev_xyz, 't', -.05, .05, 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+' - Levitus!C ZT-plot (0-1500m) '+blabla+' '+d1_d2 
     49  pltt, ts_z, 'zt', -.5, .5, INTER = .05, 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+' - Levitus!C ZT-plot (0-6000m) '+blabla+' '+d1_d2 
     53  pltt, ts_z, 'zt', -.5, .5, INTER = .05, 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/' 
     56  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ] 
     57  if KEYWORD_SET(postscript) then closeps 
    1158 
     59  IF prefix NE prefix2 THEN BEGIN 
    1260 
    13   if KEYWORD_SET(POSTSCRIPT) then begin 
    14      openps, sEXP1+'_Salinity_mean.ps', portrait = 1 
    15   endif 
    16    
    17   filename = std_iodir_climato+'data_1y_salinity_nomask.nc' 
    18   Lev = read_ncdf('vosaline', filename = filename, /nostruct ) 
    19   Lev_xyz = moyenne(Lev, 'xyz',  mask2d = masknp) 
    20   Lev_z = moyenne(Lev, 'xy', mask2d = masknp ,/KEEPBOTTOM) 
     61    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 
     62; read exp2 data 
     63    tsave = time 
     64    std_ts_read, vsal2, date1_2, date2_2, prefix2, suffix2, ts_Sal2, ts_z2 $ 
     65                 , WITHSSH = vssh2, SSHPREFIX = sshprefix2, SSHSUFFIX = sshsuffix2, LEVZ = levz 
     66    time = tsave   &   IF n_elements(time) NE jpt THEN stop 
    2167 
    22   ts_Sal = 0. 
    23   ts_Sal_2 = 0. 
    24   ts_Sal_Lev = 0. 
    25   ts_Sal_Lev_2 = 0. 
    26   ts_Time = 0. 
    27   ts_z = fltarr(jpk) 
    28   ts_z_2 = fltarr(jpk) 
     68    if KEYWORD_SET(postscript) then openps, filename+'_3.ps', portrait = 1 
    2969 
    30   masknp = read_ncdf('tmaskutil', file = 'mesh_mask_partial.nc', /nostruct) 
     70    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 
     71    pltt, ts_Sal - ts_Sal2, 't', -0.1, 0.1, date1, date2, /REMPLI $ 
     72          , COLOR = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /PORTRAIT, window = 2, _extra = ex 
    3173 
    32   for year=year1,year2 do begin 
    33       
    34      date1 = year * 10000L + 101 
    35      date2 = year * 10000L + 1231 
    36      filename = getname(s_iodir_data,sEXP1,'1Y','grid_T',date1) 
    37      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', -0.1, 0.1, INTER = .01, date1, date2, /REMPLI, style = 'so0so' $  
     76          , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 
    3877 
    39      if filename eq '' then stop  
    40      Sal   = read_ncdf('so',date1 ,date2 , filename = filename, /nostruct) 
    41      Sal_2 = read_ncdf('so',date1 ,date2 , filename = filename2, /nostruct) 
    42      SSH   =  read_ncdf('zos',date1 ,date2 , filename = filename, /nostruct) 
    43      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 
    4480 
    45      ts_Sal = [ ts_Sal, moyenne(Sal, 'xyz', mask2d = masknp, ssh = SSH) ] 
    46      ts_Sal_2 = [ ts_Sal_2, moyenne(Sal_2, 'xyz', mask2d = masknp, ssh = SSH_2) ] 
    47      ts_Sal_Lev = [ ts_Sal_Lev, moyenne(Sal, 'xyz', mask2d = masknp, ssh = SSH) - Lev_xyz ] 
    48      ts_Sal_Lev_2 = [ ts_Sal_Lev_2, moyenne(Sal_2, 'xyz', mask2d = masknp, ssh = SSH_2) - Lev_xyz ] 
    49      ts_Time = [ ts_Time, Time] ; Time = date 
     81  ENDIF 
    5082 
    51      ts_z = [ [ts_z], [moyenne(Sal, 'xy',/KEEPBOTTOM, mask2d = masknp, ssh = SSH)- Lev_z] ] 
    52      ts_z_2 = [ [ts_z_2], [moyenne(Sal_2, 'xy',/KEEPBOTTOM, mask2d = masknp, ssh = SSH_2)- Lev_z] ] 
    53  
    54   endfor 
    55    
    56   ts_Sal = ts_Sal[1:*]             ; remove first record of 0 
    57   ts_Sal_2 = ts_Sal_2[1:*]         ; remove first record of 0 
    58   ts_Sal_Lev = ts_Sal_Lev[1:*]     ; remove first record of 0 
    59   ts_Sal_Lev_2 = ts_Sal_Lev_2[1:*] ; remove first record of 0 
    60    
    61   ts_z = ts_z[*,1:*]            ; remove first record of 0 
    62   ts_z_2 = ts_z_2[*,1:*]            ; remove first record of 0 
    63    
    64   Time = ts_Time[1:*]           ; remove first record of 0 
    65   jpt = n_elements(Time) 
    66  
    67  
    68   title = sEXP1+' Global Annual Mean Salinity --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    69   pltt, ts_Sal, 't', 34.54, 34.76, INTER = .2, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $ 
    70        , small = [1,2,1], YTITLE='psu', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex 
    71  
    72   title = sEXP1+' Global Annual Mean Salinity - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    73   pltt, ts_Sal_Lev, 't', -.05, .05, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $ 
    74        , small = [1,2,2], YTITLE='psu',  TITLE = title, FORMAT = '(I3)', /NOERASE, _extra = ex 
    75  
    76 ;;;;;;;;;;;;;;;;; new page ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 
    77 ERASE 
    78   
    79   title = sEXP1+' ZT-plot Salinity Mean (0-1500m) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    80   pltt,ts_z ,'zt', -.5, .5, INTER = .025, year1 * 10000L + 101, year2 * 10000L + 1231 , /REMPLI $  
    81        , small = [1,2,1], YTITLE='psu', STYLE = 'so0so',  TITLE = title, boxzoom = 1500., /PORTRAIT, _extra = ex 
    82  
    83   title = sEXP1+' ZT-plot Salinity Mean (0-6000) - Levitus --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    84   pltt,ts_z,'zt', -.5, .5, INTER = .025, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $  
    85        , small = [1,2,2], YTITLE='psu', STYLE = 'so0so',  TITLE = title, boxzoom = 6000., zoom = 1500.,  /NOERASE, _extra = ex 
    86    
    87  
    88 if sEXP2 ne sEXP1 then begin 
    89  
    90 ERASE 
    91  
    92    title = sEXP1+' - '+sEXP2+' Global Annual Mean Salinity --- year '+strtrim(year1,1)+' - '+strtrim(year2,1) 
    93    pltt, ts_Sal - ts_Sal_2, 't', -1., 1., INTER = .2, year1 * 10000L + 101, year2 * 10000L + 1231, /REMPLI $ 
    94          , COLOR = 250, small = [1,2,1], YTITLE='psu', TITLE = title, FORMAT = '(I3)', /PORTRAIT, _extra = ex 
    95     
    96 endif 
    97  
    98    
    99   if KEYWORD_SET(POSTSCRIPT) then begin 
    100      closeps 
    101   endif 
    102    
    103    
    10483  return 
    10584end 
     85 
Note: See TracChangeset for help on using the changeset viewer.