Changeset 382

Timestamp:

11/10/10 12:14:54 (14 years ago)

Author:

brocksce

Message:

All ancillary variables are now kept from original file
Not coordinates variables that are not used, for example PFT(PFT)

File:

• trunk/libIGCM/AA_create_ts (modified) (2 diffs)

trunk/libIGCM/AA_create_ts

@@ -862 +862 @@
                     #
                     if [ ${FlagDir} = "1M" ] ; then
-                        # Compute yearly average considering calendar
                         file_before_YE=${config_UserChoices_JobName}_${DateBegin}_${DATE_FIN_JOB_B}_1Y_${var}.nc
                         file_out_YE=${config_UserChoices_JobName}_${DATE_BUILD_END}_1Y_${var}.nc
                         #
-                        # Set date of the mid year to 01/07 at 00:00
-                        # better than default where the date information in the outputfile is the date of the last contributing time step in the input file
+                        # Extract all ancillary variables except $var, they will be appended after the cdo command
+                        # Needed because cdo applies calculs on all variables
+                        IGCM_sys_ncks -Oh -x -v ${var} ${file_out} ${file_out_YE%%.nc}_tmp1.nc
+                        # This will remove the time axis because cdo rewrites severely with only considered useful variables 
+                        IGCM_sys_cdo -mulc,1 ${file_out_YE%%.nc}_tmp1.nc ${file_out_YE%%.nc}_tmp2.nc
+                        # Compute yearly average considering calendar
+                        # and set date of the mid year to 01/07 at 00:00
                         IGCM_sys_cdo -settime,00:00 -setday,1 -setmon,7 -settunits,days -divdpy -yearsum -muldpm -selvar,${var} ${file_out} ${file_out_YE}
+                        # Now append all ancillary variables
+                        IGCM_sys_ncks -h -A ${file_out_YE%%.nc}_tmp2.nc ${file_out_YE}
+                        # Put global history attribut
+                        IGCM_sys_ncatted -Oh -a history,global,o,c,"build by libIGCM/create_ts" ${file_out_YE}
                         #
                         eval IGCM_sys_Put_Out ${file_out_YE} \${R_OUT_${comp}}/Analyse/TS_MO_YE/${file_out_YE}
@@ -874 +882 @@
                     fi
                     #
-                    IGCM_sys_Rm ${file_out}
+                    IGCM_sys_Rm ${file_out} ${file_out_YE%%.nc}_tmp*.nc
                     #
                 fi