Changeset 3999 for branches/2013
- Timestamp:
- 2013-08-06T13:39:29+02:00 (11 years ago)
- Location:
- branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
- Files:
-
- 6 added
- 21 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_ArcSal.pro
r2751 r3999 26 26 if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 27 27 ; 28 domdef, 20, 380, 60, 9028 domdef, 20, 380, 50, 90 29 29 ; 30 30 varunit = S1.unit -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_GlobMeanTS.pro
r2751 r3999 25 25 title = 'Salinity Global mean!C'+std_file1_T+' - Levitus (Black)' 26 26 if std_file1_T NE std_file2_T THEN title = title+'!C'+std_file2_T+' - Levitus (Red)' 27 plt1d, S1.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -.2, MAX = . 2, boxzoom = [4., 5300.], /KEEPBOTTOM $27 plt1d, S1.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -.2, MAX = .4, boxzoom = [4., 5300.], /KEEPBOTTOM $ 28 28 , small = [1, 2, 2], XGRIDSTYLE = 2, YGRIDSTYLE = 2, TITLE = title, /NOERASE, _extra = ex 29 29 if std_file1_T NE std_file2_T then begin 30 plt1d, S2.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -2., MAX = 2., boxzoom = [4., 5300.], /KEEPBOTTOM $30 plt1d, S2.arr - SLev.arr, typein = 'z', ticklen = 1, MIN = -2., MAX = 4., boxzoom = [4., 5300.], /KEEPBOTTOM $ 31 31 , /ov1d, COLOR = 250, TITLE = title, _extra = ex 32 32 endif -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro
r3927 r3999 40 40 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 41 41 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 42 endif else begin 43 title = title+ ' - Observations' 44 plt, Ifra1.arr - Ifra_obs.arr, MIN = -1., MAX = 1., INTER = 0.1, STYLE = 'so0so', format = '(f4.1)' $ 45 , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, CELL_FILL = 2, TITLE = title $ 46 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 47 endelse 42 ;;SF finire, aggiungere lettura obs 43 ;;SF endif else begin 44 ;;SF title = title+ ' - Observations' 45 ;;SF plt, Ifra1.arr - Ifra_obs.arr, MIN = -1., MAX = 1., INTER = 0.1, STYLE = 'so0so', format = '(f4.1)' $ 46 ;;SF , small = [1, 2, 2], COAST_THICK = 2, TITLE = title, CELL_FILL = 2, TITLE = title $ 47 ;;SF , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 48 ;;SF endelse 49 endif 48 50 49 51 domdef -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceThick.pro
r3927 r3999 22 22 domdef, 20, 380, 50, 90 23 23 vmin = 0. 24 vmax = 10.24 vmax = 6. 25 25 vint = 0.5 26 26 fmt = '(i2)' -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro
r3927 r3999 55 55 if std_file1_I NE std_file2_I then begin 56 56 title = title + std_file2_I 57 plt, Isnow1.arr - I thi2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $57 plt, Isnow1.arr - Isnow2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $ 58 58 , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 59 59 , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
r3927 r3999 199 199 ENDELSE 200 200 ; 201 Isnow_1 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file1_I ) 202 help, Isnow_1.arr, jpi,jpj 203 caldat, time, mm 204 april = where(mm EQ 4, cnt) 205 Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 206 sept = where(mm EQ 9, cnt) 207 Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 208 ;undefine, Isnow_1 209 ; 210 IF std_file2_I NE std_file1_I THEN BEGIN 211 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 212 caldat, time, mm 213 april = where(mm EQ 4, cnt) 214 Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 215 sept = where(mm EQ 9, cnt) 216 Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 217 ; undefine, Isnow_2 218 ENDIF ELSE BEGIN 219 Isnow_april_2 = {arr:-1} 220 Isnow_sept_2 = {arr:-1} 221 ENDELSE 201 ;; Sf commentato perché non funziona snow 202 ;; Sf Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 203 ;; Sf help, Isnow_1.arr, jpi,jpj 204 ;; Sf caldat, time, mm 205 ;; Sf april = where(mm EQ 4, cnt) 206 ;; Sf Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 207 ;; Sf sept = where(mm EQ 9, cnt) 208 ;; Sf Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 209 ;; Sf ;undefine, Isnow_1 210 ;; Sf ; 211 ;; Sf IF std_file2_I NE std_file1_I THEN BEGIN 212 ;; Sf Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 213 ;; Sf caldat, time, mm 214 ;; Sf april = where(mm EQ 4, cnt) 215 ;; Sf Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 216 ;; Sf sept = where(mm EQ 9, cnt) 217 ;; Sf Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 218 ;; Sf ; undefine, Isnow_2 219 ;; Sf ENDIF ELSE BEGIN 220 ;; Sf Isnow_april_2 = {arr:-1} 221 ;; Sf Isnow_sept_2 = {arr:-1} 222 ;; Sf ENDELSE 222 223 ; 223 224 jpt = 1 … … 325 326 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex 326 327 327 ; ;cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH'328 ; ;IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex328 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH' 329 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 329 330 330 331 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' … … 334 335 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 335 336 336 ; ;cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH'337 ; ;IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex337 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 338 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 338 339 339 340 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' … … 354 355 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' 355 356 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 356 357 358 359 360 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT'361 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex362 363 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL'364 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex365 366 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT'367 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex368 369 cnt = cnt+1 & blabla = 'Arctic Ice Extend: '370 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC, _extra = ex371 372 cnt = cnt+1 & blabla = 'Antarctic Ice Extend: '373 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT, _extra = ex357 ;; SF non funziona, dice che deve essere una struttura e non e!!!! 358 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: APRIL' 359 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 360 ;; SF 361 ;; SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 362 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 363 ;; SF 364 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL' 365 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 366 ;; SF 367 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 368 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 369 ; 370 ; cnt = cnt+1 & blabla = 'Arctic Ice Extend: ' 371 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC, _extra = ex 372 373 ; cnt = cnt+1 & blabla = 'Antarctic Ice Extend: ' 374 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT, _extra = ex 374 375 375 376 cnt = cnt+1 & blabla = 'Meridional stream Function: Global (no Med)' -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_bsf.pro
r2751 r3999 14 14 domdef, 0, 6000 15 15 16 bb1 = bsf(U1.arr, refvalue = 0., refpoint = [25, 0] )16 bb1 = bsf(U1.arr, refvalue = 0., refpoint = [25, 0], /transport ) 17 17 title = 'Barotropic Stream Function!C'+std_file1_U 18 18 plt, bb1, min = -200., max = 200., int = 10., /portrait, FORMAT = '(I4)', STYLE = 'so0so' $ … … 20 20 21 21 if std_file1_U NE std_file2_U then BEGIN 22 bb2 = bsf(U2.arr, refvalue = 0., refpoint = [25, 0] )22 bb2 = bsf(U2.arr, refvalue = 0., refpoint = [25, 0], /transport ) 23 23 title = 'Barotropic Stream Function!C'+std_file1_U+' - '+std_file2_U 24 24 plt, bb1.arr - bb2.arr, min = -20., max = 20., int = 2., FORMAT = '(I2)', STYLE = 'so0so' $ -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_mld.pro
r2751 r3999 32 32 ENDIF ELSE BEGIN 33 33 title = titleorg+std_file1_T 34 plt, MLD1, MIN = 0., MAX = 500., INTER = 25., /NOCONTOUR, FORMAT = '(I3)' $ 35 , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 34 plt, MLD1,label=4, cb_label=[0, 15, 25, 50, 75, 100, 125, 150, 250, 400, 600, 850],/NOCONTOUR,FORMAT = '(I3)' $ 35 , small = [1, 2, 1],COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 36 ;ORI SF plt, MLD1, MIN = 0., MAX = 500., INTER = 25., /NOCONTOUR, FORMAT = '(I3)' $ 37 ;ORI SF , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, /PORTRAIT, _extra = ex 36 38 ENDELSE 37 39 ; -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_msf.pro
r2760 r3999 33 33 domdef, 0, 6000 34 34 35 mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma )35 mm1 = msf(V1.arr, msk, indexboxzoom = ind, maskout = ma, /transport ) 36 36 title = 'Meridional Stream Function, '+subname+'!C'+std_file1_V 37 37 pltz, mm1, 'yz', -20., 20., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)', /portrait $ … … 39 39 40 40 if std_file1_V NE std_file2_V then begin 41 mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma )41 mm2 = msf(V2.arr, msk, indexboxzoom = ind, maskout = ma, /transport) 42 42 title = title+' - '+std_file2_V 43 43 pltz, mm1.arr-mm2.arr, 'yz', -10., 10., int = 1., boxzoom = [ind[0:1], lat_ext, 0, 5500], /xindex, FORMAT = '(I3)' $ -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh
r3927 r3999 16 16 #===================== User PATHS ===================== 17 17 # 18 idl_command=/Applications/itt/idl64/bin/idl 18 #idl_command=/Applications/itt/idl64/bin/idl 19 idl_command=/Applications/itt/idl71/bin/idl 20 #idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl 21 #idl_command=idl71 19 22 # 20 PS_DIR=$( pwd )/ ps_plot21 PDF_DIR=$( pwd )/ pdf_plot23 PS_DIR=$( pwd )/500yfwb0_ps_plot 24 PDF_DIR=$( pwd )/500yfwb0_pdf_plot 22 25 HTML_DIR=$( pwd )/html_plot 23 26 SAXO_DIR=/Users/sflod/SAXO_DIR 24 27 # 25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada 28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 26 29 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data 27 30 DIR_MASK=/Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins) … … 34 37 #===================== DATA ===================== 35 38 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc # Levitus_T 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=data_1y_salinity_nomask.nc # Levitus_S 39 #VAR_TEMP_3D=votemper ; FILE_TEMP_3D=data_1y_potential_temperature_nomask.nc # Levitus_T 40 #VAR_SAL_3D=vosaline ; FILE_SAL_3D=data_1y_salinity_nomask.nc # Levitus_S 41 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 42 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 38 43 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 44 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux … … 52 57 #===================== EXP1 ===================== 53 58 # 54 FILE1_T=LIM2_00010101_00011231_1m_grid_T.nc # exp1 grid T input file 55 FILE1_U=LIM2_00010101_00011231_1m_grid_U.nc # exp1 grid U input file 56 FILE1_V=LIM2_00010101_00011231_1m_grid_V.nc # exp1 grid V input file 57 FILE1_I=LIM2_00010101_00011231_1m_icemod.nc # exp1 ice input file 58 VAR1_T=votemper 59 VAR1_S=vosaline 60 VAR1_QNET=sohefldo 61 VAR1_ERP=sowaflcd 62 VAR1_EMP=sowaflup 59 #FILE1_T=l3fwb178_24910101_25001231_1Y_grid_T.nc # exp1 grid T input file 60 #FILE1_U=l3fwb178_24910101_25001231_1Y_grid_U.nc # exp1 grid U input file 61 #FILE1_V=l3fwb178_24910101_25001231_1Y_grid_V.nc # exp1 grid V input file 62 #FILE1_I=l3fwb178_24910101_25001231_1M_icemod.nc # exp1 ice input file 63 FILE1_T=500yfwb0_20010101_21001231_1Y_grid_T.nc # exp1 grid T input file 64 FILE1_U=500yfwb0_20010101_21001231_1Y_grid_U.nc # exp1 grid U input file 65 FILE1_V=500yfwb0_20010101_21001231_1Y_grid_V.nc # exp1 grid V input file 66 FILE1_I=500yfwb0_20010101_21001231_1M_icemod.nc # exp1 ice input file 67 VAR1_T=votemper 68 VAR1_S=vosaline 69 VAR1_QNET=sohflxdo 70 VAR1_ERP=sowafldp 71 VAR1_EMP=sowaflup 63 72 VAR1_MLD=somxl010 64 VAR1_U=uocetr_eff 65 VAR1_V=vocetr_eff 66 VAR1_Ithick=iicethic 73 VAR1_U=vozoeftr 74 #VAR1_U=uoce 75 #VAR1_U=vozocrtx 76 VAR1_V=vomeeftr 77 #VAR1_V=voce 78 #VAR1_V=vomecrty 79 VAR1_Ithick=iicethic 67 80 VAR1_Ifrac=iiceconc 68 VAR1_ Isnow=isnowthi81 VAR1_SNOW=isnowthi 69 82 # 70 83 #===================== EXP2 ===================== 71 84 # 72 FILE2_T=LIM2_00010101_00011231_1m_grid_T.nc # exp1 grid T input file 73 FILE2_U=LIM2_00010101_00011231_1m_grid_U.nc # exp1 grid U input file 74 FILE2_V=LIM2_00010101_00011231_1m_grid_V.nc # exp1 grid V input file 75 FILE2_I=LIM2_00010101_00011231_1m_icemod.nc 85 #FILE2_T=l3fwb178_24910101_25001231_1Y_grid_T.nc # exp1 grid T input file 86 #FILE2_U=l3fwb178_24910101_25001231_1Y_grid_U.nc # exp1 grid U input file 87 #FILE2_V=l3fwb178_24910101_25001231_1Y_grid_V.nc # exp1 grid V input file 88 #FILE2_I=l3fwb178_24910101_25001231_1M_icemod.nc # exp1 ice input file 89 FILE2_T=500yfwb0_20010101_21001231_1Y_grid_T.nc # exp1 grid T input file 90 FILE2_U=500yfwb0_20010101_21001231_1Y_grid_U.nc # exp1 grid U input file 91 FILE2_V=500yfwb0_20010101_21001231_1Y_grid_V.nc # exp1 grid V input file 92 FILE2_I=500yfwb0_20010101_21001231_1M_icemod.nc # exp1 ice input file 76 93 VAR2_T=votemper 77 94 VAR2_S=vosaline 78 VAR2_QNET=soh efldo79 VAR2_ERP=sowafl cd95 VAR2_QNET=sohflxdo 96 VAR2_ERP=sowafldp 80 97 VAR2_EMP=sowaflup 81 98 VAR2_MLD=somxl010 82 VAR2_U=uocetr_eff 83 VAR2_V=vocetr_eff 99 VAR2_U=vozoeftr 100 #VAR2_U=uoce 101 #VAR2_U=vozocrtx 102 VAR2_V=vomeeftr 103 #VAR2_V=voce 104 #VAR2_V=vomecrty 84 105 VAR2_Ithick=iicethic 85 106 VAR2_Ifrac=iiceconc 86 VAR2_ Isnow=isnowthi107 VAR2_SNOW=isnowthi 87 108 # 88 109 ######################### Export Variables ############################### … … 100 121 export FILE_FLUX VAR_FLUX 101 122 export FILE_MLD VAR_MLD 102 export FILE_ICE VAR_ICE_EXT_NH VAR_ICE_EXT_SH , VAR_ICE_area_NH,VAR_ICE_area_SH123 export FILE_ICE VAR_ICE_EXT_NH VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH 103 124 export FILE_SNOW_ARC VAR_SNOW_NH 104 125 export FILE_SNOW_ANT VAR_SNOW_SH -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_EMP.pro
r2751 r3999 27 27 ts_EMP.arr = ts_EMP.arr * ( 1.E-09 * surf_oce ) & ts_EMP.unit = 'Sv' 28 28 title = prefix+' '+d1_d2+'!C'+blabla 29 pltt, ts_EMP, 't', - .001, .001, date1, date2, /REMPLI $29 pltt, ts_EMP, 't', -2., 2., date1, date2, /REMPLI $ 30 30 , small = [1, 2, 1], YTITLE = 'Sv', TITLE = title, /PORTRAIT, _extra = ex 31 31 … … 39 39 40 40 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 41 pltt, ts_EMP.arr - ts_EMP2.arr, 't', - .001, .001, date1, date2, /REMPLI $41 pltt, ts_EMP.arr - ts_EMP2.arr, 't', -2., 2., date1, date2, /REMPLI $ 42 42 , color = 250, small = [1, 2, 2], YTITLE = 'Sv', TITLE = title, /noerase, _extra = ex 43 43 -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro
r3927 r3999 88 88 ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 89 89 ; 90 title = 'Northern Hemisphere 90 title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+' Global Annual Mean Ice Area (Black SOLID simulation)'+'!C'+ 'and Extend minus 15% (Blue SOLID simulation)'+'!C'+'Observation (dashed)' 91 91 jpt=12 92 92 time=julday(1,15,1900)+30*lindgen(12) … … 127 127 ICE_N2 = {arr:ICE_N2 * 1.e-12, unit : '10^12 m^2'} 128 128 ;ICE Extent minus 15% in NORTH Hemisphere 129 ICE_N2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix _2, suffix_2, FIRSTONLY = 1 - allrec)129 ICE_N2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 130 130 msk = ICE_N2_15.arr gt 0.15 ; remove 0.15% for observations 131 131 ICE_N2_15 = grossemoyenne( msk, 'xy', /integration, mask2d = masknp) … … 137 137 ;ICE Extent in SOUTH Hemisphere 138 138 domdef, 0, jpi-1, -90, -30, /xindex 139 ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix _2, suffix_2, FIRSTONLY = 1 - allrec)139 ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 140 140 ICE_S2 = grossemoyenne(ICE_S2.arr, 'xy', /integration, mask2d = masknp) 141 141 if jpt mod 12 ne 0 then stop … … 145 145 ICE_S2 = {arr:ICE_S2 * 1.e-12, unit : '10^12 m^2'} 146 146 ;ICE Extent minus 15% in SOUTH Hemisphere 147 ICE_S2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix _2, suffix_2, FIRSTONLY = 1 - allrec)147 ICE_S2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 148 148 msk = ICE_S2_15.arr gt 0.15 ; remove 0.15% for observations 149 149 ICE_S2_15 = grossemoyenne(msk, 'xy', /integration, mask2d = masknp) … … 154 154 ICE_S2_15 = {arr:ICE_S2_15 * 1.e-12, unit : '10^12 m^2'} 155 155 ; 156 time = tsave & IF n_elements(time) NE jpt THEN stop156 ; time = tsave & IF n_elements(time) NE jpt THEN stop 157 157 158 if KEYWORD_SET(postscript) then openps, filename+' .ps', portrait = 1158 if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 159 159 160 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (Red SOLID simulation) and Extend minus 15% (Blue SOLID simulation)'161 pltt, ICE_N.arr - ICE_N2.arr, 't', -2., 2., date1, date2, /REMPLI, /PORTRAIT, window = 2 $162 , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex163 pltt, ICE_N_15.arr - ICE_N2_15.arr , 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $164 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex165 166 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.) (Red SOLID simulation) and Extend minus 15% (Blue SOLID simulation)'167 pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., date1, date2, /REMPLI, /NOERASE $168 , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex169 pltt, ICE_S_15.arr - ICE_S2_15.arr , 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $170 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex171 172 160 173 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] 161 title = 'Northern Hemisphere'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2+'!C'+' Global Annual Mean Ice Area (BLUE)'+'!C'+ 'and Extend minus 15% (RED)' 162 jpt=12 163 time=julday(1,15,1900)+30*lindgen(12) 164 pltt, ICE_N.arr - ICE_N2.arr, 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT, window = 2 $ 165 , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 166 pltt, ICE_N_15.arr - ICE_N2_15.arr , 't', -.2, .2, 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $ 167 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 168 ;SF pltt, vICE_area_NH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 169 ;SF , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 170 ;SF pltt, vICE_ext_NH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 171 ;SF , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 172 ; 173 ; title ='Southern Hemisphere'+'!C'+prefix+' - '+prefix2+' '+d1_d2_2+' - '+'!C'+'Global Annual Mean Ice Area (BLUE)'+'!C'+ 'and Extend minus 15% (RED)' 174 title ='Southern Hemisphere'+'!C' 175 pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., 19000101, 19001231, /REMPLI, /NOERASE $ 176 , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 177 pltt, ICE_S_15.arr - ICE_S2_15.arr , 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2 $ 178 , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 179 ;SF pltt, vICE_area_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 180 ;SF , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 181 ;SF pltt, vICE_ext_SH, 't', -2., 2., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ 182 ;SF , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 183 ; 184 185 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png /> ' ] 174 186 if KEYWORD_SET(postscript) then closeps 175 187 -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro
r3929 r3999 26 26 v2_Ithick = getenv('VAR2_Ithick') & prefix2 = getenv('V2It_PREF') & suffix2 = getenv('V2It_SUFF') 27 27 28 ; get ice VOLUME climatology info29 ;;30 ; std_file_ice = isafile(getenv('FILE_ICE'), title = 'ICE Extent Climatology', iodir = std_iodir_climato)31 ;32 ; time_ice = ncdf_lec( std_file_ice, VAR='time' )33 ; time_ice = (time_ice - FLOOR(time_ice) ) * 1234 ; time_ice = round(time_ice) ; round to nearest integer35 ; t1 = where(time_ice eq 0)36 ; t1 = t1[0] ; jannuary37 ; t2 = where(time_ice eq 11, count)38 ; t2 = t2[count-1] ; last day of december39 ; nyear = (t2-t1+1)/1240 ; vICE_Ithick_NH = read_arr2d(std_file_ice, getenv('VAR1_Ithick'), t1, t2 )41 ; vICE_Ithick_SH = read_arr2d(std_file_ice, getenv('VAR_ICE_EXT_SH'), t1, t2 )42 ;43 ; vICE_area_NH = read_arr2d(std_file_ice, getenv('VAR_ICE_area_NH'), t1, t2 )44 ; vICE_area_SH = read_arr2d(std_file_ice, getenv('VAR_ICE_area_SH'), t1, t2 )45 ;46 28 cdti3 = string(cnt, format = '(i3.3)') 47 29 print, cdti3 + ') ' + blabla 48 filename = cdti3 + '_ts_ICE_ '+prefix30 filename = cdti3 + '_ts_ICE_Vol'+prefix 49 31 if prefix NE prefix2 then filename = filename + '_'+prefix2 50 32 if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 33 ; 51 34 d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 35 ; 52 36 iodir = std_iodir_data 53 37 ; ICE Area(=Surface) in NORTH Hemisphere … … 58 42 ICE_vol_N = (ICE_N.arr < 1.e10 ) * ( ICE_thick.arr < 1.e10) ; limited mask value of 1.e20, because 1.e20 * 1.e20 = inf for idl 59 43 ICE_vol_N = grossemoyenne(ICE_vol_N, 'xy', /integration, mask2d = masknp) 44 ; 60 45 if jpt mod 12 ne 0 then stop 61 46 nyr=jpt/12. … … 63 48 ICE_vol_N = total(ICE_vol_N,2)/nyr 64 49 ICE_vol_N = {arr:ICE_vol_N * 1.e-9, unit : '10^9 Km^3'} 50 ; 65 51 ;ICE Area(=Surface) in SOUTH Hemisphere 66 52 domdef, 0, jpi-1, -90, -30, /xindex … … 76 62 ICE_vol_S = {arr:ICE_vol_S * 1.e-9, unit : '10^9 Km^3'} 77 63 ; 78 title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation) + Observation (dashed)'64 title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 79 65 jpt=12 80 66 time=julday(1,15,1900)+30*lindgen(12) 81 67 pltt, ICE_vol_N, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 0., MAX = 30000. $ 82 , small = [1, 2, 1], YTITLE = ' 10^9Km^3 ', TITLE = title, _extra = ex68 , small = [1, 2, 1], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 83 69 ; 84 title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation) + Observation (dashed)'70 title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 85 71 pltt, ICE_vol_S, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 11000. $ 86 , small = [1, 2, 2], YTITLE = ' 10^9Km^3 ', TITLE = title, _extra = ex72 , small = [1, 2, 2], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 87 73 ; 88 89 74 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png /> ' ] 90 75 if KEYWORD_SET(postscript) then closeps … … 96 81 domdef, 0, jpi-1, 30, 90, /xindex 97 82 ;ICE Area(=Surface) in NORTH Hemisphere 98 ICE_N2 = rseries_ncdf(vICE2, date1 , date2, prefix, suffix, FIRSTONLY = 1 - allrec)99 ICE_thick2 = rseries_ncdf(v2_Ithick, date1 , date2, prefix, suffix, FIRSTONLY = 1 - allrec)83 ICE_N2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 84 ICE_thick2 = rseries_ncdf(v2_Ithick, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 100 85 ; Volume = Area(=Surface) * Thickness 101 86 ICE_vol_N2 = (ICE_N2.arr < 1.e10 ) * ( ICE_thick2.arr < 1.e10) ; limited mask value of 1.e20, because 1.e20 * 1.e20 = inf for idl … … 109 94 ;ICE Area(=Surface) in SOUTH Hemisphere 110 95 domdef, 0, jpi-1, -90, -30, /xindex 111 ICE_S2 = rseries_ncdf(vICE2, date1 , date2, prefix, suffix, FIRSTONLY = 1 - allrec)112 ICE_thick2 = rseries_ncdf(v2_Ithick, date1 , date2, prefix, suffix, FIRSTONLY = 1 - allrec)96 ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 97 ICE_thick2 = rseries_ncdf(v2_Ithick, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 113 98 ; Volume = Area(=Surface) * Thickness 114 99 ICE_vol_S2 = (ICE_S2.arr < 1.e10 ) * ( ICE_thick.arr < 1.e10) ; limited mask value of 1.e20, because 1.e20 * 1.e20 = inf for idl … … 120 105 ICE_vol_S2 = {arr:ICE_vol_S2 * 1.e-9, unit : '10^3 Km^3'} 121 106 122 time = tsave & IF n_elements(time) NE jpt THEN stop107 ; time = tsave & IF n_elements(time) NE jpt THEN stop 123 108 124 if KEYWORD_SET(postscript) then openps, filename+' .ps', portrait = 1109 if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 125 110 126 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Volume (North. Hemisp.)' 127 pltt, ICE_vol_N.arr - ICE_vol_N2.arr, 't', -2., 2., date1, date2, /REMPLI, /PORTRAIT, window = 2 $ 128 , COLOR = 250, small = [1, 2, 1], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 129 ; 130 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)' 131 pltt, ICE_vol_S.arr - ICE_vol_S2.arr, 't', -2., 2., date1, date2, /REMPLI, /NOERASE $ 132 , COLOR = 250, small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 111 title = 'Northern Hemisphere'+'!C'+'Global Annual Mean Ice Volume'+prefix+' - '+prefix2+' '+d1_d2_2 112 jpt=12 113 time=julday(1,15,1900)+30*lindgen(12) 114 pltt, ICE_vol_N.arr - ICE_vol_N2.arr, 't', -125., -60., 19000101, 19001231, /REMPLI, /PORTRAIT, window = 2 $ 115 , small = [1, 2, 1], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 116 ; 117 title = 'Southern Hemisphere'+'!C'+'Global Annual Mean Ice Volume'+prefix+' - '+prefix2+' '+d1_d2_2 118 pltt, ICE_vol_S.arr - ICE_vol_S2.arr, 't', -620., -10., 19000101, 19001231, /REMPLI, /NOERASE $ 119 , small = [1, 2, 2], YTITLE = 'Km^3 ', TITLE = title, _extra = ex 133 120 134 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+' .png /> ' ]121 htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png /> ' ] 135 122 if KEYWORD_SET(postscript) then closeps 136 123 -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_SSH.pro
r2751 r3999 24 24 ts_SSH.arr = ts_SSH.arr * 100 & ts_SSH.unit = 'cm' 25 25 title = prefix+' '+d1_d2+'!C'+blabla 26 pltt, ts_SSH, 't', - .1, .1, date1, date2, /REMPLI $26 pltt, ts_SSH, 't', -220., 220., date1, date2, /REMPLI $ 27 27 , small = [1, 2, 1], YTITLE = 'cm', TITLE = title, /PORTRAIT, _extra = ex 28 28 … … 36 36 37 37 title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+blabla 38 pltt, ts_SSH.arr - ts_SSH2.arr, 't', - .1, .1, date1, date2, /REMPLI $38 pltt, ts_SSH.arr - ts_SSH2.arr, 't', -10., 10., date1, date2, /REMPLI $ 39 39 , color = 250, small = [1, 2, 2], YTITLE = 'cm', TITLE = title, /noerase, _extra = ex 40 40 -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_T.pro
r3610 r3999 47 47 48 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' $49 pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so' $ 50 50 , small = [1, 2, 1], YTITLE = varunit, TITLE = title, boxzoom = 1500., /portrait, window = 1, _extra = ex 51 51 52 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' $53 pltt, ts_z, 'zt', -2., 2., inter = .1, date1, date2, /REMPLI, style = 'so0so' $ 54 54 , small = [1, 2, 2], YTITLE = varunit, TITLE = title, boxzoom = 6000., /NOERASE, _extra = ex 55 55 -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro
r3927 r3999 78 78 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 79 79 80 ; cnt = cnt+1 & blabla = 'SNOW Volume'81 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SNOW_Vol, masknp, POSTSCRIPT = postscript, _extra = ex80 ;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 82 ; 83 83 IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_ts_html_body.txt', htmltxt[1:*] -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh
r3927 r3999 18 18 idl_command=/Applications/itt/idl64/bin/idl 19 19 # 20 PS_DIR=$( pwd )/ ps_lim3ada_lim3ada21 PDF_DIR=$( pwd )/ pdf_lim3ada_lim3ada20 PS_DIR=$( pwd )/ts_500yfwb0_ps 21 PDF_DIR=$( pwd )/ts_500yfwb0_pdf 22 22 HTML_DIR=$( pwd )/html 23 23 SAXO_DIR=/Users/sflod/SAXO_DIR … … 34 34 #===================== DATA ===================== 35 35 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D= data_1y_potential_temperature_nomask.nc # Levitus_T37 VAR_SAL_3D=vosaline ; FILE_SAL_3D= data_1y_salinity_nomask.nc # Levitus_S36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 38 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index .nc # Ice Extent North Emisphere42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index .nc # Ice Extent South Emisphere43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index .nc # Ice Area North Emisphere44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index .nc # Ice Area South Emisphere41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent North Emisphere 42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent South Emisphere 43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area North Emisphere 44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area South Emisphere 45 45 # 46 46 # Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing … … 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=2 010123154 DATE1=20010101 ; DATE2=21001231 55 55 #DATE1=00010101 ; DATE2=00101231 56 56 # 57 VAR1_T= thetao ; V1T_PREF=lim3ada; V1T_SUFF=_1Y_grid_T.nc58 VAR1_S= so ; V1S_PREF=lim3ada; V1S_SUFF=_1Y_grid_T.nc59 VAR1_SSH= zos ; V1SSH_PREF=lim3ada; V1SSH_SUFF=_1Y_grid_T.nc60 VAR1_Q= qt ; V1Q_PREF=lim3ada; V1Q_SUFF=_1Y_grid_T.nc61 VAR1_EMP= wfo ; V1EMP_PREF=lim3ada; V1EMP_SUFF=_1Y_grid_T.nc62 VAR1_U= uocetr_eff ; V1U_PREF=lim3ada; V1U_SUFF=_1Y_grid_U.nc63 VAR1_V=vo cetr_eff ; V1V_PREF=lim3ada; V1V_SUFF=_1Y_grid_V.nc64 VAR1_ICE=iiceconc ; V1ICE_PREF= lim3ada; V1ICE_SUFF=_1M_icemod.nc65 VAR1_Ithick=iicethic ; V1It_PREF= lim3ada; V1It_SUFF=_1M_icemod.nc66 VAR1_SNOW=isnowthi ; V1SNOW_PREF=lim3ada; V1SNOW_SUFF=_1M_icemod.nc57 VAR1_T=votemper ; V1T_PREF=500yfwb0 ; V1T_SUFF=_1Y_grid_T.nc 58 VAR1_S=vosaline ; V1S_PREF=500yfwb0 ; V1S_SUFF=_1Y_grid_T.nc 59 VAR1_SSH=sossheig ; V1SSH_PREF=500yfwb0 ; V1SSH_SUFF=_1Y_grid_T.nc 60 VAR1_Q=sohflxdo ; V1Q_PREF=500yfwb0 ; V1Q_SUFF=_1Y_grid_T.nc 61 VAR1_EMP=sowaflup ; V1EMP_PREF=500yfwb0 ; V1EMP_SUFF=_1Y_grid_T.nc 62 VAR1_U=vozocrtx ; V1U_PREF=500yfwb0 ; V1U_SUFF=_1Y_grid_U.nc 63 VAR1_V=vomecrty ; V1V_PREF=500yfwb0 ; V1V_SUFF=_1Y_grid_V.nc 64 VAR1_ICE=iiceconc ; V1ICE_PREF=500yfwb0 ; V1ICE_SUFF=_1M_icemod.nc 65 VAR1_Ithick=iicethic ; V1It_PREF=500yfwb0 ; V1It_SUFF=_1M_icemod.nc 66 VAR1_SNOW=isnowthi ; V1SNOW_PREF=500yfwb0 ; V1SNOW_SUFF=_1M_icemod.nc 67 67 # 68 68 #===================== EXP2 ===================== 69 69 # 70 DATE1_2=20010101 ; DATE2_2=20101231 70 #DATE1_2=20010101 ; DATE2_2=20101231 71 DATE1_2=20010101 ; DATE2_2=21001231 71 72 # 72 VAR2_T= thetao ; V2T_PREF=lim3ada; V2T_SUFF=_1Y_grid_T.nc73 VAR2_S= so ; V2S_PREF=lim3ada; V2S_SUFF=_1Y_grid_T.nc74 VAR2_SSH= zos ; V2SSH_PREF=lim3ada; V2SSH_SUFF=_1Y_grid_T.nc75 VAR2_Q= qt ; V2Q_PREF=lim3ada; V2Q_SUFF=_1Y_grid_T.nc76 VAR2_EMP= wfo ; V2EMP_PREF=lim3ada; V2EMP_SUFF=_1Y_grid_T.nc77 VAR2_U= uocetr_eff ; V2U_PREF=lim3ada; V2U_SUFF=_1Y_grid_U.nc78 VAR2_V=vo cetr_eff ; V2V_PREF=lim3ada; V2V_SUFF=_1Y_grid_V.nc79 VAR2_ICE=iiceconc ; V2ICE_PREF= lim3ada; V2ICE_SUFF=_1M_icemod.nc80 VAR2_Ithick=iicethic ; V2It_PREF= lim3ada; V2It_SUFF=_1M_icemod.nc81 VAR2_SNOW=isnowthi ; V2SNOW_PREF=lim3ada; V2SNOW_SUFF=_1M_icemod.nc73 VAR2_T=votemper ; V2T_PREF=500yfwb0 ; V2T_SUFF=_1Y_grid_T.nc 74 VAR2_S=vosaline ; V2S_PREF=500yfwb0 ; V2S_SUFF=_1Y_grid_T.nc 75 VAR2_SSH=sossheig ; V2SSH_PREF=500yfwb0 ; V2SSH_SUFF=_1Y_grid_T.nc 76 VAR2_Q=sohflxdo ; V2Q_PREF=500yfwb0 ; V2Q_SUFF=_1Y_grid_T.nc 77 VAR2_EMP=sowaflup ; V2EMP_PREF=500yfwb0 ; V2EMP_SUFF=_1Y_grid_T.nc 78 VAR2_U=vozocrtx ; V2U_PREF=500yfwb0 ; V2U_SUFF=_1Y_grid_U.nc 79 VAR2_V=vomecrty ; V2V_PREF=500yfwb0 ; V2V_SUFF=_1Y_grid_V.nc 80 VAR2_ICE=iiceconc ; V2ICE_PREF=500yfwb0 ; V2ICE_SUFF=_1M_icemod.nc 81 VAR2_Ithick=iicethic ; V2It_PREF=500yfwb0 ; V2It_SUFF=_1M_icemod.nc 82 VAR2_SNOW=isnowthi ; V2SNOW_PREF=500yfwb0 ; V2SNOW_SUFF=_1M_icemod.nc 82 83 # 83 84 ######################### Export Variables ############################### -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example1
r3928 r3999 18 18 idl_command=/Applications/itt/idl64/bin/idl 19 19 # 20 PS_DIR=$( pwd )/ps_ core_tmx21 PDF_DIR=$( pwd )/pdf_ core_tmx22 HTML_DIR=$( pwd )/html _core_tmx23 SAXO_DIR=/Users/s masson/SAXO_DIR20 PS_DIR=$( pwd )/ps_tag33_tag331 21 PDF_DIR=$( pwd )/pdf_tag33_tag331 22 HTML_DIR=$( pwd )/html 23 SAXO_DIR=/Users/sflod/SAXO_DIR 24 24 # 25 DIR_DATA=/Users/s masson/tmp/idlorca2/res# path of data in NetCDF format26 DIR_CLIMATO=/Users/s masson/tmp/idlorca2/clim# path of climatological data27 DIR_MASK=/Users/s masson/tmp/idlorca2/clim# path of mask files (ex: subbasins)25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 26 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data 27 DIR_MASK=/Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins) 28 28 # 29 29 #===================== Model GRID ===================== 30 30 # 31 FILE_MESH_MASK=/Users/s masson/data/vp33beta_mesh_mask_7.nc # meshmask31 FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 32 32 FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc # sub-bassin masks 33 33 # 34 34 #===================== DATA ===================== 35 35 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D= data_1y_potential_temperature_nomask.nc # Levitus_T37 VAR_SAL_3D=vosaline ; FILE_SAL_3D= data_1y_salinity_nomask.nc # Levitus_S36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 38 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index .nc # Ice Extent North Emisphere42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index .nc # Ice Extent South Emisphere43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index .nc # Ice Area North Emisphere44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index .nc # Ice Area South Emisphere41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent North Emisphere 42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent South Emisphere 43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area North Emisphere 44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area South Emisphere 45 45 # 46 46 # Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing … … 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= 10101 ; DATE2=2000123154 DATE1=20010101 ; DATE2=21001231 55 55 # 56 56 VAR1_T=thetao ; V1T_PREF=tmx_2000 ; V1T_SUFF=_1Y_grid_T.nc … … 100 100 # 101 101 export DATE1 DATE2 102 export VAR1_T V1T_PREF V1T_SUFF103 export VAR1_S V1S_PREF V1S_SUFF104 export VAR1_SSH V1SSH_PREF V1SSH_SUFF105 export VAR1_Q V1Q_PREF V1Q_SUFF106 export VAR1_EMP V1EMP_PREF V1EMP_SUFF107 export VAR1_U V1U_PREF V1U_SUFF108 export VAR1_V V1V_PREF V1V_SUFF109 export VAR1_ICE V1ICE_PREF V1ICE_SUFF110 export VAR1_Ithick V1It_PREF V1It_SUFF111 export VAR1_SNOW V1 ICE_SNOWV1SNOW_SUFF102 export VAR1_T V1T_PREF V1T_SUFF 103 export VAR1_S V1S_PREF V1S_SUFF 104 export VAR1_SSH V1SSH_PREF V1SSH_SUFF 105 export VAR1_Q V1Q_PREF V1Q_SUFF 106 export VAR1_EMP V1EMP_PREF V1EMP_SUFF 107 export VAR1_U V1U_PREF V1U_SUFF 108 export VAR1_V V1V_PREF V1V_SUFF 109 export VAR1_ICE V1ICE_PREF V1ICE_SUFF 110 export VAR1_Ithick V1It_PREF V1It_SUFF 111 export VAR1_SNOW V1SNOW_PREF V1SNOW_SUFF 112 112 113 113 #===================== EXP2 ===================== … … 122 122 export VAR2_ICE V2ICE_PREF V2ICE_SUFF 123 123 export VAR2_Ithick V2It_PREF V2It_SUFF 124 export VAR2_SNOW V2 ICE_SNOWV2SNOW_SUFF124 export VAR2_SNOW V2SNOW_PREF V2SNOW_SUFF 125 125 # -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example2
r3928 r3999 20 20 PS_DIR=$( pwd )/ps_tag33_tag331 21 21 PDF_DIR=$( pwd )/pdf_tag33_tag331 22 HTML_DIR=$( pwd )/html _tag33_tag33123 SAXO_DIR=/Users/s masson/SAXO_DIR22 HTML_DIR=$( pwd )/html 23 SAXO_DIR=/Users/sflod/SAXO_DIR 24 24 # 25 DIR_DATA=/Users/s masson/tmp/idlorca2/res# path of data in NetCDF format26 DIR_CLIMATO=/Users/s masson/tmp/idlorca2/clim# path of climatological data27 DIR_MASK=/Users/s masson/tmp/idlorca2/clim# path of mask files (ex: subbasins)25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada # path of data in NetCDF format 26 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES # path of climatological data 27 DIR_MASK=/Users/sflod/idl_PLOTS/MASK # path of mask files (ex: subbasins) 28 28 # 29 29 #===================== Model GRID ===================== 30 30 # 31 FILE_MESH_MASK=/Users/s masson/data/vp33beta_mesh_mask_7.nc # meshmask31 FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask 32 32 FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc # sub-bassin masks 33 33 # 34 34 #===================== DATA ===================== 35 35 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D= data_1y_potential_temperature_nomask.nc # Levitus_T37 VAR_SAL_3D=vosaline ; FILE_SAL_3D= data_1y_salinity_nomask.nc # Levitus_S36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 38 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index .nc # Ice Extent North Emisphere42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index .nc # Ice Extent South Emisphere43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index .nc # Ice Area North Emisphere44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index .nc # Ice Area South Emisphere41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent North Emisphere 42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent South Emisphere 43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area North Emisphere 44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area South Emisphere 45 45 # 46 46 # Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing … … 100 100 # 101 101 export DATE1 DATE2 102 export VAR1_T V1T_PREF V1T_SUFF103 export VAR1_S V1S_PREF V1S_SUFF104 export VAR1_SSH V1SSH_PREF V1SSH_SUFF105 export VAR1_Q V1Q_PREF V1Q_SUFF106 export VAR1_EMP V1EMP_PREF V1EMP_SUFF107 export VAR1_U V1U_PREF V1U_SUFF108 export VAR1_V V1V_PREF V1V_SUFF109 export VAR1_ICE V1ICE_PREF V1ICE_SUFF110 export VAR1_Ithick V1It_PREF V1It_SUFF111 export VAR1_SNOW V1 ICE_SNOWV1SNOW_SUFF102 export VAR1_T V1T_PREF V1T_SUFF 103 export VAR1_S V1S_PREF V1S_SUFF 104 export VAR1_SSH V1SSH_PREF V1SSH_SUFF 105 export VAR1_Q V1Q_PREF V1Q_SUFF 106 export VAR1_EMP V1EMP_PREF V1EMP_SUFF 107 export VAR1_U V1U_PREF V1U_SUFF 108 export VAR1_V V1V_PREF V1V_SUFF 109 export VAR1_ICE V1ICE_PREF V1ICE_SUFF 110 export VAR1_Ithick V1It_PREF V1It_SUFF 111 export VAR1_SNOW V1SNOW_PREF V1SNOW_SUFF 112 112 113 113 #===================== EXP2 ===================== … … 122 122 export VAR2_ICE V2ICE_PREF V2ICE_SUFF 123 123 export VAR2_Ithick V2It_PREF V2It_SUFF 124 export VAR2_SNOW V2 ICE_SNOWV2SNOW_SUFF124 export VAR2_SNOW V2SNOW_PREF V2SNOW_SUFF 125 125 # -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example_ipcc
r3927 r3999 34 34 #===================== DATA ===================== 35 35 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D= data_1y_potential_temperature_nomask.nc # Levitus_T37 VAR_SAL_3D=vosaline ; FILE_SAL_3D= data_1y_salinity_nomask.nc # Levitus_S36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 38 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index .nc # Ice Extent North Emisphere42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index .nc # Ice Extent South Emisphere43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index .nc # Ice Area North Emisphere44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index .nc # Ice Area South Emisphere41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent North Emisphere 42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent South Emisphere 43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area North Emisphere 44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area South Emisphere 45 45 # 46 46 # Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing … … 100 100 # 101 101 export DATE1 DATE2 102 export VAR1_T V1T_PREF V1T_SUFF103 export VAR1_S V1S_PREF V1S_SUFF104 export VAR1_SSH V1SSH_PREF V1SSH_SUFF105 export VAR1_Q V1Q_PREF V1Q_SUFF106 export VAR1_EMP V1EMP_PREF V1EMP_SUFF107 export VAR1_U V1U_PREF V1U_SUFF108 export VAR1_V V1V_PREF V1V_SUFF109 export VAR1_ICE V1ICE_PREF V1ICE_SUFF110 export VAR1_Ithick V1It_PREF V1It_SUFF111 export VAR1_SNOW V1 ICE_SNOWV1SNOW_SUFF102 export VAR1_T V1T_PREF V1T_SUFF 103 export VAR1_S V1S_PREF V1S_SUFF 104 export VAR1_SSH V1SSH_PREF V1SSH_SUFF 105 export VAR1_Q V1Q_PREF V1Q_SUFF 106 export VAR1_EMP V1EMP_PREF V1EMP_SUFF 107 export VAR1_U V1U_PREF V1U_SUFF 108 export VAR1_V V1V_PREF V1V_SUFF 109 export VAR1_ICE V1ICE_PREF V1ICE_SUFF 110 export VAR1_Ithick V1It_PREF V1It_SUFF 111 export VAR1_SNOW V1SNOW_PREF V1SNOW_SUFF 112 112 113 113 #===================== EXP2 ===================== … … 122 122 export VAR2_ICE V2ICE_PREF V2ICE_SUFF 123 123 export VAR2_Ithick V2It_PREF V2It_SUFF 124 export VAR2_SNOW V2 ICE_SNOWV2SNOW_SUFF124 export VAR2_SNOW V2SNOW_PREF V2SNOW_SUFF 125 125 # -
branches/2013/dev_r3918_CNRS_idl_plots/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example_old_names
r3927 r3999 34 34 #===================== DATA ===================== 35 35 # 36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D= data_1y_potential_temperature_nomask.nc # Levitus_T37 VAR_SAL_3D=vosaline ; FILE_SAL_3D= data_1y_salinity_nomask.nc # Levitus_S36 VAR_TEMP_3D=votemper ; FILE_TEMP_3D=potT_annual_mean.nc # PHC3 37 VAR_SAL_3D=vosaline ; FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc # PHC3 38 38 VAR_SST=sst ; FILE_SST=NewREY_ORCA2_1991_2000_1y.nc # Reynolds 39 39 VAR_FLUX=qnet ; FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc # flux 40 40 VAR_MLD=mld ; FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc # Mixed layer depth 41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index .nc # Ice Extent North Emisphere42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index .nc # Ice Extent South Emisphere43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index .nc # Ice Area North Emisphere44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index .nc # Ice Area South Emisphere41 VAR_ICE_EXT_NH=extt_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent North Emisphere 42 VAR_ICE_EXT_SH=extt_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Extent South Emisphere 43 VAR_ICE_area_NH=area_NH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area North Emisphere 44 VAR_ICE_area_SH=area_SH ; FILE_ICE=sea_ice_index_2000.nc # Ice Area South Emisphere 45 45 # 46 46 # Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing
Note: See TracChangeset
for help on using the changeset viewer.