Changeset 71 for trunk/src


Ignore:
Timestamp:
08/08/11 13:41:35 (13 years ago)
Author:
pinsard
Message:

try to interp ERA-I on climserv (not ok yet); minimize diffreneces between sources

Location:
trunk/src
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/interp_erai_dewt_1989_2009.pro

    r70 r71  
    2828;           rankdir="LR", 
    2929;           ] 
     30; 
    3031;           file_in [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/20c3m_erai_d2_TROP_1989_2009.nc"]; 
    3132;           mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/mask_oaflux_30N30S.nc"]; 
    32 ; 
    3333;           file_out [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_d2m_19890101_20091231_oafluxgrid.nc"]; 
     34; 
    3435;           interp_erai_dewt_1989_2009 [shape=box, 
    3536;           fontname=Courier, 
    … …  
    4041;           {file_in mask} -> {interp_erai_dewt_1989_2009} -> {file_out} 
    4142; 
    42 ;       } 
     43;        } 
    4344; 
    4445; SEE ALSO 
    … …  
    7172; make it work : pb on loholt1:: 
    7273; 
    73 ;   Variable is undefined: OUTMASK_IND. 
     74;   d2m[longitude,latitude,time]=d2mout:d2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr  @ globattr 
     75;   % NCDF_DIMDEF: 0 is not a valid cdfid. 
     76;   % Stop encountered: INTERP_ERAI_DEWT_1989_2009    1 
     77;      /home/pinsard/tropflux_ws/src/interp_erai_dewt_1989_2009.pro 
    7478; 
    7579; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 
    … …  
    162166; build data filename 
    163167filename='20c3m_erai_d2_TROP_1989_2009.nc' 
     168; 
     169; check if this file exists 
    164170fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST) 
    165171IF fullfilename[0] EQ '' THEN BEGIN 
    … …  
    168174   STOP 
    169175ENDIF 
    170 ; 
    171176; 
    172177; test if ${PROJECT_OD} defined 
    … …  
    234239endfor 
    235240 
    236  
    237241timein=timein & jptin=jpt 
    238242 
    … …  
    255259globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    256260 
    257  
    258261ncfields = 'd2m[longitude,latitude,time]=d2mout:d2m_attr; ' $ 
    259262                      + 'longitude[]=lon:lon_attr; ' $ 
    … …  
    264267@ncdf_quickwrite 
    265268 
    266  
    267269end 

  • trunk/src/interp_erai_lwr_1989_2009.pro

    r70 r71  
    2828;           rankdir="LR", 
    2929;           ] 
     30; 
     31;           file_in [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/20c3m_erai_str_TROP_1989_2009.nc"]; 
    3032;           mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/mask_oaflux_30N30S.nc"]; 
    31 ; 
    32 ;           file_in [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/20c3m_erai_str_TROP_1989_2009.nc"]; 
    33 ; 
    34 ;           ncfile [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_lwr_19890102_20091231_oafluxgrid.nc"]; 
     33;           file_out [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_lwr_19890102_20091231_oafluxgrid.nc"]; 
    3534; 
    3635;           interp_erai_lwr_1989_2009 [shape=box, 
    … …  
    4039;           label="${PROJECT}/src/interp_erai_lwr_1989_2009.pro"]; 
    4140; 
    42 ;           {mask file_in} -> {interp_erai_lwr_1989_2009} -> {ncfile} 
     41;           {file_in mask} -> {interp_erai_lwr_1989_2009} -> {file_out} 
    4342; 
    4443;         } 
    … …  
    7372; make it work : pb on loholt1:: 
    7473; 
    75 ;   Variable is undefined: OUTMASK_IND. 
     74;  OUTMASK_IND     UNDEFINED = <Undefined> 
     75;  SET_OUTMSKVAL   UNDEFINED = <Undefined> 
     76;   
     77;  lwr[longitude,latitude,time]=lwrout:lwr_attr; longitude[]=xlon:lon_attr; latitude[]=ylat:lat_attr; tt[*time]=tt:time_attr  @ globattr 
     78;  % NCDF_DIMDEF: 0 is not a valid cdfid. 
     79;  % Stop encountered: INTERP_ERAI_LWR_1989_2009    1 
     80;     /home/pinsard/tropflux_ws/src/interp_erai_lwr_1989_2009.pro 
    7681; 
    7782; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 
    … …  
    178183; check if this file exists 
    179184fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST) 
    180 IF fullfilename_msk[0] EQ '' THEN BEGIN 
    181    msg = 'eee : the file ' + fullfilename_msk + ' was not found.' 
     185IF fullfilename[0] EQ '' THEN BEGIN 
     186   msg = 'eee : the file ' + fullfilename + ' was not found.' 
    182187   ras = report(msg) 
    183188   STOP 
    … …  
    238243mskout=read_ncdf("msk", file=fullfilename_msk,/nostr) 
    239244 
    240  
    241245si=size(lwrin) 
    242246lwrout=fltarr(jpi,jpj,jptin) 
    … …  
    247251  tab=reform(lwrin(*,*,jt)) 
    248252  lwrout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $ 
    249      , lonout,latout,method='bilinear' $ 
    250      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout) 
     253      , lonout,latout,method='bilinear' $ 
     254      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout) 
    251255  help, OUTMASK_IND 
    252256  help, SET_OUTMSKVAL 

  • trunk/src/interp_erai_sst_1989_2009.pro

    r70 r71  
    1515; by :ref:`compute_erai_daily_region_2d.sh`. 
    1616; 
    17 ; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` containing OAFLUX grid have been p roduced by :ref:`oaflux_mask_30N30S.pro`. 
     17; :file:`${PROJECT_ID}/mask_oaflux_30N30S.nc` containing OAFLUX grid have been produced by :ref:`oaflux_mask_30N30S.pro`. 
    1818; 
    1919; Interpolated sst is written in 
    … …  
    2929;           rankdir="LR", 
    3030;           ] 
     31; 
    3132;           file_in [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/20c3m_erai_sstk_TROP_1989_2009.nc"]; 
    32 ; 
    3333;           mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/mask_oaflux_30N30S.nc"]; 
    34 ; 
    35 ;           ncfile [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc"]; 
     34;           file_out [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_sst_19890101_20091231_oafluxgrid.nc"]; 
    3635; 
    3736;           interp_erai_sst_1989_2009 [shape=box, 
    … …  
    4140;           label="${PROJECT}/src/interp_erai_sst_1989_2009.pro"]; 
    4241; 
    43 ;           {file_in mask} -> {interp_erai_sst_1989_2009} -> {ncfile} 
     42;           {file_in mask} -> {interp_erai_sst_1989_2009} -> {file_out} 
     43; 
    4444;          } 
    4545; 
    … …  
    7373; make it work : pb on loholt1:: 
    7474; 
    75 ;   Variable is undefined: OUTMASK_IND. 
     75;  sst[longitude,latitude,time]=sstkout:sstk_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[*time]=timein:time_attr  @ globattr 
     76;  % NCDF_DIMDEF: 0 is not a valid cdfid. 
     77;  % Stop encountered: INTERP_ERAI_SST_1989_2009    1 
     78;     /home/pinsard/tropflux_ws/src/interp_erai_sst_1989_2009.pro 
    7679; 
    7780; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 
    … …  
    167170; check if this file exists 
    168171fullfilename = isafile(iodirin + filename, NEW=0, /MUST_EXIST) 
    169 IF fullfilename_msk[0] EQ '' THEN BEGIN 
    170    msg = 'eee : the file ' + fullfilename_msk + ' was not found.' 
     172IF fullfilename[0] EQ '' THEN BEGIN 
     173   msg = 'eee : the file ' + fullfilename + ' was not found.' 
    171174   ras = report(msg) 
    172175   STOP 
    … …  
    233236  tab=reform(sstkin(*,*,jt)) 
    234237  sstkout(*,*,jt)=call_interp2d(tab,lonin,latin,mskin $ 
    235      , lonout,latout,method='bilinear' $ 
    236      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout) 
     238      , lonout,latout,method='bilinear' $ 
     239      , OUTMASK_IND=mskout, SET_OUTMSKVAL=mskout) 
    237240  sstkout(*,*,jt)=sstkout(*,*,jt)*mskout+(1.-mskout)*1.e20 
    238241endfor 
    239  
    240242 
    241243timein=timein & jptin=jpt 
    … …  
    257259globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    258260 
    259  
    260261ncfields = 'sst[longitude,latitude,time]=sstkout:sstk_attr; ' $ 
    261262                      + 'longitude[]=lon:lon_attr; ' $ 
    … …  
    266267@ncdf_quickwrite 
    267268 
    268  
    269269end 

  • trunk/src/interp_erai_t2m_1989_2009.pro

    r70 r71  
    2828;           rankdir="TB", 
    2929;           ] 
     30; 
    3031;           file_in [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/20c3m_erai_t2_TROP_1989_2009.nc"]; 
    3132;           mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/mask_oaflux_30N30S.nc"]; 
    32 ; 
    33 ;           ncfile [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 
     33;           file_out [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 
    3434; 
    3535;           interp_erai_t2m_1989_2009 [shape=box, 
    … …  
    3939;           label="${PROJECT}/src/interp_erai_t2m_1989_2009.pro"]; 
    4040; 
    41 ;           {file_in mask} -> {interp_erai_t2m_1989_2009} -> {ncfile} 
     41;           {file_in mask} -> {interp_erai_t2m_1989_2009} -> {file_out} 
    4242; 
    4343;          } 
    44 ; 
    4544; 
    4645; SEE ALSO 
    … …  
    7877; TODO 
    7978; ==== 
     79; 
     80; make it work : pb on loholt1:: 
     81; 
     82;  t2m[longitude,latitude,*time]=t2mout:t2m_attr; longitude[]=lon:lon_attr; latitude[]=lat:lat_attr; time[]=timein:time_attr  @globattr 
     83;  % NCDF_DIMDEF: 0 is not a valid cdfid. 
     84;  % Stop encountered: INTERP_ERAI_T2M_1989_2009    1 
     85;     /home/pinsard/tropflux_ws/src/interp_erai_t2m_1989_2009.pro 
    8086; 
    8187; strange view (lat and lon shift with ncview) : check grid init 
    … …  
    334340globattr={source:'Data are from ECMWF ERA-Interim reanalysis', timerange:cda0+' - '+cda1} 
    335341 
    336  
    337342help, t2mout 
    338343help, timein 
    … …  
    386391;++@ncdf_quickwrite 
    387392 
    388  
    389393end 

  • trunk/src/interp_erai_ws_1989_2009.pro

    r70 r71  
    3030;           rankdir="LR", 
    3131;           ] 
     32; 
    3233;           file_u10 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/20c3m_erai_u10_TROP_1989_2009.nc"]; 
    3334;           file_v10 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/20c3m_erai_v10_TROP_1989_2009.nc"]; 
    34 ; 
    3535;           mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/mask_oaflux_30N30S.nc"]; 
    36 ; 
    37 ;           ncfile [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_ws_19890101_20091231_oafluxgrid.nc"]; 
     36;           file_out [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_ws_19890101_20091231_oafluxgrid.nc"]; 
    3837; 
    3938;           interp_erai_ws_1989_2009 [shape=box, 
    … …  
    4342;           label="${PROJECT}/src/interp_erai_ws_1989_2009.pro"]; 
    4443; 
    45 ;           {file_u10 file_v10 mask} -> {interp_erai_ws_1989_2009} -> {ncfile} 
    46 ; 
    47 ;        } 
     44;           {file_u10 file_v10 mask} -> {interp_erai_ws_1989_2009} -> {file_out} 
     45; 
     46;         } 
    4847; 
    4948; SEE ALSO 
    … …  
    7675; make it work : pb on loholt1:: 
    7776; 
    78 ;   Variable is undefined: OUTMASK_IND. 
     77;  Writing fields:  
     78;  u10[longitude,latitude,time]=u10out:u10_attr; v10[longitude,latitude,time]=v10out:v10_attr; longitude[]=xlon:lon_attr; latitude[]=ylat:lat_attr; tt[*time]=tt:time_attr  @ globattr 
     79;  % NCDF_DIMDEF: 0 is not a valid cdfid. 
     80;  % Stop encountered: INTERP_ERAI_WS_1989_2009    1 
     81;     /home/pinsard/tropflux_ws/src/interp_erai_ws_1989_2009.pro 
    7982; 
    8083; check OUTMASK_IND and SET_OUTMSKVAL added to call_interp2d call 
Note: See TracChangeset for help on using the changeset viewer.