[2751] | 1 | pro std_ts_all, doplot = doplot, _extra=ex |
---|
[2709] | 2 | |
---|
[2751] | 3 | compile_opt idl2, strictarrsubs |
---|
[2709] | 4 | |
---|
[2751] | 5 | @common |
---|
| 6 | @std_common |
---|
[2709] | 7 | |
---|
| 8 | PRINT, '' |
---|
| 9 | PRINT, ' ############################################' |
---|
| 10 | PRINT, '' |
---|
| 11 | PRINT, ' plot of all TIME SERIES' |
---|
| 12 | PRINT, '' |
---|
| 13 | PRINT, ' ############################################' |
---|
| 14 | PRINT, '' |
---|
[2751] | 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) |
---|
[4320] | 26 | ; climatologies |
---|
| 27 | std_file_Levitus_T = isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato) |
---|
| 28 | std_file_Levitus_S = isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato) |
---|
| 29 | std_file_reynolds = isafile(getenv('FILE_SST'), title = 'Reynolds', iodir = std_iodir_climato) |
---|
| 30 | std_file_oaflux = isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato) |
---|
| 31 | std_file_mld = isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato) |
---|
| 32 | std_file_ice = isafile(getenv('FILE_ICE'), title = 'ICE', iodir = std_iodir_climato) |
---|
| 33 | std_file_snow_arc = isafile(getenv('FILE_SNOW_ARC'), title = 'SNOW_ARC', iodir = std_iodir_climato) |
---|
| 34 | std_file_snow_ant = isafile(getenv('FILE_SNOW_ANT'), title = 'SNOW_ANT', iodir = std_iodir_climato) |
---|
[2709] | 35 | |
---|
[4320] | 36 | |
---|
[2751] | 37 | date1 = long(getenv('DATE1')) & date2 = long(getenv('DATE2')) |
---|
| 38 | date1_2 = long(getenv('DATE1_2')) & date2_2 = long(getenv('DATE2_2')) |
---|
[2709] | 39 | |
---|
[2751] | 40 | allrec = 1 - keyword_set(long(getenv('READ_ONLY_FIRST_RECORD'))) |
---|
[2709] | 41 | |
---|
| 42 | ;######################################################################### |
---|
[2751] | 43 | ;###################### STANDARD PLOTS ################################ |
---|
[2709] | 44 | ;######################################################################### |
---|
| 45 | |
---|
[2751] | 46 | IF keyword_set(doplot) EQ 0 THEN doplot = 0 |
---|
| 47 | |
---|
[2709] | 48 | ; fixed color tabled |
---|
[2751] | 49 | lct, 64 |
---|
| 50 | cnt = 0 |
---|
| 51 | htmltxt = '' |
---|
| 52 | ; |
---|
| 53 | cnt = cnt+1 & blabla = 'Global Mean Temperature' |
---|
| 54 | IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_T, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
| 55 | ; |
---|
| 56 | cnt = cnt+1 & blabla = 'Global Mean Salinity' |
---|
| 57 | IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_S, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
| 58 | ; |
---|
| 59 | cnt = cnt+1 & blabla = 'Global Mean SSH' |
---|
| 60 | IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SSH, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
| 61 | ; |
---|
| 62 | cnt = cnt+1 & blabla = 'Global Mean Q net' |
---|
| 63 | IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Q, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
| 64 | ; |
---|
| 65 | cnt = cnt+1 & blabla = 'Global Mean EMP' |
---|
| 66 | IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_EMP, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
| 67 | ; |
---|
| 68 | cnt = cnt+1 & blabla = 'Drake Transport' |
---|
| 69 | IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Drake, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
| 70 | ; |
---|
| 71 | cnt = cnt+1 & blabla = 'Max AMOC' |
---|
| 72 | IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_AMOC, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
| 73 | ; |
---|
| 74 | cnt = cnt+1 & blabla = 'Sea-Ice cover' |
---|
| 75 | IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
[4320] | 76 | ; |
---|
| 77 | cnt = cnt+1 & blabla = 'Sea-Ice Volume' |
---|
| 78 | IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
[2709] | 79 | |
---|
[4320] | 80 | ;SF cnt = cnt+1 & blabla = 'SNOW Volume' |
---|
| 81 | ;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SNOW_Vol, masknp, POSTSCRIPT = postscript, _extra = ex |
---|
| 82 | ; |
---|
[2751] | 83 | IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_ts_html_body.txt', htmltxt[1:*] |
---|
[2709] | 84 | |
---|
| 85 | return |
---|
| 86 | END |
---|