Changeset 4518 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro
 Timestamp:
 20140225T16:20:43+01:00 (7 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro
r4517 r4518 105 105 ;climatology 106 106 QNET = read_ncdf(getenv('VAR_FLUX'), filename = std_file_oaflux ) 107 107 108 108 ; erp (evaporation damping) 109 109 ERP1 = read_ncdf(getenv('VAR1_ERP'), allrecords = allrec, direc = 't', filename = std_file1_T ) … … 186 186 april = where(mm EQ 4, cnt) 187 187 Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 188 jan = where(mm EQ 1, cnt) 188 189 sept = where(mm EQ 9, cnt) 190 Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 189 191 Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 190 192 undefine, Ithi_1 … … 195 197 april = where(mm EQ 4, cnt) 196 198 Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 199 jan = where(mm EQ 1, cnt) 197 200 sept = where(mm EQ 9, cnt) 198 201 Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 202 Ithi_jan_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 199 203 undefine, Ithi_2 200 204 ENDIF ELSE BEGIN 201 205 Ithi_april_2 = {arr:1} 202 206 Ithi_sept_2 = {arr:1} 207 Ithi_jan_2 = {arr:1} 203 208 ENDELSE 204 209 ; 205 210 Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 206 211 caldat, time, mm 212 febr = where(mm EQ 2, cnt) 213 Iage_febr_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 207 214 march = where(mm EQ 3, cnt) 208 215 Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} … … 214 221 Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 215 222 caldat, time, mm 223 febr = where(mm EQ 2, cnt) 224 Iage_febr_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 216 225 march = where(mm EQ 3, cnt) 217 226 Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_2.unit} … … 220 229 undefine, Iage_2 221 230 ENDIF ELSE BEGIN 231 Iage_febr_2 = {arr:1} 222 232 Iage_march_2 = {arr:1} 223 233 Iage_sept_2 = {arr:1} … … 225 235 ; 226 236 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 227 ;; help, Ifra_1.arr, jpi,jpj228 237 caldat, time, mm 238 febr = where(mm EQ 2, cnt) 239 Ifra_febr_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 229 240 march = where(mm EQ 3, cnt) 230 241 Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} … … 236 247 Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 237 248 caldat, time, mm 249 febr = where(mm EQ 2, cnt) 250 Ifra_febr_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 238 251 march = where(mm EQ 3, cnt) 239 252 Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} … … 242 255 undefine, Ifra_2 243 256 ENDIF ELSE BEGIN 257 Ifra_febr_2 = {arr:1} 244 258 Ifra_march_2 = {arr:1} 245 259 Ifra_sept_2 = {arr:1} … … 248 262 Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 249 263 caldat, time, mm 264 febr = where(mm EQ 2, cnt) 265 Isnow_febr_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 250 266 march = where(mm EQ 3, cnt) 251 267 Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} … … 257 273 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 258 274 caldat, time, mm 275 febr = where(mm EQ 2, cnt) 276 Isnow_febr_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 259 277 march = where(mm EQ 3, cnt) 260 278 Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} … … 263 281 ; undefine, Isnow_2 264 282 ENDIF ELSE BEGIN 283 Isnow_febr_2 = {arr:1} 265 284 Isnow_march_2 = {arr:1} 266 285 Isnow_sept_2 = {arr:1} … … 450 469 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 451 470 ; 452 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH'453 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex471 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: JAN' 472 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_jan_1, Ithi_jan_2, /ARC, /JAN, _extra = ex 454 473 ; 455 474 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 456 475 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 457 ; 458 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: SEPT' 459 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 460 ; 461 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 462 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 463 ; 476 ; 464 477 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 465 478 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex … … 475 488 ; 476 489 cnt = cnt+1 & blabla = 'Antarctic Ice Age: MARCH' 477 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_ march_1, Iage_march_2, /ANT, /MARCH, _extra = ex490 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_febr_1, Iage_febr_2, /ANT, /FEBR, _extra = ex 478 491 ; 479 492 cnt = cnt+1 & blabla = 'Antarctic Ice Age: SEPT' … … 487 500 ; 488 501 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: MARCH' 489 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_ march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex502 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_febr_1, Ifra_febr_2, /ANT, /FEBR, _extra = ex 490 503 ; 491 504 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' … … 499 512 ; 500 513 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: MARCH' 501 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_ march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex514 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_febr_1, Isnow_febr_2, /ANT, /FEBR, _extra = ex 502 515 ; 503 516 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT'
Note: See TracChangeset
for help on using the changeset viewer.