Changeset 478 for trunk

Timestamp:

01/04/12 12:27:56 (12 years ago)

Author:

lelod

Message:

adapt traite_amsuam.sh

Location:

trunk/src

Files:

• extract_amsua.pro (modified) (1 diff)
• traite_amsuab.sh (modified) (9 diffs)

trunk/src/extract_amsua.pro

@@ -266 +266 @@
 
 ; lecture de la liste des fichiers AMSU-A
-    ;list_filea = project_id_env+'list_filea'
-list_filea='list_filea'
+list_filea = project_id_env+'list_file'
+;list_filea='list_file'
     ;index_filea= 0
 nb_filea = file_lines(list_filea)

trunk/src/traite_amsuab.sh

@@ -92 +92 @@
 # If you want to use only AMSU-B::
 #
-#   $ traite_amsuab.sh -b 20060801 -e 20060802 -lonmin -25 -lonmax 25 -latmin -5 -latmax 20 -use_amsua 0
+#   $ traite_amsuab.sh -b 20060813 -e 20060814 -lonmin -25 -lonmax 25 -latmin -5 -latmax 20 -numch b5
 #
 # SEE ALSO
@@ -107 +107 @@
 #
 # Use :
-# :ref:`extract_amsuab.pro`
-# :ref:`correct_nadir_amsu.pro`
+# :ref:`extract_amsua.pro`
 # :ref:`amsu2ncdf.pro`
 #
@@ -259 +258 @@
 log_date=$(date -u +"%Y%m%dT%H%M%SZ")
 #
-usage=" Usage : ${command} [-f] -b yyyymmdd -e yyyymmdd -lonmin lonmin -lonmax lonmax -latmin latmin -latmax latmax [-use_amsua use_amsua]"
+usage=" Usage : ${command} [-f] -b yyyymmdd -e yyyymmdd -lonmin lonmin -lonmax lonmax -latmin latmin -latmax latmax -numch numch"
 #
 hostname=$(hostname)
@@ -271 +270 @@
 yyyymmddb=20060420
 yyyymmdde=20060420
-use_amsua=1
+numch='a5'
 jdeb=1
 jfin=10
+netcdf_build=0
 #
 minargcount=4
@@ -313 +313 @@
         shift
      ;;
-     -use_amsua)
-        use_amsua=${2}
+     -numch)
+        numch=${2}
         shift
      ;;
@@ -380 +380 @@
 fi
 #
-case ${use_amsua} in
-    0|1)
-    ;;
-    *)
-     echo "${command} : eee : use_amsua ${use_amsua} argument invalid"
-     exit 1
-    ;;
-esac
 #
 tool=${IDL_CMD}
@@ -435 +427 @@
 echo "latmin=${latmin}" 1>> ${log}
 echo "latmax=${latmax}" 1>> ${log}
-echo "use_amsua=${use_amsua}"  1>> ${log}
+echo "numch=${numch}"  1>> ${log}
 echo "" 1>> ${log}
 #
@@ -469 +461 @@
     mn=$(${date_cmd} -u -d "1970-01-01 ${current} sec" +%M)
     #
-    case ${use_amsua} in
-       0)
-          prefix=amsub_ch5_
-       ;;
-       1)
-          prefix=amsua_amsub_ch5_
-          # write the file ${PROJECT_ID}/list_filea containing AMSU-A file on the current date
-          ls ${PROJECT_ID}/AMSU/AMSUA*/L1C/${yyyy}/${yyyy}_${mm}_${dd}/*.L1C > ${PROJECT_ID}/list_filea
-       ;;
-    esac   
+    prefix=$(printf  "%1c" ${numch} | tr [:lower:] [:upper:])
+ 
     # write the file ${PROJECT_ID}/list_filea containing AMSU-A file on the current date
-    ls ${PROJECT_ID}/AMSU/AMSUB*/L1C/${yyyy}/${yyyy}_${mm}_${dd}/*.L1C > ${PROJECT_ID}/list_fileb
-    #
-#++more ${PROJECT_ID}/list_fileb
-#++read a 
-    cat >> ${VARAMMA}/src/traite_amsuab_${$}.pro << EOF
-.compile correct_nadir_amsu
+    ls ${PROJECT_ID}/AMSU/AMSU${prefix}*/L1C/${yyyy}/${yyyy}_${mm}_${dd}/*.L1C > ${PROJECT_ID}/list_file
+
+    cat >> ${VARAMMA}/src/traite_amsu_${$}.pro << EOF
 .compile interpol_correc
-extract_amsuab, ${yyyy},${mm}, ${dd}, ${use_amsua}, ${lonmin}, ${lonmax}, ${latmin}, ${latmax}
-correct_nadir_amsu, ${use_amsua}, ${yyyy}, ${mm},${dd}, ${lonmin}, ${lonmax}, ${latmin}, ${latmax}
+.compile extract_amsua
+resol=1
+extract_amsua, ${numch},  ${yyyy},${mm}, ${dd}, resol, ${lonmin}, ${lonmax}, ${latmin}, ${latmax}
+EOF
+    if [ ${netcdf_build} -eq 1 ]
+    then
+        cat >> ${VARAMMA}/src/traite_amsu_${$}.pro << EOF
 amsu2ncdf, ${use_amsua}, ${yyyy},${mm},${dd}, ${lonmin}, ${lonmax}, ${latmin}, ${latmax}
+EOF
+     fi
+     cat >> ${VARAMMA}/src/traite_amsu_${$}.pro << EOF
 exit
 EOF
     ${IDL_CMD} << EOF >> ${log} 2>&1
-     @traite_amsuab_${$}
+     @traite_amsu_${$}
 EOF
     status=${?}
@@ -499 +488 @@
     then
        echo "eee: ${IDL_CMD} failed : ${status}" >> ${log} 2>&1
-       cat ${VARAMMA}/src/traite_amsuab_${$}.pro >> ${log} 2>&1
+       cat ${VARAMMA}/src/traite_amsu_${$}.pro >> ${log} 2>&1
        exit 1
     fi
-    rm ${VARAMMA}/src/traite_amsuab_${$}.pro
-    rm  ${PROJECT_ID}/list_filea 2> /dev/null
-    rm  ${PROJECT_ID}/list_fileb 
+    rm ${VARAMMA}/src/traite_amsu_${$}.pro
+    rm  ${PROJECT_ID}/list_file 2> /dev/null 
     unset status
-    # concatenate yyyymmdd-hhmnss NetCDF files
-    list=${PROJECT_ID}/AMSU/${yyyy}/${mm}/${prefix}${yyyy}${mm}${dd}_??????_${geomin}_${geomax}_nadir.nc
-    fileout=${PROJECT_ID}/AMSU/${yyyy}/${mm}/${prefix}${yyyy}${mm}${dd}_${geomin}_${geomax}_nadir.nc
-    ncrcat -O ${list} ${fileout} >> ${log} 2>&1
-    status=${?}
-    if [ ${status} -ne 0 ]
+    if [ ${netcdf_build} -eq 1 ]
     then
-       echo "eee : pb with ncrcat" 1>> ${log}
-       for file in ${list}
-       do
-           ncdump -h ${file} >> ${log} 2>&1
-       done
-       exit 1
-    else
-       echo "iii : adding ${fileout}" >> ${log} 2>&1
+         # concatenate yyyymmdd-hhmnss NetCDF files
+         list=${PROJECT_ID}/AMSU/${yyyy}/${mm}/${prefix}${yyyy}${mm}${dd}_??????_${geomin}_${geomax}_nadir.nc
+         fileout=${PROJECT_ID}/AMSU/${yyyy}/${mm}/${prefix}${yyyy}${mm}${dd}_${geomin}_${geomax}_nadir.nc
+         ncrcat -O ${list} ${fileout} >> ${log} 2>&1
+         status=${?}
+         if [ ${status} -ne 0 ]
+         then
+             echo "eee : pb with ncrcat" 1>> ${log}
+             for file in ${list}
+             do
+                 ncdump -h ${file} >> ${log} 2>&1
+             done
+             exit 1
+          else
+             echo "iii : adding ${fileout}" >> ${log} 2>&1
+          fi
     fi
     unset dd