Changeset 19 for trunk/src


Ignore:
Timestamp:
12/17/10 16:18:08 (13 years ago)
Author:
pinsard
Message:

remove hard coded directory for mask_oaflux_30N30S.nc

Location:
trunk/src
Files:
10 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/TropFlux_19890101_20091231.pro

    r12 r19  
    2121;           rankdir="LR", 
    2222;           ] 
    23 ;           mask [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc"]; 
     23;           mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/mask_oaflux_30N30S.nc"]; 
    2424;           file_sst [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_sst_19890101_20091231_v20.nc"]; 
    2525;           file_wind [shape=ellipse,fontname=Courier,label="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc"]; 
    … …  
    4545; ======== 
    4646; 
     47; :ref:`tropflux_profile.sh` 
     48; 
    4749; :func:`cor30a` 
    4850; 
    … …  
    6163; coding rules 
    6264; 
     65; KNOWN ISSUES 
     66; ============ 
     67; 
     68; test of existence of fullfilename_msk not very efficient because 
     69; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 
     70; 
    6371; EVOLUTIONS 
    6472; ========== 
    6573; 
     74; - fplod 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin) 
     75; 
     76;   * remove hard coded directory for mask_oaflux_30N30S.nc 
     77; 
    6678; - fplod 20101214T112131Z aedon.locean-ipsl.upmc.fr (Darwin) 
    6779; 
    … …  
    8092pro TropFlux_19890101_20091231 
    8193@common 
     94; 
     95; check for input directory  
     96; 
     97; test if ${TROPFLUX_ID} defined 
     98tropflux_id_env=GETENV('TROPFLUX_ID') 
     99CASE tropflux_id_env OF 
     100    ''  :  BEGIN 
     101     msg = 'eee : ${TROPFLUX_ID} is not defined' 
     102     ras = report(msg) 
     103     STOP 
     104           END 
     105 ELSE: BEGIN 
     106     msg = 'iii : ${TROPFLUX_ID} is ' + tropflux_id_env 
     107     ras = report(msg) 
     108       END 
     109ENDCASE 
     110; 
     111iodirin = isadirectory(tropflux_id_env) 
     112; 
     113; existence and protection of ${TROPFLUX_ID} 
     114IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 
     115   msg = 'eee : the directory' + iodirin  + ' is not accessible.' 
     116   ras = report(msg) 
     117   STOP 
     118ENDIF 
     119; 
     120; build mask filename 
     121filename_msk='mask_oaflux_30N30S.nc' 
     122; 
     123; check if this file exists 
     124fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 
     125IF fullfilename_msk[0] EQ '' THEN BEGIN 
     126   msg = 'eee : the file ' + fullfilename_msk + ' was not found.' 
     127   ras = report(msg) 
     128   STOP 
     129ENDIF 
    82130 
    83131da1=19880101 & da2=20101231 
    84 file='/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc' 
    85 initncdf, file 
    86 msk=ncdf_lec(file,var='msk') 
     132; 
     133initncdf, fullfilename_msk  
     134msk=ncdf_lec(fullfilename_msk,var='msk') 
    87135 
    88136dir='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/' 

  • trunk/src/TropFlux_NRT_ncdf.pro

    r12 r19  
    1717;           file3 [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_19890101_20091231_v51.nc"]; 
    1818;           ncfile [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_19890101_20091231.nc"]; 
    19 ;           TropFlux_NRT_ncdf_v50  [shape=box, 
     19;           TropFlux_NRT_ncdf_v50 [shape=box, 
    2020;           fontname=Courier, 
    2121;           color=blue, 

  • trunk/src/interp_erai_dewt_1989_2009.pro

    r12 r19  
    1515;           ] 
    1616;           file_in [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/ERAI_global/20c3m_erai_d2_TROP_1989_2009.nc"]; 
    17 ;           mask [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc"]; 
     17;           mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/mask_oaflux_30N30S.nc"]; 
    1818; 
    1919;           file_out [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/TropFlux/input_uncor/erai_d2m_19890101_20091231_oafluxgrid.nc"]; 
    … …  
    3131; ======== 
    3232; 
     33; :ref:`tropflux_profile.sh` 
     34; 
    3335; EXAMPLES 
    3436; ======== 
    … …  
    4547; coding rules 
    4648; 
     49; KNOWN ISSUES 
     50; ============ 
     51; 
     52; test of existence of fullfilename_msk not very efficient because 
     53; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 
     54; 
    4755; EVOLUTIONS 
    4856; ========== 
     57; 
     58; - fplod 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin) 
     59; 
     60;   * remove hard coded directory for mask_oaflux_30N30S.nc 
    4961; 
    5062; - fplod 20101215T093710Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    6375pro interp_erai_dewt_1989_2009 
    6476@common 
    65  
     77; 
     78; check for input directory  
     79; 
     80; test if ${TROPFLUX_ID} defined 
     81tropflux_id_env=GETENV('TROPFLUX_ID') 
     82CASE tropflux_id_env OF 
     83    ''  :  BEGIN 
     84     msg = 'eee : ${TROPFLUX_ID} is not defined' 
     85     ras = report(msg) 
     86     STOP 
     87           END 
     88 ELSE: BEGIN 
     89     msg = 'iii : ${TROPFLUX_ID} is ' + tropflux_id_env 
     90     ras = report(msg) 
     91       END 
     92ENDCASE 
     93; 
     94iodirin = isadirectory(tropflux_id_env) 
     95; 
     96; existence and protection of ${TROPFLUX_ID} 
     97IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 
     98   msg = 'eee : the directory' + iodirin  + ' is not accessible.' 
     99   ras = report(msg) 
     100   STOP 
     101ENDIF 
     102; 
     103; build mask filename 
     104filename_msk='mask_oaflux_30N30S.nc' 
     105; 
     106; check if this file exists 
     107fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 
     108IF fullfilename_msk[0] EQ '' THEN BEGIN 
     109   msg = 'eee : the file ' + fullfilename_msk + ' was not found.' 
     110   ras = report(msg) 
     111   STOP 
     112ENDIF 
     113; 
    66114filein='/Volumes/PRAVEEN/ERAI_global/20c3m_erai_d2_TROP_1989_2009.nc' 
    67 gridout='/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc' 
    68115 
    69116initncdf, filein 
    … …  
    78125mskin=glamt*0.+1. 
    79126 
    80 initncdf, gridout 
     127initncdf, fullfilename_msk 
    81128domdef 
    82129latout=reform(gphit(0,*)) & lonout=reform(glamt(*,0)) 
    83130print, 'lat grid ',min(latout),max(latout),latout(1)-latout(0) 
    84131print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0) 
    85 mskout=read_ncdf("msk", file=gridout,/nostr) 
     132mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 
    86133 
    87134help, d2min,lonin,latin,mskin,lonout,latout,mskout 
    … …  
    99146timein=timein & jptin=jpt 
    100147 
    101 initncdf, gridout 
     148initncdf, fullfilename_msk 
    102149cda0=string(jul2date(timein(0)),format='(i8.8)') 
    103150cda1=string(jul2date(timein(jpt-1)),format='(i8.8)') 

  • trunk/src/interp_erai_lwr_1989_2009.pro

    r12 r19  
    11;+ 
    22; 
    3 ; .. interp_erai_lwr_1989_2009.pro: 
     3; .. _interp_erai_lwr_1989_2009.pro: 
    44; 
    55; ============================= 
    … …  
    1313;           rankdir="LR", 
    1414;           ] 
    15 ;           mask [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc"]; 
     15;           mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/mask_oaflux_30N30S.nc"]; 
    1616; 
    1717;           file_in [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/ERAI_global/20c3m_erai_str_TROP_1989_2009.nc"]; 
    … …  
    3232; ======== 
    3333; 
     34; :ref:`tropflux_profile.sh` 
     35; 
    3436; EXAMPLES 
    3537; ======== 
    … …  
    4648; coding rules 
    4749; 
     50; KNOWN ISSUES 
     51; ============ 
     52; 
     53; test of existence of fullfilename_msk not very efficient because 
     54; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 
     55; 
    4856; EVOLUTIONS 
    4957; ========== 
     58; 
     59; - fplod 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin) 
     60; 
     61;   * remove hard coded directory for mask_oaflux_30N30S.nc 
    5062; 
    5163; - fplod 20101215T094137Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    6476pro interp_erai_lwr_1989_2009 
    6577@common 
     78; 
     79; check for input directory  
     80; 
     81; test if ${TROPFLUX_ID} defined 
     82tropflux_id_env=GETENV('TROPFLUX_ID') 
     83CASE tropflux_id_env OF 
     84    ''  :  BEGIN 
     85     msg = 'eee : ${TROPFLUX_ID} is not defined' 
     86     ras = report(msg) 
     87     STOP 
     88           END 
     89 ELSE: BEGIN 
     90     msg = 'iii : ${TROPFLUX_ID} is ' + tropflux_id_env 
     91     ras = report(msg) 
     92       END 
     93ENDCASE 
     94; 
     95iodirin = isadirectory(tropflux_id_env) 
     96; 
     97; existence and protection of ${TROPFLUX_ID} 
     98IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 
     99   msg = 'eee : the directory' + iodirin  + ' is not accessible.' 
     100   ras = report(msg) 
     101   STOP 
     102ENDIF 
     103; 
     104; build mask filename 
     105filename_msk='mask_oaflux_30N30S.nc' 
     106; 
     107; check if this file exists 
     108fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 
     109IF fullfilename_msk[0] EQ '' THEN BEGIN 
     110   msg = 'eee : the file ' + fullfilename_msk + ' was not found.' 
     111   ras = report(msg) 
     112   STOP 
     113ENDIF 
    66114 
    67115st=19890101 & en=20090801 
    … …  
    72120lwrin=-1*lwrin 
    73121help, lwrin 
    74  
    75 gridout='/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc' 
    76122 
    77123initncdf, filein 
    … …  
    85131mskin=glamt*0.+1. 
    86132 
    87 initncdf, gridout 
     133initncdf, fullfilename_msk 
    88134domdef 
    89135latout=reform(gphit(0,*)) & lonout=reform(glamt(*,0)) 
    90136print, 'lat grid ',min(latout),max(latout),latout(1)-latout(0) 
    91137print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0) 
    92 mskout=read_ncdf("msk", file=gridout,/nostr) 
     138mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 
    93139 
    94140help, lwrin 
    … …  
    113159xlon=reform(glamt(*,0) ) & ylat=reform(gphit(0,*)) 
    114160 
    115 initncdf, gridout 
     161initncdf, fullfilename_msk 
    116162valmask=1.e20 
    117163 

  • trunk/src/interp_erai_sst_1989_2009.pro

    r12 r19  
    1515;           file_in [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/ERAI_global/20c3m_erai_sstk_TROP_1989_2009.nc"]; 
    1616; 
    17 ;           mask [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc"]; 
     17;           mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/mask_oaflux_30N30S.nc"]; 
    1818; 
    1919;           ncfile [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/TropFlux/input_uncor/erai_sst_19890101_20091231_oafluxgrid.nc"]; 
    … …  
    3131; ======== 
    3232; 
     33; :ref:`tropflux_profile.sh` 
     34; 
    3335; EXAMPLES 
    3436; ======== 
    … …  
    4547; coding rules 
    4648; 
     49; KNOWN ISSUES 
     50; ============ 
     51; 
     52; test of existence of fullfilename_msk not very efficient because 
     53; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 
     54; 
    4755; EVOLUTIONS 
    4856; ========== 
     57; 
     58; - fplod 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin) 
     59; 
     60;   * remove hard coded directory for mask_oaflux_30N30S.nc 
    4961; 
    5062; - fplod 20101215T100932Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    6476pro interp_erai_sst_1989_2009 
    6577@common 
     78; 
     79; check for input directory  
     80; 
     81; test if ${TROPFLUX_ID} defined 
     82tropflux_id_env=GETENV('TROPFLUX_ID') 
     83CASE tropflux_id_env OF 
     84    ''  :  BEGIN 
     85     msg = 'eee : ${TROPFLUX_ID} is not defined' 
     86     ras = report(msg) 
     87     STOP 
     88           END 
     89 ELSE: BEGIN 
     90     msg = 'iii : ${TROPFLUX_ID} is ' + tropflux_id_env 
     91     ras = report(msg) 
     92       END 
     93ENDCASE 
     94; 
     95iodirin = isadirectory(tropflux_id_env) 
     96; 
     97; existence and protection of ${TROPFLUX_ID} 
     98IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 
     99   msg = 'eee : the directory' + iodirin  + ' is not accessible.' 
     100   ras = report(msg) 
     101   STOP 
     102ENDIF 
     103; 
     104; build mask filename 
     105filename_msk='mask_oaflux_30N30S.nc' 
     106; 
     107; check if this file exists 
     108fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 
     109IF fullfilename_msk[0] EQ '' THEN BEGIN 
     110   msg = 'eee : the file ' + fullfilename_msk + ' was not found.' 
     111   ras = report(msg) 
     112   STOP 
     113ENDIF 
    66114 
    67115filein='/Volumes/PRAVEEN/ERAI_global/20c3m_erai_sstk_TROP_1989_2009.nc' 
    68 gridout='/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc' 
    69116 
    70117initncdf, filein 
    … …  
    80127;mskin=glamt*0.+1. 
    81128 
    82 initncdf, gridout 
     129initncdf, fullfilename_msk 
    83130domdef 
    84131latout=reform(gphit(0,*)) & lonout=reform(glamt(*,0)) 
    85132print, 'lat grid ',min(latout),max(latout),latout(1)-latout(0) 
    86133print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0) 
    87 mskout=read_ncdf("msk", file=gridout,/nostr) 
     134mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 
    88135;stop 
    89136help, sstkin,lonin,latin,mskin,lonout,latout,mskout 
    … …  
    101148timein=timein & jptin=jpt 
    102149 
    103 initncdf, gridout 
     150initncdf, fullfilename_msk 
    104151time=timegen(7670, units='days', start=julday(1,1,1989)) & jpt=n_elements(time) 
    105152tt=time 

  • trunk/src/interp_erai_t2m_1989_2009.pro

    r12 r19  
    11;+ 
    22; 
    3 ; .. interp_erai_t2m_1989_2009.pro: 
     3; .. _interp_erai_t2m_1989_2009.pro: 
    44; 
    55; ============================= 
    … …  
    1515;           ] 
    1616;           file_in [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/ERAI_global/20c3m_erai_t2_TROP_1989_2009.nc"]; 
    17 ;           mask [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc"]; 
     17;           mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/mask_oaflux_30N30S.nc"]; 
    1818; 
    1919;           ncfile [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 
    … …  
    3333; ======== 
    3434; 
     35; :ref:`tropflux_profile.sh` 
     36; 
    3537; EXAMPLES 
    3638; ======== 
    … …  
    4850; coding rules 
    4951; 
     52; KNOWN ISSUES 
     53; ============ 
     54; 
     55; test of existence of fullfilename_msk not very efficient because 
     56; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 
     57; 
    5058; EVOLUTIONS 
    5159; ========== 
     60; 
     61; - fplod 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin) 
     62; 
     63;   * remove hard coded directory for mask_oaflux_30N30S.nc 
    5264; 
    5365; - fplod 20101215T112140Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    6678pro interp_erai_t2m_1989_2009 
    6779@common 
     80; 
     81; check for input directory  
     82; 
     83; test if ${TROPFLUX_ID} defined 
     84tropflux_id_env=GETENV('TROPFLUX_ID') 
     85CASE tropflux_id_env OF 
     86    ''  :  BEGIN 
     87     msg = 'eee : ${TROPFLUX_ID} is not defined' 
     88     ras = report(msg) 
     89     STOP 
     90           END 
     91 ELSE: BEGIN 
     92     msg = 'iii : ${TROPFLUX_ID} is ' + tropflux_id_env 
     93     ras = report(msg) 
     94       END 
     95ENDCASE 
     96; 
     97iodirin = isadirectory(tropflux_id_env) 
     98; 
     99; existence and protection of ${TROPFLUX_ID} 
     100IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 
     101   msg = 'eee : the directory' + iodirin  + ' is not accessible.' 
     102   ras = report(msg) 
     103   STOP 
     104ENDIF 
     105; 
     106; build mask filename 
     107filename_msk='mask_oaflux_30N30S.nc' 
     108; 
     109; check if this file exists 
     110fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 
     111IF fullfilename_msk[0] EQ '' THEN BEGIN 
     112   msg = 'eee : the file ' + fullfilename_msk + ' was not found.' 
     113   ras = report(msg) 
     114   STOP 
     115ENDIF 
    68116 
    69117filein='/Volumes/PRAVEEN/ERAI_global/20c3m_erai_t2_TROP_1989_2009.nc' 
    70 gridout='/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc' 
    71118 
    72119initncdf, filein 
    … …  
    81128mskin=glamt*0.+1. 
    82129 
    83 initncdf, gridout 
     130initncdf, fullfilename_msk 
    84131domdef 
    85132latout=reform(gphit(0,*)) & lonout=reform(glamt(*,0)) 
    86133print, 'lat grid ',min(latout),max(latout),latout(1)-latout(0) 
    87134print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0) 
    88 mskout=read_ncdf("msk", file=gridout,/nostr) 
     135mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 
    89136 
    90137help, t2min,lonin,latin,mskin,lonout,latout,mskout 
    … …  
    100147timein=timein & jptin=jpt 
    101148 
    102 initncdf, gridout 
     149initncdf, fullfilename_msk 
    103150help, t2mout 
    104151tt=timein 

  • trunk/src/interp_erai_ws_1989_2009.pro

    r12 r19  
    1616;           file_v10 [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/ERAI_global/20c3m_erai_v10_TROP_1989_2009.nc"]; 
    1717; 
    18 ;           mask [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc"]; 
     18;           mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/mask_oaflux_30N30S.nc"]; 
    1919; 
    2020;           ncfile [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc"]; 
    … …  
    3333; ======== 
    3434; 
     35; :ref:`tropflux_profile.sh` 
     36; 
    3537; EXAMPLES 
    3638; ======== 
    … …  
    4850; coding rules 
    4951; 
     52; KNOWN ISSUES 
     53; ============ 
     54; 
     55; test of existence of fullfilename_msk not very efficient because 
     56; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 
     57; 
    5058; EVOLUTIONS 
    5159; ========== 
     60; 
     61; - fplod 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin) 
     62; 
     63;   * remove hard coded directory for mask_oaflux_30N30S.nc 
    5264; 
    5365; - fplod 20101215T112657Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    6678pro interp_erai_ws_1989_2009 
    6779@common 
     80; 
     81; check for input directory  
     82; 
     83; test if ${TROPFLUX_ID} defined 
     84tropflux_id_env=GETENV('TROPFLUX_ID') 
     85CASE tropflux_id_env OF 
     86    ''  :  BEGIN 
     87     msg = 'eee : ${TROPFLUX_ID} is not defined' 
     88     ras = report(msg) 
     89     STOP 
     90           END 
     91 ELSE: BEGIN 
     92     msg = 'iii : ${TROPFLUX_ID} is ' + tropflux_id_env 
     93     ras = report(msg) 
     94       END 
     95ENDCASE 
     96; 
     97iodirin = isadirectory(tropflux_id_env) 
     98; 
     99; existence and protection of ${TROPFLUX_ID} 
     100IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 
     101   msg = 'eee : the directory' + iodirin  + ' is not accessible.' 
     102   ras = report(msg) 
     103   STOP 
     104ENDIF 
     105; 
     106; build mask filename 
     107filename_msk='mask_oaflux_30N30S.nc' 
     108; 
     109; check if this file exists 
     110fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 
     111IF fullfilename_msk[0] EQ '' THEN BEGIN 
     112   msg = 'eee : the file ' + fullfilename_msk + ' was not found.' 
     113   ras = report(msg) 
     114   STOP 
     115ENDIF 
    68116 
    69117st=19880101 & en=20100930 
    … …  
    77125v10in=read_ncdf("v10",st, en,file=filein,/nostr) 
    78126 
    79 gridout='/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc' 
    80  
    81127initncdf, filein 
    82128domdef 
    … …  
    89135mskin=glamt*0.+1. 
    90136 
    91 initncdf, gridout 
     137initncdf, fullfilename_msk 
    92138domdef 
    93139latout=reform(gphit(0,*)) & lonout=reform(glamt(*,0)) 
    94140print, 'lat grid ',min(latout),max(latout),latout(1)-latout(0) 
    95141print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0) 
    96 mskout=read_ncdf("msk", file=gridout,/nostr) 
     142mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 
    97143 
    98144help, u10in, v10in,lonin,latin,mskin,lonout,latout,mskout 
    … …  
    121167xlon=reform(glamt(*,0) ) & ylat=reform(gphit(0,*)) 
    122168 
    123 initncdf, gridout 
     169initncdf, fullfilename_msk 
    124170valmask=1.e20 
    125171 

  • trunk/src/interp_olr_30n30s_1989_2009.pro

    r12 r19  
    1515;           file_in [shape=ellipse,fontname=Courier,label="/Users/pkb/data/heat_budget/NOAA_OLR/new/olr.day.mean_new.nc"]; 
    1616; 
    17 ;           mask [shape=ellipse,fontname=Courier,label="/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc"]; 
     17;           mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/mask_oaflux_30N30S.nc"]; 
    1818; 
    1919;           ncfile [shape=ellipse,fontname=Courier,label="/Users/pkb/data/heat_budget/NOAA_OLR/newolr_oafluxgrid_30n30s_19890101_20091231.nc"]; 
    … …  
    3232; ======== 
    3333; 
     34; :ref:`tropflux_profile.sh` 
     35; 
    3436; EXAMPLES 
    3537; ======== 
    … …  
    4648; coding rules 
    4749; 
     50; KNOWN ISSUES 
     51; ============ 
     52; 
     53; test of existence of fullfilename_msk not very efficient because 
     54; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 
     55; 
    4856; EVOLUTIONS 
    4957; ========== 
     58; 
     59; - fplod 20101217T140745Z aedon.locean-ipsl.upmc.fr (Darwin) 
     60; 
     61;   * remove hard coded directory for mask_oaflux_30N30S.nc 
    5062; 
    5163; - fplod 20101215T113118Z aedon.locean-ipsl.upmc.fr (Darwin) 
    … …  
    6476pro interp_olr_30n30s_1989_2009 
    6577@common 
     78; 
     79; check for input directory  
     80; 
     81; test if ${TROPFLUX_ID} defined 
     82tropflux_id_env=GETENV('TROPFLUX_ID') 
     83CASE tropflux_id_env OF 
     84    ''  :  BEGIN 
     85     msg = 'eee : ${TROPFLUX_ID} is not defined' 
     86     ras = report(msg) 
     87     STOP 
     88           END 
     89 ELSE: BEGIN 
     90     msg = 'iii : ${TROPFLUX_ID} is ' + tropflux_id_env 
     91     ras = report(msg) 
     92       END 
     93ENDCASE 
     94; 
     95iodirin = isadirectory(tropflux_id_env) 
     96; 
     97; existence and protection of ${TROPFLUX_ID} 
     98IF (FILE_TEST(iodirin, /DIRECTORY, /EXECUTABLE, /READ) EQ 0) THEN BEGIN 
     99   msg = 'eee : the directory' + iodirin  + ' is not accessible.' 
     100   ras = report(msg) 
     101   STOP 
     102ENDIF 
     103; 
     104; build mask filename 
     105filename_msk='mask_oaflux_30N30S.nc' 
     106; 
     107; check if this file exists 
     108fullfilename_msk = isafile(iodirin + filename_msk, NEW=0, /MUST_EXIST) 
     109IF fullfilename_msk[0] EQ '' THEN BEGIN 
     110   msg = 'eee : the file ' + fullfilename_msk + ' was not found.' 
     111   ras = report(msg) 
     112   STOP 
     113ENDIF 
    66114 
    67115filein='/Users/pkb/data/heat_budget/NOAA_OLR/new/olr.day.mean_new.nc' 
    68 gridout='/Volumes/PRAVEEN/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc' 
    69116 
    70117initncdf, filein 
    … …  
    79126timein=time & jptin=jpt 
    80127 
    81 initncdf, gridout 
     128initncdf, fullfilename_msk 
    82129domdef 
    83130latout=reform(gphit(0,*)) & lonout=reform(glamt(*,0)) 
    84131print, 'lat grid ',min(latout),max(latout),latout(1)-latout(0) 
    85132print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0) 
    86 mskout=read_ncdf("msk", file=gridout,/nostr) 
     133mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 
    87134 
    88135help, olrin,lonin,latin,mskin,lonout,latout,mskout 

  • trunk/src/oaflux_mask_30N30S.pro

    r18 r19  
    4444; :ref:`tropflux_profile.sh` 
    4545; 
     46; :ref:`get_oaflux.sh` 
     47; 
    4648; :func:`report <saxo:report>` 
    4749; :func:`initncdf <saxo:initncdf>` 
    … …  
    4951; :ref:`ncdf_quickwrite <saxo:ncdf_quickwrite>` 
    5052; :func:`isadirectory <saxo:isadirectory>` 
     53; :func:`isafile <saxo:isafile>` 
     54; 
     55; :ref:`TropFlux_19890101_20091231.pro` 
     56; :ref:`interp_erai_dewt_1989_2009.pro` 
     57; :ref:`interp_erai_lwr_1989_2009.pro` 
     58; :ref:`interp_erai_sst_1989_2009.pro` 
     59; :ref:`interp_erai_t2m_1989_2009.pro` 
     60; :ref:`interp_erai_ws_1989_2009.pro` 
     61; :ref:`interp_olr_30n30s_1989_2009.pro` 
    5162; 
    5263; EXAMPLES 
    … …  
    7687; 
    7788; test of existence of fullfilename_in not very efficient because 
    78 ; MUST_EXIST keyword of :func:`isafile` not yet implemented 
     89; MUST_EXIST keyword of :func:`isafile <saxo:isafile>` not yet implemented 
    7990; 
    8091; EVOLUTIONS 

  • trunk/src/time_axe_modif_TropFlux.pro

    r12 r19  
    1616;           file_inout [shape=ellipse,fontname=Courier,label="/Users/pkb/data/TropFlux/TropFlux_19890101_20091231.nc"]; 
    1717; 
    18 ;           time_axe_modif_TropFlux  [shape=box, 
     18;           time_axe_modif_TropFlux [shape=box, 
    1919;           fontname=Courier, 
    2020;           color=blue, 
Note: See TracChangeset for help on using the changeset viewer.