Changeset 3888 for branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

Timestamp:

2013-04-23T10:37:54+02:00 (11 years ago)

Author:

flavoni

Message:

update DIL_scripts to add plot of ICE extent and observations

File:

• branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro (modified) (6 diffs)

branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

@@ -28 +28 @@
   std_file_oaflux =  isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato)
   std_file_mld =  isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato)
+  std_file_ice =  isafile(getenv('FILE_ICE'), title = 'ICE', iodir = std_iodir_climato)
 
   IF strlowcase(getenv('FILE_GEOHEAT')) EQ 'no' THEN std_file_geoheat = 'no' $
@@ -160 +161 @@
   Ithi_sept_1 = {arr:1./float(cnt) * total(Ithi_1.arr[*, *, temporary(sept)], 3), unit:Ithi_1.unit} 
   undefine, Ithi_1
+;simona
+;ICE_EXT = read_ncdf(getenv('VAR_ICE'), filename = std_file_ice )
+;simona
 ;
   IF std_file2_I NE std_file1_I THEN BEGIN
@@ -175 +179 @@
 ;
   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
+  help, Ifra_1.arr, jpi,jpj
   caldat, time, mm
   march = where(mm EQ 3, cnt)
@@ -180 +185 @@
   sept = where(mm EQ 9, cnt)
   Ifra_sept_1 = {arr:1./float(cnt) * total(Ifra_1.arr[*, *, temporary(sept)], 3), unit:Ifra_1.unit}
-  undefine, Ifra_1
+  ;undefine, Ifra_1
 ;
   IF std_file2_I NE std_file1_I THEN BEGIN
@@ -189 +194 @@
     sept = where(mm EQ 9, cnt)
     Ifra_sept_2 = {arr:1./float(cnt) * total(Ifra_2.arr[*, *, temporary(sept)], 3), unit:Ifra_2.unit}
-    undefine, Ifra_2
+  ;  undefine, Ifra_2
   ENDIF ELSE BEGIN 
     Ifra_march_2 = {arr:-1}
@@ -323 +328 @@
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex
 
+;  cnt = cnt+1   &   blabla = 'Arctic Ice Extend: '
+;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex
+;  
+;  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: '
+;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex
+;
   cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex