Changeset 55

Timestamp:

08/20/14 18:45:39 (10 years ago)

Author:

lahlod

Message:

new scripts

Location:

trunk/src/scripts_Laura/ARCTIC/Travail_CEN

Files:

• daily_emis_AMSUA_AMSUB_89.py (modified) (5 diffs)
• map_emis_spec_diff_all_arctic.py (added)
• map_monthly_mean_emis_params_sea_ice.py (added)
• param_anomaly.py (modified) (2 diffs)

trunk/src/scripts_Laura/ARCTIC/Travail_CEN/daily_emis_AMSUA_AMSUB_89.py

@@ -53 +53 @@
 
 
-'''
+
 ############################################
 # time evolution (monthly) in a given zone #
@@ -62 +62 @@
 
 # select borders of zone
-yi = 1880. 
-yf = 2280.
-xi = -480.
-xf = -120.
+yi = 320.
+yf =720.
+xi = -1080. 
+xf = -720.
 
 #find corresponding index in xvec and yvec...
@@ -124 +124 @@
     diff_a23 = fichier_a23.variables['e_spec_lamb'][:]
     fichier_a23.close()
-    # AMSUB
+    # AMSUB 89
     fichier_b = Dataset('/net/argos/data/parvati/lahlod/ARCTIC/monthly_GLACE/gridded_data/cartesian_grid/res_20/cartesian_grid_res-20_monthly_data_lamb_spec_near_nadir_AMSUB89_' + month[imo] + '2009.nc', 'r', format = 'NETCDF3_CLASSIC')
     spec_b = fichier_b.variables['e_spec'][:]
@@ -207 +207 @@
 ion()
 figure()
+plot(std_year_spec_a - std_year_spec_b, 'r', label = '[std AMSUA89] - [std AMSUB89]')
+plot(np.zeros([365]), '--k')
+xlim(0, 365)
+ylim(-0.020, 0.020)
+xticks(vec_months, month, rotation = 25)
+#yticks(np.arange(0.4, 1., 0.05))
+fontP = FontProperties()
+fontP.set_size('small')
+legend(loc = 3, prop = fontP)
+grid()
+ylabel('bias of std')
+title('North Canadian Archipelago')
+plt.savefig('/usr/home/lahlod/twice_d/fig_output_ARCTIC/fig_output_sea_ice_study/ice_class_AMSUA89_AMSUB89/evolution_by_zone/monthly_evolution_bias_std_AMSUA89_AMSUB89_zone_North_Canadian_Archipelago.png')
+
+
+
+vec_months = np.array([0, 31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334])
+ion()
+figure()
 subplot(2, 1, 1)
 plot(mean_year_spec_a, '+-r', label = 'AMSUA89')
@@ -230 +249 @@
 #plot(mean_year_lamb_a - mean_year_lamb_b, 'b', label = 'lamb AMSUA - AMSUB')
 plot(np.zeros([365], float), '--k')
-plot(std_year_spec_a, '+-r', label = 'std spec AMSUA89')
+plot(std_year_spec_a , '+-r', label = 'std spec AMSUA89')
 plot(std_year_spec_b, '+-g', label = 'std spec AMSUB89')
 #plot(std_year_lamb_a, 'c', label = 'lamb AMSUA')

trunk/src/scripts_Laura/ARCTIC/Travail_CEN/param_anomaly.py

@@ -57 +57 @@
 ice_diff_anom = np.zeros([4, M, ny, nx], float)
 ice_ratio_anom = np.zeros([4, M, ny, nx], float)
-for ifr in range (0, 4):
+for ifr in range (3, 4):
     print 'frequency ' + frequ[ifr]
     for imo in range (0, M):
         print 'month ' + month[imo]
-        fichier_ice = Dataset('/net/argos/data/parvati/lahlod/ARCTIC/AMSUA_ice_class/sub_classification/cartesian_grid_map_sea_ice_extent_with-AMSUA23-and-89_' + month[imo] + '2009_AMSUA' + str(frequ[ifr]) + '_spec_thresholds.nc', 'r', format='NETCDF3_CLASSIC')
+        fichier_ice = Dataset('/net/argos/data/parvati/lahlod/ARCTIC/AMSUB_ice_class/sub_classification/cartesian_grid_map_sea_ice_extent_with-AMSUA23-and-89_' + month[imo] + '2009_AMSUB' + str(frequ[ifr]) + '_spec_thresholds.nc', 'r', format='NETCDF3_CLASSIC')
         ice_spec = fichier_ice.variables['spec_ice_area'][:]
         ice_lamb = fichier_ice.variables['lamb_ice_area'][:]
@@ -88 +88 @@
         ######################## 
         print 'stack in file month ' + str(month[imo])
-        rootgrp = Dataset('/net/argos/data/parvati/lahlod/ARCTIC/AMSUA_ice_class/sub_classification/cartesian_grid_map_sea_ice_extent_with-AMSUA23-and-89_param_anomaly_' + month[imo] + '2009_AMSUA' + str(frequ[ifr]) + '.nc', 'w', format='NETCDF3_CLASSIC')
+        rootgrp = Dataset('/net/argos/data/parvati/lahlod/ARCTIC/AMSUB_ice_class/sub_classification/cartesian_grid_map_sea_ice_extent_with-AMSUA23-and-89_param_anomaly_' + month[imo] + '2009_AMSUB' + str(frequ[ifr]) + '.nc', 'w', format='NETCDF3_CLASSIC')
         rootgrp.createDimension('longitude', nx)
         rootgrp.createDimension('latitude', ny)