Changeset 4755 for trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
- Timestamp:
- 2014-09-02T15:55:01+02:00 (10 years ago)
- Location:
- trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
- Files:
-
- 2 added
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/README
r4394 r4755 24 24 25 25 26 27 26 #---------------------------------------------------------------------- 27 # HOW TO USE 28 #---------------------------------------------------------------------- 29 ./std_main.sh -plot -pdf 30 or 31 ./std_main.sh -ts -pdf 28 32 29 33 -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro
r4559 r4755 21 21 22 22 IF keyword_set(arc) THEN BEGIN 23 domdef, 20, 380, 50, 9023 domdef, 20, 380, 40, 90 24 24 map = [90, 0, 0] 25 25 ENDIF 26 26 IF keyword_set(ant) THEN BEGIN 27 domdef, 20, 380, -90, - 5027 domdef, 20, 380, -90, -40 28 28 map = [-90, 0, 0] 29 29 ENDIF -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
r4559 r4755 172 172 Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 173 173 jan = where(mm EQ 1, cnt) 174 Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 174 175 sept = where(mm EQ 9, cnt) 175 Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}176 176 Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 177 177 undefine, Ithi_1 … … 288 288 Isal_sept_2 = {arr:-1} 289 289 ENDELSE 290 ; 290 291 291 IvelU_1 = read_ncdf(getenv('VAR1_IvelU'), allrecords = allrec, filename = std_file1_I ) 292 292 IvelV_1 = read_ncdf(getenv('VAR1_IvelV'), allrecords = allrec, filename = std_file1_I ) … … 335 335 Ivelo_sept_2 = {arr:-1} 336 336 ENDELSE 337 ; 337 338 338 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 339 339 caldat, time, mm … … 396 396 cnt = cnt+1 & blabla = 'Erp salinity damping term' 397 397 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 398 ; 398 ;; 399 399 cnt = cnt+1 & blabla = 'Evaporation - Precipitation - Runoff term' 400 400 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex … … 439 439 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD2, MLD, _extra = ex 440 440 ; 441 cnt = cnt+1 & blabla = 'Mixed layer depth ortho plan ARCTIC' 442 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_MLD_ortho, MLD1, MLD2, /ARC, _extra = ex 443 ;sf 444 cnt = cnt+1 & blabla = 'Mixed layer depth ortho plan ANTARTIC' 445 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_MLD_ortho, MLD1, MLD2, /ANT, _extra = ex 446 ;sf 441 447 cnt = cnt+1 & blabla = 'Zonal mean Mixed layer depth' 442 448 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ZonMld, MLD1, MLD2, MLD, _extra = ex … … 564 570 cnt = cnt+1 & blabla = 'Antartic Ice Velocity: SEPT' 565 571 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_sept_1, IvelU_sept_2, IvelV_sept_1, IvelV_sept_2, Ivelo_sept_1, Ivelo_sept_2, /ANT, /SEPT, _extra = ex 566 ; 572 567 573 cnt = cnt+1 & blabla = 'Arctic Ice Age: MARCH' 568 574 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh
r4562 r4755 21 21 #idl_command=idl71 22 22 # 23 PS_DIR=$( pwd )/ OR2LPv36_ps_plot24 PDF_DIR=$( pwd )/ OR2LPv36_pdf_plot23 PS_DIR=$( pwd )/p4H25a50-testht_ps_plot 24 PDF_DIR=$( pwd )/p4H25a50-testht_pdf_plot 25 25 HTML_DIR=$( pwd )/html_plot 26 26 SAXO_DIR=/Users/sflod/SAXO_DIR 27 27 # 28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/ 500y_pisces# path of data in NetCDF format28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 29 29 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data 30 30 DIR_MASK=/Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins) … … 54 54 # 55 55 #===================== EXP1 ===================== 56 # 23610101_2370123157 FILE1_T= OR2LPv36_23610101_23701231_1Y_grid_T.nc # exp1 grid T input file58 FILE1_U= OR2LPv36_23610101_23701231_1Y_grid_U.nc # exp1 grid U input file59 FILE1_V= OR2LPv36_23610101_23701231_1Y_grid_V.nc # exp1 grid V input file60 FILE1_I= OR2LPv36_23610101_23701231_1M_icemod.nc # exp1 ice input file56 # 57 FILE1_T=p4H25a50_20410101_20501231_1Y_grid_T.nc # exp1 grid T input file 58 FILE1_U=p4H25a50_20410101_20501231_1Y_grid_U.nc # exp1 grid U input file 59 FILE1_V=p4H25a50_20410101_20501231_1Y_grid_V.nc # exp1 grid V input file 60 FILE1_I=p4H25a50_20410101_20501231_1M_icemod.nc # exp1 ice input file 61 61 VAR1_T=thetao 62 62 VAR1_S=so 63 VAR1_QNET= tohfls64 VAR1_ERP= sowafldp63 VAR1_QNET=qt 64 VAR1_ERP=wfcorr 65 65 VAR1_EMP=wfo 66 66 VAR1_MLD=mldr10_1 67 67 VAR1_U=uocetr_eff 68 68 VAR1_V=vocetr_eff 69 VAR1_Ithick= iicethic70 VAR1_Ifrac=s oicecov71 VAR1_Isnow= isnowthi72 #VAR1_Isal=sisali73 #VAR1_Iage=siages74 VAR1_IvelU= iicevelu75 VAR1_IvelV= iicevelv76 #VAR1_Ivelo=sivelo69 VAR1_Ithick=sithic 70 VAR1_Ifrac=siconc 71 VAR1_Isnow=snvolu 72 VAR1_Isal=sisali 73 VAR1_Iage=siages 74 VAR1_IvelU=sivelu 75 VAR1_IvelV=sivelv 76 VAR1_Ivelo=sivelo 77 77 # 78 78 #===================== EXP2 ===================== 79 79 # 80 FILE2_T= OR2LPv36_23610101_23701231_1Y_grid_T.nc # exp1 grid T input file81 FILE2_U= OR2LPv36_23610101_23701231_1Y_grid_U.nc # exp1 grid U input file82 FILE2_V= OR2LPv36_23610101_23701231_1Y_grid_V.nc # exp1 grid V input file83 FILE2_I= OR2LPv36_23610101_23701231_1M_icemod.nc # exp1 ice input file80 FILE2_T=testht_20410101_20501231_1Y_grid_T.nc # exp1 grid T input file 81 FILE2_U=testht_20410101_20501231_1Y_grid_U.nc # exp1 grid U input file 82 FILE2_V=testht_20410101_20501231_1Y_grid_V.nc # exp1 grid V input file 83 FILE2_I=testht_20410101_20501231_1M_icemod.nc # exp1 ice input file 84 84 VAR2_T=thetao 85 85 VAR2_S=so 86 VAR2_QNET= tohfls87 VAR2_ERP= sowafldp86 VAR2_QNET=qt 87 VAR2_ERP=wfcorr 88 88 VAR2_EMP=wfo 89 89 VAR2_MLD=mldr10_1 90 90 VAR2_U=uocetr_eff 91 91 VAR2_V=vocetr_eff 92 VAR2_Ithick= iicethic93 VAR2_Ifrac=s oicecov94 VAR2_Isnow= isnowthi95 #VAR2_Isal=sisali96 #VAR2_Iage=siages97 VAR2_IvelU= iicevelu98 VAR2_IvelV= iicevelv99 #VAR2_Ivelo=sivelo92 VAR2_Ithick=sithic 93 VAR2_Ifrac=siconc 94 VAR2_Isnow=snvolu 95 VAR2_Isal=sisali 96 VAR2_Iage=siages 97 VAR2_IvelU=sivelu 98 VAR2_IvelV=sivelv 99 VAR2_Ivelo=sivelo 100 100 # 101 101 ######################### Export Variables ############################### -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names
r4559 r4755 74 74 VAR1_IvelU=iicevelu 75 75 VAR1_IvelV=iicevelv 76 VAR1_Ivelo=iicevelo76 #VAR1_Ivelo=iicevelo 77 77 # 78 78 #===================== EXP2 ===================== … … 99 99 VAR2_IvelU=iicevelu 100 100 VAR2_IvelV=iicevelv 101 VAR2_Ivelo=iicevelo101 #VAR2_Ivelo=iicevelo 102 102 # 103 103 ######################### Export Variables ############################### -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro
r4497 r4755 88 88 ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 89 89 ; 90 title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 90 ;;title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 91 title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 91 92 jpt=12 92 93 time=julday(1,15,1900)+30*lindgen(12) 93 pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $94 pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, LINESTYLE=2, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 94 95 , COLOR = 000 , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 95 pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT , LINESTYLE=2$ ;;; dashed lines is LINESTYLE=2 $96 pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT $ ;;; dashed lines is LINESTYLE=2 $ 96 97 , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 97 pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT $98 pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ 98 99 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 99 pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT , LINESTYLE=2$ ;;; dashed lines is LINESTYLE=2 $100 pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT $ ;;; dashed lines is LINESTYLE=2 $ 100 101 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 101 102 ; 102 title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area ( CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)'103 pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $103 title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 104 pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, LINESTYLE=2, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 104 105 ,COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 105 pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT , LINESTYLE=2$106 pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT $ 106 107 , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 107 pltt, vICE_area_SH, 't', /REMPLI, /PORTRAIT $108 pltt, vICE_area_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ 108 109 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 109 pltt, vICE_ext_SH, 't', /REMPLI, /PORTRAIT , LINESTYLE=2$110 pltt, vICE_ext_SH, 't', /REMPLI, /PORTRAIT $ 110 111 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 111 112 ; … … 159 160 160 161 161 title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 162 ;;title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 163 title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 162 164 jpt=12 163 165 time=julday(1,15,1900)+30*lindgen(12) 164 pltt, ICE_N, 't', MIN = 4, MAX = 16, /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $166 pltt, ICE_N, 't', MIN = 4, MAX = 16, /REMPLI, /PORTRAIT, LINESTYLE=2, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 165 167 , COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; BLACK 166 pltt, ICE_N2, 't', /REMPLI, /PORTRAIT $168 pltt, ICE_N2, 't', /REMPLI, /PORTRAIT , LINESTYLE=2 $ 167 169 , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; RED 168 pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT , LINESTYLE=2$ ; linee tratteggiate LINESTYLE=2 $170 pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $ 169 171 , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 170 pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT , LINESTYLE=2$ ; linee tratteggiate LINESTYLE=2 $172 pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $ 171 173 , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 172 pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT $174 pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ 173 175 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; light blue 174 pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT , LINESTYLE=2$176 pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT $ 175 177 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; blu scuro 176 178 ; 177 title ='Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+'Global Annual Mean Ice Area ( CONTINUS)'+'!C'+ 'and Extend minus 15% (DASHED)'179 title ='Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+'Global Annual Mean Ice Area (DASHED)'+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 178 180 ; title ='Southern Hemisphere'+'!C' 179 pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $181 pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, LINESTYLE=2, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 180 182 , COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 181 pltt, ICE_S2, 't', /REMPLI, /NOERASE $183 pltt, ICE_S2, 't', /REMPLI, /NOERASE, LINESTYLE=2 $ 182 184 , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 183 pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT , LINESTYLE=2$ ; linee tratteggiate LINESTYLE=2 $185 pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $ 184 186 , /ov1d, COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 185 pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT , LINESTYLE=2$ ; linee tratteggiate LINESTYLE=2 $187 pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $ 186 188 , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 187 pltt, vICE_area_SH, 't', /REMPLI, /PORTRAIT $189 pltt, vICE_area_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ 188 190 , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 189 pltt, vICE_ext_SH, 't', /REMPLI, /PORTRAIT , LINESTYLE=2$191 pltt, vICE_ext_SH, 't', /REMPLI, /PORTRAIT $ 190 192 , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 191 193 ; -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro
r4514 r4755 5 5 arr2d = arr2d[t1:t2] 6 6 arr2d = reform(arr2d,12,nyear) ; put in 2D array 7 arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 7 ;arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 8 arr2d = arr2d[*, nyear-1] ; select last year 8 9 9 10 return, arr2d … … 34 35 ; 35 36 d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 37 d2 = '('+strtrim(date2, 1)+')' 36 38 ; 37 39 iodir = std_iodir_data … … 74 76 ;; AREA 75 77 ICE_area_export = reform(ICE_area_export, 12, nyr) 76 ICE_area_export = total(ICE_area_export,2)/nyr 78 ;ICE_area_export = total(ICE_area_export,2)/nyr ; old version monthly mean average over all years 79 ICE_area_export = ICE_area_export[*,nyr-1] ; dim= 12, index 0-11; last year choosen 77 80 ; ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'} ; annual mean 78 81 ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'} ; monthly mean … … 80 83 ; 81 84 ICE_vol_export = reform(ICE_vol_export, 12, nyr) 82 ICE_vol_export = total(ICE_vol_export,2)/nyr 85 ;ICE_vol_export = total(ICE_vol_export,2)/nyr 86 ICE_vol_export = ICE_vol_export[*,nyr-1] 83 87 ; ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 365 , unit : '10^3 Km^3/year'} ; annual mean 84 88 ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'} ; monthly mean 85 89 86 90 ; 87 title = 'Fram Strait Areal Export'+'!C'+prefix+' '+d1_d2 91 ;title = 'Fram Strait Areal Export: LAST YEAR'+'!C'+prefix+' '+d1_d2 92 title = 'Fram Strait Areal Export: LAST YEAR'+'!C'+prefix+' '+d2 88 93 jpt=12 89 94 time=julday(1,15,1900)+30*lindgen(12) … … 99 104 xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 100 105 ; 101 title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 106 ;title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' '+d1_d2 107 title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' '+d2 102 108 pltt, ICE_vol_export, 't', /REMPLI, MIN = 60., MAX = 500. , /NOERASE, XGRIDSTYLE = 1 $ 103 109 , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex … … 118 124 119 125 d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 126 d2_2 = '('+strtrim(date2_2, 1)+')' 120 127 tsave = time 121 128 domdef, indx1, indx2, indy1, indy2, /xindex, /yindex,/memeindices … … 140 147 ;; AREA 141 148 ICE_area_export_2 = reform(ICE_area_export_2, 12, nyr) 142 ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 149 ;ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 150 ICE_area_export_2 = ICE_area_export_2[*,nyr-1] 143 151 ; ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'} ; annual mean 144 152 ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'} ; monthly mean 145 153 ; 146 154 ICE_vol_export_2 = reform(ICE_vol_export_2, 12, nyr) 147 ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 155 ;ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 156 ICE_vol_export_2 = ICE_vol_export_2[*,nyr-1] 148 157 ; ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'} ; annual mean 149 158 ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'} ; monthly mean … … 152 161 if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 153 162 154 title = 'Fram Strait Areal Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 163 ;title = 'Fram Strait Areal Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 164 title = 'Fram Strait Areal Export LAST YEAR'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d2_2 155 165 jpt=12 156 166 time=julday(1,15,1900)+30*lindgen(12) … … 170 180 xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 171 181 ; 172 title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 182 ;title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 183 title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d2_2 173 184 pltt, ICE_vol_export, 't', /REMPLI, MIN = 60., MAX = 500., /NOERASE, XGRIDSTYLE = 1 $ 174 185 , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro
r4504 r4755 70 70 ; 71 71 title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (BLACK)' 72 pltt, ICE_vol_S, 't', MIN = 0., MAX = 1 1000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $72 pltt, ICE_vol_S, 't', MIN = 0., MAX = 12000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 73 73 , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 74 74 ; … … 119 119 ; 120 120 title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Global Annual Mean Ice Volume' 121 pltt, ICE_vol_S, 't', MIN = 0., MAX = 1 1000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $121 pltt, ICE_vol_S, 't', MIN = 0., MAX = 12000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 122 122 , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 123 123 pltt, ICE_vol_S2, 't', /REMPLI, /NOERASE $ -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro
r4559 r4755 60 60 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SSH, masknp, POSTSCRIPT = postscript, _extra = ex 61 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 = ex62 ;;;; 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 64 ; 65 65 cnt = cnt+1 & blabla = 'Global Mean EMP' … … 77 77 cnt = cnt+1 & blabla = 'Sea-Ice Volume' 78 78 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 79 ; 80 cnt = cnt+1 & blabla = 'ICE Velocity '81 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex82 ; 79 ;SF; and probably useless 80 ;SF cnt = cnt+1 & blabla = 'ICE Velocity ' 81 ;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 82 ;SF; 83 83 cnt = cnt+1 & blabla = 'Sea-Ice Volume Export at Fram Strait' 84 84 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_FRAM, masknp, POSTSCRIPT = postscript, _extra = ex -
trunk/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh
r4559 r4755 18 18 idl_command=/Applications/itt/idl64/bin/idl 19 19 # 20 PS_DIR=$( pwd )/ts_ heatbug_ps21 PDF_DIR=$( pwd )/ts_ heatbug_pdf20 PS_DIR=$( pwd )/ts_p4H25a50-testht_ps 21 PDF_DIR=$( pwd )/ts_p4H25a50-testht_pdf 22 22 HTML_DIR=$( pwd )/html 23 23 SAXO_DIR=/Users/sflod/SAXO_DIR 24 24 # 25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada 25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 26 26 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data 27 27 DIR_MASK=/Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins) … … 52 52 READ_ONLY_FIRST_RECORD=0 # if 0 then read all records in files else read only the first reacord in each file 53 53 # 54 #DATE1=20010101 ; DATE2=25001231 55 DATE1=20310101 ; DATE2=20401231 56 #DATE1=00010101 ; DATE2=00101231 54 DATE1=20010101 ; DATE2=20501231 57 55 # 58 VAR1_T=thetao ; V1T_PREF=heatbug; V1T_SUFF=_1Y_grid_T.nc59 VAR1_S=so ; V1S_PREF=heatbug; V1S_SUFF=_1Y_grid_T.nc60 VAR1_SSH=zos ; V1SSH_PREF= heatbug; V1SSH_SUFF=_1Y_grid_T.nc61 VAR1_Q=tohfls ; V1Q_PREF=heatbug; V1Q_SUFF=_1Y_grid_T.nc62 VAR1_EMP=wfo ; V1EMP_PREF= heatbug; V1EMP_SUFF=_1Y_grid_T.nc63 VAR1_U=uocetr_eff ; V1U_PREF=heatbug; V1U_SUFF=_1Y_grid_U.nc64 VAR1_V=vocetr_eff ; V1V_PREF=heatbug; V1V_SUFF=_1Y_grid_V.nc65 VAR1_ICE=siconc ; V1ICE_PREF= heatbug; V1ICE_SUFF=_1M_icemod.nc66 VAR1_Ithick=sithic ; V1It_PREF=heatbug; V1It_SUFF=_1M_icemod.nc67 VAR1_SNOW=sn dept ; V1SNOW_PREF=heatbug; V1SNOW_SUFF=_1M_icemod.nc68 VAR1_IvelV=sivelv ; V1IvV_PREF= heatbug ; V1IvV_SUFF=_1M_icemod.nc69 VAR1_Ivel=sivelo ; V1Iv_PREF=heatbug; V1Iv_SUFF=_1M_icemod.nc56 VAR1_T=thetao ; V1T_PREF=p4H25a50 ; V1T_SUFF=_1Y_grid_T.nc 57 VAR1_S=so ; V1S_PREF=p4H25a50 ; V1S_SUFF=_1Y_grid_T.nc 58 VAR1_SSH=zos ; V1SSH_PREF=p4H25a50 ; V1SSH_SUFF=_1Y_grid_T.nc 59 ##VAR1_Q=qt ; V1Q_PREF=p4H25a50 ; V1Q_SUFF=_1Y_grid_T.nc 60 VAR1_EMP=wfo ; V1EMP_PREF=p4H25a50 ; V1EMP_SUFF=_1Y_grid_T.nc 61 VAR1_U=uocetr_eff ; V1U_PREF=p4H25a50 ; V1U_SUFF=_1Y_grid_U.nc 62 VAR1_V=vocetr_eff ; V1V_PREF=p4H25a50 ; V1V_SUFF=_1Y_grid_V.nc 63 VAR1_ICE=siconc ; V1ICE_PREF=p4H25a50 ; V1ICE_SUFF=_1M_icemod.nc 64 VAR1_Ithick=sithic ; V1It_PREF=p4H25a50 ; V1It_SUFF=_1M_icemod.nc 65 VAR1_SNOW=snthic ; V1SNOW_PREF=p4H25a50 ; V1SNOW_SUFF=_1M_icemod.nc 66 VAR1_IvelV=sivelv ; V1IvV_PREF=p4H25a50 ; V1IvV_SUFF=_1M_icemod.nc 67 VAR1_Ivel=sivelo ; V1Iv_PREF=p4H25a50 ; V1Iv_SUFF=_1M_icemod.nc 70 68 # 71 69 #===================== EXP2 ===================== 72 70 # 73 #DATE1_2=20010101 ; DATE2_2=20101231 74 DATE1_2=20310101 ; DATE2_2=20401231 75 #DATE1_2=20010101 ; DATE2_2=25001231 71 DATE1_2=20010101 ; DATE2_2=20501231 76 72 # 77 VAR2_T=thetao ; V2T_PREF=heatbug; V2T_SUFF=_1Y_grid_T.nc78 VAR2_S=so ; V2S_PREF=heatbug; V2S_SUFF=_1Y_grid_T.nc79 VAR2_SSH=zos ; V2SSH_PREF=heatbug; V2SSH_SUFF=_1Y_grid_T.nc80 VAR2_Q=tohfls ; V2Q_PREF=heatbug; V2Q_SUFF=_1Y_grid_T.nc81 VAR2_EMP=wfo ; V2EMP_PREF=heatbug; V2EMP_SUFF=_1Y_grid_T.nc82 VAR2_U=uocetr_eff ; V2U_PREF=heatbug; V2U_SUFF=_1Y_grid_U.nc83 VAR2_V=vocetr_eff ; V2V_PREF=heatbug; V2V_SUFF=_1Y_grid_V.nc84 VAR2_ICE=siconc ; V2ICE_PREF=heatbug; V2ICE_SUFF=_1M_icemod.nc85 VAR2_Ithick=sithic ; V2It_PREF=heatbug; V2It_SUFF=_1M_icemod.nc86 VAR2_SNOW=sn dept ; V2SNOW_PREF=heatbug; V2SNOW_SUFF=_1M_icemod.nc87 VAR2_IvelV=sivelv ; V2IvV_PREF=heatbug ; V2IvV_SUFF=_1M_icemod.nc88 VAR2_Ivel=sivelo ; V2Iv_PREF=heatbug ; V2Iv_SUFF=_1M_icemod.nc73 VAR2_T=thetao ; V2T_PREF=testht ; V2T_SUFF=_1Y_grid_T.nc 74 VAR2_S=so ; V2S_PREF=testht ; V2S_SUFF=_1Y_grid_T.nc 75 VAR2_SSH=zos ; V2SSH_PREF=testht ; V2SSH_SUFF=_1Y_grid_T.nc 76 ##VAR2_Q=qt ; V2Q_PREF=testht ; V2Q_SUFF=_1Y_grid_T.nc 77 VAR2_EMP=wfo ; V2EMP_PREF=testht ; V2EMP_SUFF=_1Y_grid_T.nc 78 VAR2_U=uocetr_eff ; V2U_PREF=testht ; V2U_SUFF=_1Y_grid_U.nc 79 VAR2_V=vocetr_eff ; V2V_PREF=testht ; V2V_SUFF=_1Y_grid_V.nc 80 VAR2_ICE=siconc ; V2ICE_PREF=testht ; V2ICE_SUFF=_1M_icemod.nc 81 VAR2_Ithick=sithic ; V2It_PREF=testht ; V2It_SUFF=_1M_icemod.nc 82 VAR2_SNOW=snthic ; V2SNOW_PREF=testht ; V2SNOW_SUFF=_1M_icemod.nc 83 VAR2_IvelV=sivelv ; V2IvV_PREF=testht ; V2IvV_SUFF=_1M_icemod.nc 84 VAR2_Ivel=sivelo ; V2IvV_PREF=testht ; V2IvV_SUFF=_1M_icemod.nc 89 85 # 90 86 ######################### Export Variables ###############################
Note: See TracChangeset
for help on using the changeset viewer.