New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3888 for branches/2012 – NEMO

Changeset 3888 for branches/2012


Ignore:
Timestamp:
2013-04-23T10:37:54+02:00 (11 years ago)
Author:
flavoni
Message:

update DIL_scripts to add plot of ICE extent and observations

Location:
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_main.sh

    r2751 r3888  
    173173       pdflist=${pdflist}' '${PDF_DIR}/${file}pdf 
    174174   done 
    175    texexec --pdfarrange --result=$PDF_DIR/$filepdf $pdflist 
     175#SF   texexec --pdfarrange --result=$PDF_DIR/$filepdf $pdflist 
     176# needed pdfjam installed 
     177   pdfjoin ${pdflist} --outfile $PDF_DIR/${filepdf} 
    176178   rm -f $PDF_DIR/$( basename $filepdf .pdf ).aux $PDF_DIR/$( basename $filepdf .pdf ).log 
    177    echo ${PDF_DIR}/$filepdf done 
     179   if [ ! -e ${PDF_DIR}/$filepdf ] 
     180       then 
     181       echo "${PDF_DIR}/$filepdf do not exist" 
     182       else 
     183       echo "${PDF_DIR}/$filepdf done" 
     184   fi 
     185#if [ ! -f ${PDF_DIR}/$filepdf ] echo ${PDF_DIR}/$filepdf not exist 
     186#SF   echo ${PDF_DIR}/$filepdf done 
    178187#commented because convert gives images of bad quality 
    179188#  convert ${PS_DIR}/*.ps /tmp/all_${PLOTTYPE}.pdf 
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r2755 r3888  
    2828  std_file_oaflux =  isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato) 
    2929  std_file_mld =  isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato) 
     30  std_file_ice =  isafile(getenv('FILE_ICE'), title = 'ICE', iodir = std_iodir_climato) 
    3031 
    3132  IF strlowcase(getenv('FILE_GEOHEAT')) EQ 'no' THEN std_file_geoheat = 'no' $ 
     
    160161  Ithi_sept_1 = {arr:1./float(cnt) * total(Ithi_1.arr[*, *, temporary(sept)], 3), unit:Ithi_1.unit}  
    161162  undefine, Ithi_1 
     163;simona 
     164;ICE_EXT = read_ncdf(getenv('VAR_ICE'), filename = std_file_ice ) 
     165;simona 
    162166; 
    163167  IF std_file2_I NE std_file1_I THEN BEGIN 
     
    175179; 
    176180  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
     181  help, Ifra_1.arr, jpi,jpj 
    177182  caldat, time, mm 
    178183  march = where(mm EQ 3, cnt) 
     
    180185  sept = where(mm EQ 9, cnt) 
    181186  Ifra_sept_1 = {arr:1./float(cnt) * total(Ifra_1.arr[*, *, temporary(sept)], 3), unit:Ifra_1.unit} 
    182   undefine, Ifra_1 
     187  ;undefine, Ifra_1 
    183188; 
    184189  IF std_file2_I NE std_file1_I THEN BEGIN 
     
    189194    sept = where(mm EQ 9, cnt) 
    190195    Ifra_sept_2 = {arr:1./float(cnt) * total(Ifra_2.arr[*, *, temporary(sept)], 3), unit:Ifra_2.unit} 
    191     undefine, Ifra_2 
     196  ;  undefine, Ifra_2 
    192197  ENDIF ELSE BEGIN  
    193198    Ifra_march_2 = {arr:-1} 
     
    323328  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
    324329 
     330;  cnt = cnt+1   &   blabla = 'Arctic Ice Extend: ' 
     331;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex 
     332 
     333;  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: ' 
     334;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex 
     335; 
    325336  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)' 
    326337  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex 
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro

    r2751 r3888  
    1 pro std_ts_ICE, masknp, s_iodir_data, POSTSCRIPT = postscript, _extra = ex 
     1pro std_ts_ICE, masknp, s_iodir_data,  POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    1010; get exp2 info 
    1111  vICE2 = getenv('VAR2_ICE')   &   prefix2 = getenv('V2ICE_PREF')   &   suffix2 = getenv('V2ICE_SUFF') 
     12; get ice climatology info 
     13  std_file_ice =  isafile(getenv('FILE_ICE'), title = 'ICE Extent Climatology', iodir = std_iodir_climato) 
    1214; 
     15  time_ice = ncdf_lec( std_file_ice, VAR='time' ) 
     16  time_ice = (time_ice - FLOOR(time_ice) ) * 12 
     17  time_ice = round(time_ice) ; arrondir a intero + vicino 
     18  t1 = where(time_ice eq 0) 
     19  t1 = t1[0] ;  gennaio 
     20  t2 = where(time_ice eq 11, count) 
     21  t2 = t2[count-1] ; ultimo giorno di dicembre 
     22  nyear = (t2-t1+1)/12 
     23  vICE_ext_NH = ncdf_lec(std_file_ice, VAR=getenv('VAR_ICE_EXT_NH') ) 
     24  vICE_ext_NH = vICE_ext_NH[t1:t2] 
     25  vICE_ext_NH = reform(vICE_ext_NH,12,nyear); putin 2d array 
     26  vICE_ext_NH = total(vICE_ext_NH,2)/nyear ; faccio il totale sulla seconda dimensione (gli anni) 
     27;  
     28  vICE_ext_SH = ncdf_lec(std_file_ice, VAR=getenv('VAR_ICE_EXT_SH') ) 
     29  vICE_ext_SH = vICE_ext_SH[t1:t2] 
     30  vICE_ext_SH = reform(vICE_ext_SH,12,nyear); putin 2d array 
     31  vICE_ext_SH = total(vICE_ext_SH,2)/nyear ; faccio il totale sulla seconda dimensione (gli anni) 
     32; 
     33  vICE_area_NH = ncdf_lec(std_file_ice, VAR=getenv('VAR_ICE_area_NH') ) 
     34  vICE_area_NH = vICE_area_NH[t1:t2] 
     35  vICE_area_NH = reform(vICE_area_NH,12,nyear); putin 2d array 
     36  vICE_area_NH = total(vICE_area_NH,2)/nyear ; faccio il totale sulla seconda dimensione (gli anni) 
     37; 
     38  vICE_area_SH = ncdf_lec(std_file_ice, VAR=getenv('VAR_ICE_area_SH') ) 
     39  vICE_area_SH = vICE_area_SH[t1:t2] 
     40  vICE_area_SH = reform(vICE_area_SH,12,nyear); putin 2d array 
     41  vICE_area_SH = total(vICE_area_SH,2)/nyear ; faccio il totale sulla seconda dimensione (gli anni) 
     42 
    1343  cdti3 = string(cnt, format = '(i3.3)') 
    1444  print, cdti3 + ') ' + blabla 
    15   filename = cdti3 + '_ts_AMOC_'+prefix 
     45  filename = cdti3 + '_ts_ICE_'+prefix 
    1646  if prefix NE prefix2 then filename = filename + '_'+prefix2 
    17   if KEYWORD_SET(postscript) then openps, filename+'_1.ps', portrait = 1 
     47  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
    1848; 
    1949  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 
    2050; 
    2151  iodir = std_iodir_data 
    22  
     52  ; ICE Extent and Ice Extent minus 15% in NORTH Hemisphere 
    2353  domdef, 0, jpi-1, 30, 90, /xindex 
    24   ICE_N = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 
    25   ICE_N.arr = ICE_N.arr * 1.e-12   &   ICE_N.unit = '10^12 m^2' 
    26   domdef, 0, jpi-1, -90, -30, /xindex 
    27   ICE_S = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 
    28   ICE_S.arr = ICE_S.arr * 1.e-12   &   ICE_S.unit = '10^12 m^2' 
    29  
    30   title = prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.)' 
    31   pltt, ICE_N, 't', 0., 15., date1, date2, /REMPLI, /PORTRAIT $ 
     54  ICE_N = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     55  ICE_N = grossemoyenne(ICE_N.arr, 'xy', /integration, mask2d = masknp) 
     56  if jpt mod 12 ne 0 then stop 
     57  nyr=jpt/12. 
     58  ice_N = reform(ice_n, 12,nyr) 
     59  ice_n = total(ice_n,2)/nyr 
     60  ICE_N = {arr:ICE_N * 1.e-12, unit : '10^12 m^2'} 
     61  ; extend di ghiaccio meno il 15% 
     62  ICE_N_15 = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     63  msk = ICE_N_15.arr gt 0.15 ; remove 0.15% for observations 
     64  ICE_N_15 = grossemoyenne( msk, 'xy', /integration, mask2d = masknp) 
     65  if jpt mod 12 ne 0 then stop 
     66  nyr=jpt/12. 
     67  ice_N_15 = reform(ice_n_15, 12,nyr) 
     68  ice_n_15 = total(ice_n_15,2)/nyr 
     69  ICE_N_15 = {arr:ICE_N_15 * 1.e-12, unit : '10^12 m^2'} 
     70  ; 
     71  ;ICE Extent and Ice Extent minus 15% in NORTH Hemisphere 
     72  domdef, 0, jpi-1, -90, -30, /xindex  
     73  ICE_S = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     74  ICE_S = grossemoyenne(ICE_S.arr, 'xy', /integration, mask2d = masknp) 
     75  if jpt mod 12 ne 0 then stop 
     76  nyr=jpt/12. 
     77  ice_S = reform(ice_S, 12,nyr) 
     78  ice_S = total(ice_S,2)/nyr 
     79  ICE_S = {arr:ICE_S * 1.e-12, unit : '10^12 m^2'} 
     80  ; 
     81  ;extend di ghiaccio meno il 15% 
     82  ICE_S_15 = rseries_ncdf(vICE1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec) 
     83  msk = ICE_S_15.arr gt 0.15 ; remove 0.15% for observations 
     84  ICE_S_15 = grossemoyenne(msk, 'xy', /integration, mask2d = masknp) 
     85  if jpt mod 12 ne 0 then stop 
     86  nyr=jpt/12. 
     87  ice_S_15 = reform(ice_S_15, 12,nyr) 
     88  ice_S_15 = total(ice_S_15,2)/nyr 
     89  ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 
     90  ; 
     91  title = prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.) (Black SOLID simulation) and Extend minus 15% (Blue SOLID simulation)'+'!C'+'Observation (dashed)' 
     92  jpt=12 
     93  time=julday(1,15,1900)+30*lindgen(12) 
     94  pltt, ICE_N, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT,MIN = 4., MAX = 16. $ 
    3295        , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    33    
    34   title = prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)' 
    35   pltt, ICE_S, 't', 0., 15., date1, date2, /REMPLI, /NOERASE $ 
     96  pltt, ICE_N_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
     97        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
     98  pltt,  vICE_area_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate  $ 
     99         , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     100  pltt,  vICE_ext_NH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate  $ 
     101        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     102;;;DA RIMETTERE? return  
     103  title = prefix+' '+d1_d2+' - '+'Global Annual Mean Ice Area (South. Hemisp.) (Black) and Extend minus 15% (Blue)'+'!C'+'Observation (dashed)' 
     104  pltt, ICE_S, 't', 0., 15.,19000101 ,19001231 , /REMPLI, /NOERASE , MIN = 0., MAX = 20. $ 
    36105        , small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    37    
    38   htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_1.png  />  ' ] 
     106  pltt, ICE_S_15, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
     107        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
     108  pltt,  vICE_area_SH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate  $ 
     109         , /ov1d, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     110  pltt,  vICE_ext_SH, 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate  $ 
     111        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     112;;return 
     113  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
    39114  if KEYWORD_SET(postscript) then closeps 
    40115 
     
    45120    tsave = time 
    46121    domdef, 0, jpi-1, 30, 90, /xindex 
    47     ICE_N2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 
    48     ICE_N2.arr = ICE_N2.arr * 1.e-12   &   ICE_N2.unit = '10^12 m^2' 
    49     domdef, 0, jpi-1, -90, -30, /xindex 
    50     ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec, direc = 'xy', /integration, mask2d = masknp) 
    51     ICE_S2.arr = ICE_S2.arr * 1.e-12   &   ICE_S2.unit = '10^12 m^2' 
     122    ICE_N2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec) 
     123    ICE_N2 = grossemoyenne(ICE_N2.arr, 'xy', /integration, mask2d = masknp) 
     124    if jpt mod 12 ne 0 then stop 
     125    nyr=jpt/12. 
     126    ICE_N2 = reform(ICE_N2, 12,nyr) 
     127    ICE_N2 = total(ICE_N2,2)/nyr 
     128    ICE_N2 = {arr:ICE_N2 * 1.e-12, unit : '10^12 m^2'} 
     129    ; extend di ghiaccio meno il 15% 
     130    ICE_N2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix_2, suffix_2, FIRSTONLY = 1 - allrec) 
     131    msk = ICE_N2_15.arr gt 0.15  ; remove 0.15% for observations 
     132    ICE_N2_15 = grossemoyenne( msk, 'xy', /integration, mask2d = masknp) 
     133    if jpt mod 12 ne 0 then stop 
     134    nyr=jpt/12. 
     135    ICE_N2_15 = reform(ICE_N2_15, 12,nyr) 
     136    ICE_N2_15 = total(ICE_N2_15,2)/nyr 
     137    ICE_N2_15 = {arr:ICE_N2_15 * 1.e-12, unit : '10^12 m^2'} 
     138     ;ICE Extent and Ice Extent minus 15% in NORTH Hemisphere 
     139    domdef, 0, jpi-1, -90, -30, /xindex  
     140    ICE_S2 = rseries_ncdf(vICE2, date1_2, date2_2, prefix_2, suffix_2, FIRSTONLY = 1 - allrec) 
     141    ICE_S2 = grossemoyenne(ICE_S2.arr, 'xy', /integration, mask2d = masknp) 
     142    if jpt mod 12 ne 0 then stop 
     143    nyr=jpt/12. 
     144    ICE_S2 = reform(ICE_S2, 12,nyr) 
     145    ICE_S2 = total(ICE_S2,2)/nyr 
     146    ICE_S2 = {arr:ICE_S2 * 1.e-12, unit : '10^12 m^2'} 
     147    ; 
     148    ;extend di ghiaccio meno il 15% 
     149    ICE_S2_15 = rseries_ncdf(vICE2, date1_2, date2_2, prefix_2, suffix_2, FIRSTONLY = 1 - allrec) 
     150    msk = ICE_S2_15.arr gt 0.15  ; remove 0.15% for observations 
     151    ICE_S2_15 = grossemoyenne(msk, 'xy', /integration, mask2d = masknp) 
     152    if jpt mod 12 ne 0 then stop 
     153    nyr=jpt/12. 
     154    ice_S2_15 = reform(ICE_S2_15, 12,nyr) 
     155    ice_S2_15 = total(ICE_S2_15,2)/nyr 
     156    ICE_S2_15 = {arr:ICE_S2_15 * 1.e-12, unit : '10^12 m^2'} 
     157  ; 
    52158    time = tsave   &   IF n_elements(time) NE jpt THEN stop 
    53159 
    54     if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
     160    if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1 
    55161 
    56     title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.)' 
     162    title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (North. Hemisp.) (Red SOLID simulation)  and Extend minus 15% (Blue SOLID simulation)' 
    57163    pltt, ICE_N.arr - ICE_N2.arr, 't', -2., 2., date1, date2, /REMPLI, /PORTRAIT, window = 2 $ 
    58164          , COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     165    pltt, ICE_N_15.arr -  ICE_N2_15.arr , 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
     166      , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    59167     
    60     title = prefix+' '+d1_d2+' - '+prefix2+' '+d1_d2_2+'!C'+'Global Annual Mean Ice Area (South. Hemisp.)' 
     168    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)' 
    61169    pltt, ICE_S.arr - ICE_S2.arr, 't', -2., 2., date1, date2, /REMPLI, /NOERASE $ 
    62170          , COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     171    pltt, ICE_S_15.arr -  ICE_S2_15.arr , 't', 0., 15., 19000101, 19001231, /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
     172      , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
     173     
    63174 
    64     htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ] 
     175    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
    65176    if KEYWORD_SET(postscript) then closeps 
    66177 
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro

    r2751 r3888  
    6262  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_AMOC, masknp, POSTSCRIPT = postscript, _extra = ex 
    6363; 
    64   cnt = cnt+1   &   blabla = 'Sea-Ice cover' 
     64  cnt = cnt+1   &   blabla = 'Sea-Ice cover and Extent' 
    6565  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE, masknp, POSTSCRIPT = postscript, _extra = ex 
    66  
     66; 
     67; 
     68   
    6769  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_ts_html_body.txt', htmltxt[1:*] 
    6870 
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example1

    r2753 r3888  
    4848READ_ONLY_FIRST_RECORD=0   # if 0 then read all records in files else read only the first reacord in each file 
    4949# 
    50 DATE1=10101   ;   DATE2=20001231 
     50DATE1=10101   ;   DATE2=101231 
    5151# 
    5252VAR1_T=thetao     ;   V1T_PREF=tmx_2000     ;   V1T_SUFF=_1Y_grid_T.nc 
     
    5757VAR1_U=uoce_eff   ;   V1U_PREF=tmx_2000     ;   V1U_SUFF=_1Y_grid_U.nc 
    5858VAR1_V=voce_eff   ;   V1V_PREF=tmx_2000     ;   V1V_SUFF=_1Y_grid_V.nc 
    59 VAR1_ICE=sic      ;   V1ICE_PREF=tmx_2000   ;   V1ICE_SUFF=_1Y_icemod.nc 
     59VAR1_ICE=sic      ;   V1ICE_PREF=tmx_2000   ;   V1ICE_SUFF=_1M_icemod.nc 
    6060# 
    6161#===================== EXP2 ===================== 
    6262# 
    63 DATE1_2=10101   ;   DATE2_2=20001231 
     63DATE1_2=10101   ;   DATE2_2=101231 
    6464# 
    6565VAR2_T=thetao     ;   V2T_PREF=core2000     ;   V2T_SUFF=_1Y_grid_T.nc 
     
    7070VAR2_U=uoce_eff   ;   V2U_PREF=core2000     ;   V2U_SUFF=_1Y_grid_U.nc 
    7171VAR2_V=voce_eff   ;   V2V_PREF=core2000     ;   V2V_SUFF=_1Y_grid_V.nc 
    72 VAR2_ICE=sic      ;   V2ICE_PREF=core2000   ;   V2ICE_SUFF=_1Y_icemod.nc 
     72VAR2_ICE=sic      ;   V2ICE_PREF=core2000   ;   V2ICE_SUFF=_1M_icemod.nc 
    7373# 
    7474######################### Export Variables ############################### 
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh_example2

    r2753 r3888  
    5757VAR1_U=uocetr_eff ;   V1U_PREF=tag33     ;   V1U_SUFF=_1Y_grid_U.nc 
    5858VAR1_V=vocetr_eff ;   V1V_PREF=tag33     ;   V1V_SUFF=_1Y_grid_V.nc 
    59 VAR1_ICE=sic      ;   V1ICE_PREF=tag33   ;   V1ICE_SUFF=_1Y_icemod.nc 
     59VAR1_ICE=sic      ;   V1ICE_PREF=tag33   ;   V1ICE_SUFF=_1M_icemod.nc 
    6060# 
    6161#===================== EXP2 ===================== 
     
    7070VAR2_U=uocetr_eff ;   V2U_PREF=trunk331     ;   V2U_SUFF=_1Y_grid_U.nc 
    7171VAR2_V=vocetr_eff ;   V2V_PREF=trunk331     ;   V2V_SUFF=_1Y_grid_V.nc 
    72 VAR2_ICE=sic      ;   V2ICE_PREF=trunk331   ;   V2ICE_SUFF=_1Y_icemod.nc 
     72VAR2_ICE=sic      ;   V2ICE_PREF=trunk331   ;   V2ICE_SUFF=_1M_icemod.nc 
    7373# 
    7474######################### Export Variables ############################### 
Note: See TracChangeset for help on using the changeset viewer.