Changeset 482 for trunk


Ignore:
Timestamp:
01/04/12 15:45:40 (12 years ago)
Author:
lelod
Message:

deplacement corrections land/sea

Location:
trunk/src
Files:
5 added
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/extract_amsua.pro

    r481 r482  
    1414; generaliser pour amsub aussi) 
    1515; 
    16 ; decode les noms des fichiers dans la date choisie, puis 
    17 ; appelle le prgm de lecture 
     16; decode les noms des fichiers dans la liste de fichiers, puis 
     17; appelle le prgm de lecture L1C 
    1818; 
    1919; appelle :ref:`interpol_correc`, qui fournit les fonctions de correction au 
    … …  
    2222; 
    2323; applique une interpolation tenant compte de la dimension de la tache 
    24 ; au sol selon la position dans la fauchee 
     24; au sol selon la position dans la fauchee : appel :ref:`interpolswath` 
    2525; 
    2626; ecrit les donnees pour la zone (lat/long) choisie 
    27 ; 
    28 ; ++ le fichier de sortie contient dans l'ordre mm, jj (jour de 
    29 ; l'annee), le temps (identique pour les points d'une meme fauchee), le 
    30 ; no satellite, le no de fov (1 a 40 ou 1 a 90), long, lat, tb 
    31 ; 
    32 ; Les valeurs manquantes sont codées par NaN. 
     27; la structure de sortie contient un en-tete (date, zone, resol) 
     28; puis les donnees du jour et du canal traite 
     29; 
    3330; 
    3431; .. only:: man 
    … …  
    7168; 
    7269;   $ mkdir -p ${PROJECT_ID}/AMSU/2006/08 
    73 ;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0110.E0255.B4271112.GC.L1C" > ${PROJECT_ID}/list_filea 
    74 ;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0250.E0437.B4271213.GC.L1C" >> ${PROJECT_ID}/list_filea 
     70;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0110.E0255.B4271112.GC.L1C" > ${PROJECT_ID}/list_file 
     71;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0250.E0437.B4271213.GC.L1C" >> ${PROJECT_ID}/list_file 
    7572; 
    7673; Using AMSU-A channel 5:: 
    … …  
    192189; called by :ref:`traite_amsuab.sh` 
    193190; 
    194 ; use :ref:`read_amsua1c.pro`, :ref:`read_amsub1c.pro`, 
     191; use :ref:`read_amsua1c`, :ref:`read_amsub1c.pro`, 
    195192; :func:`geolocation_to_string_idl`, :func:`mem_to_file_amsu_t2` 
    196 ; 
    197 ; next step :  ++ 
     193; use :ref:`interpol_correc`  
     194; use :ref:`interpolswath` 
     195; 
     196; next step :  au choix: amsu2ncdf ou cresamsu pour transformation en netcdf 
    198197; 
    199198; TODO 
    … …  
    203202; interpolswath ne gere pas cette situation 
    204203; 
    205 ; comprendre les points bizarres autour de -40,30 visibles avec l'exemple 
    206 ; de 20060801 
    207 ; 
    208 ; lever le doute sur le contenu du fichier écrit par le printf (je (fplod) 
    209 ; crains des zesors louches!!).louche 
    210 ; 
    211204; vérifier que pas NaN dans les lignes écrites 
    212205; 
    213206; faire un test avec un zone géographique dans le pacifique pour voir comment 
    214 ; le système réagi (devrait dire : rien à faire !) 
     207; le système réagit 
    215208; 
    216209; changer le nom des fichiers de listes pour éviter des conflits lors de runs 
    … …  
    221214; eviter la creation d'un fichier vide (si pas de points dans la zone) 
    222215; 
    223 ; tester la cohérence entre les fichiers présents dans list_filea 
     216; tester la cohérence entre les fichiers présents dans list_file 
    224217; en terme de date avec yyyy, mm et dd donnés en paramètre 
    225218; 

  • trunk/src/get_amsu.sh

    r396 r482  
    4545#           amsu_ref_t1 [shape=ellipse,fontname=Courier,label="/bdd/AMSU-1C/AMSU*/L1C/yyyy/yyyy_mm_dd/*.L1C"]; 
    4646# 
    47 #           amsu_ref_a_cor_land [shape=ellipse,fontname=Courier,label="/homedata/eymard/TRAIT_AMSU/CORR_LAND_AMSUA_JUIL2006.DAT"]; 
    48 #           amsu_ref_a_cor_sea [shape=ellipse,fontname=Courier,label="/homedata/eymard/TRAIT_AMSU/CORR_SEA_AMSUA_JUIL2006.DAT"]; 
    49 #           amsu_ref_b_cor_land [shape=ellipse,fontname=Courier,label="/homedata/eymard/TRAIT_AMSU/CORR_LAND_AMSUB_JUIL2006.DAT"]; 
    50 #           amsu_ref_b_cor_sea [shape=ellipse,fontname=Courier,label="/homedata/eymard/TRAIT_AMSU/AMSU/CORR_SEA_AMSUB_JUIL2006.DAT"]; 
    5147# 
    5248#           amsu_t1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/AMSU/AMSU*/L1C/yyyy/yyyy_mm_dd/*.L1C"]; 
    5349# 
    54 #           amsu_a_cor_land [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_LAND_AMSUA_JUIL2006.DAT"]; 
    55 #           amsu_a_cor_sea [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_SEA_AMSUA_JUIL2006.DAT"]; 
    56 #           amsu_b_cor_land [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_LAND_AMSUB_JUIL2006.DAT"]; 
    57 #           amsu_b_cor_sea [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_SEA_AMSUB_JUIL2006.DAT"]; 
    5850# 
    5951#           log [shape=ellipse,fontname=Courier,label="${PROJECT_LOG}/get_amsu.log.YYYYMMDDTHHMMSSZ"]; 
    … …  
    357349           if [ ! -d ${loc} ] 
    358350           then 
    359               mkdir -p ${loc} 
     351              mkdir -p ${loc} 2> /dev/null 
    360352              echo "iii : adding directory ${loc}" >> ${log} 
    361353           fi 
    … …  
    370362fi 
    371363# 
    372 # looking for correction files 
    373 list_url="" 
    374 dirref_base=/homedata/eymard/TRAIT_AMSU/ 
    375 for bfile in CORR_LAND_AMSUA_JUIL2006.DAT CORR_SEA_AMSUA_JUIL2006.DAT CORR_LAND_AMSUB_JUIL2006.DAT CORR_SEA_AMSUB_JUIL2006.DAT 
    376 do 
    377     file=${dirref_base}/${bfile} 
    378     if [ ! -f ${file} ] 
    379     then 
    380         echo "iii : ${file} not found" >> ${log} 2>&1 
    381     else 
    382         echo "iii : ${file} found" >> ${log} 2>&1 
    383         list_url="${list_url} ${file}" 
    384     fi 
    385 done 
    386 # 
    387 if [ "${list_url}" = "" ] 
    388 then 
    389    echo "iii : no files founds" >> ${log} 2>&1 
    390 else 
    391     # 
    392     # if file is not already in ${PROJECT_ID}, get it 
    393     for url in ${list_url} 
    394     do 
    395         file=$(echo ${url} |  sed -e "s+${dirref_base}+${PROJECT_ID}/AMSU/+") 
    396         if [ ! -f ${file} ] 
    397         then 
    398            loc=$(dirname ${file}) 
    399            if [ ! -d ${loc} ] 
    400            then 
    401               mkdir -p ${loc} 
    402               echo "iii : adding directory ${loc}" >> ${log} 
    403            fi 
    404            ln -s ${url} ${file} 
    405            echo "iii : adding file ${file}" >> ${log} 
    406         else 
    407            echo "iii : ${file} exists" >> ${log} 
    408         fi 
    409     done 
    410     unset file 
    411     unset url 
    412 fi 
    413 unset list_url 
     364 
    414365# 
    415366unset command 

  • trunk/src/interpol_correc.pro

    r475 r482  
    2626;        digraph interpol_correc { 
    2727; 
    28 ;           amsu_b_cor_land [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_LAND_AMSUB_JUIL2006.DAT"]; 
    29 ;           amsu_b_cor_sea [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_SEA_AMSUB_JUIL2006.DAT"]; 
     28;           amsu_b_cor_land [shape=diamond,fontname=Courier,label="${PROJECT}/src/dataref_amsu/CORR_LAND_AMSUB_JUIL2006.DAT"]; 
     29;           amsu_b_cor_sea [shape=diamond,fontname=Courier,label="${PROJECT}/src/dataref_amsu/CORR_SEA_AMSUB_JUIL2006.DAT"]; 
    3030; 
    3131;           interpol_correc [shape=box, 
    … …  
    160160 
    161161 
    162 path=project_id_env 
     162path=project_env + 'src/dataref_amsu/' 
    163163; lecture des fichiers d'ajustement au nadir pour AMSUA et B 
    164164; cas amsua 
    165165if nomcanal eq 'a' then begin 
    166    cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA_JUIL2006.DAT' 
    167    cor_seaa=project_id_env+'AMSU/CORR_SEA_AMSUA_JUIL2006.DAT' 
     166   cor_landa=path + '/CORR_LAND_AMSUA_JUIL2006.DAT' 
     167   cor_seaa=path + '/CORR_SEA_AMSUA_JUIL2006.DAT' 
    168168  ;test de l'autre fichier de correction: quel est le bon??? 
    169169  ; cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA.DAT' 
    … …  
    213213;cas amsub 
    214214if nomcanal eq 'b' then begin 
    215    cor_landb=path+'AMSU/CORR_LAND_AMSUB_JUIL2006.DAT' 
    216    cor_seab=path+'AMSU/CORR_SEA_AMSUB_JUIL2006.DAT' 
     215   cor_landb=path+'CORR_LAND_AMSUB_JUIL2006.DAT' 
     216   cor_seab=path+'CORR_SEA_AMSUB_JUIL2006.DAT' 
    217217   ;cor_landb=path+'AMSU/CORR_LAND_AMSUB.DAT' 
    218218   ;cor_seab=path+'AMSU/CORR_SEA_AMSUB.DAT' 

  • trunk/src/traite_amsuab.sh

    r478 r482  
    5555# 
    5656#           amsu_t1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/AMSU/AMSU*/L1C/yyyy/yyyy_mm_dd/*.L1C"]; 
    57 #           amsu_a_cor_land [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_LAND_AMSUA_JUIL2006.DAT"]; 
    58 #           amsu_a_cor_sea [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_SEA_AMSUA_JUIL2006.DAT"]; 
    59 #           amsu_b_cor_land [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_LAND_AMSUB_JUIL2006.DAT"]; 
    60 #           amsu_b_cor_sea [shape=diamond,fontname=Courier,label="${PROJECT_ID}/AMSU/CORR_SEA_AMSUB_JUIL2006.DAT"]; 
    61 #           amsu_t2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/yyyy/mm/amsua_amsub_yyyymmdd_geomin_geomax.dat"]; 
    62 #           amsu_t3 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/yyyy/mm/amsua_amsub_yyyymmdd_geomin_geomax_nadir.dat"]; 
    63 #           amsu_t4 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/yyyy/mm/amsua_amsub_yyyymmdd_geomin_geomax_nadir.nc"]; 
     57#           amsu_a_cor_land [shape=diamond,fontname=Courier,label="${PROJECT}/src/dataref_amsu/CORR_LAND_AMSUA_JUIL2006.DAT"]; 
     58#           amsu_a_cor_sea [shape=diamond,fontname=Courier,label="${PROJECT}/src/dataref_amsu/CORR_SEA_AMSUA_JUIL2006.DAT"]; 
     59#           amsu_b_cor_land [shape=diamond,fontname=Courier,label="${PROJECT}/src/dataref_amsu/CORR_LAND_AMSUB_JUIL2006.DAT"]; 
     60#           amsu_b_cor_sea [shape=diamond,fontname=Courier,label="${PROJECT}/src/dataref_amsu/CORR_SEA_AMSUB_JUIL2006.DAT"]; 
     61#           amsu_t2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/yyyy/mm/numch_yyyymmdd_geomin_geomax.dat"]; 
     62#           amsu_t3 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/yyyy/mm/numch_yyyymmdd_geomin_geomax_nadir.dat"]; 
     63#           amsu_t4 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/yyyy/mm/numch_yyyymmdd_geomin_geomax_nadir.nc"]; 
    6464#           mask [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/MASK/ETOPO1_Ice_g_gmt4.nc"]; 
    6565# 
    … …  
    9292# If you want to use only AMSU-B:: 
    9393# 
    94 #   $ traite_amsuab.sh -b 20060813 -e 20060814 -lonmin -25 -lonmax 25 -latmin -5 -latmax 20 -numch b5 
     94#   $ traite_amsuab.sh -b 20060813 -e 20060815 -lonmin -25 -lonmax 25 -latmin -5 -latmax 20 -numch a5 
    9595# 
    9696# SEE ALSO 
    … …  
    117117# ==== 
    118118# 
    119 # still not working for use_amsua = 1 because of pb inside correct_nadir_amsu 
     119# test si list_file non vide pour lancer IDL 
    120120# 
    121121# tri chrono dans les list_filea ou apres la production du fichier .dat  
    … …  
    129129# /homedata/pinsard/varamma_d/AMSU/AMSUBN17/L1C/2006/2006_08_01/NSS.AMBX.NM.D06213.S1029.E1224.B2132223.WI.L1C 
    130130# avec use_amsua 0 
    131 # 
    132 # why need .compile ?! 
    133131# 
    134132# build/complete catalog 
    … …  
    465463    # write the file ${PROJECT_ID}/list_filea containing AMSU-A file on the current date 
    466464    ls ${PROJECT_ID}/AMSU/AMSU${prefix}*/L1C/${yyyy}/${yyyy}_${mm}_${dd}/*.L1C > ${PROJECT_ID}/list_file 
    467  
     465    # 
     466    # build output directory if needed 
     467    mkdir -p ${PROJECT_ID}/AMSU/${yyyy}/${mm} 2> /dev/null 
     468    # 
     469    # build IDL commands file 
    468470    cat >> ${VARAMMA}/src/traite_amsu_${$}.pro << EOF 
    469471.compile interpol_correc 
    … …  
    492494    fi 
    493495    rm ${VARAMMA}/src/traite_amsu_${$}.pro 
    494     rm  ${PROJECT_ID}/list_file 2> /dev/null  
     496    rm ${PROJECT_ID}/list_file 2> /dev/null  
    495497    unset status 
    496498    if [ ${netcdf_build} -eq 1 ] 
Note: See TracChangeset for help on using the changeset viewer.