Changeset 4902 for branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG
- Timestamp:
- 2014-11-27T17:13:38+01:00 (9 years ago)
- Location:
- branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG
- Files:
-
- 15 edited
- 2 copied
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg
r4493 r4902 14 14 cn_ocerst_in = "amm12_restart_oce" ! suffix of ocean restart name (input) 15 15 cn_ocerst_out = "restart_oce_out" ! suffix of ocean restart name (input) 16 nn_istate = 1 ! output the initial state (1) or not (0)17 16 nn_stock = 1296 ! frequency of creation of a restart file (modulo referenced to 1) 18 17 nn_write = 144 ! frequency of write in the output file (modulo referenced to nit000) -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/README
r4394 r4902 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 -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro
r4559 r4902 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 -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
r4559 r4902 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 -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh
r4562 r4902 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 ############################### -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names
r4559 r4902 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 ############################### -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro
r4497 r4902 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 ; -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro
r4514 r4902 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 -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro
r4504 r4902 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 $ -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro
r4559 r4902 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 -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh
r4559 r4902 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 ############################### -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/SHARED/1_namelist_ref
r4900 r4902 31 31 nn_leapy = 0 ! Leap year calendar (1) or not (0) 32 32 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 33 nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=.true. 33 34 nn_rstctl = 0 ! restart control => activated only if ln_rstart = T 34 35 ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist … … 119 120 ! 120 121 rn_rdt = 5760. ! time step for the dynamics (and tracer if nn_acc=0) 121 nn_baro = 64 ! number of barotropic time step ("key_dynspg_ts")122 122 rn_atfp = 0.1 ! asselin time filter parameter 123 123 nn_acc = 0 ! acceleration of convergence : =1 used, rdt < rdttra(k) … … 152 152 / 153 153 !----------------------------------------------------------------------- 154 &namsplit ! time splitting parameters ("key_dynspg_ts") 155 !----------------------------------------------------------------------- 156 ln_bt_fw = .TRUE. ! Forward integration of barotropic equations 157 ln_bt_av = .TRUE. ! Time filtering of barotropic variables 158 ln_bt_nn_auto = .TRUE. ! Set nn_baro automatically to be just below 159 ! a user defined maximum courant number (rn_bt_cmax) 160 nn_baro = 30 ! Number of iterations of barotropic mode 161 ! during rn_rdt seconds. Only used if ln_bt_nn_auto=F 162 rn_bt_cmax = 0.8 ! Maximum courant number allowed if ln_bt_nn_auto=T 163 nn_bt_flt = 1 ! Time filter choice 164 ! = 0 None 165 ! = 1 Boxcar over nn_baro barotropic steps 166 ! = 2 Boxcar over 2*nn_baro " " 167 / 168 !----------------------------------------------------------------------- 154 169 &namcrs ! Grid coarsening for dynamics output and/or 155 170 ! passive tracer coarsened online simulations … … 168 183 / 169 184 !----------------------------------------------------------------------- 185 &namc1d ! 1D configuration options ("key_c1d") 186 !----------------------------------------------------------------------- 187 rn_lat1d = 50 ! Column latitude (default at PAPA station) 188 rn_lon1d = -145 ! Column longitude (default at PAPA station) 189 ln_c1d_locpt= .true. ! Localization of 1D config in a grid (T) or independant point (F) 190 / 191 !----------------------------------------------------------------------- 170 192 &namtsd ! data : Temperature & Salinity 171 193 !----------------------------------------------------------------------- … … 173 195 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 174 196 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 175 sn_tem = 'data_1m_potential_temperature_nomask', -1 ,'votemper' , .true. , .true. , 'yearly' , '' , '', ''176 sn_sal = 'data_1m_salinity_nomask' , -1 ,'vosaline' , .true. , .true. , 'yearly' , '' , '', ''197 sn_tem = 'data_1m_potential_temperature_nomask', -1 ,'votemper' , .true. , .true. , 'yearly' , '' , '' , '' 198 sn_sal = 'data_1m_salinity_nomask' , -1 ,'vosaline' , .true. , .true. , 'yearly' , '' , '' , '' 177 199 ! 178 200 cn_dir = './' ! root directory for the location of the runoff files 179 201 ln_tsd_init = .true. ! Initialisation of ocean T & S with T &S input data (T) or not (F) 180 ln_tsd_tradmp = . false. ! damping of ocean T & S toward T &S input data (T) or not (F)202 ln_tsd_tradmp = .true. ! damping of ocean T & S toward T &S input data (T) or not (F) 181 203 / 182 204 !!====================================================================== … … 270 292 &namsbc_core ! namsbc_core CORE bulk formulae 271 293 !----------------------------------------------------------------------- 272 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !273 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !274 sn_wndi = 'u_10.15JUNE2009_fill' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , 'weights_core2_agrif_bicubic.nc' , 'Uwnd' , ''275 sn_wndj = 'v_10.15JUNE2009_fill' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , 'weights_core2_agrif_bicubic.nc' , 'Vwnd' , ''294 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 295 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 296 sn_wndi = 'u_10.15JUNE2009_fill' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , 'weights_core2_agrif_bicubic.nc' , 'Uwnd' , '' 297 sn_wndj = 'v_10.15JUNE2009_fill' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , 'weights_core2_agrif_bicubic.nc' , 'Vwnd' , '' 276 298 sn_qsr = 'ncar_rad.15JUNE2009_fill' , 24 , 'SWDN_MOD', .false. , .true. , 'yearly' , 'weights_core2_agrif_bilinear.nc' , '' , '' 277 299 sn_qlw = 'ncar_rad.15JUNE2009_fill' , 24 , 'LWDN_MOD', .false. , .true. , 'yearly' , 'weights_core2_agrif_bilinear.nc' , '' , '' … … 348 370 &namtra_qsr ! penetrative solar radiation 349 371 !----------------------------------------------------------------------- 350 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights 351 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename 352 sn_chl ='chlorophyll', -1 , 'CHLA' , .true. , .true. , 'yearly' , 'weights_bilin' 372 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 373 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 374 sn_chl ='chlorophyll', -1 , 'CHLA' , .true. , .true. , 'yearly' , 'weights_bilin' , '' , '' 353 375 354 376 cn_dir = './' ! root directory for the location of the runoff files … … 546 568 nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state 547 569 ! = 1, bdy data are read in 'bdydata .nc' files 570 cn_ice_lim = 'none' ! 571 nn_ice_lim_dta = 0 ! = 0, bdy data are equal to the initial state 572 ! = 1, bdy data are read in 'bdydata .nc' files 573 rn_ice_tem = 270. ! lim3 only: arbitrary temperature of incoming sea ice 574 rn_ice_sal = 10. ! lim3 only: -- salinity -- 575 rn_ice_age = 30. ! lim3 only: -- age -- 576 548 577 ln_tra_dmp =.false. ! open boudaries conditions for tracers 549 578 ln_dyn3d_dmp =.false. ! open boundary condition for baroclinic velocities … … 566 595 bn_tem = 'amm12_bdyT_tra' , 24 , 'votemper' , .true. , .false. , 'daily' , '' , '' , '' 567 596 bn_sal = 'amm12_bdyT_tra' , 24 , 'vosaline' , .true. , .false. , 'daily' , '' , '' , '' 597 ! for lim2 598 ! bn_frld = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true. , .false. , 'daily' , '' , '' , '' 599 ! bn_hicif = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true. , .false. , 'daily' , '' , '' , '' 600 ! bn_hsnif = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true. , .false. , 'daily' , '' , '' , '' 601 ! for lim3 602 ! bn_a_i = 'amm12_bdyT_ice' , 24 , 'ileadfra' , .true. , .false. , 'daily' , '' , '' , '' 603 ! bn_ht_i = 'amm12_bdyT_ice' , 24 , 'iicethic' , .true. , .false. , 'daily' , '' , '' , '' 604 ! bn_ht_s = 'amm12_bdyT_ice' , 24 , 'isnowthi' , .true. , .false. , 'daily' , '' , '' , '' 568 605 cn_dir = 'bdydta/' 569 606 ln_full_vel = .false. … … 590 627 ! = 2 : nonlinear friction 591 628 rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 592 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) 629 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 630 rn_bfri2_max = 1.e-1 ! max. bottom drag coefficient (non linear case and ln_loglayer=T) 593 631 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) 594 rn_bfrz0 = 3.e-3 ! bottom roughness for loglayer bfr coeff 632 rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T 633 ln_loglayer = .false. ! logarithmic formulation (non linear case) 595 634 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 596 635 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T) … … 726 765 ln_vvl_layer = .false. ! full layer vertical coordinate 727 766 ln_vvl_ztilde_as_zstar = .false. ! ztilde vertical coordinate emulating zstar 767 ln_vvl_zstar_at_eqtor = .false. ! ztilde near the equator 728 768 rn_ahe3 = 0.0e0 ! thickness diffusion coefficient 729 769 rn_rst_e3t = 30.e0 ! ztilde to zstar restoration timescale [days] … … 941 981 ! (no physical validity of the results) 942 982 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 943 /944 !-----------------------------------------------------------------------945 &namc1d ! 1D configuration options ("key_c1d")946 !-----------------------------------------------------------------------947 rn_lat = 50 ! Column latitude948 rn_lon = -145 ! Column longitude949 983 / 950 984 !----------------------------------------------------------------------- -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/SHARED/field_def.xml
r4901 r4902 225 225 <field id="emp_x_sst" long_name="Concentration/Dilution term on SST" unit="kgC/m2/s" /> 226 226 <field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kgPSU/m2/s" /> 227 228 227 229 228 <field id="iceconc" long_name="ice concentration" unit="%" /> 230 <field id="icebopr" long_name="daily bottom thermo ice prod." unit="km3/day" />231 <field id="icedypr" long_name="daily dynamic ice prod." unit="km3/day" />232 <field id="ioceflxb" long_name="Oceanic flux at the ice base" unit="W/m2" />233 229 <field id="uice_ipa" long_name="Ice velocity along i-axis at I-point (ice presence average)" unit="m/s" /> 234 230 <field id="vice_ipa" long_name="Ice velocity along j-axis at I-point (ice presence average)" unit="m/s" /> … … 238 234 <field id="qsr_oce" long_name="solar heat flux at ocean surface" unit="W/m2" /> 239 235 <field id="qns_oce" long_name="non-solar heat flux at ocean surface" unit="W/m2" /> 240 <field id="hfbri" long_name="heat flux due to brine release" unit="W/m2" /> 236 <field id="qt_ice" long_name="total flux at ice surface" unit="W/m2" /> 237 <field id="qsr_ice" long_name="solar heat flux at ice surface" unit="W/m2" /> 238 <field id="qns_ice" long_name="non-solar heat flux at ice surface" unit="W/m2" /> 239 <field id="qtr_ice" long_name="solar heat flux transmitted thru the ice" unit="W/m2" /> 241 240 <field id="utau_ice" long_name="Wind stress along i-axis over the ice at i-point" unit="N/m2" /> 242 241 <field id="vtau_ice" long_name="Wind stress along j-axis over the ice at i-point" unit="N/m2" /> 243 <field id="qsr_io" long_name="Ice-Oce downward solar heat flux" unit="W/m2" />244 <field id="qns_io" long_name="Ice-Oce downward non-solar heat flux" unit="W/m2" />245 242 <field id="micesalt" long_name="Mean ice salinity" unit="psu" /> 246 243 <field id="miceage" long_name="Mean ice age" unit="years" /> 247 <field id="icelapr" long_name="daily lateral thermo ice prod." unit="km3/day" /> 248 <field id="icesipr" long_name="daily snowice ice prod." unit="km3/day" /> 244 245 <field id="iceage_cat" long_name="Ice age for categories" unit="days" axis_ref="ncatice" /> 246 <field id="iceconc_cat" long_name="Ice concentration for categories" unit="%" axis_ref="ncatice" /> 247 <field id="icethic_cat" long_name="Ice thickness for categories" unit="m" axis_ref="ncatice" /> 248 <field id="snowthic_cat" long_name="Snow thicknessi for categories" unit="m" axis_ref="ncatice" /> 249 <field id="salinity_cat" long_name="Sea-Ice Bulk salinity for categories" unit="ppt" axis_ref="ncatice" /> 250 <field id="brinevol_cat" long_name="Brine volume for categories" unit="%" axis_ref="ncatice" /> 251 249 252 <field id="micet" long_name="Mean ice temperature" unit="degC" /> 250 <field id="icehc" long_name="ice total heat content" unit="10^9 253 <field id="icehc" long_name="ice total heat content" unit="10^9J" /> 251 254 <field id="isnowhc" long_name="snow total heat content" unit="10^9J" /> 252 255 <field id="icest" long_name="ice surface temperature" unit="degC" /> 253 <field id="sfxbri" long_name="brine salt flux" unit="psu*kg/m2/day" />254 <field id="sfxthd" long_name="equivalent FW salt flux" unit="psu*kg/m2/day" />255 256 <field id="ibrinv" long_name="brine volume" unit="%" /> 256 257 <field id="icecolf" long_name="frazil ice collection thickness" unit="m" /> 257 258 <field id="icestr" long_name="ice strength" unit="N/m" /> 258 259 <field id="icevel" long_name="ice velocity" unit="m/s" /> 259 <field id="isume" long_name="surface melt" unit="km3/day" />260 <field id="ibome" long_name="bottom melt" unit="km3/day" />261 260 <field id="idive" long_name="divergence" unit="10-8s-1" /> 262 261 <field id="ishear" long_name="shear" unit="10-8s-1" /> 263 <field id="icerepr" long_name="daily resultant ice prod./melting from limupdate" unit="km3/day" /> 264 <field id="icevolu" long_name="ice volume" unit="km3" /> 265 <field id="snowvol" long_name="snow volume" unit="km3" /> 266 <field id="sfxmec" long_name="salt flux from ridging rafting" unit="psu*kg/m2/day" /> 262 <field id="icevolu" long_name="ice volume" unit="m" /> 263 <field id="snowvol" long_name="snow volume" unit="m" /> 264 265 <field id="icetrp" long_name="ice volume transport" unit="m/day" /> 266 <field id="snwtrp" long_name="snw volume transport" unit="m/day" /> 267 <field id="deitrp" long_name="advected ice enhalpy" unit="W/m2" /> 268 <field id="destrp" long_name="advected snw enhalpy" unit="W/m2" /> 269 270 <field id="sfxbri" long_name="brine salt flux" unit="psu*kg/m2/day" /> 271 <field id="sfxdyn" long_name="salt flux from ridging rafting" unit="psu*kg/m2/day" /> 267 272 <field id="sfxres" long_name="salt flux from lipupdate (resultant)" unit="psu*kg/m2/day" /> 268 <field id="icetrp" long_name="ice volume transport" unit="km3/day" /> 269 270 273 <field id="sfxbog" long_name="salt flux from bot growth" unit="psu*kg/m2/day" /> 274 <field id="sfxbom" long_name="salt flux from bot melt" unit="psu*kg/m2/day" /> 275 <field id="sfxsum" long_name="salt flux from surf melt" unit="psu*kg/m2/day" /> 276 <field id="sfxsni" long_name="salt flux from snow-ice formation" unit="psu*kg/m2/day" /> 277 <field id="sfxopw" long_name="salt flux from open water ice formation" unit="psu*kg/m2/day" /> 278 <field id="sfx" long_name="salt flux total" unit="psu*kg/m2/day" /> 279 280 <field id="vfxbog" long_name="daily bottom thermo ice prod." unit="m/day" /> 281 <field id="vfxdyn" long_name="daily dynamic ice prod." unit="m/day" /> 282 <field id="vfxopw" long_name="daily lateral thermo ice prod." unit="m/day" /> 283 <field id="vfxsni" long_name="daily snowice ice prod." unit="m/day" /> 284 <field id="vfxsum" long_name="surface melt" unit="m/day" /> 285 <field id="vfxbom" long_name="bottom melt" unit="m/day" /> 286 <field id="vfxres" long_name="daily resultant ice prod./melting from limupdate" unit="m/day" /> 287 <field id="vfxice" long_name="ice melt/growth" unit="m/day" /> 288 <field id="vfxsnw" long_name="snw melt/growth" unit="m/day" /> 289 <field id="vfxsub" long_name="snw sublimation" unit="m/day" /> 290 <field id="vfxspr" long_name="snw precipitation on ice" unit="m/day" /> 291 292 <field id="hfxsum" long_name="heat fluxes causing surface ice melt" unit="W/m2" /> 293 <field id="hfxbom" long_name="heat fluxes causing bottom ice melt" unit="W/m2" /> 294 <field id="hfxbog" long_name="heat fluxes causing bottom ice growth" unit="W/m2" /> 295 <field id="hfxdif" long_name="heat fluxes causing ice temperature change" unit="W/m2" /> 296 <field id="hfxopw" long_name="heat fluxes causing open water ice formation" unit="W/m2" /> 297 <field id="hfxsnw" long_name="heat fluxes causing snow melt" unit="W/m2" /> 298 <field id="hfxerr" long_name="heat fluxes error after heat diffusion" unit="W/m2" /> 299 <field id="hfxerr_rem" long_name="heat fluxes error after remapping" unit="W/m2" /> 300 <field id="hfxout" long_name="total heat fluxes received by the ocean" unit="W/m2" /> 301 <field id="hfxin" long_name="total heat fluxes at the ice/ocean surface" unit="W/m2" /> 302 303 <!-- heat flux associated with mass exchange --> 304 <field id="hfxthd" long_name="heat fluxes from ice-ocean mass exchange during thermo" unit="W/m2" /> 305 <field id="hfxdyn" long_name="heat fluxes from ice-ocean mass exchange during dynamic" unit="W/m2" /> 306 <field id="hfxres" long_name="heat fluxes from ice-ocean mass exchange during resultant" unit="W/m2" /> 307 <field id="hfxsub" long_name="heat fluxes from ice-atm. mass exchange during sublimation" unit="W/m2" /> 308 <field id="hfxspr" long_name="heat fluxes from ice-atm. mass exchange during snow precip" unit="W/m2" /> 309 310 <!-- diags --> 311 <field id="hfxdhc" long_name="Heat content variation in snow and ice" unit="W/m2" /> 312 <field id="hfxtur" long_name="turbulent heat flux at the ice base" unit="W/m2" /> 313 314 271 315 </field_group> 272 316 -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref
r4900 r4902 52 52 ecc = 2.0 ! eccentricity of the elliptical yield curve 53 53 ahi0 = 350.e0 ! horizontal eddy diffusivity coefficient for sea-ice [m2/s] 54 nevp = 300 ! number of iterations for subcycling in EVP, SB, 40055 telast =9600.0 ! timescale for elastic waves, SB, 720.054 nevp = 120 ! number of iterations for subcycling in EVP 55 relast = 0.333 ! ratio of elastic timescale over ice time step (1/3 if nevp=120 ; 1/9 if nevp=300) 56 56 alphaevp = 1.0 ! coefficient for the solution of internal ice stresses 57 57 hminrhg = 0.001 ! ice volume (a*h in m) below which ice velocity equal ocean velocity -
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/SHARED/namelist_ref
r4901 r4902 244 244 ! =1 global mean of e-p-r set to zero at each time step 245 245 ! =2 annual global mean of e-p-r set to zero 246 ! =3 global emp set to zero and spread out over erp area247 246 ln_wave = .false. ! Activate coupling with wave (either Stokes Drift or Drag coefficient, or both) (T => fill namsbc_wave) 248 247 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => fill namsbc_wave)
Note: See TracChangeset
for help on using the changeset viewer.