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

    r2740 r2751  
    1 pro std_ts, std_EXP1, std_EXP2, year1, year2, _extra=ex 
    2  
    3  
    4 PRINT, 'usage: postscript ' 
    5 PRINT, '========================================================================' 
    6 PRINT, 'std_ts, ''std_EXP1'', ''std_EXP2'', ''year_start'', ''year_end'', /postscript' 
    7 PRINT, '========================================================================' 
    8 PRINT, 'std_ts, ''core2000'', ''core2000'', ''1'', ''2000'', /postscript' 
    9 PRINT, '========================================================================' 
    10  
     1pro std_ts_all, doplot = doplot, _extra=ex 
    112 
    123  compile_opt idl2, strictarrsubs 
     4 
     5@common 
     6@std_common 
    137 
    148  PRINT, '' 
    … …  
    1913  PRINT, '  ############################################' 
    2014  PRINT, '' 
     15; 
     16  std_iodir_data    = isadirectory(getenv('DIR_DATA'), title = 'path of data in NetCdf format') 
     17  std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'), title = 'path of climatological data') 
     18  std_iodir_mask    = isadirectory(getenv('DIR_MASK'), title = 'path of mask files (ex: subbasins)') 
     19; meshmask 
     20  std_file_mesh = isafile(getenv('FILE_MESH_MASK'), title = 'mesh_mask', iodir = std_iodir_mask) 
     21  std_file_msksub = isafile(getenv('FILE_MASK_SUBDOMAIN'), title = 'sub-bassin masks', iodir = std_iodir_mask) 
     22; load the grid 
     23  load_orca, std_file_mesh   
     24; reading variables 
     25  masknp = read_ncdf('tmaskutil', file = std_file_mesh, /nostruct, /cont_nofill) 
    2126 
     27  date1   = long(getenv('DATE1'))     &   date2   = long(getenv('DATE2')) 
     28  date1_2 = long(getenv('DATE1_2'))   &   date2_2 = long(getenv('DATE2_2')) 
    2229 
    23 @common 
    24 @initorca2 
    25 @std_com 
    26  
    27  
    28 ;  std_iodir_data    = '/Users/sflod/idl_PLOTS/DATA_STORE/RUN_INTERAN/mercure/' ; path of data in NetCdf format 
    29   std_iodir_data    = '/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/brodie/' ; path of data in NetCdf format 
    30  
     30  allrec =  1 - keyword_set(long(getenv('READ_ONLY_FIRST_RECORD'))) 
    3131 
    3232;######################################################################### 
    33 ;######################  TIME SERIES   ################################ 
     33;######################  STANDARD PLOTS   ################################ 
    3434;######################################################################### 
    3535 
     36  IF keyword_set(doplot) EQ 0 THEN doplot = 0 
     37 
    3638; fixed color tabled 
    37   lct,64 
     39  lct, 64 
     40  cnt = 0 
     41  htmltxt = '' 
     42; 
     43  cnt = cnt+1   &   blabla = 'Global Mean Temperature' 
     44  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_T, masknp, POSTSCRIPT = postscript, _extra = ex 
     45; 
     46  cnt = cnt+1   &   blabla = 'Global Mean Salinity' 
     47  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_S, masknp, POSTSCRIPT = postscript, _extra = ex 
     48; 
     49  cnt = cnt+1   &   blabla = 'Global Mean SSH' 
     50  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SSH, masknp, POSTSCRIPT = postscript, _extra = ex 
     51; 
     52  cnt = cnt+1   &   blabla = 'Global Mean Q net' 
     53  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Q, masknp, POSTSCRIPT = postscript, _extra = ex 
     54; 
     55  cnt = cnt+1   &   blabla = 'Global Mean EMP' 
     56  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_EMP, masknp, POSTSCRIPT = postscript, _extra = ex 
     57; 
     58  cnt = cnt+1   &   blabla = 'Drake Transport' 
     59  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Drake, masknp, POSTSCRIPT = postscript, _extra = ex 
     60; 
     61  cnt = cnt+1   &   blabla = 'Max AMOC' 
     62  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_AMOC, masknp, POSTSCRIPT = postscript, _extra = ex 
     63; 
     64  cnt = cnt+1   &   blabla = 'Sea-Ice cover' 
     65  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE, masknp, POSTSCRIPT = postscript, _extra = ex 
    3866 
    39   
    40  
    41   ;; Temperature Mean 
    42   ;; ======================================= 
    43   PRINT, ' ' 
    44   PRINT, '         Temperature Mean' 
    45   PRINT, '1 ' 
    46   ts_mean_T, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex 
    47  
    48   ;; Salinity Mean 
    49   ;; ======================================= 
    50   PRINT, ' ' 
    51   PRINT, '         Salinity Mean' 
    52   PRINT, '2 ' 
    53   ts_mean_S, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex 
    54    
    55   ;; Sea Surface Height 
    56   ;; ======================================= 
    57   PRINT, ' ' 
    58   PRINT, '         SSH' 
    59   PRINT, '3 ' 
    60   ts_mean_SSH, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex 
    61  
    62    ;; Heat flux Mean 
    63   ;; ======================================= 
    64   PRINT, ' ' 
    65   PRINT, '         Heat Mean' 
    66   ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 
    67   ; commented because not good in iodef.xml file for 33beta!!!!!!!!!!! 
    68   ;PRINT, '4 ' 
    69   ts_mean_Q, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex 
    70  
    71   ;; EMP Mean 
    72   ;; ======================================= 
    73   PRINT, ' ' 
    74   PRINT, '         EMP Mean' 
    75   ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 
    76   ; commented because not good in iodef.xml file for 33beta!!!!!!!!!!! 
    77   ;PRINT, '5 ' 
    78   ts_mean_EMP, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex 
    79  
    80  
    81   ;; Drake transport 
    82   ;; ======================================= 
    83   PRINT, ' ' 
    84   PRINT, '              Drake Transport' 
    85   ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 
    86   ; commented because not good in iodef.xml file for 33beta!!!!!!!!!!! 
    87   ;PRINT, '6 ' 
    88   ts_mean_U, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex 
    89  
    90   ;; Max AMOC 
    91   ;; ======================================= 
    92   PRINT, ' ' 
    93   PRINT, '              Max AMOC' 
    94   ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 
    95   ; commented because not good in iodef.xml file for 33beta!!!!!!!!!!! 
    96   ;PRINT, '7 ' 
    97   ts_mean_V, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex 
    98  
    99   ;; Sea-Ice Area  
    100   ;; ======================================= 
    101   PRINT, ' ' 
    102   PRINT, '              Sea-Ice Area' 
    103   PRINT, '8 ' 
    104   ts_mean_ICE, std_EXP1, std_EXP2, year1, year2, std_iodir_data, _extra = ex 
     67  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_ts_html_body.txt', htmltxt[1:*] 
    10568 
    10669  return 
Note: See TracChangeset for help on using the changeset viewer.