- Timestamp:
- 04/17/12 10:27:26 (12 years ago)
- Location:
- trunk/src
- Files:
-
- 14 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/amsu2ncdf.pro
r523 r525 28 28 ; :raise dd: required 29 29 ; 30 ; :param lon _min: longitude min, W < 031 ; :units lon _min: deg32 ; :type lon _min: double33 ; :raise lon _min: required34 ; 35 ; :param lon _max: longitude max, W < 036 ; :units lon _max: deg37 ; :type lon _max: double38 ; :raise lon _max: required39 ; 40 ; :param lat _min: latitude min, N > 041 ; :units lat _min: deg42 ; :type lat _min: double43 ; :raise lat _min: required44 ; 45 ; :param lat _max: latitude max, N > 046 ; :units lat _max: deg47 ; :type lat _max: double48 ; :raise lat _max: required30 ; :param lonmin: longitude min, W < 0 31 ; :units lonmin: deg 32 ; :type lonmin: double 33 ; :raise lonmin: required 34 ; 35 ; :param lonmax: longitude max, W < 0 36 ; :units lonmax: deg 37 ; :type lonmax: double 38 ; :raise lonmax: required 39 ; 40 ; :param latmin: latitude min, N > 0 41 ; :units latmin: deg 42 ; :type latmin: double 43 ; :raise latmin: required 44 ; 45 ; :param latmax: latitude max, N > 0 46 ; :units latmax: deg 47 ; :type latmax: double 48 ; :raise latmax: required 49 49 ; 50 50 ; read :file:`${PROJECT_ID}/AMSU/yyyy/mm/dd/amsua[_amsub]_{yyyymmdd}_{geomin}_{geomax}_nadir.dat` … … 89 89 ; mm=8 90 90 ; dd=01 91 ; lon _min=-60.92 ; lon _max=50.93 ; lat _min=-30.94 ; lat _max=45.95 ; amsu2ncdf, use_amsua, yyyy, mm, dd, lon _min, lon_max, lat_min, lat_max91 ; lonmin=-60. 92 ; lonmax=50. 93 ; latmin=-30. 94 ; latmax=45. 95 ; amsu2ncdf, use_amsua, yyyy, mm, dd, lonmin, lonmax, latmin, latmax 96 96 ; 97 97 ; or:: … … 101 101 ; mm=8 102 102 ; dd=01 103 ; lon _min=-60.104 ; lon _max=50.105 ; lat _min=-30.106 ; lat _max=45.107 ; amsu2ncdf, use_amsua, yyyy, mm, dd, lon _min, lon_max, lat_min, lat_max103 ; lonmin=-60. 104 ; lonmax=50. 105 ; latmin=-30. 106 ; latmax=45. 107 ; amsu2ncdf, use_amsua, yyyy, mm, dd, lonmin, lonmax, latmin, latmax 108 108 ; 109 109 ; … … 160 160 ; 161 161 ; * complete description with parameters 162 ; * add lon _min, lon_max, lat_min, lat_max parameter162 ; * add lonmin, lonmax, latmin, latmax parameter 163 163 ; * align parameters order to correct_nadir_amsu 164 164 ; * change output filename … … 212 212 ; 213 213 ;- 214 PRO amsu2ncdf, use_amsua, yyyy, mm, dd, lon _min, lon_max, lat_min, lat_max214 PRO amsu2ncdf, use_amsua, yyyy, mm, dd, lonmin, lonmax, latmin, latmax 215 215 ; 216 216 compile_opt idl2, strictarrsubs … … 231 231 look='filename' 232 232 scale = 1. 233 geomin=geolocation_to_string_idl(lon _min,lat_min, look, scale)234 geomax=geolocation_to_string_idl(lon _max,lat_max, look, scale)233 geomin=geolocation_to_string_idl(lonmin,latmin, look, scale) 234 geomax=geolocation_to_string_idl(lonmax,latmax, look, scale) 235 235 filename=project_id_env+ 'AMSU/' $ 236 236 + string(yyyy,format='(I4.4)') + '/' $ -
trunk/src/cal_hcl.m
r368 r525 1 function [tps,alt,lat,lon,tpt,ep_cl]=cal_hcl(tps_min,tps_max,lat _min,lat_max,lon_min,lon_max)1 function [tps,alt,lat,lon,tpt,ep_cl]=cal_hcl(tps_min,tps_max,latmin,latmax,lonmin,lonmax) 2 2 3 3 %CAL_HCL Calcul de la hauteur de la couche limite ECMWF opera … … 23 23 % tps_max 24 24 % borne max pour la date au format matlab 25 % lat _min25 % latmin 26 26 % borne min pour la latitude 27 % lat _max27 % latmax 28 28 % borne max pour la latitude 29 % lon _min29 % lonmin 30 30 % borne min pour la longitude 31 % lon _max31 % lonmax 32 32 % borne max pour la longitude 33 33 % … … 54 54 % >> tps_min=datenum([2006 07 01]) 55 55 % >> tps_max=datenum([2006 07 10]); 56 % >> lat _min=-5; lat_max=5;57 % >> lon _min=-10; lon_max=10;58 % >> [tps,alt,lat,lon,tpt,ep_cl]=cal_hcl(tps_min,tps_max,lat _min,lat_max,lon_min,lon_max);56 % >> latmin=-5; latmax=5; 57 % >> lonmin=-10; lonmax=10; 58 % >> [tps,alt,lat,lon,tpt,ep_cl]=cal_hcl(tps_min,tps_max,latmin,latmax,lonmin,lonmax); 59 59 % >> moytps_epcl=squeeze(nanmean(ep_cl,1)); 60 60 % >> % Graphique … … 101 101 % 102 102 %- 103 usage='[tps,alt,lat,lon,tpt,ep_cl]=cal_hcl(tps_min,tps_max,lat _min,lat_max,lon_min,lon_max)';103 usage='[tps,alt,lat,lon,tpt,ep_cl]=cal_hcl(tps_min,tps_max,latmin,latmax,lonmin,lonmax)'; 104 104 if nargin~=nargin(mfilename) 105 105 disp(['Incorrect number of input arguments ' num2str(nargin)]); … … 122 122 if m_min == m_max 123 123 124 [tps,level,lat,lon,ta]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/ta.2006%02d.amh.WESTAFR_025.nc',m_min),{'ta'},'time',tps_min_mod,tps_max_mod,'level',61.,91.,'lat',lat _min,lat_max,'lon',lon_min,lon_max);125 [tps,level,lat,lon,q]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/q.2006%02d.amh.WESTAFR_025.nc',m_min),{'q'},'time',tps_min_mod,tps_max_mod,'level',61.,91.,'lat',lat _min,lat_max,'lon',lon_min,lon_max);126 [tps,lat,lon,msl]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_SF/2006/msl.2006%02d.ash.WESTAFR_025.nc',m_min),{'msl'},'time',tps_min_mod,tps_max_mod,'lat',lat _min,lat_max,'lon',lon_min,lon_max);124 [tps,level,lat,lon,ta]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/ta.2006%02d.amh.WESTAFR_025.nc',m_min),{'ta'},'time',tps_min_mod,tps_max_mod,'level',61.,91.,'lat',latmin,latmax,'lon',lonmin,lonmax); 125 [tps,level,lat,lon,q]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/q.2006%02d.amh.WESTAFR_025.nc',m_min),{'q'},'time',tps_min_mod,tps_max_mod,'level',61.,91.,'lat',latmin,latmax,'lon',lonmin,lonmax); 126 [tps,lat,lon,msl]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_SF/2006/msl.2006%02d.ash.WESTAFR_025.nc',m_min),{'msl'},'time',tps_min_mod,tps_max_mod,'lat',latmin,latmax,'lon',lonmin,lonmax); 127 127 128 128 else 129 129 130 [tps,level,lat,lon,ta]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/ta.2006%02d.amh.WESTAFR_025.nc',m_min),{'ta'},'time',tps_min_mod,NaN,'level',61.,91.,'lat',lat _min,lat_max,'lon',lon_min,lon_max);131 [tps,level,lat,lon,q]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/q.2006%02d.amh.WESTAFR_025.nc',m_min),{'q'},'time',tps_min_mod,NaN,'level',61.,91.,'lat',lat _min,lat_max,'lon',lon_min,lon_max);132 [tps,lat,lon,msl]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_SF/2006/msl.2006%02d.ash.WESTAFR_025.nc',m_min),{'msl'},'time',tps_min_mod,NaN,'lat',lat _min,lat_max,'lon',lon_min,lon_max);130 [tps,level,lat,lon,ta]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/ta.2006%02d.amh.WESTAFR_025.nc',m_min),{'ta'},'time',tps_min_mod,NaN,'level',61.,91.,'lat',latmin,latmax,'lon',lonmin,lonmax); 131 [tps,level,lat,lon,q]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/q.2006%02d.amh.WESTAFR_025.nc',m_min),{'q'},'time',tps_min_mod,NaN,'level',61.,91.,'lat',latmin,latmax,'lon',lonmin,lonmax); 132 [tps,lat,lon,msl]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_SF/2006/msl.2006%02d.ash.WESTAFR_025.nc',m_min),{'msl'},'time',tps_min_mod,NaN,'lat',latmin,latmax,'lon',lonmin,lonmax); 133 133 134 134 i=1; 135 135 while m_min+i ~= m_max 136 136 137 [tps_tmp,level,lat,lon,ta_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/ta.2006%02d.amh.WESTAFR_025.nc',m_min+i),{'ta'},'time',NaN,NaN,'level',61.,91.,'lat',lat _min,lat_max,'lon',lon_min,lon_max);138 [tps_tmp,level,lat,lon,q_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/q.2006%02d.amh.WESTAFR_025.nc',m_min+i),{'q'},'time',NaN,NaN,'level',61.,91.,'lat',lat _min,lat_max,'lon',lon_min,lon_max);139 [tps_tmp,lat,lon,msl_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_SF/2006/msl.2006%02d.ash.WESTAFR_025.nc',m_min+i),{'msl'},'time',NaN,NaN,'lat',lat _min,lat_max,'lon',lon_min,lon_max);137 [tps_tmp,level,lat,lon,ta_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/ta.2006%02d.amh.WESTAFR_025.nc',m_min+i),{'ta'},'time',NaN,NaN,'level',61.,91.,'lat',latmin,latmax,'lon',lonmin,lonmax); 138 [tps_tmp,level,lat,lon,q_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/q.2006%02d.amh.WESTAFR_025.nc',m_min+i),{'q'},'time',NaN,NaN,'level',61.,91.,'lat',latmin,latmax,'lon',lonmin,lonmax); 139 [tps_tmp,lat,lon,msl_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_SF/2006/msl.2006%02d.ash.WESTAFR_025.nc',m_min+i),{'msl'},'time',NaN,NaN,'lat',latmin,latmax,'lon',lonmin,lonmax); 140 140 141 141 tps=[tps;tps_tmp]; … … 147 147 end 148 148 149 [tps_tmp,level,lat,lon,ta_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/ta.2006%02d.amh.WESTAFR_025.nc',m_min+i),{'ta'},'time',NaN,tps_max_mod,'level',61.,91.,'lat',lat _min,lat_max,'lon',lon_min,lon_max);150 [tps_tmp,level,lat,lon,q_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/q.2006%02d.amh.WESTAFR_025.nc',m_min+i),{'q'},'time',NaN,tps_max_mod,'level',61.,91.,'lat',lat _min,lat_max,'lon',lon_min,lon_max);151 [tps_tmp,lat,lon,msl_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_SF/2006/msl.2006%02d.ash.WESTAFR_025.nc',m_min+i),{'msl'},'time',NaN,tps_max_mod,'lat',lat _min,lat_max,'lon',lon_min,lon_max);149 [tps_tmp,level,lat,lon,ta_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/ta.2006%02d.amh.WESTAFR_025.nc',m_min+i),{'ta'},'time',NaN,tps_max_mod,'level',61.,91.,'lat',latmin,latmax,'lon',lonmin,lonmax); 150 [tps_tmp,level,lat,lon,q_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_ML/2006/q.2006%02d.amh.WESTAFR_025.nc',m_min+i),{'q'},'time',NaN,tps_max_mod,'level',61.,91.,'lat',latmin,latmax,'lon',lonmin,lonmax); 151 [tps_tmp,lat,lon,msl_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_SF/2006/msl.2006%02d.ash.WESTAFR_025.nc',m_min+i),{'msl'},'time',NaN,tps_max_mod,'lat',latmin,latmax,'lon',lonmin,lonmax); 152 152 153 153 -
trunk/src/correct_nadir_amsu-v2.pro
r513 r525 60 60 ; canal (string) 61 61 ; resol 62 ; lon _min63 ; lon _max64 ; lat _min65 ; lat _max62 ; lonmin 63 ; lonmax 64 ; latmin 65 ; latmax 66 66 ; 67 67 ; bornes domaines necessaires pour le masque terre-mer … … 97 97 ; canal='a5' 98 98 ; resol=1 99 ; lon _min=-60.100 ; lon _max=50.101 ; lat _min=-30.102 ; lat _max=45.103 ; correct_nadir_amsu, canal,resol, lon _min, lon_max, lat_min, lat_max99 ; lonmin=-60. 100 ; lonmax=50. 101 ; latmin=-30. 102 ; latmax=45. 103 ; correct_nadir_amsu, canal,resol, lonmin, lonmax, latmin, latmax 104 104 ; 105 105 ; :History: … … 176 176 ; 177 177 ;- 178 PRO correct_nadir_amsu, canal, resol, lon _min, lon_max, lat_min, lat_max178 PRO correct_nadir_amsu, canal, resol, lonmin, lonmax, latmin, latmax 179 179 ; 180 180 compile_opt idl2, strictarrsubs … … 190 190 ; lecture fichier land - sea (S. Masson) 191 191 file=project_id_env+'/MASK/ETOPO1_Ice_g_gmt4.nc' 192 domdef,lon _min, lon_max, lat_min, lat_max192 domdef,lonmin, lonmax, latmin, latmax 193 193 194 194 bate = ncdf_lec(file, var = 'z') GT 0 … … 199 199 200 200 ; extraction zone donnees AMSU 201 domdef,lon _min,lon_max,lat_min,lat_max201 domdef,lonmin,lonmax,latmin,latmax 202 202 bate = ncdf_lec(file, var = 'z') GT 0 203 203 xxe = float(ncdf_lec(file, var = 'lon')) … … 286 286 if (ifov eq nbpix) do begin 287 287 interpolswath,tb,lat,lon,resol,nbgrid,tbgrid,latgrid,longrid 288 zone=where((longrid ge lon _min) and (longrid le lon_max) $289 and (latgrid ge lat _min) and (latgrid le lat_max),npt)288 zone=where((longrid ge lonmin) and (longrid le lonmax) $ 289 and (latgrid ge latmin) and (latgrid le latmax),npt) 290 290 if npt ne 0 then begin 291 291 temps[zone]=hour[i] -
trunk/src/correct_nadir_amsu.pro
r513 r525 59 59 ; yyyy 60 60 ; mm 61 ; lon _min62 ; lon _max63 ; lat _min64 ; lat _max61 ; lonmin 62 ; lonmax 63 ; latmin 64 ; latmax 65 65 ; 66 66 ; bornes domaines necessaires pour le masque terre-mer … … 98 98 ; mm=8 99 99 ; dd=1 100 ; lon _min=-60.101 ; lon _max=50.102 ; lat _min=-30.103 ; lat _max=45.104 ; correct_nadir_amsu, use_amsua, yyyy, mm, dd, lon _min, lon_max, lat_min, lat_max100 ; lonmin=-60. 101 ; lonmax=50. 102 ; latmin=-30. 103 ; latmax=45. 104 ; correct_nadir_amsu, use_amsua, yyyy, mm, dd, lonmin, lonmax, latmin, latmax 105 105 ; 106 106 ; Using AMSU-B only :: … … 110 110 ; mm=8 111 111 ; dd=1 112 ; lon _min=-60.113 ; lon _max=50.114 ; lat _min=-30.115 ; lat _max=45.116 ; correct_nadir_amsu, use_amsua, yyyy, mm, dd, lon _min, lon_max, lat_min, lat_max112 ; lonmin=-60. 113 ; lonmax=50. 114 ; latmin=-30. 115 ; latmax=45. 116 ; correct_nadir_amsu, use_amsua, yyyy, mm, dd, lonmin, lonmax, latmin, latmax 117 117 ; 118 118 ; :History: … … 182 182 ; 183 183 ;- 184 PRO correct_nadir_amsu, use_amsua, yyyy, mm, dd, lon _min, lon_max, lat_min, lat_max184 PRO correct_nadir_amsu, use_amsua, yyyy, mm, dd, lonmin, lonmax, latmin, latmax 185 185 ; 186 186 compile_opt idl2, strictarrsubs … … 304 304 ; fichier land - sea (S. Masson) 305 305 file=project_id_env+'/MASK/ETOPO1_Ice_g_gmt4.nc' 306 domdef,lon _min, lon_max, lat_min, lat_max306 domdef,lonmin, lonmax, latmin, latmax 307 307 308 308 bate = ncdf_lec(file, var = 'z') GT 0 … … 313 313 314 314 ; extraction zone donnees AMSU :passage en parametres 315 domdef,lon _min,lon_max,lat_min,lat_max315 domdef,lonmin,lonmax,latmin,latmax 316 316 bate = ncdf_lec(file, var = 'z') GT 0 317 317 xxe = float(ncdf_lec(file, var = 'lon')) … … 333 333 look = 'filename' 334 334 scale = 1. 335 geomin = geolocation_to_string_idl(lon _min, lat_min, look, scale)336 geomax = geolocation_to_string_idl(lon _max, lat_max, look, scale)335 geomin = geolocation_to_string_idl(lonmin, latmin, look, scale) 336 geomax = geolocation_to_string_idl(lonmax, latmax, look, scale) 337 337 fichier1=project_id_env+ 'AMSU/' $ 338 338 + string(yyyy,format='(I4.4)') + '/' $ -
trunk/src/cresamsu.pro
r519 r525 31 31 ; ddb=1 32 32 ; dde=3 33 ; lon _min=-25.34 ; lon _max=25.35 ; lat _min=-5.36 ; lat _max=20.33 ; lonmin=-25. 34 ; lonmax=25. 35 ; latmin=-5. 36 ; latmax=20. 37 37 ; temp='m' 38 38 ; rtemp=0.5 … … 40 40 ; cresamsu, numch, yyyyb, mmb, ddb $ 41 41 ; , yyyye, mme, dde $ 42 ; , lon _min, lon_max, lat_min, lat_max,temp,rtemp,rhdeg42 ; , lonmin, lonmax, latmin, latmax,temp,rtemp,rhdeg 43 43 ; 44 44 ; SEE ALSO … … 130 130 ;- 131 131 132 pro cresamsu, numch, yyyyb, mmb, ddb, yyyye,mme,dde,lon _min,lon_max,lat_min,lat_max,temp,rtemp,rhdeg132 pro cresamsu, numch, yyyyb, mmb, ddb, yyyye,mme,dde,lonmin,lonmax,latmin,latmax,temp,rtemp,rhdeg 133 133 134 134 ; les parametres de l'interpolation … … 197 197 198 198 ; grille des donnees en sortie 199 lo1=lon _min+rhdeg200 lo2=lon _max-rhdeg201 la1=lat _min+rhdeg202 la2=lat _max-rhdeg199 lo1=lonmin+rhdeg 200 lo2=lonmax-rhdeg 201 la1=latmin+rhdeg 202 la2=latmax-rhdeg 203 203 print,"bornes de la grille de sortie",lo1,lo2,la1,la2,jdeb,jfin 204 204 kmax=fix((lo2-lo1)/pash) ; pour valeurs de 0 a kmax … … 241 241 cont2=120L ; ????? 242 242 look = 'filename' 243 geomin = geolocation_to_string_idl(lon _min, lat_min, look,1)244 geomax = geolocation_to_string_idl(lon _max, lat_max, look,1)243 geomin = geolocation_to_string_idl(lonmin, latmin, look,1) 244 geomax = geolocation_to_string_idl(lonmax, latmax, look,1) 245 245 246 246 testfilename='' … … 262 262 ; lecture du fichier journalier 263 263 result=file_amsu_t2_to_mem( yyyy, mm, dd, numch $ 264 , lon _min, lon_max, lat_min, lat_max $264 , lonmin, lonmax, latmin, latmax $ 265 265 , testfilename) 266 266 if (size(result,/TYPE) EQ 3) THEN BEGIN -
trunk/src/cresamsu_cdf.pro
r503 r525 21 21 ; 22 22 ; :: 23 ; cresamsu_cdf,use_amsua,yyyyb,mmb,ddb,yyyye,mme,dde,lon _min,lon_max, lat_min, lat_max,listcanal23 ; cresamsu_cdf,use_amsua,yyyyb,mmb,ddb,yyyye,mme,dde,lonmin,lonmax, latmin, latmax,listcanal 24 24 ; 25 25 ; :Examples: … … 35 35 ; ddb=1 36 36 ; dde=3 37 ; lon _min=-2038 ; lon _max=4039 ; lat _min=040 ; lat _max=3041 ; cresamsu_cdf,use_amsua,yyyyb,mmb,ddb,yyyye,mme,dde,lon _min,lon_max, lat_min, lat_max,listcanal37 ; lonmin=-20 38 ; lonmax=40 39 ; latmin=0 40 ; latmax=30 41 ; cresamsu_cdf,use_amsua,yyyyb,mmb,ddb,yyyye,mme,dde,lonmin,lonmax, latmin, latmax,listcanal 42 42 ; 43 43 ; :Restrictions: … … 122 122 ;- 123 123 124 pro cresamsu_cdf,use_amsua,yyyyb,mmb,ddb,yyyye,mme,dde,lon _min, lon_max, lat_min, lat_max,listcanal124 pro cresamsu_cdf,use_amsua,yyyyb,mmb,ddb,yyyye,mme,dde,lonmin, lonmax, latmin, latmax,listcanal 125 125 126 126 ; les parametres de l'interpolation … … 177 177 178 178 ; grille des donnees en sortie 179 lo1=lon _min+rhdeg180 lo2=lon _max-rhdeg181 la1=lat _min+rhdeg182 la2=lat _max-rhdeg179 lo1=lonmin+rhdeg 180 lo2=lonmax-rhdeg 181 la1=latmin+rhdeg 182 la2=latmax-rhdeg 183 183 print,"bornes de la grille de sortie",lo1,lo2,la1,la2,jdeb,jfin 184 184 kmax=fix((lo2-lo1)/pash) ; pour valeurs de 0 a kmax … … 238 238 cont2=120L 239 239 look = 'filename' 240 geomin = geolocation_to_string_idl(lon _min, lat_min, look,1)241 geomax = geolocation_to_string_idl(lon _max, lat_max, look,1)240 geomin = geolocation_to_string_idl(lonmin, latmin, look,1) 241 geomax = geolocation_to_string_idl(lonmax, latmax, look,1) 242 242 243 243 for yyyy=yyyyb,yyyye do begin -
trunk/src/extract_amsua.pro
r521 r525 61 61 ; :param mm: 62 62 ; :param resol: 63 ; :param lon _min:64 ; :param lon _max:65 ; :param lat _min:66 ; :param lat _max:63 ; :param lonmin: 64 ; :param lonmax: 65 ; :param latmin: 66 ; :param latmax: 67 67 ; 68 68 ; EXAMPLES … … 94 94 ; dd=13 95 95 ; resol=1 96 ; lon _min=-25.97 ; lon _max=25.98 ; lat _min=-5.99 ; lat _max=20.96 ; lonmin=-25. 97 ; lonmax=25. 98 ; latmin=-5. 99 ; latmax=20. 100 100 ; files_list = search_amsufiles(numch, yyyy, mm, dd) 101 ; extract_amsua, numch, files_list, yyyy, mm, dd, resol, lon _min, lon_max, lat_min, lat_max101 ; extract_amsua, numch, files_list, yyyy, mm, dd, resol, lonmin, lonmax, latmin, latmax 102 102 ; 103 103 ; :file:`${PROJECT_ID}/AMSU/2006/08/a5_20060813_025w05s_025e20n.dat` … … 275 275 ;- 276 276 PRO extract_amsua, numch, files_list, yyyy, mm, dd, resol $ 277 , lon _min, lon_max, lat_min, lat_max277 , lonmin, lonmax, latmin, latmax 278 278 ; 279 279 compile_opt idl2, strictarrsubs … … 355 355 amzeni=REFORM(ama_scan.angles[0,*,*]/100.) 356 356 ttt=REFORM(ama_scan.scnlintime/3600000.) 357 jnd=where(amalong[midpix,*] gt lon _min-15 and amalong[midpix,*] lt lon_max+15 and amalati[midpix,*] gt lat_min-15 and amalati[midpix,*] lt lat_max+15 ,nzon)357 jnd=where(amalong[midpix,*] gt lonmin-15 and amalong[midpix,*] lt lonmax+15 and amalati[midpix,*] gt latmin-15 and amalati[midpix,*] lt latmax+15 ,nzon) 358 358 print,"iii : nb de points du fichier dans le domaine geographique +15deg ",nzon, SYSTIME() 359 359 if nzon gt 1 then begin … … 444 444 ;print,'wwwtest: nb points fov interpole ',n_elements(tbgrid) 445 445 ; selection des taches au sol situees dans la zone d'interet 446 zone=where((longrid ge lon _min) and (longrid le lon_max) $447 and (latgrid ge lat _min) and (latgrid le lat_max) and (tbgrid gt tbmin) and (tbgrid lt tbmax), npt)446 zone=where((longrid ge lonmin) and (longrid le lonmax) $ 447 and (latgrid ge latmin) and (latgrid le latmax) and (tbgrid gt tbmin) and (tbgrid lt tbmax), npt) 448 448 ;print,'wwwtest: nb points zone d interet ',npt 449 449 if npt ne 0 then begin … … 511 511 , resol : resol $ 512 512 } 513 header2= { lon _min: lon_min $514 , lon _max: lon_max $515 , lat _min: lat_min $516 , lat _max: lat_max $513 header2= { lonmin: lonmin $ 514 , lonmax: lonmax $ 515 , latmin: latmin $ 516 , latmax: latmax $ 517 517 } 518 518 -
trunk/src/extract_amsuab.pro
r517 r525 57 57 ; yyyy 58 58 ; mm 59 ; lon _min60 ; lon _max61 ; lat _min62 ; lat _max59 ; lonmin 60 ; lonmax 61 ; latmin 62 ; latmax 63 63 ; 64 64 ; :Examples: … … 77 77 ; mm=8 78 78 ; dd=1 79 ; lon _min=-60.80 ; lon _max=50.81 ; lat _min=-30.82 ; lat _max=45.83 ; extract_amsuab, yyyy, mm, dd, use_amsua, lon _min, lon_max, lat_min, lat_max79 ; lonmin=-60. 80 ; lonmax=50. 81 ; latmin=-30. 82 ; latmax=45. 83 ; extract_amsuab, yyyy, mm, dd, use_amsua, lonmin, lonmax, latmin, latmax 84 84 ; 85 85 ; ${PROJECT_ID}/AMSU/2006/08/amsub_20060801_060w30s_050e45n.dat must have been created. … … 91 91 ; mm=8 92 92 ; dd=1 93 ; lon _min=-60.94 ; lon _max=50.95 ; lat _min=-30.96 ; lat _max=45.97 ; extract_amsuab, use_amsua, yyyy, mm, dd, lon _min, lon_max, lat_min, lat_max93 ; lonmin=-60. 94 ; lonmax=50. 95 ; latmin=-30. 96 ; latmax=45. 97 ; extract_amsuab, use_amsua, yyyy, mm, dd, lonmin, lonmax, latmin, latmax 98 98 ; 99 99 ; ${PROJECT_ID}/AMSU/2006/08/amsua_amsub_20060801_060w30s_050e45n.dat must have been created. … … 191 191 ;- 192 192 PRO extract_amsuab, yyyy, mm, dd, use_amsua $ 193 , lon _min, lon_max, lat_min, lat_max193 , lonmin, lonmax, latmin, latmax 194 194 ; 195 195 compile_opt idl2, strictarrsubs … … 209 209 look = 'filename' 210 210 scale = 1. 211 geomin = geolocation_to_string_idl(lon _min, lat_min, look, scale)212 geomax = geolocation_to_string_idl(lon _max, lat_max, look, scale)211 geomin = geolocation_to_string_idl(lonmin, latmin, look, scale) 212 geomax = geolocation_to_string_idl(lonmax, latmax, look, scale) 213 213 + string(mm,format='(I2.2)') + '/' 214 214 ; … … 451 451 ; ici on selectionne la zone definie 452 452 ;----------------- 453 zone=where((amalon ge lon _min) and (amalon le lon_max) $454 and (amalat ge lat _min) and (amalat le lat_max) )453 zone=where((amalon ge lonmin) and (amalon le lonmax) $ 454 and (amalat ge latmin) and (amalat le latmax) ) 455 455 456 456 nn=n_elements(zone) -
trunk/src/extract_nc.m
r368 r525 34 34 35 35 % [tps,lon,lat,extract1,extract2]=extract_nc('nomdufichier',{'data1' 'data2'}, 36 % ... 'tps',tps_min,tps_max,'lon',lon _min,lon_max,'lat',lat_min,lat_max);36 % ... 'tps',tps_min,tps_max,'lon',lonmin,lonmax,'lat',latmin,latmax); 37 37 % 38 38 % voir :ref:`cal_hcl.m` pour exemple d'utilisation -
trunk/src/file_amsu_t2_to_mem.pro
r512 r525 5 5 ; ======================= 6 6 ; 7 ; .. function:: file_amsu_t2_to_mem(yyyy, mm, dd, numch, lon _min, lon_max, lat_min, lat_max, testfilename)7 ; .. function:: file_amsu_t2_to_mem(yyyy, mm, dd, numch, lonmin, lonmax, latmin, latmax, testfilename) 8 8 ; 9 9 ; DESCRIPTION … … 30 30 ; mm=8 31 31 ; dd=13 32 ; lon _min=-25.33 ; lon _max=25.34 ; lat _min=-5.35 ; lat _max=20.32 ; lonmin=-25. 33 ; lonmax=25. 34 ; latmin=-5. 35 ; latmax=20. 36 36 ; testfilename='' 37 ; result = file_amsu_t2_to_mem(yyyy, mm, dd, numch, lon _min, lon_max, lat_min, lat_max, testfilename)37 ; result = file_amsu_t2_to_mem(yyyy, mm, dd, numch, lonmin, lonmax, latmin, latmax, testfilename) 38 38 ; 39 39 ; impression de controle : structure et 1re ligne de données:: … … 106 106 ;- 107 107 FUNCTION file_amsu_t2_to_mem, yyyy, mm, dd, numch $ 108 , lon _min, lon_max, lat_min, lat_max $108 , lonmin, lonmax, latmin, latmax $ 109 109 , testfilename 110 110 ; … … 130 130 look = 'filename' 131 131 scale = 1. 132 geomin = geolocation_to_string_idl(lon _min, lat_min, look, scale)133 geomax = geolocation_to_string_idl(lon _max, lat_max, look, scale)132 geomin = geolocation_to_string_idl(lonmin, latmin, look, scale) 133 geomax = geolocation_to_string_idl(lonmax, latmax, look, scale) 134 134 fullfilename = project_id_env+ 'AMSU/' $ 135 135 + string(yyyy,format='(I4.4)') + '/' $ … … 167 167 , resol : 0.0 $ 168 168 } 169 header2= { lon _min: 0.0 $170 , lon _max: 0.0 $171 , lat _min: 0.0 $172 , lat _max: 0.0 $169 header2= { lonmin: 0.0 $ 170 , lonmax: 0.0 $ 171 , latmin: 0.0 $ 172 , latmax: 0.0 $ 173 173 } 174 174 -
trunk/src/file_bathy_to_mem.pro
r521 r525 2 2 ;+ 3 3 ; 4 ; .. function:: file_bathy_to_mem(lon _min, lon_max, lat_min, lat_max)4 ; .. function:: file_bathy_to_mem(lonmin, lonmax, latmin, latmax) 5 5 ; 6 6 ; .. _file_bathy_to_mem.pro: … … 37 37 ; } 38 38 ; 39 ; :param lon _min:40 ; :param lon _max:41 ; :param lat _min:42 ; :param lat _max:39 ; :param lonmin: 40 ; :param lonmax: 41 ; :param latmin: 42 ; :param latmax: 43 43 ; 44 44 ; EXAMPLES … … 47 47 ; ++:: 48 48 ; 49 ; lon _min=-25.50 ; lon _max=25.51 ; lat _min=-5.52 ; lat _max=20.53 ; result = file_bathy_to_mem(lon _min, lon_max, lat_min, lat_max)49 ; lonmin=-25. 50 ; lonmax=25. 51 ; latmin=-5. 52 ; latmax=20. 53 ; result = file_bathy_to_mem(lonmin, lonmax, latmin, latmax) 54 54 ; 55 55 ; SEE ALSO … … 91 91 ;- 92 92 function file_bathy_to_mem $ 93 , lon _min, lon_max, lat_min, lat_max93 , lonmin, lonmax, latmin, latmax 94 94 ; 95 95 compile_opt idl2, strictarrsubs … … 107 107 file=project_id_env+'/MASK/ETOPO1_Ice_g_gmt4.nc' 108 108 ;attention: zone prise en compte limitee par le choix des nos en x et 109 ;y dans initncdf. La region traitee ne doit pas exceder lon _min=-25 &110 ;lon _max=40 & lat_min=-5 & lat_max=40109 ;y dans initncdf. La region traitee ne doit pas exceder lonmin=-25 & 110 ;lonmax=40 & latmin=-5 & latmax=40 111 111 112 112 initncdf,file, zaxisname='toto',xaxisname='lon',yaxisname='lat',XMINMESH=8000,XMAXMESH=15000, YMINMESH=4500,YMAXMESH=8500 113 domdef,lon _min-15., lon_max+15., lat_min-15., lat_max+15.113 domdef,lonmin-15., lonmax+15., latmin-15., latmax+15. 114 114 bathy=read_ncdf('z',file=file) 115 115 bate=bathy.arr gt 0 116 116 xxe1=reform(glamt[*,0]) 117 zonx=where(xxe1 ge lon _min-15. and xxe1 le lon_max+15.)117 zonx=where(xxe1 ge lonmin-15. and xxe1 le lonmax+15.) 118 118 xxe=xxe1[zonx] 119 119 yye1=reform(gphit[0,*]) 120 zony=where(yye1 ge lat _min-15. and yye1 le lat_max+15.)120 zony=where(yye1 ge latmin-15. and yye1 le latmax+15.) 121 121 yye=yye1[zony] 122 122 -
trunk/src/map_amab.pro
r513 r525 13 13 ; jour_ama 14 14 ; flagg 15 ; lon _min16 ; lon _max17 ; lat _min18 ; lat _max15 ; lonmin 16 ; lonmax 17 ; latmin 18 ; latmax 19 19 ; 20 20 ; :Uses: … … 42 42 ; amafile=project_id_env+"/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0110.E0255.B4271112.GC.L1C" 43 43 ; ambfile=project_id_env+"/AMSUBN15/L1C/2006/2006_08_01/NSS.AMBX.NK.D06213.S0110.E0255.B4271112.GC.L1C" 44 ; lon _min=-60.45 ; lon _max=5046 ; lat _min=-3047 ; lat _max=4548 ; map_amab, use_amsua, use_amsub, amafile,ambfile,dam,jour_ama,flagg,lon _min,lon_max,lat_min,lat_max44 ; lonmin=-60. 45 ; lonmax=50 46 ; latmin=-30 47 ; latmax=45 48 ; map_amab, use_amsua, use_amsub, amafile,ambfile,dam,jour_ama,flagg,lonmin,lonmax,latmin,latmax 49 49 ; print, dam 50 50 ; print, jour_ama … … 64 64 ; 65 65 ;- 66 PRO map_amab,use_amsua, use_amsub, amafile,ambfile,dam,jour_ama,flagg,lon _min,lon_max,lat_min,lat_max66 PRO map_amab,use_amsua, use_amsub, amafile,ambfile,dam,jour_ama,flagg,lonmin,lonmax,latmin,latmax 67 67 ; 68 68 ;++compile_opt idl2, strictarrsubs … … 254 254 ;;on definit ici par precaution contre les donnees erronees seulement la zone geographique 255 255 256 zone=where((amalon ge lon _min) and (amalon le lon_max) $257 and (amalat ge lat _min) and (amalat le lat_max) )256 zone=where((amalon ge lonmin) and (amalon le lonmax) $ 257 and (amalat ge latmin) and (amalat le latmax) ) 258 258 259 259 nn=n_elements(zone) -
trunk/src/mem_to_file_amsu_t2.pro
r517 r525 31 31 ; dd=1 32 32 ; testfilename='' 33 ; lon _min=-25.34 ; lon _max=25.35 ; lat _min=-5.36 ; lat _max=20.37 ; amsu_t2 = file_amsu_t2_to_mem(yyyy, mm, dd, numch, lon _min, lon_max, lat_min, lat_max, testfilename)33 ; lonmin=-25. 34 ; lonmax=25. 35 ; latmin=-5. 36 ; latmax=20. 37 ; amsu_t2 = file_amsu_t2_to_mem(yyyy, mm, dd, numch, lonmin, lonmax, latmin, latmax, testfilename) 38 38 ; testfilename='./a' 39 39 ; result = mem_to_file_amsu_t2(amsu_t2, testfilename) … … 130 130 look = 'filename' 131 131 scale = 1. 132 geomin = geolocation_to_string_idl(amsu_t2.header2.lon _min, amsu_t2.header2.lat_min, look, scale)133 geomax = geolocation_to_string_idl(amsu_t2.header2.lon _max, amsu_t2.header2.lat_max, look, scale)132 geomin = geolocation_to_string_idl(amsu_t2.header2.lonmin, amsu_t2.header2.latmin, look, scale) 133 geomax = geolocation_to_string_idl(amsu_t2.header2.lonmax, amsu_t2.header2.latmax, look, scale) 134 134 fullfilename = project_id_env+ 'AMSU/' $ 135 135 + string(amsu_t2.header1.yyyy,format='(I4.4)') + '/' $ -
trunk/src/plot_amsu.pro
r502 r525 15 15 ; :param yyyy: 16 16 ; :param mm: 17 ; :param lon _min:18 ; :param lon _max:19 ; :param lat _min:20 ; :param lat _max:17 ; :param lonmin: 18 ; :param lonmax: 19 ; :param latmin: 20 ; :param latmax: 21 21 ; 22 22 ; EXAMPLES … … 30 30 ; dd=13 31 31 ; resol=1 32 ; lon _min=-60.33 ; lon _max=50.34 ; lat _min=-30.35 ; lat _max=45.36 ; plot_amsu, numch, yyyy, mm, dd, lon _min, lon_max, lat_min, lat_max32 ; lonmin=-60. 33 ; lonmax=50. 34 ; latmin=-30. 35 ; latmax=45. 36 ; plot_amsu, numch, yyyy, mm, dd, lonmin, lonmax, latmin, latmax 37 37 ; 38 38 ; :file:`${PROJECT_ID}/AMSU/2006/08/a5_20060813_060w30s_050e45n.dat` … … 55 55 ; 56 56 ;- 57 pro plot_amsu, numch,yyyy,mm,dd,lon _min,lon_max,lat_min,lat_max57 pro plot_amsu, numch,yyyy,mm,dd,lonmin,lonmax,latmin,latmax 58 58 59 59 compile_opt idl2, strictarrsubs … … 64 64 look = 'filename' 65 65 66 geomin = geolocation_to_string_idl(lon _min, lat_min, look,1)67 geomax = geolocation_to_string_idl(lon _max, lat_max, look,1)66 geomin = geolocation_to_string_idl(lonmin, latmin, look,1) 67 geomax = geolocation_to_string_idl(lonmax, latmax, look,1) 68 68 69 result=file_amsu_t2_to_mem( yyyy,mm,dd,numch,lon _min,lon_max,lat_min,lat_max,testfilename)69 result=file_amsu_t2_to_mem( yyyy,mm,dd,numch,lonmin,lonmax,latmin,latmax,testfilename) 70 70 desc=result.data.desc 71 71 hour=result.data.hour … … 88 88 window,0 89 89 loadct,39 90 plot,lon,lat,xrange=[lon _min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,/nodata90 plot,lon,lat,xrange=[lonmin,lonmax],yrange=[latmin,latmax],xstyle=1,ystyle=1,/nodata 91 91 colscale=intarr(n_elements(tb)) 92 92 for j=0,n_elements(tb)-1 do colscale[j]=(tb[j]-min(tb))*255/(max(tb)-min(tb)) … … 101 101 window,1 102 102 loadct,39 103 plot,lon,lat,xrange=[lon _min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,/nodata103 plot,lon,lat,xrange=[lonmin,lonmax],yrange=[latmin,latmax],xstyle=1,ystyle=1,/nodata 104 104 colscale=intarr(n_elements(tb)) 105 105 for j=0,n_elements(tb)-1 do colscale[j]=(tb[j]-min(tb))*255/(max(tb)-min(tb))
Note: See TracChangeset
for help on using the changeset viewer.